Changeset 3234

Timestamp:

Feb 23, 2024, 3:30:17 PM (9 months ago)

Author:

emillour

Message:

Mars PCM:
Minor fix in topmons_setup: handle case when mount is "resolved" for high
resolution runs.
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/topmons_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4523 +4523 @@
 Fixing a bug in 1D: atm_wat_profile should evolve freely when considering
 relaxation, and should not be set to qsat in cold conditions. 
+
+== 23/02/2024 == EM
+Minor fix in topmons_setup: handle case when mount is "resolved" for high
+resolution runs.

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

@@ -1056 +1056 @@
             alpha_hmons(ig)= 0.5*(hmons(ig)-hmin)/(hmax-hmin)
             ! Sanity check
-            if (alpha_hmons(ig).le.0) then
+            if (alpha_hmons(ig).lt.0) then
+             write(*,*) "topmons_setup: ig=",ig," alpha_hmons(ig)=",alpha_hmons(ig)
+             write(*,*) "topmons_setup: hmin=",hmin," hmax=",hmax
               call abort_physic("topmons_setup","ERROR: alpha_hmons cannot be <0 "//&
                                 "if the mesh contains a mountain. Please check your "//&
                                 "formula or your start files.",1)
+            endif
+            if (alpha_hmons(ig).eq.0) then
+              write(*,*) "topmons_setup: warning! For ig=",ig, &
+                         " alpha_mons(ig)=0; ", &
+                         " i.e. mount resolved at this resolution"
+              contains_mons(ig)=.false.
             endif
           else
@@ -1080 +1088 @@
             hsummit(:)=14000.
       ENDIF ! (ngrid.gt.1)
-      
+
       end subroutine topmons_setup
 !********************************************************************************