Index: /issm/trunk-jpl/test/NightlyRun/test4003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 26738)
+++ /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 26739)
@@ -4,10 +4,10 @@
 %Script control parameters
 steps=1:11;
-steps=[1:3 6:11];
+%steps=[1:3 6:11];
 
 %To download and recompile MITgcm from scratch:
-%!rm -rf ${ISSM_DIR}/test/MITgcm/install
-%!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
-%!rm -rf Models
+!rm -rf ${ISSM_DIR}/test/MITgcm/install
+!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
+!rm -rf Models
 
 %Organizer
@@ -273,8 +273,8 @@
 
     %level sets:
-    md.mask.ocean_levelset=iceshelf_mask(:);
-    md.mask.ice_levelset=-ones(md.mesh.numberofvertices,1);
-
-    pos=find(~ice_mask(:) & thickness(:)==0); md.mask.ice_levelset(pos)=1;
+    md.mask.ice_levelset=iceshelf_mask(:);
+    pos=find(md.mask.ice_levelset==1); md.mask.ice_levelset(pos)=-1;
+    pos=find(md.mask.ice_levelset==0); md.mask.ice_levelset(pos)=1;
+    md.mask.ocean_levelset=-ones(md.mesh.numberofvertices,1);
 
     savemodel(org,md);
@@ -324,10 +324,10 @@
 	pos=find(md.mesh.lat==min(md.mesh.lat));
 	md.masstransport.spcthickness(pos)=md.geometry.thickness(pos);
-	md.stressbalance.spcvx(pos)=-800;
-	md.stressbalance.spcvy(pos)=0;
+	md.stressbalance.spcvy(pos)=800;
+	md.stressbalance.spcvx(pos)=0;
 
 	%deal with boundaries, excluding icefront:
 	pos=find(md.mesh.long==min(md.mesh.long) | md.mesh.long==max(md.mesh.long));
-	md.stressbalance.spcvy(pos)=0;
+	md.stressbalance.spcvx(pos)=0;
 
 	point1=find(md.mesh.y==min(md.mesh.y)); point2=find(md.mesh.x==max(md.mesh.x));
@@ -362,5 +362,5 @@
 	md.groundingline.melt_interpolation='NoMeltOnPartiallyFloating';
 	md.groundingline.friction_interpolation='SubelementFriction2';
-	md.masstransport.stabilization=0;
+	md.masstransport.stabilization=1;
 
 	md.cluster=generic('name',oshostname(),'np',2);