Changeset 2297 for trunk/LMDZ.GENERIC/libf/phystd/callkeys_mod.F90

Timestamp:

Apr 26, 2020, 7:45:46 PM (5 years ago)

Author:

jvatant

Message:

Add a generic n-layer aerosol scheme to replace the former buggy 2-layer scheme as well as the hard-coded NH3 cloud.

It can be called using 'aeronlay=.true.' in callphys.def, and set the number of layers (up to 4) with 'nlayaero'.
Then, the following parameters are read as arrays of size nlayaero in callphys.def (separated by blank space)

*aeronlay_tauref (Optical depth of aerosol layer at ref wavelenght)
*aeronlay_lamref (Ref wavelenght (m))
*aeronlay_choice (Choice for vertical profile - 1:tau follows atm scale height btwn top and bottom - 2:tau follows it own scale height)
*aeronlay_pbot (Bottom pressure (Pa))
*aeronlay_ptop (Top pressure (Pa) - useful only if choice=1)
*aeronlay_sclhght (Ratio of aerosol layer scale height / atmospheric scale height - useful only if choice=2 )
*aeronlay_size (Particle size (m))
*optprop_aeronlay_vis File for VIS opt properties.
*optprop_aeronlay_ir File for IR opt properties.

+Extra info :

+ In addition of solving the bug from 2-layer it enables different optical properties.
+ The former scheme are left for retrocompatibility (for now) but you should use the new one.
+ See aeropacity.F90 for the calculations

+ Each layer can have different optical properties, size of particle ...
+ You have different choices for vertical profile of the aerosol layers :

• aeronlay_choice = 1 : Layer tau is spread between ptop and pbot following atm scale height.
• aeronlay_choice = 2 : Layer tau follows its own scale height above cloud deck (pbot).

In this case ptop is dummy and sclhght gives the ratio H_cl/H_atm.

+ The reference wavelenght for input optical depth is now read as input (aeronlay_lamref)
+ Following the last point some comment is added in suaer_corrk about the 'not-really-dummy'ness of IR lamref..

File:

• trunk/LMDZ.GENERIC/libf/phystd/callkeys_mod.F90 (modified) (4 diffs)

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

@@ -43 +43 @@
       logical,save :: aeroco2,aeroh2o,aeroh2so4,aeroback2lay
 !$OMP THREADPRIVATE(aeroco2,aeroh2o,aeroh2so4,aeroback2lay)
-      logical,save :: aeronh3, aeroaurora
-!$OMP THREADPRIVATE(aeronh3,aeroaurora)
+      logical,save :: aeronh3, aeronlay, aeroaurora
+!$OMP THREADPRIVATE(aeronh3,aeronlay,aeroaurora)
       logical,save :: aerofixco2,aerofixh2o
 !$OMP THREADPRIVATE(aerofixco2,aerofixh2o)
@@ -62 +62 @@
 !$OMP THREADPRIVATE(photochem)
 
-      integer,save :: versH2H2cia
       integer,save :: iddist
       integer,save :: iaervar
       integer,save :: iradia
       integer,save :: startype
-!$OMP THREADPRIVATE(versH2H2cia,iddist,iaervar,iradia,startype)
+      integer,save :: versH2H2cia
+      integer,save :: nlayaero
+!$OMP THREADPRIVATE(iddist,iaervar,iradia,startype,versH2H2cia,nlayaero)
+      integer,dimension(:),allocatable,save :: aeronlay_choice
+!$OMP THREADPRIVATE(aeronlay_choice)
+
+      character(64),save :: optprop_back2lay_vis
+      character(64),save :: optprop_back2lay_ir
+      character(64),dimension(:),allocatable,save :: optprop_aeronlay_vis
+      character(64),dimension(:),allocatable,save :: optprop_aeronlay_ir
+!$OMP THREADPRIVATE(optprop_back2lay_vis,optprop_back2lay_ir,optprop_aeronlay_vis,optprop_aeronlay_ir)
 
       real,save :: tplanckmin
@@ -84 +93 @@
       real,save :: obs_tau_col_strato
 !$OMP THREADPRIVATE(Tstrat,tplanet,obs_tau_col_tropo,obs_tau_col_strato)
-      character(64),save :: optprop_back2lay_vis
-      character(64),save :: optprop_back2lay_ir
       real,save :: pres_bottom_tropo
       real,save :: pres_top_tropo
@@ -101 +108 @@
       real,save :: tau_nh3_cloud
 !$OMP THREADPRIVATE(size_nh3_cloud, pres_nh3_cloud, tau_nh3_cloud)
+      real,dimension(:),allocatable,save :: aeronlay_tauref
+      real,dimension(:),allocatable,save :: aeronlay_lamref
+      real,dimension(:),allocatable,save :: aeronlay_ptop
+      real,dimension(:),allocatable,save :: aeronlay_pbot
+      real,dimension(:),allocatable,save :: aeronlay_sclhght
+      real,dimension(:),allocatable,save :: aeronlay_size
+!$OMP THREADPRIVATE(aeronlay_tauref,aeronlay_lamref,aeronlay_ptop,aeronlay_pbot,aeronlay_sclhght,aeronlay_size)
       real,save :: co2supsat
       real,save :: pceil