Changeset 2830 for trunk/LMDZ.MARS/libf

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/libf/phymars/physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -2985 +2985 @@
 c        call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",3,zplay)
 c        call WRITEDIAGFI(ngrid,"pplev","Pressure","Pa",3,zplev)
+         call WRITEDIAGFI(ngrid,"zzlev","Interlayer altitude",
+     &                    "m",3,zzlev(:,1:nlayer))
          call writediagfi(ngrid,"pphi","Geopotential","m2s-2",3,
      &                    pphi)