Changeset 2247 for trunk/LMDZ.MARS/libf

Timestamp:

Mar 4, 2020, 12:35:33 PM (5 years ago)

Author:

abierjon

Message:

Mars GCM:
Partially resolved ticket #14 : when calling writediagfi for all kinds of dust Density-Scaled Opacity (dsodust,dsords,dsotop,dso), an information about the IR wavelength used by the GCM is added in the "title" attribute of the netcdf variable ("tes" = 9.3micron ; "mcs" = 21.6micron).
NB : these dust DSO calculated by the GCM are extinction opacities, hence not always corresponding to the actual opacity/optical depth measured by the instruments (ex: TES = 9.3micron absorption optical depth).
AB

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (2 diffs)

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

@@ -2924 +2924 @@
              call WRITEDIAGFI(ngrid,'dustN','Dust number',
      &                        'part/kg',3,ndust)
-             call WRITEDIAGFI(ngrid,'dsodust',
-     &                        'density scaled optical depth in the IR',
-     &                        'm2.kg-1',3,dsodust)
-             call WRITEDIAGFI(ngrid,'dso',
-     &                        'density scaled optical depth in the IR',
-     &                        'm2.kg-1',3,dsodust+dsords+dsotop)
-             
+             
+             select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsodust',
+     &  'density scaled extinction opacity of std dust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsodust)
+               call WRITEDIAGFI(ngrid,'dso',
+     &  'density scaled extinction opacity of all dust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsodust+dsords+dsotop)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsodust',
+     &  'density scaled extinction opacity of std dust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsodust)
+               call WRITEDIAGFI(ngrid,'dso',
+     &  'density scaled extinction opacity of all dust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsodust+dsords+dsotop)
+             end select
            else ! (doubleq=.false.)
              do iq=1,dustbin
@@ -2987 +2997 @@
              call WRITEDIAGFI(ngrid,'totaldustq','total dust mass',
      &                        'kg/kg',3,qdusttotal)
-             call WRITEDIAGFI(ngrid,'dsords',
-     &                       'density scaled opacity of stormdust',
-     &                       'm2.kg-1',3,dsords)
+             
+             select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsords',
+     &  'density scaled extinction opacity of stormdust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsords)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsords',
+     &  'density scaled extinction opacity of stormdust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsords)
+             end select
            endif ! (rdstorm)
 
            if (slpwind) then
-             call WRITEDIAGFI(ngrid,'dsotop',
-     &                       'density scaled opacity of topdust',
-     &                       'm2.kg-1',3,dsotop)
+             select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsotop',
+     &  'density scaled extinction opacity of topdust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsotop)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsotop',
+     &  'density scaled extinction opacity of topdust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsotop)
+             end select
            endif ! (slpwind)