Changeset 3103 for palm

Timestamp:

Jul 4, 2018 5:30:52 PM (6 years ago)

Author:

suehring

Message:

Reference lowest terrain height to zero level; new check for negative terrain heights in static input file

Location:

palm/trunk/SOURCE

Files:

• init_grid.f90 (modified) (3 diffs)
• netcdf_data_input_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Reference lowest terrain height to zero level
+! 
+! 3068 2018-06-12 14:49:41Z Giersch
 ! New warning message concerning grid stretching has been introduced
 ! 
@@ -1276 +1279 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, land_surface, ocean, urban_surface
+        ONLY:  bc_lr_cyc, bc_ns_cyc, land_surface, message_string, ocean,      &
+               urban_surface
 
     USE indices,                                                               &
@@ -1313 +1317 @@
 
     REAL(wp)                            ::  ocean_offset        !< offset to consider inverse vertical coordinate at topography definition
+    REAL(wp)                            ::  oro_min = 0.0_wp    !< minimum terrain height in entire model domain, used to reference terrain to zero
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  oro_max             !< maximum terrain height occupied by an building with certain id
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  oro_max_l           !< maximum terrain height occupied by an building with certain id, on local subdomain
 
+
+!
+!-- Reference lowest terrain height to zero. In case the minimum terrain height
+!-- is non-zero, all grid points of the lower vertical grid levels might be 
+!-- entirely below the surface, meaning a waste of computational resources.
+!-- In order to avoid this, remove the lowest terrain height. Please note,
+!-- in case of a nested run, the global minimum from all parent and childs 
+!-- need to be remove to avoid steep edges at the child-domain boundaries.
+    IF ( input_pids_static )  THEN
+
+       CALL MPI_ALLREDUCE( MINVAL( terrain_height_f%var ), oro_min, 1,         &
+                           MPI_REAL, MPI_MIN, MPI_COMM_WORLD, ierr )
+                           
+       terrain_height_f%var = terrain_height_f%var - oro_min
+!                           
+!--    Give an informative message that terrain height is referenced to zero    
+       IF ( oro_min > 0.0_wp )  THEN
+          WRITE( message_string, * ) 'Terrain height was internally shifted '//&
+                          'downwards by ', oro_min, 'meter(s) to save ' //     &
+                          'computational resources.'
+          CALL message( 'init_grid', 'PA0505', 0, 0, 0, 6, 0 )
+       ENDIF
+    ENDIF    
+    
 !
 !-- In the following, buildings and orography are further preprocessed 

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! New check for negative terrain heights
+! 
+! 3089 2018-06-27 13:20:38Z suehring
 ! Revise call for message routine in case of local data inconsistencies.
 ! 
@@ -2883 +2886 @@
           message_string = 'NetCDF variable zt is not ' //                     &
                            'allowed to have missing data'
+          CALL message( 'netcdf_data_input_mod', 'NDI013', 2, 2, myid, 6, 0 )
+       ENDIF
+!
+!--    Check for negative terrain heights
+       IF ( ANY( terrain_height_f%var < 0.0_wp ) )  THEN
+          message_string = 'NetCDF variable zt is not ' //                     &
+                           'allowed to have negative values'
           CALL message( 'netcdf_data_input_mod', 'NDI013', 2, 2, myid, 6, 0 )
        ENDIF