Index: /issm/trunk/examples/SquareIceShelf/Square.par
===================================================================
--- /issm/trunk/examples/SquareIceShelf/Square.par	(revision 8311)
+++ /issm/trunk/examples/SquareIceShelf/Square.par	(revision 8312)
@@ -13,5 +13,5 @@
 md.drag_type=2; %0 none 1 plastic 2 viscous
 pos=find(md.elementoniceshelf);
-md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
+md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 md.drag_coefficient(md.elements(pos,:))=0;
 md.drag_p=ones(md.numberofelements,1);
@@ -19,5 +19,11 @@
 
 disp('      creating temperature');
-md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
+
+disp('      initial velocity');
+md.vx=zeros(md.numberofnodes,1);
+md.vy=zeros(md.numberofnodes,1);
+md.vz=zeros(md.numberofnodes,1);
+md.vel=zeros(md.numberofnodes,1);
 
 disp('      creating flow law paramter');
Index: /issm/trunk/examples/SquareIceShelf/runme.m
===================================================================
--- /issm/trunk/examples/SquareIceShelf/runme.m	(revision 8312)
+++ /issm/trunk/examples/SquareIceShelf/runme.m	(revision 8312)
@@ -0,0 +1,10 @@
+md=model;
+md=mesh(md,'DomainOutline.exp',100000);
+md=geography(md,'all','');
+md=parameterize(md,'Square.par');
+md=setelementstype(md,'macayeal','all');
+md.dt=1;
+md.ndt=3;
+md.forcing_accumulation=[ones(md.numberofnodes,1),2*ones(md.numberofnodes,1),3*ones(md.numberofnodes,1)];
+md.forcing_accumulation(end+1,:)=[0.5,1.5,2.5]*md.yts;
+md=solve(md,Transient2dSolutionEnum);