Index: /issm/workshop/2012/Talks/07_Ice_flow_models/Code/eismint.m
===================================================================
--- /issm/workshop/2012/Talks/07_Ice_flow_models/Code/eismint.m	(revision 14129)
+++ /issm/workshop/2012/Talks/07_Ice_flow_models/Code/eismint.m	(revision 14130)
@@ -1,21 +1,24 @@
 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