Changeset 5606

Timestamp:

08/27/10 08:42:22 (15 years ago)

Author:

seroussi

Message:

do not print the images in validation

Location:

issm/trunk/test/NightlyRun

Files:

• test1101.m (modified) (1 diff)
• test1102.m (modified) (3 diffs)
• test1103.m (modified) (1 diff)
• test1105.m (modified) (1 diff)
• test1107.m (modified) (1 diff)
• test1109.m (modified) (1 diff)
• test1201.m (modified) (1 diff)
• test1202.m (modified) (1 diff)
• test1203.m (modified) (1 diff)
• test1205.m (modified) (1 diff)
• test1207.m (modified) (1 diff)
• test1208.m (modified) (2 diffs)
• test1301.m (modified) (1 diff)
• test1302.m (modified) (1 diff)
• test1303.m (modified) (1 diff)
• test1304.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test1101.m

@@ -1 +1 @@
 %This test is a test from the ISMP-HOM Intercomparison project
 %Pattyn and Payne 2006
-printingflag=true;
+printingflag=false;
 
 L_list={5000,10000,20000,40000,80000,160000};

issm/trunk/test/NightlyRun/test1102.m

@@ -1 +1 @@
 %This test is a test from the ISMP-HOM Intercomparison project
 %Pattyn and Payne 2006
+printingflag=false;
 
 L_list={5000,10000,20000,40000,80000,160000};
+L_list={80000};
 results={};
+minvx=[];
+maxvx=[];
 
 for i=1:length(L_list),
@@ -12 +16 @@
         md=squaremesh(md,L,L,nx,ny);
         md=geography(md,'',''); %ice sheet test
+
+%       %Find elements at the corner and extract model
+%       posnodes=find((md.x==0 | md.x==max(md.x)) & (md.y==0 | md.y==max(md.y)));
+%       [a,b]=find(ismember(md.elements,posnodes));
+%       elements=ones(md.numberofelements,1);
+%       elements(a)=0;
+%       md=modelextract(md,elements);
+
         md=parameterize(md,'../Par/ISMIPA.par');
-        md=extrude(md,9,1);
-
-        md=setelementstype(md,'pattyn','all','stokes','all');
+        md=extrude(md,10,1);
+        md=setelementstype(md,'stokes','all');
 
         %Create dirichlet on the bed only
-        md.spcvelocity=zeros(md.numberofgrids,6);
+        %md.spcvelocity=zeros(md.numberofgrids,6);
         pos=find(md.gridonbed);
-        md.spcvelocity(pos,1:2)=1;
+        md.spcvelocity(pos,1:3)=1;
 
         %Create MPCs to have periodic boundary conditions
-        posx=find(md.x==0);
-        posx2=find(md.x==max(md.x));
+        %posx=find(md.x==0);
+        %posx2=find(md.x==max(md.x));
+        %posx=find(md.x==0 & md.y~=0 & md.y~=max(md.y) & ~md.gridonbed);
+        %posx2=find(md.x==max(md.x) &  md.y~=0 & md.y~=max(md.y) & ~md.gridonbed);
 
-        posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x)); %Don't take the same grids two times
-        posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x));
+        %posy=find(md.y==0 & md.x~=0 & md.x~=max(md.x) & ~md.gridonbed); %Don't take the same grids two times
+        %posy2=find(md.y==max(md.y) & md.x~=0 & md.x~=max(md.x) & ~md.gridonbed);
 
-        md.penalties=[posx,posx2;posy,posy2];
+        %md.penalties=[posx,posx2;posy,posy2];
 
         %Compute the diagnostic
+        md.verbose=3;
+        md.eps_abs=NaN;
+        md.eps_rel=NaN;
+        md.eps_res=1;
         md=solve(md,'analysis_type',DiagnosticSolutionEnum);
 
@@ -38 +55 @@
         vy=PatchToVec(md.results.DiagnosticSolution.Vy);
         vz=PatchToVec(md.results.DiagnosticSolution.Vz);
+        pressure=PatchToVec(md.results.DiagnosticSolution.Pressure);
         results{i}=md.results.DiagnosticSolution;
+        minvx(i)=min(vx(end-md.numberofgrids2d+1:end));
+        maxvx(i)=max(vx(end-md.numberofgrids2d+1:end));
 
