Changeset 9684 for issm/trunk/test/Par/SquareThermal.par

Timestamp:

09/08/11 10:57:56 (14 years ago)

Author:

seroussi

Message:

added initialization class

File:

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

issm/trunk/test/Par/SquareThermal.par

@@ -10 +10 @@
 
 disp('      creating velocities');
-md.vx=zeros(md.numberofnodes,1);
-md.vy=zeros(md.numberofnodes,1);
-md.vz=zeros(md.numberofnodes,1);
+md.initialization.vx=zeros(md.numberofnodes,1);
+md.initialization.vy=zeros(md.numberofnodes,1);
+md.initialization.vz=zeros(md.numberofnodes,1);
 
 disp('      creating drag');
@@ -23 +23 @@
 
 disp('      creating temperatures');
-md.temperature=(273-20)*ones(md.numberofnodes,1);
+md.initialization.temperature=(273-20)*ones(md.numberofnodes,1);
 
 disp('      creating flow law paramter');
-md.materials.rheology_B=paterson(md.temperature);
+md.materials.rheology_B=paterson(md.initialization.temperature);
 md.materials.rheology_n=3*ones(md.numberofelements,1);
 
@@ -39 +39 @@
 
 disp('      boundary conditions for thermal model');
-md.thermal.spctemperature(:)=md.temperature;
+md.thermal.spctemperature(:)=md.initialization.temperature;
 md.basalforcings.geothermalflux=zeros(md.numberofnodes,1); 
 pos=find(md.mask.elementongroundedice);md.basalforcings.geothermalflux(md.elements(pos,:))=1*10^-3; %1 mW/m^2