Index: /issm/trunk-jpl/test/NightlyRun/test4003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 23885)
+++ /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 23886)
@@ -269,12 +269,12 @@
 
 	%friction: 
-	md.friction.coefficient=30*ones(md.mesh.numberofvertices,1);
-	pos=find(md.mask.groundedice_levelset<=0);
-	md.friction.coefficient(pos)=0;
+	md.friction.coefficient=0*ones(md.mesh.numberofvertices,1);
+	pos=find(md.mask.groundedice_levelset>0);
+	md.friction.coefficient(pos)=5;
 	md.friction.p=ones(md.mesh.numberofelements,1);
 	md.friction.q=ones(md.mesh.numberofelements,1);
 
 	%temperatures and surface mass balance:
-	md.initialization.temperature=(273.15-20)*ones(md.mesh.numberofvertices,1);
+	md.initialization.temperature=(273.15-22)*ones(md.mesh.numberofvertices,1);
 	md.initialization.pressure=md.materials.rho_ice*md.constants.g*(md.geometry.surface-md.geometry.base);
 	md.smb.mass_balance = [1*ones(md.mesh.numberofvertices,1); 1];
@@ -302,6 +302,8 @@
 
 	%get some flux at the ice divide: 
-	pos=find(md.mesh.lat==min(md.mesh.lat));
-	md.stressbalance.spcvy(pos)=200;
+	pos=find(md.mesh.y==min(md.mesh.y));
+	md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
+	md.stressbalance.vx(pos)=0;
+	md.stressbalance.vy(pos)=1500;
 
 	%deal with boundaries, excluding icefront: 
@@ -309,8 +311,5 @@
 	vertex_on_boundary(md.mesh.segments(:,1:2))=1;
 	pos=find(vertex_on_boundary & md.mask.groundedice_levelset<=0);
-	md.stressbalance.spcvx(pos)=md.initialization.vx(pos);
-	md.stressbalance.spcvy(pos)=md.initialization.vy(pos);
-	md.stressbalance.spcvz(pos)=md.initialization.vz(pos);
-	md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
+	md.stressbalance.spcvx(pos)=0;
 
 	md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
@@ -334,8 +333,8 @@
 	md.timestepping.final_time=final_time;
 	md.timestepping.time_step=time_step;
-	md.transient.isgroundingline=1;
+	md.transient.isgroundingline=0;
 	md.transient.isthermal=0;
 	md.groundingline.migration='SubelementMigration';
-	md.groundingline.melt_interpolation='SubelementMelt2';
+	md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
 	md.groundingline.friction_interpolation='SubelementFriction2';
 
@@ -449,5 +448,5 @@
 	md.transient.isthermal=0;
 	md.groundingline.migration='SubelementMigration';
-	md.groundingline.melt_interpolation='SubelementMelt2';
+	md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
 	md.groundingline.friction_interpolation='SubelementFriction2';