Index: /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/DomainOutline.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/DomainOutline.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Pattyn.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Pattyn.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Pattyn.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-10 -10
+1000010 -10
+1000010 1000010
+-10 1000010
+-10 -10
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Square.par
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Square.par	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/Square.par	(revision 761)
@@ -0,0 +1,35 @@
+%Start defining model parameters here
+
+disp('      creating thickness');
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+disp('      creating drag');
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+disp('      creating temperature');
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+disp('      creating flow law paramter');
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceSheetBC(md);
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=1;
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/NoFront/runme.m	(revision 761)
@@ -0,0 +1,23 @@
+% This file can be run to compare an icesheet in hydrostatic equilibrium with an iceshelf.
+% This test deals with a  mesh includind a iceshelf/icesheet with no ice front. The geometry 
+% is square. Just run this file in Matlab, with a properly setup Ice code. 
+% The result should be exactly the same since the ice sheet is in hydrostatic equilibrium
+% so the friction at the base should be zero.
+
+md=model;
+md=mesh(md,'DomainOutline.exp',50000);
+md=geography(md,'all','');
+md=parameterize(md,'Square.par');
+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 3d model
+md=geography(md,'','');
+md=solve(md,'analysis_type','diagnostic');
+vel_sheet=md.results.diagnostic.vel;
+
+%Plot of the velocity from the ice sheet and he ice shelf model
+plotmodel(md,'data',vel_shelf,'title','Velocity of ice shelf','data',vel_sheet,'title','Velocity of ice sheet','data',vel_sheet-vel_shelf,'title','Absolute difference between ice sheet and ice shelf velocities','data',abs((vel_sheet-vel_shelf)./vel_shelf)*100,'title','Relative difference between ice sheet and ice shelf velocities [in %]')
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/DomainOutline.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/DomainOutline.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Front.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Front.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Front.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Pattyn.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Pattyn.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Pattyn.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-10 -10
+1000010 -10
+1000010 1000010
+-10 1000010
+-10 -10
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Shelf.exp
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Shelf.exp	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Shelf.exp	(revision 761)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 800000
+-1000 1100000
+1100000 1100000
+1100000 800000
+-1000 800000
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Square.par
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Square.par	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/Square.par	(revision 761)
@@ -0,0 +1,35 @@
+%Start defining model parameters here
+
+disp('      creating thickness');
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+disp('      creating drag');
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+disp('      creating temperature');
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+disp('      creating flow law paramter');
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=1;
Index: /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m
===================================================================
--- /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m	(revision 761)
+++ /issm/trunk/test/Validation/HydrostaticIceSheet/WithFront/runme.m	(revision 761)
@@ -0,0 +1,23 @@
+% This file can be run to compare an icesheet in hydrostatic equilibrium with an iceshelf.
+% This test deals with a  mesh includind a iceshelf/icesheet therminated with a floating iceshelf. The geometry 
+% is square. Just run this file in Matlab, with a properly setup Ice code. 
+% The result should be exactly the same since the ice sheet is in hydrostatic equilibrium
+% so the friction at the base should be zero.
+
+md=model;
+md=mesh(md,'DomainOutline.exp',50000);
+md=geography(md,'all','');
+md=parameterize(md,'Square.par');
+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 3d model
+md=geography(md,'Shelf.exp','');
+md=solve(md,'analysis_type','diagnostic');
+vel_sheet=md.results.diagnostic.vel;
+
+%Plot of the velocity from the ice sheet and he ice shelf model
+plotmodel(md,'data',vel_shelf,'title','Velocity of ice shelf','data',vel_sheet,'title','Velocity of ice sheet','data',vel_sheet-vel_shelf,'title','Absolute difference between ice sheet and ice shelf velocities','data',abs((vel_sheet-vel_shelf)./vel_shelf)*100,'title','Relative difference between ice sheet and ice shelf velocities [in %]')