Changeset 2810 for trunk/LMDZ.GENERIC

Timestamp:

Oct 29, 2022, 8:09:31 PM (2 years ago)

Author:

emillour

Message:

Generic GCM:
Minor fix related to the RGCS scheme: when not using it, "iaero_generic(1)"
was tested although not defined; using tests on "aerogeneric" (number
of generic aerosols) works better.
EM

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/aeropacity.F90 (modified) (4 diffs)
• libf/phystd/callkeys_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/README

@@ -1734 +1734 @@
 Tracer names are extracted from traceur.def via initracer.F90 and no
 longer transfered from the dynamics to the physics
+
+== 19/10/2022 == LT
+Initialisation of Radiative Generic Condensable Aerosols
+We can activate the scheme by putting aerogeneric = # of aerosols in callphys.def. This is the only needed thing for activating the radiative effects. They must be tracer in modern tracer.def
+Commented out the abort if we use more than 4 aerosols
+Added reading of optprop files for Radiative Generic Condensable Aerosols
+We use the same file for IR and VI channel. For now, only MnS, Cr,Fe and Mg2SiO4 can be read. If you want to add another specie, check the code, it is explained how to quickly do that (right above the Initializations)
+Added radii calculation for Radiative Generic Condensable aerosols
+Changed the hardcoded size of the totalemit array
+The hardcoded size is now 1900 instead of 100 so we don't exceed the array size when working at high spectral resolution (very rare case)
+Added opacity computation for Radiative Generic Condensable aerosols
+We do this computation in the same fashion as what's performed on water and dust.
+switch iniaerosol and initracer order, to prepare for the RGCS scheme
+Needed to switch the order of initialization so we can use the RGCS scheme without the assumption that ice and vap tracer of the same specie are following each other in the traceur.def fil
+The RGCS scheme: Now we can use this scheme to take into account radiative effect of a Generic Condensable Specie. One needs to add the keyword is_rgcs=1 to the ModernTraceur.def to activate the scheme + specify the number of aerogeneric specie in callphys.def. Also, one needs to go into suaer_corrk.F90 to add the Mie scattering file into the code for the specie you want. Currently, MnS, Fe, Cr, Mg2SiO4 are available. This has been tested on a Hot Jupiter case, might have issue in different configuration. Also, the numerical stability of the scheme is not guaranted when using strong scatterers such as silicates. Physical oscillations in the evaporation/condensation can occur, and need fixing. This will be for another commit.
+
+== 29/10/2022 == EM
+Minor fix related to the RGCS scheme: when not using it, "iaero_generic(1)"
+was tested although not defined; using tests on "aerogeneric" (number
+of generic aerosols) works better.

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

@@ -3 +3 @@
 
        use radinc_h, only : L_TAUMAX,naerkind
-       use aerosol_mod
+       use aerosol_mod, only: iaero_nlay, iaero_generic, &
+                              iaero_aurora, iaero_back2lay, iaero_co2, &
+                              iaero_dust, iaero_h2o, iaero_h2so4, &
+                              iaero_nh3, i_rgcs_ice, noaero
        USE tracer_h, only: noms,rho_co2,rho_ice,rho_q,mmol
        use comcstfi_mod, only: g, pi, mugaz, avocado
@@ -145 +148 @@
           print*,'iaero_aurora= ',iaero_aurora
         endif
-        if (iaero_generic(1) .ne. 0) then 
+        if (aerogeneric .ne. 0) then 
           print*,"iaero_generic= ",iaero_generic(:)
         endif
@@ -658 +661 @@
 !    Lucas Teinturier (2022)
 !===========================================================================
-      if (iaero_generic(1) .ne. 0) then ! we enter the scheme
+      if (aerogeneric .ne. 0) then ! we enter the scheme
         do ia=1,aerogeneric
           iaer = iaero_generic(ia)
@@ -674 +677 @@
           enddo !ig=1,ngrid
         enddo !ia=1,aerogeneric
-      endif !iaergo_generic(1) .ne. 0
+      endif !aerogeneric .ne. 0
 
 ! --------------------------------------------------------------------------

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

@@ -8 +8 @@
       logical,save :: season,diurnal,tlocked,rings_shadow,lwrite
 !$OMP THREADPRIVATE(season,diurnal,tlocked,rings_shadow,lwrite)
-      logical,save :: callstats,calleofdump
-!$OMP THREADPRIVATE(callstats,calleofdump)
+      logical,save :: callstats ! .true. to output stats.nc files
+!$OMP THREADPRIVATE(callstats)
       logical,save :: callgasvis,continuum,graybody
 !$OMP THREADPRIVATE(callgasvis,continuum,graybody)
@@ -50 +50 @@
       logical,save :: aerofixco2,aerofixh2o
 !$OMP THREADPRIVATE(aerofixco2,aerofixh2o)
-      integer,save :: aerogeneric
+      integer,save :: aerogeneric ! number of aerosols of "generic" kind
 !$OMP THREADPRIVATE(aerogeneric)
       logical,save :: hydrology