-        plotmodel(md,'data',vx,'data',vy,'data',vz,'layer#all',md.numlayers)
+        %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)
+        if printingflag, 
+                set(gcf,'Color','w')
+                printmodel(['ismipastokesvx' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+                system(['mv ismipastokesvx' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+        end
+        plotmodel(md,'data',vy,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',3)
+        if printingflag, 
+                set(gcf,'Color','w')
+                printmodel(['ismipastokesvy' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+                system(['mv ismipastokesvy' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+        end
+        plotmodel(md,'data',vz,'layer#all',md.numlayers,'xlim',[0 L/10^3],'ylim',[0 L/10^3],'unit','km','figure',4)
+        if printingflag, 
+                set(gcf,'Color','w')
+                printmodel(['ismipastokesvz' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+                system(['mv ismipastokesvz' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+        end
 
+        if(L==5000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP5000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[10 18],'xlim',[0 5000],'title','','xlabel','')
+        elseif(L==10000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP10000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[10 30],'xlim',[0 10000],'title','','xlabel','')
+        elseif(L==20000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP20000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[0 50],'xlim',[0 20000],'title','','xlabel','')
+        elseif(L==40000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP40000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[0 80],'xlim',[0 40000],'title','','xlabel','')
+        elseif(L==80000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP80000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[0 100],'xlim',[0 80000],'title','','xlabel','')
+        elseif(L==160000),
+                plotmodel(md,'data',vx,'sectionvalue','../Exp/ISMIP160000.exp','layer',md.numlayers,...
+                        'resolution',[10 10],'ylim',[0 120],'xlim',[0 160000],'title','','xlabel','')
+        end
+        if printingflag, 
+                set(gcf,'Color','w')
+                printmodel(['ismipastokesvxsec' num2str(L)],'png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+                system(['mv ismipastokesvxsec' num2str(L) '.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+        end
 end
 
+%Now plot the min and max values of vx for each size of the square
+plot([5 10 20 40 80 160],minvx);ylim([0 18])
+if printingflag, 
+        set(gcf,'Color','w')
+        printmodel('ismipastokesminvx','png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+        system(['mv ismipastokesminvx.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+end
+plot([5 10 20 40 80 160],maxvx);ylim([0 120])
+if printingflag, 
+        set(gcf,'Color','w')
+        printmodel('ismipastokesmaxvx','png','margin','on','marginsize',25,'frame','off','resolution',1.5,'hardcopy','off');
+        system(['mv ismipastokesmaxvx.png ' ISSM_DIR '/../website/doc_pdf/validation/Images/ISMIP/TestA ']);
+end
 %Fields and tolerances to track changes
 field_names     ={...

issm/trunk/test/NightlyRun/test1103.m

@@ -1 +1 @@
 %This test is a test from the ISMP-HOM Intercomparison project
 %Pattyn and Payne 2006
-printingflag=true;
+printingflag=false;
 
 L_list={5000,10000,20000,40000,80000,160000};

issm/trunk/test/NightlyRun/test1105.m

@@ -1 +1 @@
 %This test is a test from the ISMP-HOM Intercomparison project
 %Pattyn and Payne 2006
-printingflag=true;
+printingflag=false;
 
 L_list={5000,10000,20000,40000,80000,160000};

issm/trunk/test/NightlyRun/test1107.m

@@ -1 +1 @@
 %This test is a test from the ISMP-HOM Intercomparison project
 %Pattyn and Payne 2006
-printingflag=true;
+printingflag=false;
 
 L_list={5000,10000,20000,40000,80000,160000};

issm/trunk/test/NightlyRun/test1109.m

@@ -5 +5 @@
 %                   - Pattyn with some sliding
 %                   - Stokes with some sliding
-printingflag=true;
+printingflag=false;
 results={};
 

issm/trunk/test/NightlyRun/test1201.m

@@ -1 +1 @@
 %This test is a test from the EISMINT for Ice shelves Vincent Rommelaere 1996
-printingflag=true;
+printingflag=false;
 
 results={};

issm/trunk/test/NightlyRun/test1202.m

@@ -1 +1 @@
 %Test on the diagnostic model and the prognostic in 2d
-printingflag=true;
+printingflag=false;
 
 %tests 3 and 4: using Glen's flow law

issm/trunk/test/NightlyRun/test1203.m

@@ -1 +1 @@
 %Test on the diagnostic model and the prognostic in 2d
-printingflag=true;
+printingflag=false;
 
 %test 5 and 6 : 

issm/trunk/test/NightlyRun/test1205.m

@@ -44 +44 @@
 %Plot of the velocity from the exact and calculated solutions
 figure(1)
+set(gcf,'Position',[1 1 1580 1150])
 subplot(2,2,1)
 p=patch('Faces',md.elements2d,'Vertices',[md.x2d md.y2d],'FaceVertexCData',...

issm/trunk/test/NightlyRun/test1207.m

@@ -17 +17 @@
 %We extrude the model to have a 3d model
 md=extrude(md,numlayers,1);
-md=setelementstype(md,'pattyn','all','stokes','all');
+md=setelementstype(md,'stokes','all');
 
 %Spc the grids on the bed

issm/trunk/test/NightlyRun/test1208.m

@@ -1 +1 @@
 %EISMINT benchmark experiment A
 numlayers=10;
-resolution=30000;
+resolution=40000;
 
 %To begin with the numerical model
@@ -18 +18 @@
 
 %Adapt the time steps to the resolution
-md.dt=25;
+md.dt=20;
+md.verbose=10;
 md.output_frequency=20;
 

issm/trunk/test/NightlyRun/test1301.m

@@ -3 +3 @@
 % at upper surface and an impose flux at its base. The result must be a linear temperature from the upper to the lower
 % surface with an imposed slope (Geothermal flux). if it is not the case, something is thermal modeling has been changed...
-printingflag=true;
+printingflag=false;
 
 md=model;

issm/trunk/test/NightlyRun/test1302.m

@@ -2 +2 @@
 % There is u=v=0 and w=cst everywhere the only thermal boundary conditions are an imposed temperature
 % at upper surface and an impose flux at its base.
-printingflag=true;
+printingflag=false;
 
 md=model;

issm/trunk/test/NightlyRun/test1303.m

@@ -3 +3 @@
 % at the lower and upper surface. The result must be a linear temperature from the upper to the lower
 % surface. if it is not the case, something is thermal modeling has been changed...
-printingflag=true;
+printingflag=false;
 
 md=model;

issm/trunk/test/NightlyRun/test1304.m

@@ -3 +3 @@
 % at upper surface and an impose flux at its base. The result must be a linear temperature from the upper to the lower
 % surface with an imposed slope (Geothermal flux). if it is not the case, something is thermal modeling has been changed...
-printingflag=true;
+printingflag=false;
 
 md=model;