Changeset 2817 for trunk/LMDZ.MARS/util/aeroptical.def

Timestamp:

Nov 14, 2022, 6:39:11 PM (2 years ago)

Author:

abierjon

Message:

Mars GCM:
Non-retrocompatible changes of the program util/aeroptical.F90 :

• make the program more generic, as the user now specify in aeroptical.def the aerosols to be computed, as well as the necessary information for each aerosol, that are : aerosol_name,aerosol_mmr(GCM variable name),aerosol_reff(GCM variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
• creation of a separate module aeropt_mod.F90 that can be used by other programs (e.g. for the MCD), with one subroutine reading the ASCII optical properties file, one subroutine interpolating optical properties at a given (wvl,reff) tuple, and one subroutine deallocating the module variables.
• modify the compile program to compile this module aeropt_mod.F90 along with the aeroptical.F90 program. There is no change for the user, which just has to run the command 'compile aeroptical' like before.
• correct the handling of user's wavelength or particle sizes that do not lie within the boundaries of the optical properties file.

AB

File:

• trunk/LMDZ.MARS/util/aeroptical.def (modified) (1 diff)

trunk/LMDZ.MARS/util/aeroptical.def

@@ -1 +1 @@
 diagfi.nc
-0
+~/datadir
 5.e-7
-~/datadir
 ext
-
+dust,dustq,reffdust,2500,optprop_dustvis_TM_n50.dat
+h2o_ice,h2o_ice,5.e-6,920,optprop_icevis_n30.dat
+stormdust,rdsdustq,reffstormdust,2500,optprop_dustvis_TM_n50.dat
+topdust,topdustq,refftopdust,2500,optprop_dustvis_TM_n50.dat
 
 
 ##----------------------------------------------------------
 1) Name of the GCM input file
-2) Value of the water ice effective radius (in meters)
-  (if 0, reffwice is read in the input file ;
-   usually between 3.e-6m and 5.e-6m )
+2) Directory where the optical properties files are
 3) The wavelength at which the output opacities are calibrated (in meters)
-4) Directory where the optical properties files are
-5) Opacity type : extinction (ext) or absorption (abs)
+4) Opacity type : extinction (ext) or absorption (abs)
+5) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
+6) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
+...
+N) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
 
+##----------------------------------------------------------
+USUAL GCM OPTPROPFILES :
+- for water ice : "optprop_icevis_n30.dat" for visible (wvl<=5.e-6) ; "optprop_iceir_n30.dat" for infrared (wvl>5.e-6)
+- for dust (any type) : "optprop_dustvis_TM_n50.dat" for visible (wvl<=5.e-6) ; "optprop_dustir_n50.dat" for infrared (wvl>5.e-6)