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

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

@@ -23 +23 @@
       use rocketduststorm_mod, only: rocketduststorm, dustliftday
       use calcstormfract_mod, only: calcstormfract
-      use topmons_mod, only: topmons,alpha_hmons
+      use topmons_mod, only: topmons,topmons_setup
       use tracer_mod, only: noms, mmol, igcm_co2, igcm_n2, igcm_co2_ice,
      &                      igcm_co, igcm_o, igcm_h2o_vap, igcm_h2o_ice,
@@ -308 +308 @@
                                     !            - in a mesh with stormdust and background dust (false)
                                     !            - in a mesh with background dust only (true)
-c     entrainment by slope winds
+c     entrainment by mountain top dust flows
       logical nohmons               ! nohmons used to calculate twice the radiative 
-                                    ! transfer when slpwind is active : 
+                                    ! transfer when topflows is active : 
                                     !            - in a mesh with topdust and background dust (false)
                                     !            - in a mesh with background dust only (true)
@@ -507 +507 @@
       REAL co2totB
 
-c entrainment by slope winds above sb-grid scale topography
+c entrainment by mountain top dust flows above sub-grid scale topography
       REAL pdqtop(ngrid,nlayer,nq) ! tendency for dust after topmons
-      REAL hmax,hmin
-      REAL hsummit(ngrid)
 
 c when no startfi file is asked for init
@@ -718 +716 @@
          endif
 
-c        Initialize mountain mesh fraction for the entrainment by slope wind param.
+c        Initialize mountain mesh fraction for the entrainment by top flows param.
 c        ~~~~~~~~~~~~~~~
-        if (slpwind) then
-          !! alpha_hmons calculation
-          if (ngrid.gt.1) then
-            call planetwide_maxval(hmons,hmax )
-            call planetwide_minval(hmons,hmin )
-            do ig=1,ngrid
-              alpha_hmons(ig)= 0.5*(hmons(ig)-hmin)/(hmax-hmin)
-            enddo
-          else
-            hmin=0.
-            hmax=23162.1 !set here the height of the sub-grid scaled topography
-            do ig=1,ngrid               
-              alpha_hmons(ig)= (hmons(ig)-hmin)/(hmax-hmin) !0.1*(hmons(ig)-hmin)/(hmax-hmin)
-              print*,"1D, hmons=",hmons(ig),"alpha=",alpha_hmons(ig)
-            enddo
-          endif ! (ngrid.gt.1)
-        endif ! if (slpwind)
+        if (topflows) call topmons_setup(ngrid)
 
 #endif
@@ -972 +954 @@
 c          Transfer through CO2 (except NIR CO2 absorption)
 c            and aerosols (dust and water ice)
-           ! callradite for background dust
+           ! callradite for background dust (out of the rdstorm fraction)
            clearatm=.true.
-           !! callradite for background dust in the case of slope wind entrainment
+           !! callradite for background dust (out of the topflows fraction)
            nohmons=.true.
 
@@ -988 +970 @@
      &     taucloudtes,rdust,rice,nuice,riceco2,nuiceco2,co2ice,
      &     rstormdust,rtopdust,totstormfract,clearatm,dsords,dsotop,
-     &     alpha_hmons,nohmons,clearsky,totcloudfrac)
+     &     nohmons,clearsky,totcloudfrac)
 
            ! case of sub-grid water ice clouds: callradite for the clear case
@@ -1006 +988 @@
      &              rice,nuice,riceco2, nuiceco2,co2ice,rstormdust,
      &              rtopdust,totstormfract,
-     &              clearatm,dsords,dsotop,alpha_hmons,nohmons,
-     &              clearsky,totcloudfrac)
+     &              clearatm,dsords,dsotop,
+     &              nohmons,clearsky,totcloudfrac)
                clearsky = .false.  ! just in case.
                ! Sum the fluxes and heating rates from cloudy/clear
@@ -1200 +1182 @@
      &                      clearsky,totcloudfrac,
 c               input sub-grid scale topography
-     &                      nohmons,alpha_hmons,
+     &                      nohmons,
 c               output
      &                      pdqrds,wspeed,dsodust,dsords,dsotop,
@@ -1250 +1232 @@
 c     3.2 Dust entrained from the PBL up to the top of sub-grid scale topography
 c    -------------------------------------------
-      IF (slpwind) THEN
-         if (ngrid.gt.1) then
-           hsummit(:)=summit(:)-phisfi(:)/g
-         else
-           hsummit(:)=14000.
-         endif
+      IF (topflows) THEN
          clearatm=.true. ! stormdust is not accounted in the extra heating on top of the mountains
          nohmons=.false.
@@ -1266 +1243 @@
      &                totstormfract,clearatm,
      &                clearsky,totcloudfrac,
-     &                nohmons,hsummit,
+     &                nohmons,
      &                pdqtop,wtop,dsodust,dsords,dsotop,
      &                tau_pref_scenario,tau_pref_gcm)
@@ -1287 +1264 @@
          ENDDO
 
-      ENDIF ! end of if (slpwind)
+      ENDIF ! end of if (topflows)
 
 c     3.3 Dust injection from the surface
@@ -3387 +3364 @@
            endif ! (rdstorm)
 
-           if (slpwind) then
+           if (topflows) then
              call WRITEDIAGFI(ngrid,'refftopdust','refftopdust',
      &                        'm',3,rtopdust*ref_r0)
@@ -3404 +3381 @@
      &         'm2.kg-1',3,dsotop)
              end select
-           endif ! (slpwind)
+           endif ! (topflows)
 
            if (dustscaling_mode==2) then