Changeset 2199 for trunk/LMDZ.MARS/libf/phymars/conf_phys.F

Timestamp:

Dec 13, 2019, 2:06:35 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Implementation of a new parametrization of the dust entrainment by slope winds above the sub-grid scale topography. The parametrization is activated with the flag slpwind=.true. (set to "false" by
default) in callphys.def. The new parametrization involves the new tracers topdust_mass and topdust_number.
MV

File:

• trunk/LMDZ.MARS/libf/phymars/conf_phys.F (modified) (3 diffs)

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

@@ -298 +298 @@
         write(*,*)" coeff_detrainment = ",coeff_detrainment
 
+! entrainment by slope wind scheme
+         write(*,*)"call slope wind lifting parametrization"
+         slpwind=.false. ! default value
+         call getin("slpwind",slpwind)
+         write(*,*)" slpwind = ",slpwind
+
 ! latent heat release from ground water ice sublimation/condensation
          write(*,*)"latent heat release during sublimation", 
@@ -430 +436 @@
            endif
          endif
-! rocket dust storm
+! rocket dust storm and entrainment by slope wind
 ! Test of incompatibility:
-! if rdstorm is used, then doubleq should be true
-         if (rdstorm.and..not.doubleq) then
-           print*,'if rdstorm is used, then doubleq should be used !'
-           stop
-         endif
-         if (rdstorm.and..not.active) then
-           print*,'if rdstorm is used, then active should be used !'
+! 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 
+     &            should be used !'
+           stop
+         endif
+         if ((rdstorm.or.slpwind).and..not.active) then
+           print*,'if rdstorm or slpwind is used, then active 
+     &            should be used !'
            stop
          endif
          if (rdstorm.and..not.lifting) then
-           print*,'if rdstorm is used, then lifting should be used !'
-           stop
-         endif
-         if (rdstorm.and..not.freedust) then
-           print*,'if rdstorm is used, then freedust should be used !'
+           print*,'if rdstorm is used, then lifting 
+     &            should be used !'
+           stop
+         endif
+         if ((rdstorm.or.slpwind).and..not.freedust) then
+           print*,'if rdstorm or slpwind is used, then freedust 
+     &            should be used !'
            stop
          endif
          if (rdstorm.and.(dustinjection.eq.0)) then
-           print*,'if rdstorm is used, then dustinjection should 
-     &             be used !'
+           print*,'if rdstorm is used, then dustinjection
+     &            should be used !'
            stop
          endif
@@ -765 +775 @@
           ! and picky compilers who know name_iaer(2) is out of bounds
           j=2
-         IF (rdstorm.AND..NOT.activice) name_iaer(2) = 
+         IF (rdstorm.AND..NOT.activice.AND..NOT.slpwind) name_iaer(2) = 
      &                       "stormdust_doubleq" !! storm dust two-moment scheme
-         IF (rdstorm.AND.water.AND.activice) name_iaer(3) = 
-     &                                             "stormdust_doubleq" 
+         IF (rdstorm.AND.water.AND.activice.AND..NOT.slpwind)
+     &              name_iaer(3) = "stormdust_doubleq" 
+         IF (slpwind.AND..NOT.activice.AND..NOT.rdstorm) name_iaer(2) = 
+     &                       "topdust_doubleq" !! storm dust two-moment scheme
+         IF (slpwind.AND.water.AND.activice.AND..NOT.rdstorm)
+     &              name_iaer(3) =  "topdust_doubleq"
+         IF (rdstorm.AND.slpwind.AND..NOT.activice) THEN 
+             name_iaer(2) = "stormdust_doubleq"
+             name_iaer(3) = "topdust_doubleq"
+         ENDIF
+         IF (rdstorm.AND.slpwind.AND.water.AND.activice) THEN 
+             name_iaer(3) = "stormdust_doubleq"
+             name_iaer(4) = "topdust_doubleq"
+         ENDIF
          IF (water.AND.activice) name_iaer(j) = "h2o_ice"      !! radiatively-active clouds
          IF (submicron.AND.active) name_iaer(j) = "dust_submicron" !! JBM experimental stuff