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

Timestamp:

Sep 22, 2020 2:38:38 PM (4 years ago)

Author:

suehring

Message:

Reference model topography to the lowest grid point also in ASCII input case

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Reference model topography to the lowest grid point also in ASCII input case
+!
+! 4653 2020-08-27 08:54:43Z pavelkrc
 ! file re-formatted to follow the PALM coding standard
 !
@@ -950 +953 @@
 !--    Update reference height used within output files
        init_model%origin_z = init_model%origin_z + oro_min
-
+!
+!-- ASCII topography branch. In this case, in contrast to the static driver input, topography is
+!-- assumed to be building; thus the minimum building height is substracted from the building array.
+    ELSE
+#if defined( __parallel )
+       CALL MPI_ALLREDUCE( MINVAL( buildings_f%var_2d ), oro_min, 1, MPI_REAL, MPI_MIN,            &
+                           MPI_COMM_WORLD, ierr )
+#else
+       oro_min = MINVAL( buildings_f%var_2d )
+#endif
+       buildings_f%var_2d = buildings_f%var_2d - oro_min
+!
+!--    Update reference height used within output files
+       init_model%origin_z = init_model%origin_z + oro_min
     ENDIF