Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/util/README

    r2817 r2830  
    192192--------------------------------------------------------------------
    193193
    194 Program to compute opacities [1/km] of aerosols at a wavelength given by the user.
     194Program to compute opacity profiles opa_[aer] (1/km) of aerosols
     195at a wavelength given by the user.
    195196Computation is made from the aerosol mass mixing ratios and effective radius,
    196197associated with air density (rho) and files containing the aerosol
    197198optical properties (generally present in the GCM datadir/ directory).
    198199The user have to precise the type of opacity : extinction or absorption.
     200The user can also ask to compute the column-integrated optical depth tau_[aer].
    199201
    200202NB : this program requires to compile the module aeropt_mod.F90
Note: See TracChangeset for help on using the changeset viewer.