Changeset 2628 for trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90

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/surfdat_h.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/surfdat_h.F90

@@ -35 +35 @@
 !$OMP                TESice_Scoef,iceradius,dtemisice,                           &
 !$OMP                zmea,zstd,zsig,zgam,zthe,hmons,summit,base,z0,z0_default )
+
+  !! mountain top dust flows
+  REAL,SAVE,ALLOCATABLE :: alpha_hmons(:) ! sub-grid scale mountain mesh fraction
+  REAL,SAVE,ALLOCATABLE :: hsummit(:) ! mountain height above the GCM surface
+  LOGICAL,SAVE,ALLOCATABLE :: contains_mons(:) ! is there a mountain in the grid mesh ?
+          
+!$OMP THREADPRIVATE(alpha_hmons,hsummit,contains_mons)
 
   !! variables
@@ -75 +82 @@
     allocate(summit(ngrid))
     allocate(base(ngrid))
-    
+    allocate(alpha_hmons(ngrid))
+    allocate(hsummit(ngrid))
+    allocate(contains_mons(ngrid))
+       
   end subroutine ini_surfdat_h
 
@@ -83 +93 @@
   implicit none
 
-    if (allocated(albedodat))   deallocate(albedodat)
-    if (allocated(phisfi))      deallocate(phisfi)
-    if (allocated(watercaptag)) deallocate(watercaptag)
-    if (allocated(dryness))     deallocate(dryness)
-    if (allocated(zmea))        deallocate(zmea)
-    if (allocated(zstd))        deallocate(zstd)
-    if (allocated(zsig))        deallocate(zsig)
-    if (allocated(zgam))        deallocate(zgam)
-    if (allocated(zthe))        deallocate(zthe)
-    if (allocated(z0))          deallocate(z0)
-    if (allocated(qsurf))       deallocate(qsurf)
-    if (allocated(tsurf))       deallocate(tsurf)
-    if (allocated(co2ice))      deallocate(co2ice)
-    if (allocated(watercap))    deallocate(watercap)
-    if (allocated(emis))        deallocate(emis)
-    if (allocated(capcal))      deallocate(capcal)
-    if (allocated(fluxgrd))     deallocate(fluxgrd)
-    if (allocated(hmons))       deallocate(hmons)
-    if (allocated(summit))      deallocate(summit)
-    if (allocated(base))        deallocate(base)
+    if (allocated(albedodat))     deallocate(albedodat)
+    if (allocated(phisfi))        deallocate(phisfi)
+    if (allocated(watercaptag))   deallocate(watercaptag)
+    if (allocated(dryness))       deallocate(dryness)
+    if (allocated(zmea))          deallocate(zmea)
+    if (allocated(zstd))          deallocate(zstd)
+    if (allocated(zsig))          deallocate(zsig)
+    if (allocated(zgam))          deallocate(zgam)
+    if (allocated(zthe))          deallocate(zthe)
+    if (allocated(z0))            deallocate(z0)
+    if (allocated(qsurf))         deallocate(qsurf)
+    if (allocated(tsurf))         deallocate(tsurf)
+    if (allocated(co2ice))        deallocate(co2ice)
+    if (allocated(watercap))      deallocate(watercap)
+    if (allocated(emis))          deallocate(emis)
+    if (allocated(capcal))        deallocate(capcal)
+    if (allocated(fluxgrd))       deallocate(fluxgrd)
+    if (allocated(hmons))         deallocate(hmons)
+    if (allocated(summit))        deallocate(summit)
+    if (allocated(base))          deallocate(base)
+    if (allocated(alpha_hmons))   deallocate(alpha_hmons)
+    if (allocated(hsummit))       deallocate(hsummit)
+    if (allocated(contains_mons)) deallocate(contains_mons)
     
   end subroutine end_surfdat_h