Changeset 14130

Timestamp:

12/11/12 10:22:48 (12 years ago)

Author:

seroussi

Message:

Added solution

File:

• issm/workshop/2012/Talks/07_Ice_flow_models/Code/eismint.m (modified) (1 diff)

issm/workshop/2012/Talks/07_Ice_flow_models/Code/eismint.m

@@ -1 +1 @@
 md=model();
 
-%Create mesh with roundmesh and remove 'return;'
-return;
+%Create mesh with roundmesh
+md=roundmesh(md,750000,30000);
 
 %Set mask
-return;
+md=setmask(md,'','');
 
 %Parameterize model
-return;
+md=parameterize(md,'EISMINT.par');
 
 %We extrude the model to have a 3d model
-return;
+md=extrude(md,10,1);
 
 %Set ice flow approximation
-return;
+md=setflowequation(md,’hutter’,'all');
 
 %Create boundary conditions: zero velocity on the bed
-return;
+pos=find(md.mesh.vertexonbed);
+md.diagnostic.spcvx(pos)=0;
+md.diagnostic.spcvy(pos)=0;
+md.diagnostic.spcvz(pos)=0;
 
 %Go Solve