Changeset 2711 for trunk/LMDZ.GENERIC/libf

Timestamp:

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

Author:

aslmd

Message:

now data for generic tracers are called only one time !

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (2 diffs)
• initracer.F90 (modified) (5 diffs)

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

@@ -7 +7 @@
                 pt, pq, pdt, pdq, pdtlsc, pdqvaplsc, pdqliqlsc)
         use ioipsl_getin_p_mod, only: getin_p !-> to get the metallicity 
-        use generic_cloud_common_h, only : RLVTT, cpp, &
-        Psat_generic,Lcpdqsat_generic,specie_parameters,specie_parameters_table
+        use generic_cloud_common_h
         USE tracer_h
         IMPLICIT none
@@ -107 +106 @@
                         end if
                         if (igcm_generic_ice .eq. -1) then
-                                write(*,*) "ERROR : You set a vap traceur but you forgot to set the corresponding ice traceur"
+                                write(*,*) "ERROR : You set a vap traceur but you forgot to set the corresponding ice traceur, &
+                                or the pair vap/ice is not written one after another in traceur.def"
                         endif
 
-                        ! Need to call specie_parameters of the considered specie 
-                        !call specie_parameters(noms(iq)(9:len(trim(noms(iq)))-4))
-                        write(*,*) 'looking specie parameters for : ',noms(iq)(1:len(trim(noms(iq)))-4)
-                        call specie_parameters_table(noms(iq)(1:len(trim(noms(iq)))-4))
-                        
+                        m=constants_mass(iq)
+                        delta_vapH=constants_delta_vapH(iq)
+                        Tref=constants_Tref(iq)
+                        Pref=constants_Pref(iq)
+                        epsi=constants_epsi(iq)
+                        RLVTT=constants_RLVTT(iq)
+                        metallicity_coeff=constants_metallicity_coeff(iq)
+
                         Lcp=RLVTT/cpp ! need to be init here
 

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

@@ -6 +6 @@
       USE recombin_corrk_mod, ONLY: ini_recombin
       USE mod_phys_lmdz_para, only: is_master, bcast
+      use generic_cloud_common_h
       IMPLICIT NONE
 !=======================================================================
@@ -117 +118 @@
        ENDIF
        IF (.NOT. allocated(is_condensable)) allocate(is_condensable(nq)) !LT
+       IF (.NOT. allocated(constants_mass)) allocate(constants_mass(nq))
+       IF (.NOT. allocated(constants_delta_vapH)) allocate(constants_delta_vapH(nq))
+       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_metallicity_coeff)) allocate(constants_metallicity_coeff(nq))
+
        !! initialization
        alpha_lift(:)     = 0.
@@ -133 +142 @@
        radius(:)=0.
        qext(:)=0.
-       is_condensable(:)= 0 !LT
+
+       ! For condensable tracers, by Lucas Teinturier and Noé Clément (2022)
+
+       is_condensable(:)= 0
+
+       constants_mass(:)=0
+       constants_delta_vapH(:)=0
+       constants_Tref(:)=0
+       constants_Pref(:)=0
+       constants_epsi(:)=0
+       constants_RLVTT(:)=0
+       constants_metallicity_coeff(:)=0
+
        rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
 
@@ -192 +213 @@
       igcm_ch3co=0
       igcm_hcaer=0
-      ! ! LT generic tracers 
-      ! igcm_generic_Fe_vap = 0
-      ! igcm_generic_Fe_ice = 0 
-      ! igcm_generic_Cr_vap = 0
-      ! igcm_generic_Cr_ice = 0
-
 
       ! 1. find dust tracers
@@ -435 +450 @@
 
       if (is_master) close(407)
+
+      ! Get specific data of condensable tracers
+      do iq=1,nq
+        if((is_condensable(iq)==1)) then
+                write(*,*) "There is a specie which is condensable, for generic condensation : ", noms(iq)
+                write(*,*) 'looking specie parameters for : ',noms(iq)(1:len(trim(noms(iq)))-4)
+                call specie_parameters_table(noms(iq)(1:len(trim(noms(iq)))-4))
+                constants_mass(iq)=m
+                constants_delta_vapH(iq)=delta_vapH
+                constants_Tref(iq)=Tref
+                constants_Pref(iq)=Pref
+                constants_epsi(iq)=epsi
+                constants_RLVTT(iq)=RLVTT
+                constants_metallicity_coeff(iq)=metallicity_coeff
+        else
+                write(*,*) "This tracer is not condensable, for generic condensation :  : ", noms(iq)
+                write(*,*) "We keep condensable constants at zero"
+        endif !(is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
+      enddo ! iq=1,nq
 
       ! Calculate number of species in the chemistry