Changeset 2690 for trunk/LMDZ.GENERIC/libf

Timestamp:

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

Author:

aslmd

Message:

1st commit of generic clouds scheme

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• evap_generic.F90 (added)
• genericcommon_h.F90 (added)
• initracer.F90 (added)
• largescale_generic_mod.F90 (added)
• physiq_mod.F90 (modified) (5 diffs)
• tracer_h.F90 (modified) (4 diffs)

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

@@ -32 +32 @@
                           alpha_lift, alpha_devil, qextrhor, &
                           igcm_h2o_ice, igcm_h2o_vap, igcm_dustbin, &
-                          igcm_co2_ice, nesp, is_chim
+                          igcm_co2_ice, nesp, is_chim, is_generic
       use time_phylmdz_mod, only: ecritphy, iphysiq, nday
       use phyetat0_mod, only: phyetat0
@@ -70 +70 @@
       use turb_mod, only : q2,sensibFlux,turb_resolved
       use mass_redistribution_mod, only: mass_redistribution
+      use largescale_generic_mod, only: largescale_generic 
 #ifndef MESOSCALE
       use vertical_layers_mod, only: presnivs, pseudoalt
@@ -295 +296 @@
       ! For Atmospheric Temperatures : (K/s)    
       real dtlscale(ngrid,nlayer)                             ! Largescale routine.
+      real dtlscale_generic(ngrid,nlayer)                     ! largescale_generic routine. 
       real zdtc(ngrid,nlayer)                                 ! Condense_co2 routine.
       real zdtdif(ngrid,nlayer)                               ! Turbdiff/vdifc routines.
@@ -327 +329 @@
       real dqvaplscale(ngrid,nlayer)  ! Largescale routine.
       real dqcldlscale(ngrid,nlayer)  ! Largescale routine.
+      real dqvaplscale_generic(ngrid,nlayer) !largescale_generic routine.
+      real dqcldlscale_generic(ngrid,nlayer) !largescale_generic routine.
       REAL,allocatable,save :: zdqchim(:,:,:) ! Calchim_asis routine
       REAL,allocatable,save :: zdqschim(:,:)  ! Calchim_asis routine
@@ -1542 +1546 @@
   !   VI.3. Aerosol particles
   ! -------------------------
