Index: /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m	(revision 4867)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m	(revision 4868)
@@ -11,12 +11,15 @@
 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
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m	(revision 4867)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m	(revision 4868)
@@ -11,12 +11,15 @@
 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