Changeset 4868

Timestamp:

07/29/10 08:48:18 (15 years ago)

Author:

seroussi

Message:

Hydrostatic ice sheet test

Location:

issm/trunk/test/Validation/HydrostaticIceSheet

Files:

• NoFront/runme.m (modified) (1 diff)
• WithFront/runme.m (modified) (1 diff)

issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m

@@ -11 +11 @@
 md=setelementstype(md,'Macayeal','all');
 
-%Compute solution for a 2d model
-md=solve(md,'analysis_type','diagnostic');
-vel_shelf=md.results.diagnostic.vel;
+%Compute solution for a 2d model of ice shelf
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+vel_shelf=zeros(md.numberofgrids,1);
+vel_shelf(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 
-%Compute solution for a 3d model
+%Compute solution for a 2d model of ice sheet
 md=geography(md,'','');
-md=solve(md,'analysis_type','diagnostic');
-vel_sheet=md.results.diagnostic.vel;
+md=parameterize(md,'Square.par');
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+vel_sheet=zeros(md.numberofgrids,1);
+vel_sheet(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 
 %Plot of the velocity from the ice sheet and he ice shelf model

issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m

@@ -11 +11 @@
 md=setelementstype(md,'Macayeal','all');
 
-%Compute solution for a 2d model
-md=solve(md,'analysis_type','diagnostic');
-vel_shelf=md.results.diagnostic.vel;
+%Compute solution for a 2d model of ice shelf
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+vel_shelf=zeros(md.numberofgrids,1);
+vel_shelf(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 
-%Compute solution for a 3d model
+%Compute solution for a 2d model of ice sheet
 md=geography(md,'Shelf.exp','');
-md=solve(md,'analysis_type','diagnostic');
-vel_sheet=md.results.diagnostic.vel;
+md=parameterize(md,'Square.par');
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
+vel_sheet=zeros(md.numberofgrids,1);
+vel_sheet(md.results.DiagnosticSolution.Vel.index)=md.results.DiagnosticSolution.Vel.value;
 
 %Plot of the velocity from the ice sheet and he ice shelf model