Changeset 4146 for trunk/LMDZ.GENERIC/libf/phygeneric/thermcell_env.F90

Timestamp:

Mar 19, 2026, 2:35:46 PM (12 days ago)

Author:

gmilcareck

Message:

Thermodynamics update on LMDZ.GENERIC

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/thermcell_env.F90 (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/thermcell_env.F90

@@ -26 +26 @@
       USE watercommon_h, ONLY: RLvCp, RETV, Psat_water
       USE tracer_h, ONLY: igcm_h2o_vap, igcm_h2o_ice
-      USE callkeys_mod, ONLY: water, generic_condensation
-      USE comcstfi_mod, ONLY: r, cpp, mugaz
+      USE callkeys_mod, ONLY: water, generic_condensation, metallicity
+      USE comcstfi_mod, ONLY: rd_ref, cppd_ref, mugaz_ref
       USE generic_cloud_common_h, ONLY: Psat_generic, epsi_generic, RLVTT_generic
       USE generic_tracer_index_mod, ONLY: generic_tracer_index
@@ -75 +75 @@
       INTEGER :: igcm_generic_vap, igcm_generic_ice! index of the vap and ice of generic_tracer
       LOGICAL :: call_ice_vap_generic ! to call only one time the ice/vap pair of a tracer
-      REAL, SAVE :: metallicity ! metallicity of planet --- is not used here, but necessary to call function Psat_generic
-      !$OMP THREADPRIVATE(metallicity)
       REAL :: RETV_generic, RV_generic, RLvCp_generic
 
@@ -90 +88 @@
       zqt(:,:) = 0.
       zql(:,:) = 0.
-      
-      metallicity=0.0 ! default value --- is not used here but necessary to call function Psat_generic
-       
+         
 !===============================================================================
 ! Condensation and latent heat release
@@ -132 +128 @@
                ENDDO
 
-               RV_generic = (8.314463*1000.)/(epsi_generic*mugaz)
-               RETV_generic = RV_generic/r-1.
-               RLvCp_generic = RLVTT_generic/cpp
+               RV_generic = (8.314463*1000.)/(epsi_generic*mugaz_ref)
+               RETV_generic = RV_generic/rd_ref-1.
+               RLvCp_generic = RLVTT_generic/cppd_ref
 
                DO k = 1,nlay