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

Last change on this file since 3493 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
RevLine 
[2317]1diagfi.nc
[2817]2~/datadir
[2317]35.e-7
[2443]4ext
[2830]5yes
[2817]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
[2317]10
11
12##----------------------------------------------------------
131) Name of the GCM input file
[2817]142) Directory where the optical properties files are
[2317]153) The wavelength at which the output opacities are calibrated (in meters)
[2817]164) Opacity type : extinction (ext) or absorption (abs)
[2830]175) Computation of the column-integrated optical depth? (yes/no)
18
[2817]196) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
[2830]207) aerosol_name,aerosol_mmr(variable name),aerosol_reff(variable name or single value in m),aerosol_rho(value in kg/m3),optpropfile_name
[2817]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
[2317]23
[2817]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.