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

Timestamp:

Nov 23, 2022, 4:09:32 PM (2 years ago)

Author:

abierjon

Message:

Mars GCM:
The program util/aeroptical.F90 can now compute column-integrated optical depths
of the aerosols, in addition to the opacity profiles. The user has to specify it
('yes' or 'no') in aeroptical.def
Detailed changes :

• update of the aeroptical.def file to ask the user if the column optical depth should be computed (non-retrocompatible change)
• add in the init2 subroutine a computation of the layers' height delta_z, with adaptable method depending on the availability of some variables in the input file and on wether or not the input file has been zrecasted before. Preliminary validation shows that these different methods yield a +/-3% error on the output column optical depth tau_[aer]. The delta_z variable is also written in the output file.
• add a log file for warnings in the interpolation subroutine in aeropt_mod.F90 + add some comments in the code
• add the ouput "zzlev" (= interlayer altitudes) in libf/phymars/physiq_mod.F so that it can be directly used by aeroptical.F90 to compute delta_z

AB

File:

• trunk/LMDZ.MARS/util/aeroptical.def (modified) (2 diffs)

trunk/LMDZ.MARS/util/aeroptical.def

@@ -3 +3 @@
 5.e-7
 ext
+yes
 dust,dustq,reffdust,2500,optprop_dustvis_TM_n50.dat
 h2o_ice,h2o_ice,5.e-6,920,optprop_icevis_n30.dat
@@ -14 +15 @@
 3) The wavelength at which the output opacities are calibrated (in meters)
 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
+5) Computation of the column-integrated optical depth? (yes/no)
+
 6) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
+7) 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