Changeset 15771 for issm/trunk-jpl/test/NightlyRun/test1101.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/test1101.m (modified) (2 diffs)

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

@@ -21 +21 @@
 
         %Create dirichlet on the bed only
-        md.diagnostic.spcvx=NaN*ones(md.mesh.numberofvertices,1);
-        md.diagnostic.spcvy=NaN*ones(md.mesh.numberofvertices,1);
-        md.diagnostic.spcvz=NaN*ones(md.mesh.numberofvertices,1);
+        md.stressbalance.spcvx=NaN*ones(md.mesh.numberofvertices,1);
+        md.stressbalance.spcvy=NaN*ones(md.mesh.numberofvertices,1);
+        md.stressbalance.spcvz=NaN*ones(md.mesh.numberofvertices,1);
 
         pos=find(md.mesh.vertexonbed);
-        md.diagnostic.spcvx(pos)=0.;
-        md.diagnostic.spcvy(pos)=0.;
+        md.stressbalance.spcvx(pos)=0.;
+        md.stressbalance.spcvy(pos)=0.;
 
         %Create MPCs to have periodic boundary conditions
@@ -36 +36 @@
         posy2=find(md.mesh.y==max(md.mesh.y) & md.mesh.x~=0. & md.mesh.x~=max(md.mesh.x));
 
-        md.diagnostic.vertex_pairing=[posx,posx2;posy,posy2];
+        md.stressbalance.vertex_pairing=[posx,posx2;posy,posy2];
 
-        %Compute the diagnostic
+        %Compute the stressbalance
         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);
-        results{i}=md.results.DiagnosticSolution;
+        vx=(md.results.StressbalanceSolution.Vx);
+        vy=(md.results.StressbalanceSolution.Vy);
+        vz=(md.results.StressbalanceSolution.Vz);
+        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));