Changeset 9679 for issm/trunk/test/NightlyRun/test236.m

Timestamp:

09/08/11 08:07:18 (14 years ago)

Author:

seroussi

Message:

added diagnostic class

File:

• issm/trunk/test/NightlyRun/test236.m (modified) (2 diffs)

issm/trunk/test/NightlyRun/test236.m

@@ -34 +34 @@
 
 %Boundary conditions:
-md.spcvx=NaN*ones(md.numberofnodes,1);
-md.spcvy=NaN*ones(md.numberofnodes,1);
-md.spcvz=NaN*ones(md.numberofnodes,1);
+md.diagnostic.spcvx=NaN*ones(md.numberofnodes,1);
+md.diagnostic.spcvy=NaN*ones(md.numberofnodes,1);
+md.diagnostic.spcvz=NaN*ones(md.numberofnodes,1);
 
 %constrain flanks to 0 normal velocity
 pos=find(md.x==xmin | md.x==xmax);
-md.spcvx(pos)=0;
-md.spcvz(pos)=NaN;
+md.diagnostic.spcvx(pos)=0;
+md.diagnostic.spcvz(pos)=NaN;
 
 %constrain grounding line to 0 velocity
 pos=find(md.y==ymin);
-md.spcvx(pos)=0;
-md.spcvy(pos)=0;
+md.diagnostic.spcvx(pos)=0;
+md.diagnostic.spcvy(pos)=0;
 
 %icefront
 nodeonicefront=zeros(md.numberofnodes,1);
 pos=find(md.y==ymax); nodeonicefront(pos)=1;
-pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); pressureload=md.segments(pos,:);
-pressureload=[pressureload 1*md.mask.elementonfloatingice(pressureload(:,end))];
-md.pressureload=pressureload;
+pos=find(nodeonicefront(md.segments(:,1)) | nodeonicefront(md.segments(:,2))); diagnostic.icefront=md.segments(pos,:);
+diagnostic.icefront=[diagnostic.icefront 1*md.mask.elementonfloatingice(diagnostic.icefront(:,end))];
+md.diagnostic.icefront=diagnostic.icefront;
 
 %partitioning
@@ -77 +77 @@
 
 %imperative! 
-md.eps_rel=10^-10; %tighten for qmu analysese
+md.diagnostic.reltol=10^-10; %tighten for qmu analysese
 md.qmu.isdakota=1;