Changeset 2803 for trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90

Timestamp:

Oct 19, 2022, 11:24:43 AM (2 years ago)

Author:

aslmd

Message:

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 file

File:

• trunk/LMDZ.GENERIC/libf/phystd/aeropacity.F90 (modified) (7 diffs)

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

@@ -4 +4 @@
        use radinc_h, only : L_TAUMAX,naerkind
        use aerosol_mod
-       USE tracer_h, only: noms,rho_co2,rho_ice
-       use comcstfi_mod, only: g, pi
+       USE tracer_h, only: noms,rho_co2,rho_ice,rho_q,mmol
+       use comcstfi_mod, only: g, pi, mugaz, avocado
        use geometry_mod, only: latitude
        use callkeys_mod, only: aerofixco2,aerofixh2o,kastprof,cloudlvl, &
@@ -13 +13 @@
                 tau_nh3_cloud, pres_nh3_cloud,                          &
                 nlayaero, aeronlay_tauref, aeronlay_choice,             & 
-                aeronlay_pbot, aeronlay_ptop, aeronlay_sclhght
-                  
+                aeronlay_pbot, aeronlay_ptop, aeronlay_sclhght,         &
+                aerogeneric
+        use generic_tracer_index_mod, only: generic_tracer_index
        implicit none
 
@@ -30 +31 @@
 !     dust removal, simplification by Robin Wordsworth (2009)
 !     Generic n-layer aerosol - J. Vatant d'Ollone (2020)
+!     Radiative Generic Condensable aerosols - Lucas Teinturier (2022)
 !
 !     Input
@@ -90 +92 @@
 
       real CLFtot
-
+      integer igen_ice,igen_vap ! to store the index of generic tracer
+      logical dummy_bool ! dummy boolean just in case we need one
       ! identify tracers
       IF (firstcall) THEN
@@ -141 +144 @@
           print*,'iaero_aurora= ',iaero_aurora
         endif
-
+        if (iaero_generic(1) .ne. 0) then 
+          print*,"iaero_generic= ",iaero_generic(:)
+        endif
         firstcall=.false.
       ENDIF ! of IF (firstcall)
@@ -362 +367 @@
             ENDDO
          ENDDO
-
+         
 ! 1/700. is assuming a "sulfurtau" of 1
 ! Sulfur aerosol routine to be improved.
@@ -645 +650 @@
 
       end if ! if Auroral aerosols  
-
+!===========================================================================
+!    Radiative Generic Condensable aerosols scheme
+!    Only used when we give aerogeneric != 0 in callphys.def
+!    Computes the generic aerosols' opacity in the same fashion as water of 
+!    dust, using the QREFvis3d of the concerned specie
+!    Lucas Teinturier (2022)
+!===========================================================================
+      if (iaero_generic(1) .ne. 0) then ! we enter the scheme
+        do ia=1,aerogeneric
+          iaer = iaero_generic(ia)
+          ! Initialization
+          aerosol(1:ngrid,1:nlayer,iaer) = 0.D0
+          igen_ice = -1
+          igen_vap = -1
+          ! Get index of the ice tracer (for pq, mol, rho_q)
+          call generic_tracer_index(nq,2*iaer,igen_vap,igen_ice,dummy_bool) 
+          if (igen_ice .eq. -1) then 
+            !if igen_ice not changed by generic_tracer_index, it means it's an ice tracer => igen_ice=2*iaer
+            igen_ice = 2*iaer 
+          endif
+          ! Let's loop on the horizontal and vertical grid
+          do ig=1,ngrid
+            do l=1,nlayer
+              aerosol(ig,l,iaer) = ( 0.75*QREFvis3d(ig,l,iaer)  / &
+                                  (rho_q(igen_ice) * reffrad(ig,l,iaer)) ) * &
+                                  (pq(ig,l,igen_ice)+1E-9 ) *                &
+                                  (pplev(ig,l) - pplev(ig,l+1)) /g 
+            enddo !l=1,nlayer
+          enddo !ig=1,ngrid
+        enddo !ia=1,aerogeneric
+      endif !iaergo_generic(1) .ne. 0
 
 ! --------------------------------------------------------------------------