Changeset 2064 for trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

Timestamp:

Jan 10, 2019, 3:44:36 PM (6 years ago)

Author:

aboissinot

Message:

Fix a bug in thermcell_closure arguments which implies wrong fluxes normalization.

File:

• trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

@@ -39 +39 @@
 ! AB : linf is used to set the lowest possible first level because we allow it
 !      to begin higher than the surface. It is set to 2 in order to remove the
-!      first layer for gas giant. If there is a surface, it has to be set to 1.
+!      first layer for gas giant.
+!      If there is a surface, it has to be set to 1.
+!      If someone want to call more than once the thermal plume model in one
+!      or more grid point, this variable must become a saved array of INTEGER
+!      with size ngrid. 
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 INTEGER,PARAMETER :: linf                       = 2
@@ -59 +63 @@
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 ! AB : Parameters needed only for a loop in thermcell_alp (diagnoses).
-!      Maybe to remove.
+!      Maybe to be removed.
 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 INTEGER,PARAMETER :: nbsrf                      = 1
@@ -65 +69 @@
 
 CONTAINS
-
-   SUBROUTINE init_thermcell_mod(g, rcp, r, pi, T_h2o_ice_liq, RV)
-   
-   IMPLICIT NONE
-   
-   REAL g
-   REAL rcp
-   REAL r
-   REAL pi
-   REAL T_h2o_ice_liq
-   REAL RV
-   
-   RTT = T_h2o_ice_liq
-   RG = g
-   RKAPPA = rcp
-   RPI = pi
-   RD = r
-   
-   RETURN
-   END SUBROUTINE init_thermcell_mod
+      
+      SUBROUTINE init_thermcell_mod(g, rcp, r, pi, T_h2o_ice_liq, RV)
+      
+      IMPLICIT NONE
+      
+      REAL g
+      REAL rcp
+      REAL r
+      REAL pi
+      REAL T_h2o_ice_liq
+      REAL RV
+      
+      RTT = T_h2o_ice_liq
+      RG = g
+      RKAPPA = rcp
+      RPI = pi
+      RD = r
+      
+      RETURN
+      END SUBROUTINE init_thermcell_mod
 
 END MODULE thermcell_mod