Changeset 2297 for trunk/LMDZ.GENERIC/README

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/README (modified) (2 diffs)

trunk/LMDZ.GENERIC/README

@@ -1512 +1512 @@
   Everybody should be using AU, that's all.
 
-== 24/02/2020 (r2244) == JVO
+== 24/02/2020 (r2245) == JVO
 + Add a 'versH2H2cia' int key in callphys that allows two values (2011 or 2018) to
   deal with updated HITRAN file (for interpolateH2H2.F90) from 2018 that includes the
@@ -1533 +1533 @@
   + Add a sanity check in callcorrk instead of leaving out-of-bounds planckir indexes.
 
-
+== 26/04/2020 (r2297) == JVO
+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..