Changeset 2915 for LMDZ5/trunk/libf/phylmd/soil.F90

Timestamp:

Jun 16, 2017, 6:03:16 PM (7 years ago)

Author:

jbmadeleine

Message:

Improved the way thermal inertia works in simulations using imposed SST and sea ice.
Removed inertie_ice which was confusing; Added inertie_sic for sea-ice and inertie_lic for
land-ice. When snow is present, inertie_sno is used instead of inertie_sol, inertie_sic or
inertie_lic. For sea-ice, added a flag called iflag_sic:
If iflag_sic=0, thermal inertia is changed over sea-ice when snow is present
If iflag_sic=1, thermal inertia is kept constant and equal to inertie_sic

File:

• LMDZ5/trunk/libf/phylmd/soil.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/soil.F90

@@ -152 +152 @@
 !-----------------------------------------------------------------------
 !   Calcul de l'inertie thermique a partir de la variable rnat.
-!   on initialise a inertie_ice meme au-dessus d'un point de mer au cas 
+!   on initialise a inertie_sic meme au-dessus d'un point de mer au cas 
 !   ou le point de mer devienne point de glace au pas suivant
 !   on corrige si on a un point de terre avec ou sans glace
 !
-!-----------------------------------------------------------------------
+!   iophys can be used to write the ztherm_i variable in a phys.nc file
+!   and check the results; to do so, add "CALL iophys_ini" in physiq_mod
+!   and add knindex to the list of inputs in all the calls to soil.F90
+!   (and to soil.F90 itself !)
+!-----------------------------------------------------------------------
+
   IF (indice == is_sic) THEN
      DO ig = 1, knon
-        ztherm_i(ig)   = inertie_ice
+        ztherm_i(ig)   = inertie_sic
+     ENDDO
+     IF (iflag_sic == 0) THEN
+       DO ig = 1, knon
+         IF (snow(ig) > 0.0) ztherm_i(ig)   = inertie_sno
+       ENDDO
+!      Otherwise sea-ice keeps the same inertia, even when covered by snow
+     ENDIF
+!    CALL iophys_ecrit_index('ztherm_sic', 1, 'ztherm_sic', 'USI', &
+!      knon, knindex, ztherm_i)
+  ELSE IF (indice == is_lic) THEN
+     DO ig = 1, knon
+        ztherm_i(ig)   = inertie_lic
         IF (snow(ig) > 0.0) ztherm_i(ig)   = inertie_sno
      ENDDO
-  ELSE IF (indice == is_lic) THEN
-     DO ig = 1, knon
-        ztherm_i(ig)   = inertie_ice
-        IF (snow(ig) > 0.0) ztherm_i(ig)   = inertie_sno
-     ENDDO
+!    CALL iophys_ecrit_index('ztherm_lic', 1, 'ztherm_lic', 'USI', &
+!      knon, knindex, ztherm_i)
   ELSE IF (indice == is_ter) THEN
      DO ig = 1, knon
@@ -172 +186 @@
         IF (snow(ig) > 0.0) ztherm_i(ig)   = inertie_sno
      ENDDO
+!    CALL iophys_ecrit_index('ztherm_ter', 1, 'ztherm_ter', 'USI', &
+!      knon, knindex, ztherm_i)
   ELSE IF (indice == is_oce) THEN
      DO ig = 1, knon
-        ztherm_i(ig)   = inertie_ice
-     ENDDO
+!       This is just in case, but SST should be used by the model anyway
+        ztherm_i(ig)   = inertie_sic
+     ENDDO
+!    CALL iophys_ecrit_index('ztherm_oce', 1, 'ztherm_oce', 'USI', &
+!      knon, knindex, ztherm_i)
   ELSE
      WRITE(lunout,*) "valeur d indice non prevue", indice