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

Timestamp:

Apr 22, 2020 6:21:45 PM (4 years ago)

Author:

gronemeier

Message:

netcdf_data_input_mod:

• bugfix: check terrain height for fill values directly after reading (PA0550)
• changes:
  • remove check for negative zt (PA0551)
  • add reference height from input file to PALM reference height (origin_z)

init_grid:

• update origin_z with shifting height of orography (oro_min)
• remove message PA0505

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! update origin_z with shifting height of orography (oro_min)
+! 
+! 4457 2020-03-11 14:20:43Z raasch
 ! use statement for exchange horiz added,
 ! bugfix for call of exchange horiz 2d
@@ -936 +939 @@
 
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, message_string, ocean_mode
+        ONLY:  bc_lr_cyc, bc_ns_cyc, ocean_mode
 
     USE exchange_horiz_mod,                                                    &
@@ -946 +949 @@
 
     USE netcdf_data_input_mod,                                                 &
-        ONLY:  buildings_f, building_id_f, building_type_f, input_pids_static, &
+        ONLY:  buildings_f, building_id_f, building_type_f,                    &
+               init_model,                                                     &
+               input_pids_static,                                              &
                terrain_height_f
 
@@ -983 +988 @@
 
 !
-!-- 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.
+!-- Reference lowest terrain height to zero. This ensures that first, 
+!-- non-required gird levels (those which lie entirely below the minimum 
+!-- orography) are avoided, and second, that also negative orography can be used
+!-- within the input file.
+!-- 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
     
@@ -998 +1003 @@
 #endif
        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
+!
+!--    Update reference height used within output files
+       init_model%origin_z = init_model%origin_z + oro_min
+
     ENDIF    
     
@@ -2604 +2605 @@
 
  END SUBROUTINE set_topo_flags
-
-
-
-