Changeset 8305 for issm/trunk/test/Par/SquareEISMINT.par

Timestamp:

05/16/11 15:50:32 (14 years ago)

Author:

seroussi

Message:

grid in par files

File:

• issm/trunk/test/Par/SquareEISMINT.par (modified) (4 diffs)

issm/trunk/test/Par/SquareEISMINT.par ¶

@@ -4 +4 @@
 ymin=min(md.y);
 ymax=max(md.y);
-md.thickness=500*ones(md.numberofgrids,1);
-md.firn_layer=0*ones(md.numberofgrids,1);
+md.thickness=500*ones(md.numberofnodes,1);
+md.firn_layer=0*ones(md.numberofnodes,1);
 md.bed=-md.rho_ice/md.rho_water*md.thickness;
 md.surface=md.bed+md.thickness;
@@ -11 +11 @@
 disp('      creating drag');
 md.drag_type=2; %0 none 1 plastic 2 viscous
-md.drag_coefficient=200*ones(md.numberofgrids,1); %q=1.
+md.drag_coefficient=200*ones(md.numberofnodes,1); %q=1.
 %Take care of iceshelves: no basal drag
 pos=find(md.elementoniceshelf);
@@ -19 +19 @@
 
 disp('      creating temperature');
-md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.observed_temperature=(273-20)*ones(md.numberofnodes,1);
 
 disp('      creating flow law paramter');
 %md.B=paterson(md.observed_temperature);
-md.rheology_B=1.7687*10^8*ones(md.numberofgrids,1);
+md.rheology_B=1.7687*10^8*ones(md.numberofnodes,1);
 md.rheology_n=3*ones(md.numberofelements,1);
 
 disp('      creating accumulation rates');
-md.accumulation_rate=0.2*ones(md.numberofgrids,1); %0m/a
-md.melting_rate=0*ones(md.numberofgrids,1); %0m/a
+md.accumulation_rate=0.2*ones(md.numberofnodes,1); %0m/a
+md.melting_rate=0*ones(md.numberofnodes,1); %0m/a
 
 disp('      boundary conditions ');
@@ -37 +37 @@
 
 %Evolution of the ice shelf
-pos=find(md.y==200000); %grids on the upper boundary condition
+pos=find(md.y==200000); %nodes on the upper boundary condition
 md.spcthickness(pos,1)=1;
 md.spcthickness(pos,2)=500;