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

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

@@ -181 +181 @@
         enddo
       endif ! of if (rdstorm)
-       if (slpwind) then
+       if (topflows) then
         do iq=1,nq
           if (noms(iq).eq."topdust_mass") then
@@ -192 +192 @@
           endif
         enddo
-      endif ! of if (slpwind)    
+      endif ! of if (topflows)    
       ! 2. find chemistry and water tracers
       do iq=1,nq
@@ -613 +613 @@
         end if !(rdstorm)
 !c ----------------------------------------------------------------------
-!c slope wind scheme
+!c mountain top dust flows scheme
 !c you need a radius value for topdust to active its sedimentation
 !c we take the same value as for the normal dust
-        if (slpwind) then
+        if (topflows) then
           rho_q(igcm_topdust_mass)=rho_dust
           rho_q(igcm_topdust_number)=rho_dust
           radius(igcm_topdust_mass) = 3.e-6
           radius(igcm_topdust_number) = 3.e-6
-        end if !(slpwind)
+        end if !(topflows)
 !c ----------------------------------------------------------------------
       
@@ -859 +859 @@
        endif
 
-       if (slpwind) then 
+       if (topflows) then 
        ! verify that we indeed have topdust_mass and topdust_number tracers 
          if (igcm_topdust_mass.eq.0) then
            write(*,*) "initracer: error !!"
-           write(*,*) "  cannot use slpwind option without ",
+           write(*,*) "  cannot use topflows option without ",
      &                "a topdust_mass tracer !"
-           call abort_physic("initracer","slpwind issue",1)
+           call abort_physic("initracer","topflows issue",1)
          endif
          if (igcm_topdust_number.eq.0) then
            write(*,*) "initracer: error !!"
-           write(*,*) "  cannot use slpwind option without ",
+           write(*,*) "  cannot use topflows option without ",
      &                "a topdust_number tracer !"
-           call abort_physic("initracer","slpwind issue",1)
+           call abort_physic("initracer","topflows issue",1)
          endif
        endif