Changeset 26882

Timestamp:

02/10/22 15:50:36 (3 years ago)

Author:

dmenemen

Message:

M MITgcm/code_4003/shelfice_step_icemass.F
M MITgcm/coupling_rules.txt
A MITgcm/tools/lookat_4003.m

Location:

issm/trunk-jpl/test/MITgcm

Files:

• code_4003/shelfice_step_icemass.F (modified) (2 diffs)
• coupling_rules.txt (modified) (1 diff)
• tools/lookat_4003.m (added)

issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F

@@ -45 +45 @@
 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 +73 @@
           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

issm/trunk-jpl/test/MITgcm/coupling_rules.txt

@@ -10 +10 @@
 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