Changeset 9684 for issm/trunk/test/NightlyRun/test1302.m

Timestamp:

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

Author:

seroussi

Message:

added initialization class

File:

• issm/trunk/test/NightlyRun/test1302.m (modified) (3 diffs)

issm/trunk/test/NightlyRun/test1302.m

@@ -14 +14 @@
 pos1=find(md.elementonbed);     md.thermal.spctemperature(md.elements(pos1,1:3))=10;
 pos2=find(md.elementonsurface); md.thermal.spctemperature(md.elements(pos2,4:6))=0;
-md.vz=0.1*ones(md.numberofnodes,1);
-md.vel=sqrt( md.vx.^2+ md.vy.^2+ md.vz.^2);
-md.pressure=zeros(md.numberofnodes,1);
+md.initialization.vz=0.1*ones(md.numberofnodes,1);
+md.initialization.vel=sqrt( md.initialization.vx.^2+ md.initialization.vy.^2+ md.initialization.vz.^2);
+md.initialization.pressure=zeros(md.numberofnodes,1);
 
 %analytical results
@@ -23 +23 @@
 A=10/(exp(alpha*(-1000))-1);    %A=T(bed)/(exp(alpha*bed)-1)  with bed=-1000 T(bed)=10
 B=-A;
-md.temperature=A*exp(alpha*md.z)+B;
+md.initialization.temperature=A*exp(alpha*md.z)+B;
 
 %modeled  results
@@ -32 +32 @@
 comp_temp=zeros(md.numberofnodes,1);
 comp_temp(md.results.ThermalSolution.Temperature.index)=md.results.ThermalSolution.Temperature.value;
-relative=abs((comp_temp-md.temperature)./md.temperature)*100;
-relative(find(comp_temp==md.temperature))=0;
+relative=abs((comp_temp-md.initialization.temperature)./md.initialization.temperature)*100;
+relative(find(comp_temp==md.initialization.temperature))=0;
 plotmodel(md,'data',comp_temp,'title','Modeled temperature [K]','data','temperature','view',3,...
-        'title','Analytical temperature [K]','view',3,'data',comp_temp-md.temperature,...
+        'title','Analytical temperature [K]','view',3,'data',comp_temp-md.initialization.temperature,...
         'title','Absolute error [K]','view',3,'data',relative,'title','Relative error [%]','view',3,...
         'figposition','mathieu','FontSize#all',20)