Changeset 16454

Timestamp:

10/18/13 08:17:30 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixed tutorials with new mask and solution names

Location:

issm/trunk/examples

Files:

• Jakobshavn/Jks.par (modified) (2 diffs)
• Jakobshavn/runme.m (modified) (3 diffs)
• SquareIceShelf/Square.par (modified) (1 diff)
• SquareIceShelf/runme.m (modified) (1 diff)

issm/trunk/examples/Jakobshavn/Jks.par

@@ -44 +44 @@
 disp('   Construct basal friction parameters');
 md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
-pos=find(md.mask.vertexonfloatingice);
-md.friction.coefficient(pos)=0; %no friction applied on floating ice
+md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);
@@ -58 +57 @@
 disp('   Set other boundary conditions');
 md=SetMarineIceSheetBC(md,'./Front.exp');
-md.diagnostic.icefront(:,end)=1; %Apply water pressure at the ice front
 md.basalforcings.melting_rate=zeros(md.mesh.numberofvertices,1);

issm/trunk/examples/Jakobshavn/runme.m

@@ -42 +42 @@
         md=loadmodel('JksPar.mdl');
 
-        md=setflowequation(md,'macayeal','all');
+        md=setflowequation(md,'SSA','all');
 
         %Control general
@@ -64 +64 @@
 
         %Additional parameters
-        md.diagnostic.restol=0.01;
-        md.diagnostic.reltol=0.1;
-        md.diagnostic.abstol=NaN;
+        md.stressbalance.restol=0.01;
+        md.stressbalance.reltol=0.1;
+        md.stressbalance.abstol=NaN;
 
         %Go solve
         md.cluster=generic('name',oshostname,'np',4);
-        md.solver=addoptions(md.solver,NoneAnalysisEnum,asmoptions);
-        md.solver=addoptions(md.solver,DiagnosticVertAnalysisEnum,jacobiasmoptions);
-        md=solve(md,DiagnosticSolutionEnum);
+        md=solve(md,StressbalanceSolutionEnum);
         
         save JksControl.mdl md
@@ -83 +81 @@
                 'FontSize#all',12,...
                 'data',md.inversion.vel_obs,'title','Observed velocity',...
-                'data',md.results.DiagnosticSolution.Vel,'title','Modeled Velocity',...
+                'data',md.results.StressbalanceSolution.Vel,'title','Modeled Velocity',...
                 'colorbar#1','off','colorbar#2','on','colorbartitle#2','[m/yr]',...
                 'caxis#1-2',[0,7000],...
                 'data',md.geometry.bed,'title','Bed elevation',...
-                'data',md.results.DiagnosticSolution.FrictionCoefficient,...
+                'data',md.results.StressbalanceSolution.FrictionCoefficient,...
                 'title','Friction Coefficient',...
                 'colorbar#3','on','colorbartitle#3','[m]','colorbar#4','on');

issm/trunk/examples/SquareIceShelf/Square.par

@@ -11 +11 @@
 
 disp('      creating drag');
-pos=find(md.mask.vertexonfloatingice);
 md.friction.coefficient=200*ones(md.mesh.numberofvertices,1);
-md.friction.coefficient(pos)=0;
+md.friction.coefficient(find(md.mask.groundedice_levelset<0.))=0.;
 md.friction.p=ones(md.mesh.numberofelements,1);
 md.friction.q=ones(md.mesh.numberofelements,1);

issm/trunk/examples/SquareIceShelf/runme.m

@@ -3 +3 @@
 md=setmask(md,'all','');
 md=parameterize(md,'Square.par');
-md=setflowequation(md,'macayeal','all');
+md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname,'np',2);
-md=solve(md,DiagnosticSolutionEnum);
+md=solve(md,StressbalanceSolutionEnum);