source: trunk/LMDZ.MARS/util/aeroptical.def @ 3026

Last change on this file since 3026 was 2830, checked in by abierjon, 2 years ago

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 size: 1.3 KB
Line 
1diagfi.nc
2~/datadir
35.e-7
4ext
5yes
6dust,dustq,reffdust,2500,optprop_dustvis_TM_n50.dat
7h2o_ice,h2o_ice,5.e-6,920,optprop_icevis_n30.dat
8stormdust,rdsdustq,reffstormdust,2500,optprop_dustvis_TM_n50.dat
9topdust,topdustq,refftopdust,2500,optprop_dustvis_TM_n50.dat
10
11
12##----------------------------------------------------------
131) Name of the GCM input file
142) Directory where the optical properties files are
153) The wavelength at which the output opacities are calibrated (in meters)
164) Opacity type : extinction (ext) or absorption (abs)
175) Computation of the column-integrated optical depth? (yes/no)
18
196) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
207) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
21...
22N) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
23
24##----------------------------------------------------------
25USUAL GCM OPTPROPFILES :
26- for water ice : "optprop_icevis_n30.dat" for visible (wvl<=5.e-6) ; "optprop_iceir_n30.dat" for infrared (wvl>5.e-6)
27- 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)
Note: See TracBrowser for help on using the repository browser.