Changeset 2725 for trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90

Timestamp:

Jun 21, 2022, 11:05:47 AM (2 years ago)

Author:

aslmd

Message:

switch to epsi_generic & RLVTT_generic

File:

• trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90 (modified) (13 diffs)

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

@@ -19 +19 @@
     real,save :: Tref ! Ref temperature for Clausis-Clapeyron (K)
     real,save :: Pref ! Reference pressure for Clausius Clapeyron (Pa)
-    real, save :: epsi ! fractionnal molecular mass (m/mugaz)
-    real,save :: RLVTT !Latent heat of vaporization (J/kg)
+    real, save :: epsi_generic ! fractionnal molecular mass (m/mugaz)
+    real,save :: RLVTT_generic !Latent heat of vaporization (J/kg)
     real,save :: metallicity_coeff ! Coefficient to take into account the metallicity
 
@@ -35 +35 @@
     !   Values are based on Vissher (2010) and Morley (2012).
     !
-    !   Also set up few others useful parameters, such as epsi=m/mugaz, RLVTT and
+    !   Also set up few others useful parameters, such as epsi_generic=m/mugaz, RLVTT_generic and
     !   the metallicity_coeff.
     !   Authors 
@@ -50 +50 @@
             Pref =   1.0E5 !in Pa
             m = 140.6931
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
             metallicity_coeff=1.0*log(10.0)
             
@@ -59 +59 @@
             Pref =  1.0E5
             m = 78.04454 !Na2S
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
             metallicity_coeff=0.5*log(10.0)
             
@@ -68 +68 @@
             Pref  = 1.0E5
             m = 55.8450
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
             metallicity_coeff=0.0*log(10.0)
             
@@ -77 +77 @@
             Pref  = 1.0E5
             m = 51.9961
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
             metallicity_coeff=0.0*log(10.0)
             
@@ -87 +87 @@
             metallicity_coeff=0.0*log(10.0)
             m = 74.5498
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
         else if (trim(specname) .eq. "Mn") then 
             print*,"Loading data for Mn"
@@ -95 +95 @@
             metallicity_coeff=1.0*log(10.0)
             m = 87.003049
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
         else if (trim(specname) .eq. "Zn") then 
             print*,"Loading data for Zn"
@@ -103 +103 @@
             metallicity_coeff=1.0*log(10.0)
             m = 97.445
-            RLVTT = delta_vapH/(m/1000.)
+            RLVTT_generic = delta_vapH/(m/1000.)
         else
             print*,"Unknow species (not in Mg, Fe, Na, KCl, Cr, Mn or Zn)"
         endif
-        epsi = m/mugaz
+        epsi_generic = m/mugaz
     end subroutine specie_parameters
 
@@ -118 +118 @@
     !   From a table of traceurs
     !
-    !   Also set up few others useful parameters, such as epsi=m/mugaz, RLVTT and
+    !   Also set up few others useful parameters, such as epsi_generic=m/mugaz, RLVTT_generic and
     !   the metallicity_coeff.
     !   Authors 
@@ -178 +178 @@
         close(117)
 
-        RLVTT=delta_vapH/(m/1000.)
-
-        write(*,*) 'RLVTT', RLVTT
-
-        epsi = m/mugaz
+        RLVTT_generic=delta_vapH/(m/1000.)
+
+        write(*,*) 'RLVTT_generic', RLVTT_generic
+
+        epsi_generic = m/mugaz
 
     end subroutine specie_parameters_table
@@ -208 +208 @@
             qsat = 1. ! is the maximum amount of vapor that a parcel can hold without condensation, in specific concentration.
         else 
-            qsat = epsi *psat/(p-(1-epsi)*psat)
+            qsat = epsi_generic *psat/(p-(1-epsi_generic)*psat)
         endif
     end subroutine Psat_generic 
@@ -240 +240 @@
             dqsat =0.
         else 
-            dqsat = (RLVTT/cpp) *qsat*(p/(p-(1-epsi)*psat))*dummy
-            dlnpsat = (RLVTT/cpp) * dummy 
+            dqsat = (RLVTT_generic/cpp) *qsat*(p/(p-(1-epsi_generic)*psat))*dummy
+            dlnpsat = (RLVTT_generic/cpp) * dummy 
         endif