Changeset 4895 for palm/trunk/SOURCE/init_masks.f90

Timestamp:

Mar 3, 2021 3:39:08 PM (3 years ago)

Author:

suehring

Message:

Remove offset in terrain-following masked output and allow only mask_k_over_surface >= 1

File:

• palm/trunk/SOURCE/init_masks.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_masks.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Allow only mask_k_over_surface >= 1 for terrain-following mask output
+!
+! 4828 2021-01-05 11:21:41Z Giersch
 ! Implement snow and graupel (bulk microphysics)
 !
@@ -521 +524 @@
           DO  WHILE ( mask_k_over_surface(mid, count+1) >= 0 )
              m = mask_k_over_surface(mid, count+1)
+             IF ( m < 1 )  THEN
+                WRITE ( message_string, '(I3,A,I3,A)' )  m, ' in mask ', mid, ' must be > 0.'
+                CALL message( 'init_masks', 'PA0743', 1, 2, 0, 6, 0 )
+             ENDIF
              IF ( m > nz+1 )  THEN
                 WRITE ( message_string, '(I3,A,I3,A,I1,3A,I3)' )  m, ' in mask ', mid,             &