Changeset 4245 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Sep 30, 2019 8:40:37 AM (5 years ago)

Author:

pavelkrc

Message:

Merge branch resler into trunk

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/init_grid.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Store oro_max (building z-offset) in 2D for building surfaces
+! 
+! 4189 2019-08-26 16:19:38Z suehring
 ! - Add check for proper setting of namelist parameter topography
 ! - Set flag to indicate land surfaces in case no topography is provided
@@ -1312 +1315 @@
        ENDIF
 !
+!--    Allocate array for storing terrain height under buildings
+       IF ( buildings_f%from_file )  THEN
+          ALLOCATE( buildings_f%oro_max(nysg:nyng,nxlg:nxrg) )
+          buildings_f%oro_max = buildings_f%fill1
+       END IF
+!
 !--    Map orography as well as buildings onto grid. 
        DO  i = nxl, nxr
@@ -1325 +1334 @@
                    nr = MINLOC( ABS( build_ids_final -                         &
                                      building_id_f%var(j,i) ), DIM = 1 )
+!
+!--                Save grid-indexed oro_max
+                   buildings_f%oro_max(j,i) = oro_max(nr)
                 ENDIF
              ENDIF
@@ -1433 +1445 @@
           ENDDO
        ENDDO
+!
+!--    Horizontal exchange the oro_max array
+       CALL exchange_horiz_2d( buildings_f%oro_max(:,:), nbgp )
 !
 !--    Deallocate temporary arrays required for processing and reading data