Changeset 3234 for trunk/LMDZ.MARS/libf

Timestamp:

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

Author:

emillour

Message:

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

File:

• trunk/LMDZ.MARS/libf/phymars/topmons_mod.F90 (modified) (2 diffs)

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
 !********************************************************************************