Changeset 2915 for LMDZ5/trunk/libf

Timestamp:

Jun 16, 2017, 6:03:16 PM (8 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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• comsoil.h (modified) (1 diff)
• conf_phys_m.F90 (modified) (4 diffs)
• sisvat/sisvat.F (modified) (1 diff)
• soil.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/comsoil.h

@@ -3 +3 @@
 !
 
-      common /comsoil/inertie_sol,inertie_sno,inertie_ice
-      real inertie_sol,inertie_sno,inertie_ice
+      common /comsoil/inertie_sol,inertie_sno,inertie_sic,inertie_lic,  &
+     &                iflag_sic
+      real inertie_sol,inertie_sno,inertie_sic,inertie_lic
+      integer iflag_sic
 !$OMP THREADPRIVATE(/comsoil/)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -173 +173 @@
     REAL,SAVE :: exposant_glace_omp
     REAL,SAVE :: rei_min_omp, rei_max_omp
-    REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
+    INTEGER,SAVE :: iflag_sic_omp
+    REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_sic_omp
+    REAL,SAVE :: inertie_lic_omp
     REAL,SAVE :: qsol0_omp
     REAL,SAVE :: evap0_omp
@@ -1133 +1135 @@
     !-----------------------------------------------------------------------
     !
-    !Config Key  = inertie_ice
+    !Config Key  = iflag_sic
+    !Config Desc =  
+    !Config Def  = 0
+    !Config Help = 
+    !
+    iflag_sic_omp = 0
+    CALL getin('iflag_sic',iflag_sic_omp)
+    !
+    !Config Key  = inertie_sic
     !Config Desc =  
     !Config Def  = 2000.
     !Config Help = 
     !
-    inertie_ice_omp = 2000.
-    CALL getin('inertie_ice',inertie_ice_omp)
+    inertie_sic_omp = 2000.
+    CALL getin('inertie_sic',inertie_sic_omp)
+    !
+    !Config Key  = inertie_lic
+    !Config Desc =  
+    !Config Def  = 2000.
+    !Config Help = 
+    !
+    inertie_lic_omp = 2000.
+    CALL getin('inertie_lic',inertie_lic_omp)
     !
     !Config Key  = inertie_sno
@@ -2132 +2150 @@
     evap0 = evap0_omp
     albsno0 = albsno0_omp
+    iflag_sic = iflag_sic_omp
     inertie_sol = inertie_sol_omp
-    inertie_ice = inertie_ice_omp
+    inertie_sic = inertie_sic_omp
+    inertie_lic = inertie_lic_omp
     inertie_sno = inertie_sno_omp
     rad_froid = rad_froid_omp
@@ -2545 +2565 @@
     write(lunout,*)' evap0 = ', evap0
     write(lunout,*)' albsno0 = ', albsno0
+    write(lunout,*)' iflag_sic = ', iflag_sic
     write(lunout,*)' inertie_sol = ', inertie_sol
-    write(lunout,*)' inertie_ice = ', inertie_ice
+    write(lunout,*)' inertie_sic = ', inertie_sic
+    write(lunout,*)' inertie_lic = ', inertie_lic
     write(lunout,*)' inertie_sno = ', inertie_sno
     write(lunout,*)' f_cdrag_ter = ',f_cdrag_ter

LMDZ5/trunk/libf/phylmd/sisvat/sisvat.F

@@ -7070 +7070 @@
 
       DO ig = 1,knonv
-        ztherm_i(ig)   = inertie_ice
+        ztherm_i(ig)   = inertie_lic
         IF (isnoSV(ig) > 0) ztherm_i(ig)   = inertie_sno
       ENDDO

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