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

Timestamp:

Apr 13, 2018 11:52:24 AM (6 years ago)

Author:

suehring

Message:

Bugfixes in initalization of land-surface model as well as timeseries data output in case of elevated model surfaces.

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in initialization in case of elevated model surface. Introduce
+! index for minimum topography-top. 
+! 
+! 2955 2018-04-09 15:14:01Z suehring
 ! Improve topography filter routine and add ghost-point exchange for building 
 ! ID and building type. 
@@ -332 +336 @@
                nzb_max, nzb_s_inner, nzb_s_outer, nzb_u_inner,                 &
                nzb_u_outer, nzb_v_inner, nzb_v_outer, nzb_w_inner,             &
-               nzb_w_outer, nzt
+               nzb_w_outer, nzt, topo_min_level
     
     USE kinds
@@ -348 +352 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i             !< index variable along x 
-    INTEGER(iwp) ::  j             !< index variable along y
-    INTEGER(iwp) ::  k             !< index variable along z
-    INTEGER(iwp) ::  k_top         !< topography top index on local PE
-    INTEGER(iwp) ::  l             !< loop variable 
-    INTEGER(iwp) ::  nzb_local_max !< vertical grid index of maximum topography height
-    INTEGER(iwp) ::  nzb_local_min !< vertical grid index of minimum topography height
+    INTEGER(iwp) ::  i                !< index variable along x 
+    INTEGER(iwp) ::  j                !< index variable along y
+    INTEGER(iwp) ::  k                !< index variable along z
+    INTEGER(iwp) ::  k_top            !< topography top index on local PE
+    INTEGER(iwp) ::  l                !< loop variable 
+    INTEGER(iwp) ::  nzb_local_max    !< vertical grid index of maximum topography height
+    INTEGER(iwp) ::  nzb_local_min    !< vertical grid index of minimum topography height
                                       
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local      !< index for topography top at cell-center
@@ -553 +557 @@
 !   
 !-- Finally, if topography extents up to the model top, limit nzb_max to nzt.
-    nzb_max = MIN( nzb_max, nzt ) 
-
+    nzb_max = MIN( nzb_max, nzt )
+!
+!-- Determine minimum index of topography. Usually, this will be nzb. In case
+!-- there is elevated topography, however, the lowest topography will be higher. 
+!-- This index is e.g. used to calculate mean first-grid point atmosphere 
+!-- temperature, surface pressure and density, etc. .
+    topo_min_level   = 0
+#if defined( __parallel )
+    CALL MPI_ALLREDUCE( MINVAL( get_topography_top_index( 's' ) ),             &
+                        topo_min_level, 1, MPI_INTEGER, MPI_MIN, comm2d, ierr )
+#else
+    topo_min_level = MINVAL( get_topography_top_index( 's' ) )
+#endif
 !
 !-- Initialize boundary conditions via surface type
     CALL init_bc
-
 !
 !-- Allocate and set topography height arrays required for data output
@@ -621 +635 @@
     nzb_local(nys:nyn,nxl:nxr) = get_topography_top_index( 's' )
     CALL exchange_horiz_2d_int( nzb_local, nys, nyn, nxl, nxr, nbgp )
-
+!
+!-- Check topography for consistency with model domain. Therefore, use
+!-- maximum and minium topography-top indices. Note, minimum topography top
+!-- index is already calculated.  
     IF ( TRIM( topography ) /= 'flat' )  THEN
 #if defined( __parallel )
        CALL MPI_ALLREDUCE( MAXVAL( get_topography_top_index( 's' ) ),          &
-                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )
-       CALL MPI_ALLREDUCE( MINVAL( get_topography_top_index( 's' ) ),          &
-                           nzb_local_min, 1, MPI_INTEGER, MPI_MIN, comm2d, ierr )                   
+                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )             
 #else
        nzb_local_max = MAXVAL( get_topography_top_index( 's' ) )
-       nzb_local_min = MINVAL( get_topography_top_index( 's' ) )
 #endif
+       nzb_local_min = topo_min_level
 !
 !--    Consistency checks