Changeset 6330

Timestamp:

10/18/10 09:00:01 (14 years ago)

Author:

seroussi

Message:

test1207 tolerance

Location:

issm/trunk/test/NightlyRun

Files:

• test1110.m (modified) (1 diff)
• test1207.m (modified) (1 diff)
• test1208.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test1110.m

@@ -5 +5 @@
 ny=30;
 md=model;
-%md=squaremesh(md,L,L,nx,ny);
-md=mesh(md,'../Exp/SquareISMIP.exp',5500);
+md=squaremesh(md,L,L,nx,ny);
+%md=mesh(md,'../Exp/SquareISMIP.exp',5500);
 md=geography(md,'',''); %ice sheet test
 md=parameterize(md,'../Par/ISMIPF.par');
-md=extrude(md,5,1);
+md=extrude(md,2,1);
 
 md=setelementstype(md,'pattyn','all');
 %md=setelementstype(md,'pattyn','all','stokes','all'); 
 
+md.spcvelocity(:,:)=0;
 %Create dirichlet on the bed if no slip
 pos=find(md.gridonbed);
 md.spcvelocity(pos,1:2)=1;
-md.spcvelocity(pos,4:5)=0;
+md.spctemperature(:,2)=255;
+md.spctemperature(:,1)=255;
 
-%%Create MPCs to have periodic boundary conditions
-%posx=find(md.x==0);
-%posx2=find(md.x==max(md.x));
-%
-%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));
-%
-%md.penalties=[posx,posx2;posy,posy2];
+%Create MPCs to have periodic boundary conditions
+posx=find(md.x==0);
+posx2=find(md.x==max(md.x));
+
+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));
+
+md.penalties=[posx,posx2;posy,posy2];
+
+md.dt=5;
+md.ndt=2000;
 
 %Compute the diagnostic
-md=solve(md,'analysis_type',Transient3DSolutionEnum);
+md.cluster=wilkes;
+md.verbose=10;
+md.rheology_B=md.rheology_B/2;
+%md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+%md=solve(md,'analysis_type',ThermalSolutionEnum);
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
 
-vx=PatchToVec(md.results.DiagnosticSolution.Vx);
-vy=PatchToVec(md.results.DiagnosticSolution.Vy);
-vz=PatchToVec(md.results.DiagnosticSolution.Vz);
+%vx=PatchToVec(md.results.DiagnosticSolution.Vx);
+%vy=PatchToVec(md.results.DiagnosticSolution.Vy);
+%vz=PatchToVec(md.results.DiagnosticSolution.Vz);
 
-plotmodel(md,'data',vx,'data',vy,'data',vz,'layer#all',md.numlayers,'contourlevels',{9.2,18.4,27.6,46,73.6,82.8},'contourlevels',{0,-1.2,-0.9,-0.6,-0.3,1.2},'contourlevels',{-5.2,-3.9,-1.3,0,1.3,3.9,5.2})
+plotmodel(md,'data',vx,'data',vy,'data',vz,'layer#all',md.numlayers)

issm/trunk/test/NightlyRun/test1207.m

@@ -89 +89 @@
 };
 field_tolerances={...
-        1e-13,1e-13,1e-13, ...
+        1e-12,1e-12,1e-12, ...
 };
 field_values={

issm/trunk/test/NightlyRun/test1208.m

@@ -19 +19 @@
 md.dt=20;
 md.output_frequency=20;
-md.ndt=200000;
-md.artificial_diffusivity=2;
+md.ndt=50000;
+md.artificial_diffusivity=1;
 
 %Now we can solve the problem