Changeset 26452 for issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F

Timestamp:

09/21/21 15:59:26 (3 years ago)

Author:

dmenemen

Message:

first attempt to modify test4003 for remeshing
has not been tested and debugged yet

File:

• issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F (modified) (2 diffs)

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
+      _RL  issm_density
+      issm_density = 917.
 
       IF ( SHELFICEMassStepping ) THEN
@@ -65 +67 @@
         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