Changeset 3925

Timestamp:

Apr 23, 2019 11:30:16 AM (6 years ago)

Author:

suehring

Message:

Minor bugfix in building mapping when all building IDs in the model domain are missing

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor bugfix in building mapping when all building IDs in the model domain
+! are missing
+! 
+! 3857 2019-04-03 13:00:16Z knoop
 ! In projection of non-building 3D objects onto numerical grid remove 
 ! dependency on building_type
@@ -1360 +1364 @@
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  oro_max_l           !< maximum terrain height occupied by an building with certain id, on local subdomain
 
-
+building_id_f%var = building_id_f%fill
+building_type_f%var = building_type_f%fill
 !
 !-- Reference lowest terrain height to zero. In case the minimum terrain height
@@ -1424 +1429 @@
           num_buildings   = 0
 !
-!--       Allocate at least one element for building ids, 
+!--       Allocate at least one element for building ids and give it an inital
+!--       negative value that will be overwritten later. This, however, is 
+!--       necessary in case there all IDs in the model domain are fill values. 
           ALLOCATE( build_ids_l(1) )
+          build_ids_l = -1 
           DO  i = nxl, nxr
              DO  j = nys, nyn