Changeset 8312

Timestamp:

05/16/11 16:35:47 (14 years ago)

Author:

seroussi

Message:

forcing test

Location:

issm/trunk/examples/SquareIceShelf

Files:

• Square.par (modified) (2 diffs)
• runme.m (added)

issm/trunk/examples/SquareIceShelf/Square.par

@@ -13 +13 @@
 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 +19 @@
 
 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');