Changeset 2628 for trunk/LMDZ.MARS/libf/phymars/callradite_mod.F

Timestamp:

Feb 28, 2022, 6:46:07 PM (3 years ago)

Author:

abierjon

Message:

Mars GCM:
Big changes on mountain top dust flows for GCM6:

• the scheme now activates only in grid meshes that contain a mountain among a hard-written list, instead of every meshes. This is done to prevent strong artificial reinjections of dust in places that don't present huge converging slopes enabling the concentration of dust (ex: Valles Marineris, Hellas). Topdust is now also detrained as soon as it leaves the column it originated from.
• the list of the 19 allowed mountains is used by the subroutine topmons_setup in module topmons_mod, to compute a logical array contains_mons(ngrid). alpha_hmons and hsummit are also set up once and for all by this subroutine, which is called in physiq_mod's firstcall.
• contains_mons, alpha_hmons and hsummit are now saved variables of the module surfdat_h, and are called as such and not as arguments in the subroutines using them.
• the logical flag "slpwind" and the comments in the code have also been updated to the new terminology "mountain top dust flows", accordingly to ticket #71. The new flag read in callphys.def is "topflows".

AB

File:

• trunk/LMDZ.MARS/libf/phymars/callradite_mod.F (modified) (4 diffs)

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

@@ -12 +12 @@
      $     taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,
      $     rstormdust,rtopdust,totstormfract,clearatm,dsords,dsotop,
-     $     alpha_hmons,nohmons,clearsky,totcloudfrac)
+     $     nohmons,clearsky,totcloudfrac)
 
       use aeropacity_mod, only: aeropacity
@@ -213 +213 @@
       REAL,INTENT(OUT) :: dsords(ngrid,nlayer) ! density scaled opacity for rocket dust storm dust
       
-c     entrainment by slope wind
-      LOGICAL, INTENT(IN) :: nohmons ! true for background dust 
-      REAL, INTENT(IN) :: alpha_hmons(ngrid) ! sub-grid scale topography mesh fraction
+c     entrainment by mountain top dust flows
+      LOGICAL, INTENT(IN) :: nohmons ! true for background dust
       REAL,INTENT(OUT) :: rtopdust(ngrid,nlayer)  ! Topdust geometric mean radius (m)
       REAL,INTENT(OUT) :: dsotop(ngrid,nlayer) ! density scaled opacity for topmons dust
@@ -376 +375 @@
            enddo
          end if
-         if (slpwind.AND.active) then
+         if (topflows.AND.active) then
            do iaer=1,naerkind
              if (name_iaer(iaer).eq."topdust_doubleq") then
@@ -448 +447 @@
      &    QREFvis3d,QREFir3d,omegaREFir3d,
      &    totstormfract,clearatm,dsords,dsotop,
-     &    alpha_hmons,nohmons,
-     &    clearsky,totcloudfrac)
+     &    nohmons,clearsky,totcloudfrac)
+     
 c     Starting loop on sub-domain
 c     ----------------------------