Changeset 2628 for trunk/LMDZ.MARS/libf/phymars/callsedim_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/callsedim_mod.F (modified) (5 diffs)

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

@@ -372 +372 @@
        ENDIF !of if (rdstorm)
 
-       IF (slpwind) THEN ! identifying topdust tracers for sedimentation
+       IF (topflows) THEN ! identifying topdust tracers for sedimentation
            itopdust_mass=0      ! dummy initialization
            itopdust_number=0    ! dummy initialization
@@ -392 +392 @@
              write(*,*) 'callsedim: error! could not identify'
              write(*,*) ' tracers for topdust mass and number mixing'
-             write(*,*) ' ratio and slpwind is activated!'
+             write(*,*) ' ratio and topflows is activated!'
              call abort_physic(modname,"missing topdust tracers",1)
            endif
-       ENDIF !of if (slpwind)
+       ENDIF !of if (topflows)
 
         firstcall=.false.
@@ -452 +452 @@
         end do
       endif
-c     entrainment by slope wind
-      if (slpwind) then
+c     entrainment by mountain top dust flows
+      if (topflows) then
         do l=1,nlay
           do ig=1, ngrid
@@ -710 +710 @@
       endif ! of if (rdstorm)
 
-      if (slpwind) then
+      if (topflows) then
        DO l = 1, nlay
         DO ig=1,ngrid
@@ -718 +718 @@
         ENDDO
        ENDDO
-      endif ! of if (slpwind)
+      endif ! of if (topflows)
  
 c     Update the ice particle size "rice"