Index: /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F
===================================================================
--- /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F	(revision 26881)
+++ /issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F	(revision 26882)
@@ -45,8 +45,4 @@
 C     i,j, bi,bj  :: loop indices
       INTEGER bi,bj,i,j
-#ifdef ALLOW_CPL_ISSM
-      _RL  issm_density
-      issm_density = 917.
-#endif /* ALLOW_CPL_ISSM */
 
       IF ( SHELFICEMassStepping ) THEN
@@ -77,6 +73,5 @@
           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
+             shelficeMass(i,j,bi,bj) = R_shelfIce(i,j,bi,bj) * rhoConst
            ENDDO
           ENDDO
Index: /issm/trunk-jpl/test/MITgcm/coupling_rules.txt
===================================================================
--- /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 26881)
+++ /issm/trunk-jpl/test/MITgcm/coupling_rules.txt	(revision 26882)
@@ -10,8 +10,8 @@
 xC                     Nx*Ny       Real*8  10001005
 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
+OceanMelt              Nx*Ny       Real*8  10001007
+    melt_mesh[i]=-melt_mesh[i]/rho_ice; //heat flux provided by ocean is in kg/m^2/s
+IceBase                Nx*Ny       Real*8  10001008
+    in m, 0 is sea level, negative is below sea level, 9999 is no ice
 
 >>>> INITIALIZATION FROM SCRATCH
Index: /issm/trunk-jpl/test/MITgcm/tools/lookat_4003.m
===================================================================
--- /issm/trunk-jpl/test/MITgcm/tools/lookat_4003.m	(revision 26882)
+++ /issm/trunk-jpl/test/MITgcm/tools/lookat_4003.m	(revision 26882)
@@ -0,0 +1,37 @@
+pn='~/mitgcm/slr/components/issm/trunk-jpl/test/NightlyRun/';
+p1=[pn 'RunUncoupledMITgcm/'];
+p2=[pn 'RunCoupledMITgcmISSM/'];
+
+fld='Eta';
+ts=1;
+v1=rdmds([p1 fld],ts);
+v2=rdmds([p2 fld],ts);
+clf
+subplot(311), mypcolor(v1); title(fld), colorbar
+subplot(312), mypcolor(v2); title('coupled'), colorbar
+subplot(313), mypcolor(v2-v1); title('coupled-uncoupled'), colorbar
+
+ts=8;
+fld='R_shelfIce1_';
+v1=readbin([p2 fld myint2str(ts,10) '.data'],[3 200]);
+fld='R_shelfIce2_';
+v2=readbin([p2 fld myint2str(ts,10) '.data'],[3 200]);
+clf
+subplot(311), mypcolor(v1); title(fld), colorbar
+subplot(312), mypcolor(v2); title('coupled'), colorbar
+subplot(313), mypcolor(v2-v1); title('coupled-uncoupled'), colorbar
+
+fld='surfDiag';
+ts=2;
+v1=rdmds([p1 fld],ts);
+v2=rdmds([p2 fld],ts);
+clf
+fld={'ETAN','RSURF','oceTAUX','oceTAUY','oceQnet','oceFWflx', ...
+     'MXLDEPTH','SHIfwFlx','SHIhtFlx','SHIgammT','SHIgammS', ...
+     'SHI_mass','SHIuStar'};
+for i=1:length(fld)
+    subplot(311), mypcolor(v1(:,:,i)); title(fld{i}), colorbar
+    subplot(312), mypcolor(v2(:,:,i)); title('coupled'), colorbar
+    subplot(313), mypcolor(v2(:,:,i)-v1(:,:,i)); title('coupled-uncoupled'), colorbar
+    pause
+end