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

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

@@ -314 +314 @@
         write(*,*)" coeff_detrainment = ",coeff_detrainment
 
-! entrainment by slope wind scheme
-         write(*,*)"call slope wind lifting parametrization"
-         slpwind=.false. ! default value
-         call getin_p("slpwind",slpwind)
-         write(*,*)" slpwind = ",slpwind
+! entrainment by mountain top dust flows scheme
+         write(*,*)"call mountain top dust flows parametrization"
+         topflows=.false. ! default value
+         call getin_p("topflows",topflows)
+         write(*,*)" topflows = ",topflows
 
 ! latent heat release from ground water ice sublimation/condensation
@@ -471 +471 @@
            endif
          endif
-! rocket dust storm and entrainment by slope wind
+! rocket dust storm and entrainment by top flows
 ! Test of incompatibility:
-! if rdstorm or slpwind is used, then doubleq should be true
-         if ((rdstorm.or.slpwind).and..not.doubleq) then
-           print*,'if rdstorm or slpwind is used, then doubleq 
+! if rdstorm or topflows is used, then doubleq should be true
+         if ((rdstorm.or.topflows).and..not.doubleq) then
+           print*,'if rdstorm or topflows is used, then doubleq 
      &            should be used !'
            call abort_physic(modname,
-     &          "rdstorm or slpwind requires doubleq",1)
-         endif
-         if ((rdstorm.or.slpwind).and..not.active) then
-           print*,'if rdstorm or slpwind is used, then active 
+     &          "rdstorm or topflows requires doubleq",1)
+         endif
+         if ((rdstorm.or.topflows).and..not.active) then
+           print*,'if rdstorm or topflows is used, then active 
      &            should be used !'
            call abort_physic(modname,
-     &          "rdstorm or slpwind requires activ",1)
+     &          "rdstorm or topflows requires activ",1)
          endif
          if (rdstorm.and..not.lifting) then
@@ -492 +492 @@
      &          "rdstorm requires lifting",1)
          endif
-         if ((rdstorm.or.slpwind).and..not.freedust) then
-           print*,'if rdstorm or slpwind is used, then freedust 
+         if ((rdstorm.or.topflows).and..not.freedust) then
+           print*,'if rdstorm or topflows is used, then freedust 
      &            should be used !'
            call abort_physic(modname,
-     &          "rdstorm or slpwind requires freedust",1)
+     &          "rdstorm or topflows requires freedust",1)
          endif
          if (rdstorm.and.(dustinjection.eq.0)) then
@@ -900 +900 @@
            ! and picky compilers who know name_iaer(2) is out of bounds
            j=2
-           IF (rdstorm.AND..NOT.activice.AND..NOT.slpwind) then
+           IF (rdstorm.AND..NOT.activice.AND..NOT.topflows) then
              name_iaer(j) = "stormdust_doubleq" !! storm dust two-moment scheme
              j = j+1
            END IF
 
-           IF (rdstorm.AND.water.AND.activice.AND..NOT.slpwind) then
+           IF (rdstorm.AND.water.AND.activice.AND..NOT.topflows) then
              name_iaer(j) = "stormdust_doubleq" 
              j = j+1
            END IF
 
-           IF (slpwind.AND..NOT.activice.AND..NOT.rdstorm) then
+           IF (topflows.AND..NOT.activice.AND..NOT.rdstorm) then
              name_iaer(j) = "topdust_doubleq" !! storm dust two-moment scheme
              j = j+1
            END IF
  
-           IF (slpwind.AND.water.AND.activice.AND..NOT.rdstorm) then
+           IF (topflows.AND.water.AND.activice.AND..NOT.rdstorm) then
              name_iaer(j) =  "topdust_doubleq"
              j = j+1
            END IF
 
-           IF (rdstorm.AND.slpwind.AND..NOT.activice) THEN 
+           IF (rdstorm.AND.topflows.AND..NOT.activice) THEN 
              name_iaer(j) = "stormdust_doubleq"
              name_iaer(j+1) = "topdust_doubleq"
@@ -926 +926 @@
            ENDIF
 
-           IF (rdstorm.AND.slpwind.AND.water.AND.activice) THEN 
+           IF (rdstorm.AND.topflows.AND.water.AND.activice) THEN 
              name_iaer(j) = "stormdust_doubleq"
              name_iaer(j+1) = "topdust_doubleq"