Index: /issm/trunk/test/Verification/test26_Balanced2d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test26_Balanced2d/DomainOutline.exp	(revision 3531)
+++ /issm/trunk/test/Verification/test26_Balanced2d/DomainOutline.exp	(revision 3531)
@@ -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/Verification/test26_Balanced2d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test26_Balanced2d/Front.exp	(revision 3531)
+++ /issm/trunk/test/Verification/test26_Balanced2d/Front.exp	(revision 3531)
@@ -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/Verification/test26_Balanced2d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test26_Balanced2d/Square.par	(revision 3531)
+++ /issm/trunk/test/Verification/test26_Balanced2d/Square.par	(revision 3531)
@@ -0,0 +1,42 @@
+%Start defining model parameters here
+
+%dynamics
+md.verbose=10;
+md.dt=1; %1 year
+md.ndt=md.dt*10; 
+md.artificial_diffusivity=1;
+
+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;
+
+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;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceShelfBC(md,'Front.exp');
+
+%Spc the thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test26_Balanced2d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test26_Balanced2d/configuration.m	(revision 3531)
+++ /issm/trunk/test/Verification/test26_Balanced2d/configuration.m	(revision 3531)
@@ -0,0 +1,27 @@
+%name
+testname='test01_IceShelfIceFrontM2d';
+
+%meshing
+resolution=50000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=0;
+extrusionexponent=NaN;
+
+%sequences         analysis       sub_analysis  qmu   control  control_fit parallel
+sequences={  {'balancedthickness',  'none',      0 ,     0,      'nan'        0    };...
+				 {'balancedthickness',  'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test26_Balanced2d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test26_Balanced2d/testpresolve.m	(revision 3531)
+++ /issm/trunk/test/Verification/test26_Balanced2d/testpresolve.m	(revision 3531)
@@ -0,0 +1,3 @@
+if md.numberofgrids==340
+	load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
+end