Index: /issm/trunk-jpl/test/MITgcm/code_4003/SHELFICE_OPTIONS.h
===================================================================
--- /issm/trunk-jpl/test/MITgcm/code_4003/SHELFICE_OPTIONS.h	(revision 26451)
+++ /issm/trunk-jpl/test/MITgcm/code_4003/SHELFICE_OPTIONS.h	(revision 26452)
@@ -16,16 +16,15 @@
 
 C     allow code for simple ISOMIP thermodynamics
-!#define ALLOW_ISOMIP_TD
+#undef ALLOW_ISOMIP_TD
 
 C     allow friction velocity-dependent transfer coefficient
 C     following Holland and Jenkins, JPO, 1999
 #define SHI_ALLOW_GAMMAFRICT
-C     in uStar expression, use wet-point method to average velocity
-C     at grid-cell center
-!#define SHI_USTAR_WETPOINT
 
+C     allow (vertical) remeshing whenever ocean top thickness factor
+C     exceeds thresholds
 #define ALLOW_SHELFICE_REMESHING
-#define SHI_USTAR_TOPDR
-#undef ALLOW_SHELFICE_GROUNDED_ICE
+C     and allow to print message to STDOUT when this happens
+#define SHELFICE_REMESH_PRINT
 
 #endif /* ALLOW_SHELFICE */
Index: /issm/trunk-jpl/test/MITgcm/code_4003/packages.conf
===================================================================
--- /issm/trunk-jpl/test/MITgcm/code_4003/packages.conf	(revision 26451)
+++ /issm/trunk-jpl/test/MITgcm/code_4003/packages.conf	(revision 26452)
@@ -1,5 +1,4 @@
 gfd
 shelfice
-streamice
 obcs
 diagnostics
Index: /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F
===================================================================
--- /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F	(revision 26451)
+++ /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F	(revision 26452)
@@ -45,4 +45,6 @@
 C     i,j, bi,bj  :: loop indices
       INTEGER bi,bj,i,j
+      _RL  issm_density
+      issm_density = 917.
 
       IF ( SHELFICEMassStepping ) THEN
@@ -65,4 +67,19 @@
         ENDDO
 #endif /* ALLOW_STREAMICE */
+
+       ELSEIF (useCoupler)
+
+#ifdef ALLOW_CPL_ISSM
+        DO bj = myByLo(myThid), myByHi(myThid)
+         DO bi = myBxLo(myThid), myBxHi(myThid)
+          DO j=1-OLy,sNy+OLy-1
+           DO i=1-OLx+1,sNx+OLx-1
+             shelficeMass(i,j,bi,bj) =
+     &        R_shelfIce(i,j,bi,bj) * issm_density
+           ENDDO
+          ENDDO
+         ENDDO
+        ENDDO
+#endif /* ALLOW_CPL_ISSM */
 
        ELSE
Index: /issm/trunk-jpl/test/MITgcm/coupling_rules.txt
===================================================================
--- /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 26451)
+++ /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 26452)
@@ -11,4 +11,6 @@
 yC                     Nx*Ny       Real*8  10001006
 IceBase                Nx*Ny       Real*8  10001007
+    (in m, 0 is sea level, negative is below sea level, 9999 is no ice,
+     use density 917 kg/m^3)
 OceanMelt              Nx*Ny       Real*8  10001008
 
Index: /issm/trunk-jpl/test/MITgcm/input_4003/data.shelfice
===================================================================
--- /issm/trunk-jpl/test/MITgcm/input_4003/data.shelfice	(revision 26452)
+++ /issm/trunk-jpl/test/MITgcm/input_4003/data.shelfice	(revision 26452)
@@ -0,0 +1,20 @@
+# ===================================
+# | Parameters for SHELFICE package |
+# ===================================
+ &SHELFICE_PARM01
+ SHELFICEconserve = .TRUE.,
+ SHELFICEuseGammaFrict = .TRUE.,
+ SHELFICEDragQuadratic = 0.0015,
+ shiCdrag = 0.0015,
+ SHELFICEselectDragQuadr= 2,
+ SHELFICEMassStepping = .TRUE.,
+#--
+ SHELFICEremeshFrequency = 600.0,
+#- need to satisfy: splitThrs > 1 + mergeThrs / Sdz
+#       with: Sdz = min{ delR(k+1)/delR(k) }_[k=1:Nr-1]
+ SHELFICEsplitThreshold = 1.12,
+ SHELFICEmergeThreshold = 0.10,
+#--
+ SHELFICEtopoFile='icetopo.bin',
+ SHELFICEwriteState = .TRUE.,
+ &
Index: /issm/trunk-jpl/test/NightlyRun/test4003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 26451)
+++ /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 26452)
@@ -4,4 +4,5 @@
 %Script control parameters
 steps=1:11;
+steps=1:5;
 final_time=1/365;