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

Timestamp:

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

Author:

seroussi

Message:

grid in par files

File:

• issm/trunk/test/Par/ISMIPE.par (modified) (3 diffs)

issm/trunk/test/Par/ISMIPE.par

@@ -4 +4 @@
 data=load('../Data/ISMIPE.data','-mat');
 data=data.data;
-md.surface=zeros(md.numberofgrids,1);
-md.bed=zeros(md.numberofgrids,1);
-for i=1:md.numberofgrids
+md.surface=zeros(md.numberofnodes,1);
+md.bed=zeros(md.numberofnodes,1);
+for i=1:md.numberofnodes
         y=md.y(i);
         point1=floor(y/100)+1;
@@ -15 +15 @@
 end
 md.thickness=md.surface-md.bed;
-md.firn_layer=0*ones(md.numberofgrids,1);
+md.firn_layer=0*ones(md.numberofnodes,1);
 md.thickness(find(~md.thickness))=0.01;
 md.bed=md.surface-md.thickness;
@@ -21 +21 @@
 disp('      creating drag');
 md.drag_type=2; %0 none 1 plastic 2 viscous
-md.drag_coefficient=zeros(md.numberofgrids,1);
+md.drag_coefficient=zeros(md.numberofnodes,1);
 md.drag_p=ones(md.numberofelements,1);
 md.drag_q=ones(md.numberofelements,1);
 
 disp('      creating flow law paramter');
-md.rheology_B=6.8067*10^7*ones(md.numberofgrids,1);
+md.rheology_B=6.8067*10^7*ones(md.numberofnodes,1);
 md.rheology_n=3*ones(md.numberofelements,1);