Changeset 4302 for palm/trunk/SOURCE

Timestamp:

Nov 22, 2019 1:15:56 PM (5 years ago)

Author:

suehring

Message:

Omit tall canopy mapped on top of buildings, which may happen due to topography filtering

File:

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

palm/trunk/SOURCE/plant_canopy_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Omit tall canopy mapped on top of buildings
+! 
+! 4279 2019-10-29 08:48:17Z scharf
 ! unused variables removed
 ! 
@@ -932 +935 @@
            ONLY: message_string, ocean_mode
 
+       USE indices,                                                            &
+           ONLY:  wall_flags_0
+
        USE netcdf_data_input_mod,                                              &
            ONLY:  leaf_area_density_f
+
+       USE pegrid
 
        USE surface_mod,                                                        &
@@ -1046 +1054 @@
           DO k = 1, pch_index
              lad(k) = 0.5 * ( pre_lad(k-1) + pre_lad(k) )
-          ENDDO          
+          ENDDO
 
        ENDIF
@@ -1053 +1061 @@
 !--    Allocate 3D-array for the leaf area density (lad_s).
        ALLOCATE( lad_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
-
 !
 !--    Initialize canopy parameters cdc (canopy drag coefficient), 
@@ -1096 +1103 @@
                          lad_s(k,j,i) = leaf_area_density_f%var(k,j,i)
                       ENDDO
-                   ENDDO
-                ENDDO
-                
+!
+!--                   Check if resolved vegetation is mapped onto buildings.
+!--                   Even if the user defines LAD only for non-building covered
+!--                   grid points, the filter process of topography may cause
+!--                   that LAD overlaps with buildings. This case, an 
+!--                   informative message is given and the LAD field is set to
+!--                   zero at these 
+                      IF ( ANY( lad_s(:,j,i) /= 0.0_wp )  .AND.                &
+                           ANY( BTEST( wall_flags_0(:,j,i), 6 ) ) )  THEN
+                         lad_s(:,j,i) = 0.0_wp
+                         WRITE( message_string, * )                            &
+                                          'Resolved plant-canopy is ' //       &
+                                          'is defined on top of a building ' //&
+                                          '- lad is omitted at this grid '   //&
+                                          'point: (i,j) = ', i, j
+                         CALL message( 'pcm_init', 'PA0313', 0, 0, myid, 6, 0 )
+                      ENDIF
+                   ENDDO
+                ENDDO
                 CALL exchange_horiz( lad_s, nbgp )
 !