Changeset 2720 for trunk/LMDZ.GENERIC

Timestamp:

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

Author:

aslmd

Message:

new module to get back indexes of cond tracers

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (5 diffs)
• generic_tracer_index_mod.F90 (added)

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

@@ -10 +10 @@
         USE tracer_h
         USE mod_phys_lmdz_para, only: is_master
+        use generic_tracer_index_mod, only: generic_tracer_index
         IMPLICIT none
 
@@ -45 +46 @@
 !       Options : 
         real, save :: metallicity !metallicity of planet
-!$OMP THREADPRIVATE(metallicity)
+        REAL, SAVE :: qvap_deep   ! deep mixing ratio of water vapor when simulating bottom less planets
+!$OMP THREADPRIVATE(metallicity, qvap_deep)
 
 !       Local variables
+
+        ! to call only one time the ice/vap pair of a tracer
+        logical one_call_ice_vap_generic
+
         INTEGER i, k , nn, iq
         INTEGER,PARAMETER :: nitermax=5000
+        INTEGER,PARAMETER :: tau=14400.
         DOUBLE PRECISION,PARAMETER :: alpha=.1,qthreshold=1.d-8
         ! JL13: if "careful, T<Tmin in psat water" appears often, you may want to stabilise the model by
@@ -72 +79 @@
                 call getin_p("metallicity",metallicity)
                 write(*,*) " metallicity = ",metallicity
+
+                write(*,*) "Deep water vapor mixing ratio ? (no effect if negative) "
+                qvap_deep=-1. ! default value
+                call getin_p("qvap_deep",qvap_deep)
+                write(*,*) " qvap_deep = ",qvap_deep   
+
                 firstcall = .false.
         ENDIF
@@ -84 +97 @@
         ! Let's loop on tracers 
         do iq=1,nq
-                if((is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
-                        if (is_master) then
-                                write(*,*) "There is a specie which is condensable, for generic condensation : ", noms(iq)
-                        endif
-!                       Let's get the index of our tracers (we look for igcm _generic_vap and igcm_generic_ice)
-                        igcm_generic_vap=iq
-                        igcm_generic_ice = -1
+                one_call_ice_vap_generic = .false.
 
-                        ! We look for the corresponding ice traceur (before or after in the list of traceurs, maybe could be generalised to the whole list)
-                        if (iq .ne. nq) then
-                                if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq+1)(1:len(trim(noms(iq+1)))-4)) .and. (index(noms(iq+1),"ice") .ne. 0)) then
-                                        igcm_generic_ice = iq+1
-                                end if
-                        end if
-                        if ((iq .gt. 1)) then
-                                if ((noms(iq)(1:len(trim(noms(iq)))-4) .eq. noms(iq-1)(1:len(trim(noms(iq-1)))-4)) .and. (index(noms(iq-1),"ice") .ne. 0)) then
-                                        igcm_generic_ice = iq-1
-                                end if
-                        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, &
-                                or the pair vap/ice is not written one after another in traceur.def"
-                        endif
+                call generic_tracer_index(nq,iq,igcm_generic_vap,igcm_generic_ice,one_call_ice_vap_generic)
 
+                if(one_call_ice_vap_generic) then ! to call only one time the ice/vap pair of a tracer
                         m=constants_mass(iq)
                         delta_vapH=constants_delta_vapH(iq)
@@ -154 +148 @@
 
                         Enddo ! k= nlayer, 1, -1
-                endif !(is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
+
+                        if (qvap_deep >= 0.) then
+                                !brings lower generic vapor ratio to a fixed value.
+                                ! tau=3600. seems too fast
+                                pdqvaplsc(1:ngrid,1,igcm_generic_vap) = (qvap_deep - pq(1:ngrid,1,igcm_generic_vap))/tau - pdq(1:ngrid,1,igcm_generic_vap)
+                        endif
+                endif !if(one_call_ice_vap_generic)
         enddo ! iq=1,nq