Changeset 26452

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

Location:

issm/trunk-jpl/test

Files:

• MITgcm/code_4003/SHELFICE_OPTIONS.h (modified) (1 diff)
• MITgcm/code_4003/packages.conf (modified) (1 diff)
• MITgcm/code_4003/shelfice_step_icemass.F (modified) (2 diffs)
• MITgcm/coupling_rules.txt (modified) (1 diff)
• MITgcm/input_4003/data.shelfice (added)
• NightlyRun/test4003.m (modified) (1 diff)

issm/trunk-jpl/test/MITgcm/code_4003/SHELFICE_OPTIONS.h ¶

@@ -16 +16 @@
 
 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 */

issm/trunk-jpl/test/MITgcm/code_4003/packages.conf ¶

@@ -1 +1 @@
 gfd
 shelfice
-streamice
 obcs
 diagnostics

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

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

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

issm/trunk-jpl/test/NightlyRun/test4003.m ¶

@@ -4 +4 @@
 %Script control parameters
 steps=1:11;
+steps=1:5;
 final_time=1/365;