Changeset 2725 for trunk/LMDZ.GENERIC/libf

Timestamp:

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

Author:

aslmd

Message:

switch to epsi_generic & RLVTT_generic

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (1 diff)
• evap_generic.F90 (modified) (1 diff)
• generic_cloud_common_h.F90 (modified) (13 diffs)
• initracer.F90 (modified) (3 diffs)
• rain_generic.F90 (modified) (3 diffs)
• tracer_h.F90 (modified) (1 diff)

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

@@ -107 +107 @@
                         Tref=constants_Tref(iq)
                         Pref=constants_Pref(iq)
-                        epsi=constants_epsi(iq)
-                        RLVTT=constants_RLVTT(iq)
+                        epsi_generic=constants_epsi_generic(iq)
+                        RLVTT_generic=constants_RLVTT_generic(iq)
                         metallicity_coeff=constants_metallicity_coeff(iq)
 
-                        Lcp=RLVTT/cpp ! need to be init here
+                        Lcp=RLVTT_generic/cpp ! need to be init here
 
                         !  Vertical loop (from top to bottom)

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

@@ -38 +38 @@
 
     !   Evaporate all the ice from the tracer
-    zlvdcp = RLVTT/cpp ! RLVTT is the latent heat of vaporization (comes from generic_cloud_common_h attention)
+    zlvdcp = RLVTT_generic/cpp ! RLVTT_generic is the latent heat of vaporization (comes from generic_cloud_common_h attention)
  
     DO l=1,nlayer

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
 

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

@@ -122 +122 @@
        IF (.NOT. allocated(constants_Tref)) allocate(constants_Tref(nq))
        IF (.NOT. allocated(constants_Pref)) allocate(constants_Pref(nq))
-       IF (.NOT. allocated(constants_epsi)) allocate(constants_epsi(nq))
-       IF (.NOT. allocated(constants_RLVTT)) allocate(constants_RLVTT(nq))
+       IF (.NOT. allocated(constants_epsi_generic)) allocate(constants_epsi_generic(nq))
+       IF (.NOT. allocated(constants_RLVTT_generic)) allocate(constants_RLVTT_generic(nq))
        IF (.NOT. allocated(constants_metallicity_coeff)) allocate(constants_metallicity_coeff(nq))
 
@@ -151 +151 @@
        constants_Tref(:)=0
        constants_Pref(:)=0
-       constants_epsi(:)=0
-       constants_RLVTT(:)=0
+       constants_epsi_generic(:)=0
+       constants_RLVTT_generic(:)=0
        constants_metallicity_coeff(:)=0
 
@@ -461 +461 @@
                 constants_Tref(iq)=Tref
                 constants_Pref(iq)=Pref
-                constants_epsi(iq)=epsi
-                constants_RLVTT(iq)=RLVTT
+                constants_epsi_generic(iq)=epsi_generic
+                constants_RLVTT_generic(iq)=RLVTT_generic
                 constants_metallicity_coeff(iq)=metallicity_coeff
         else

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

@@ -133 +133 @@
          Tref=constants_Tref(iq)
          Pref=constants_Pref(iq)
-         epsi=constants_epsi(iq)
-         RLVTT=constants_RLVTT(iq)
+         epsi_generic=constants_epsi_generic(iq)
+         RLVTT_generic=constants_RLVTT_generic(iq)
 
          RCPD = cpp
@@ -265 +265 @@
                      if(zt(i,k).gt.Tsat(i,k))then ! if temperature of the layer box is greater than Tsat
                         ! treat the case where all liquid water should boil
-                        zqev=MIN((zt(i,k)-Tsat(i,k))*RCPD*dmass(i,k)/RLVTT/ptimestep,precip_rate(i)) ! on évapore
+                        zqev=MIN((zt(i,k)-Tsat(i,k))*RCPD*dmass(i,k)/RLVTT_generic/ptimestep,precip_rate(i)) ! on évapore
                         precip_rate(i)=MAX(precip_rate(i)-zqev,0.) ! we withdraw from precip_rate the evaporated part
                         d_q(i,k)=zqev/dmass(i,k)*ptimestep ! quantité évaporée
-                        d_t(i,k) = - d_q(i,k) * RLVTT/RCPD
+                        d_t(i,k) = - d_q(i,k) * RLVTT_generic/RCPD
                      else
                         zqev = MAX (0.0, (zqs(i,k)-q(i,k)))*dmass(i,k)/(ptimestep*(1.d0+dqsat(i,k)))
@@ -280 +280 @@
 
                         d_q(i,k) = - (precip_rate_tmp(i)-precip_rate(i))/dmass(i,k)*ptimestep
-                        d_t(i,k) = - d_q(i,k) * RLVTT/RCPD
+                        d_t(i,k) = - d_q(i,k) * RLVTT_generic/RCPD
                         precip_rate(i)  = precip_rate_tmp(i)
                      end if

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

@@ -47 +47 @@
        integer, save, allocatable :: constants_Tref(:)                 ! Ref temperature for Clausis-Clapeyron (K)
        integer, save, allocatable :: constants_Pref(:)                 ! Reference pressure for Clausius Clapeyron (Pa)
-       integer, save, allocatable :: constants_epsi(:)                 ! fractionnal molecular mass (m/mugaz)
-       integer, save, allocatable :: constants_RLVTT(:)                ! Latent heat of vaporization (J/kg)
+       integer, save, allocatable :: constants_epsi_generic(:)                 ! fractionnal molecular mass (m/mugaz)
+       integer, save, allocatable :: constants_RLVTT_generic(:)                ! Latent heat of vaporization (J/kg)
        integer, save, allocatable :: constants_metallicity_coeff(:)    ! Coefficient to take into account the metallicity