+         ! Generic condensation and evaporation
+         ! i need the tracers and to check wether it ends with _vap
+         do iq=1,nq
+            ! print*,"starting the generic looping"
+            ! write(*,*)"for iq = ",iq, "is_generic = ",is_generic(iq)
+            ! write(*,*)"noms(iq) = ",noms(iq)
+            ! write(*,*),"noms(iq)(len(noms(iq)-2:)=",noms(iq)(len(noms(iq))-2:)
+            !loop on tracers and check if is_generic and _vap in noms
+            if((is_generic(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
+               ! call to largescale_generic
+               print*,"we're in the GENERIC condition"
+               call largescale_generic(ngrid,nlayer,nq,ptimestep,pplev,pplay,   &
+                                       pt,pq,pdt,pdq,noms(iq),dtlscale_generic, &
+                                       dqvaplscale_generic,dqcldlscale_generic)
+               !question sur la tendance de temperature qui sort. On l'ajoute dans une variable temp 
+               ! et apres la fin de la boucle sur les tracers on l'ajouteà la tendance du pas de temps physique ?
+               !temp_dt = temp_dt+dtlscale_generic ?
+               pdq(1:ngrid,1:nlayer,ig) = pdq(1:ngrid,1:nlayer,ig)+dtlscale_generic(1:ngrid,1:nlayer)
+               pdq(1:ngrid,1:nlayer,ig+1) = pdq(1:ngrid,1:nlayer,ig+1)+dqcldlscale_generic(1:ngrid,1:nlayer)
+               !!!!huge assumption in the last line : ice tracer has index +1 of vap tracer
+            endif ! (is_generic(iq)==1 .and. (index(noms(iq),"vap")==1))
+         enddo ! iq=1,nq
+         !should I add the pdt = pdt + temp_dt ?
 
          ! Sedimentation.

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

@@ -6 +6 @@
        integer, save :: nqtot ! total number of tracers
        integer, save :: nesp  ! number of species in the chemistry
+       integer, save :: ngt   ! number of generic tracers
 !$OMP THREADPRIVATE(nqtot,nesp)
 
@@ -12 +13 @@
 
        character*30, save, allocatable :: noms(:)   ! name of the tracer
-       real, save, allocatable :: mmol(:)     ! mole mass of tracer (g/mol-1) 
+       real, save, allocatable :: mmol(:)     ! mole mass of tracer (g/mol) 
        real, save, allocatable :: aki(:)      ! to compute coefficient of thermal concduction if photochem
        real, save, allocatable :: cpi(:)      ! to compute cpnew in concentration.F if photochem
@@ -39 +40 @@
        integer, save, allocatable :: is_recomb_qotf(:) ! 1 if tracer recombination is done on-the-fly, else 0 (if 1, must have is_recomb_qset=0)
 !$OMP THREADPRIVATE(is_recomb,is_recomb_qset,is_recomb_qotf)
-
+       integer, save, allocatable :: is_generic(:)      ! 1 if tracer is generic, else 0 (added LT)
+!$OMP THREADPRIVATE(is_generic)   !also added by LT
 ! tracer indexes: these are initialized in initracer and should be 0 if the
 !                 corresponding tracer does not exist
@@ -88 +90 @@
        integer,save :: igcm_hcaer
 
-
+       ! generic tracer (added LT 2022)
+       integer,save :: igcm_generic_Fe_vap
+       integer,save :: igcm_generic_Fe_ice
+       integer,save :: igcm_generic_Cr_vap
+       integer,save :: igcm_generic_Cr_ice
+       integer,save :: igcm_generic_Mg_vap
+       integer,save :: igcm_generic_Mg_ice
+       integer,save :: igcm_generic_Na_vap
+       integer,save :: igcm_generic_Na_ice
+       integer,save :: igcm_generic_KCl_vap
+       integer,save :: igcm_generic_KCl_ice
+       integer,save :: igcm_generic_Mn_vap
+       integer,save :: igcm_generic_Mn_ice
+       integer,save :: igcm_generic_Zn_vap
+       integer,save :: igcm_generic_Zn_ice
 
        ! other tracers
        integer,save :: igcm_ar_n2 ! for simulations using co2 +neutral gaz
        integer,save :: igcm_co2_ice ! CO2 ice 
-!$OMP THREADPRIVATE(igcm_dustbin,igcm_dust_mass,igcm_dust_number,igcm_h2o_vap,igcm_h2o_ice, &
-        !$OMP igcm_co2,igcm_co,igcm_o,igcm_o1d,igcm_o2,igcm_o3,igcm_h,igcm_h2,igcm_oh,      &
-        !$OMP igcm_ho2,igcm_h2o2,igcm_n2,igcm_ar,igcm_ar_n2,igcm_co2_ice,                   &
-        !$OMP igcm_n,igcm_no,igcm_no2,igcm_n2d,igcm_ch4,igcm_ch3,igcm_ch,igcm_3ch2,         &
-        !$OMP igcm_1ch2,igcm_cho,igcm_ch2o,igcm_ch3o,igcm_c,igcm_c2,igcm_c2h,igcm_c2h2,     &
-        !$OMP igcm_c2h3,igcm_c2h4,igcm_c2h6,igcm_ch2co,igcm_ch3co,igcm_hcaer)
+!$OMP THREADPRIVATE(igcm_dustbin,igcm_dust_mass,igcm_dust_number,igcm_h2o_vap,igcm_h2o_ice,      &
+        !$OMP igcm_co2,igcm_co,igcm_o,igcm_o1d,igcm_o2,igcm_o3,igcm_h,igcm_h2,igcm_oh,        &
+        !$OMP igcm_ho2,igcm_h2o2,igcm_n2,igcm_ar,igcm_ar_n2,igcm_co2_ice,                         &
+       !$OMP igcm_n,igcm_no,igcm_no2,igcm_n2d,igcm_ch4,igcm_ch3,igcm_ch,igcm_3ch2,               &
+       !$OMP igcm_1ch2,igcm_cho,igcm_ch2o,igcm_ch3o,igcm_c,igcm_c2,igcm_c2h,igcm_c2h2,           &
+       !$OMP igcm_c2h3,igcm_c2h4,igcm_c2h6,igcm_ch2co,igcm_ch3co,igcm_hcaer,igcm_generic_Fe_vap, &
+       !$OMP igcm_generic_Fe_ice, igcm_generic_Cr_vap,igcn_generic_Cr_ice,igcm_generic_Mg_vap,   &
+       !$OMP igcm_generic_Mg_ice, igcm_generic_Na_vap,igcm_generic_Na_ice,igcm_generic_KCl_vap,  &
+       !$OMP igcm_generic_KCl_ice,igcm_generic_Mn_vap,igcm_generic_Mn_ice,igcm_generic_Zn_vap,   &
+       !$OMP igcm_generic_Zn_ice)
 
        end module tracer_h