Changeset 15771 for issm/trunk-jpl/test/NightlyRun/test1102.m

Timestamp:

08/09/13 10:59:47 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: Diagnostic is now Stressbalance

File:

• issm/trunk-jpl/test/NightlyRun/test1102.m (modified) (2 diffs)

issm/trunk-jpl/test/NightlyRun/test1102.m

@@ -29 +29 @@
         %Create dirichlet on the bed only
         pos=find(md.mesh.vertexonbed);
-        md.diagnostic.spcvx(pos)=0.;
-        md.diagnostic.spcvy(pos)=0.;
-        md.diagnostic.spcvz(pos)=0.;
+        md.stressbalance.spcvx(pos)=0.;
+        md.stressbalance.spcvy(pos)=0.;
+        md.stressbalance.spcvz(pos)=0.;
 
 %       %Create MPCs to have periodic boundary conditions
@@ -42 +42 @@
 %       posy2=find(md.mesh.y==max(md.mesh.y) & md.mesh.x~=0. & md.mesh.x~=max(md.mesh.x) & ~md.mesh.vertexonbed);
 
-%       md.diagnostic.vertex_pairing=[posx,posx2;posy,posy2];
+%       md.stressbalance.vertex_pairing=[posx,posx2;posy,posy2];
 
-        %Compute the diagnostic
-        md.diagnostic.abstol=NaN;
-        md.diagnostic.reltol=NaN;
-        md.diagnostic.restol=1.;
+        %Compute the stressbalance
+        md.stressbalance.abstol=NaN;
+        md.stressbalance.reltol=NaN;
+        md.stressbalance.restol=1.;
         md.cluster=generic('name',oshostname(),'np',8);
-        md=solve(md,DiagnosticSolutionEnum());
+        md=solve(md,StressbalanceSolutionEnum());
 
         %Plot the results and save them
-        vx=(md.results.DiagnosticSolution.Vx);
-        vy=(md.results.DiagnosticSolution.Vy);
-        vz=(md.results.DiagnosticSolution.Vz);
-        pressure=(md.results.DiagnosticSolution.Pressure);
-        results{i}=md.results.DiagnosticSolution;
+        vx=(md.results.StressbalanceSolution.Vx);
+        vy=(md.results.StressbalanceSolution.Vy);
+        vz=(md.results.StressbalanceSolution.Vz);
+        pressure=(md.results.StressbalanceSolution.Pressure);
+        results{i}=md.results.StressbalanceSolution;
         minvx(i)=min(vx(end-md.mesh.numberofvertices2d+1:end));
         maxvx(i)=max(vx(end-md.mesh.numberofvertices2d+1:end));