Changeset 2248

Timestamp:

Mar 4, 2020, 6:31:04 PM (5 years ago)

Author:

abierjon

Message:

Mars GCM:
Bug fixes in the slope winds parametrization and its interactions with rocketduststorm : now, it is considered that the sub-grid dust storm is decorrelated from the dust entrainment above the sub-grid scale topography, so that : 1) the background dust available for slope winds is updated after the rocketduststorm parametrization ; 2) the extra-heating on top of the mountains is calculated regarding topdust and background dust tracers' optical depths only (and not stormdust optical depth).
AB

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/physiq_mod.F (modified) (1 diff)
• libf/phymars/topmons_mod.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -2866 +2866 @@
 + some cleaning and commenting of the code
 
-== 02/03/2020 == AB
-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).
+== 04/03/2020 == AB
+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).
+
+== 04/03/2020 == AB
+Bug fixes in the slope winds parametrization and its interactions with rocketduststorm : now, it is considered that the sub-grid dust storm is decorrelated from the dust entrainment above the sub-grid scale topography, so that : 1) the background dust available for slope winds is updated after the rocketduststorm parametrization ; 2) the extra-heating on top of the mountains is calculated regarding topdust and background dust tracers' optical depths only (and not stormdust optical depth).

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

@@ -1124 +1124 @@
            hsummit(:)=14000.
          endif
+         clearatm=.true.
          nohmons=.false.
          pdqtop(:,:,:)=0.  

trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90

@@ -257 +257 @@
 
       !! Update the initial mixing ratios
-      zq0(1:ngrid,1:nlayer,1:nq)=pq(1:ngrid,1:nlayer,1:nq) !+pdq(1:ngrid,1:nlayer,1:nq)*ptimestep update after rocket dust storm scheme ?
+      zq0(1:ngrid,1:nlayer,1:nq)=pq(1:ngrid,1:nlayer,1:nq)+pdq(1:ngrid,1:nlayer,1:nq)*ptimestep ! update of the background dust after rocket dust storm scheme
       zq(1:ngrid,1:nlayer,1:nq)=zq0(1:ngrid,1:nlayer,1:nq)
       zq_dust_mass(1:ngrid,1:nlayer)=zq(1:ngrid,1:nlayer,igcm_dust_mass)