Changeset 9725 for issm/trunk/test/NightlyRun/test1107.m

Timestamp:

09/09/11 08:35:15 (14 years ago)

Author:

Mathieu Morlighem

Message:

Some more objects in mesh

File:

• issm/trunk/test/NightlyRun/test1107.m (modified) (4 diffs)

issm/trunk/test/NightlyRun/test1107.m

@@ -21 +21 @@
 
         %We need one grd on dirichlet: the 4 corners are set to zero
-        md.diagnostic.spcvx=NaN*ones(md.numberofnodes,1);
-        md.diagnostic.spcvy=NaN*ones(md.numberofnodes,1);
-        md.diagnostic.spcvz=NaN*ones(md.numberofnodes,1);
+        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);
 
         %Create MPCs to have periodic boundary conditions
@@ -65 +65 @@
         vz=PatchToVec(md.results.DiagnosticSolution.Vz);
         results{i}=md.results.DiagnosticSolution;
-        minvx(i)=min(vx(end-md.numberofnodes2d+1:end));
-        maxvx(i)=max(vx(end-md.numberofnodes2d+1:end));
+        minvx(i)=min(vx(end-md.mesh.numberofvertices2d+1:end));
+        maxvx(i)=max(vx(end-md.mesh.numberofvertices2d+1:end));
 
         %Now plot vx, vy, vz and vx on a cross section
-        plotmodel(md,'data',vx,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',2)
+        plotmodel(md,'data',vx,'layer#all',md.mesh.numberoflayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',2)
         if printingflag, 
                 set(gcf,'Color','w')
@@ -75 +75 @@
                 system(['mv ismipdpattynvx' num2str(L) '.png ' ISSM_DIR '/website/doc_pdf/validation/Images/ISMIP/TestD ']);
         end
-        plotmodel(md,'data',vz,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',3)
+        plotmodel(md,'data',vz,'layer#all',md.mesh.numberoflayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',3)
         if printingflag, 
                 set(gcf,'Color','w')
@@ -83 +83 @@
 
         if(L==5000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP5000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP5000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[0 20],'xlim',[0 5000],'title','','xlabel','','figure',4)
         elseif(L==10000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP10000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP10000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[0 20],'xlim',[0 10000],'title','','xlabel','','figure',4)
         elseif(L==20000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP20000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP20000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[0 30],'xlim',[0 20000],'title','','xlabel','','figure',4)
         elseif(L==40000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP40000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP40000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[10 60],'xlim',[0 40000],'title','','xlabel','','figure',4)
         elseif(L==80000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP80000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP80000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[0 200],'xlim',[0 80000],'title','','xlabel','','figure',4)
         elseif(L==160000),
-                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP160000.exp','layer',md.numlayers,...
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP160000.exp','layer',md.mesh.numberoflayers,...
                         'resolution',[10 10],'ylim',[0 400],'xlim',[0 160000],'title','','xlabel','','figure',4)
         end