Changeset 3101 for trunk

Timestamp:

Oct 25, 2023, 1:46:14 PM (13 months ago)

Author:

emillour

Message:

correction integer to real + added RCPV_generic

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• generic_cloud_common_h.F90 (modified) (2 diffs)
• initracer.F90 (modified) (3 diffs)
• tracer_h.F90 (modified) (1 diff)

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

@@ -22 +22 @@
     real,save :: RLVTT_generic !Latent heat of vaporization (J/kg)
     real,save :: metallicity_coeff ! Coefficient to take into account the metallicity
-
+    real,save :: RCPV_generic ! specific heat capacity of the tracer vapor at Tref
     
     contains 
@@ -158 +158 @@
                    print*, 'metallicity_coeff', metallicity_coeff
                 end if
-
+                if (index(table_traceurs_line,'RCPV_generic='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'RCPV_generic=')+len('RCPV_generic='):),*) RCPV_generic
+                    print*, 'RCPV_generic', RCPV_generic
+                end if
+                         
                 ios=1
             end if

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

@@ -128 +128 @@
        IF (.NOT. allocated(constants_RLVTT_generic)) allocate(constants_RLVTT_generic(nq))
        IF (.NOT. allocated(constants_metallicity_coeff)) allocate(constants_metallicity_coeff(nq))
+       IF (.NOT. allocated(constants_RCPV_generic)) allocate(constants_RCPV_generic(nq))
 
        !! initialization
@@ -157 +158 @@
        constants_RLVTT_generic(:)=0
        constants_metallicity_coeff(:)=0
+       constants_RCPV_generic(:)=0
 
        rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
@@ -477 +479 @@
                 constants_RLVTT_generic(iq)=RLVTT_generic
                 constants_metallicity_coeff(iq)=metallicity_coeff
+                constants_RCPV_generic(iq)=RCPV_generic
         else
                 write(*,*) "This tracer is not condensable, for generic condensation :  : ", noms(iq)

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

@@ -44 +44 @@
        integer,save,allocatable :: is_rgcs(:)               ! 1 if tracer is a radiative generic condensable specie, else 0 (added LT 2022)
        ! Lists of constants for condensable tracers
-       integer, save, allocatable :: constants_mass(:)                    ! molecular mass of the specie (g/mol)
-       integer, save, allocatable :: constants_delta_vapH(:)           ! Enthalpy of vaporization (J/mol)
-       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_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
+       real, save, allocatable :: constants_mass(:)                    ! molecular mass of the specie (g/mol)
+       real, save, allocatable :: constants_delta_vapH(:)           ! Enthalpy of vaporization (J/mol)
+       real, save, allocatable :: constants_Tref(:)                 ! Ref temperature for Clausis-Clapeyron (K)
+       real, save, allocatable :: constants_Pref(:)                 ! Reference pressure for Clausius Clapeyron (Pa)
+       real, save, allocatable :: constants_epsi_generic(:)                 ! fractionnal molecular mass (m/mugaz)
+       real, save, allocatable :: constants_RLVTT_generic(:)                ! Latent heat of vaporization (J/kg)
+       real, save, allocatable :: constants_metallicity_coeff(:)    ! Coefficient to take into account the metallicity
+       real, save, allocatable :: constants_RCPV_generic(:)                   ! specific heat capacity of the tracer vapor at Tref
 !$OMP THREADPRIVATE(constants_mass,constants_delta_vapH,constants_Tref)
 !$OMP THREADPRIVATE(constants_Pref,constants_epsi_generic)
-!$OMP THREADPRIVATE(constants_RLVTT_generic,constants_metallicity_coeff)
+!$OMP THREADPRIVATE(constants_RLVTT_generic,constants_metallicity_coeff,constants_RCPV_generic)
 
 !$OMP THREADPRIVATE(is_condensable,is_rgcs)   !also added by LT