Changeset 4381 for palm/trunk/SOURCE/plant_canopy_model_mod.f90

Timestamp:

Jan 20, 2020 1:51:46 PM (5 years ago)

Author:

suehring

Message:

Land-surface model: Bugfix in nested soil initialization in case no dynamic input file is present; give local error messages only onces

File:

• palm/trunk/SOURCE/plant_canopy_model_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Give error message 313 only once
+! 
+! 4363 2020-01-07 18:11:28Z suehring
 ! Fix for last commit
 ! 
@@ -973 +976 @@
        INTEGER(iwp) ::  m   !< running index
 
+       LOGICAL      ::  lad_on_top = .FALSE.  !< dummy flag to indicate that LAD is defined on a building roof
+
        REAL(wp) ::  canopy_height   !< canopy height for lad-profile construction
        REAL(wp) ::  gradient        !< gradient for lad-profile construction
@@ -1227 +1232 @@
                            ANY( BTEST( wall_flags_total_0(:,j,i), 4 ) ) )  THEN
                          lad_s(:,j,i) = 0.0_wp
-                         WRITE( message_string, * )                            &
+                         lad_on_top   = .TRUE.
+                      ENDIF
+                   ENDDO
+                ENDDO
+#if defined( __parallel )
+               CALL MPI_ALLREDUCE( MPI_IN_PLACE, lad_on_top, 1, MPI_LOGICAL,  &
+                                   MPI_LOR, comm2d, ierr)
+#endif
+                IF ( lad_on_top )  THEN
+                   WRITE( message_string, * )                                  &
                                         'Resolved plant-canopy is ' //         &
                                         'defined on top of an artificially '// &
-                                        'created building grid point ' //      &
+                                        'created building grid point(s) ' //   &
                                         '(emerged from the filtering) - ' //   &
-                                        'LAD profile is omitted at this ' //   &
-                                        'grid point: (i,j) = ', i, j
-                         CALL message( 'pcm_init', 'PA0313', 0, 0, myid, 6, 0 )
-                      ENDIF
-                   ENDDO
-                ENDDO
+                                        'LAD profile is omitted at this / ' // &
+                                        'these grid point(s).'
+                   CALL message( 'pcm_init', 'PA0313', 0, 0, 0, 6, 0 )
+                ENDIF
                 CALL exchange_horiz( lad_s, nbgp )
 !