Changeset 2787 for palm

Timestamp:

Feb 5, 2018 8:06:52 PM (6 years ago)

Author:

suehring

Message:

Check if 3D building input is consistent to numeric grid.

File:

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

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Check if 3D building input is consistent to numeric grid.
+! 
+! 2773 2018-01-30 14:12:54Z suehring
 ! - Enable initialization with 3D topography.
 ! - Move check for correct initialization in nesting mode to check_parameters.
@@ -230 +233 @@
        INTEGER(iwp)    ::  nz                                    !< number of vertical layers in file
        INTEGER(KIND=1), DIMENSION(:,:,:), ALLOCATABLE ::  var_3d !< 3d variable (lod = 2)
+
+       REAL(wp), DIMENSION(:), ALLOCATABLE ::  z                 !< vertical coordinate for 3D building, used for consistency check
 
        LOGICAL ::  from_file = .FALSE.  !< flag indicating whether an input variable is available and read from file or default values are used  
@@ -1715 +1720 @@
  
                    IF ( buildings_f%lod == 2 )  THEN
-                      ALLOCATE( buildings_f%var_3d(nzb:buildings_f%nz,         &
+                      ALLOCATE( buildings_f%z(nzb:buildings_f%nz-1) )
+                      CALL get_variable( id_topo, 'z', buildings_f%z )
+
+                      ALLOCATE( buildings_f%var_3d(nzb:buildings_f%nz-1,       &
                                                    nys:nyn,nxl:nxr) )
                       buildings_f%var_3d = 0
@@ -2724 +2732 @@
     SUBROUTINE netcdf_data_input_check_static
 
+       USE arrays_3d,                                                          &
+           ONLY:  zu
+
        USE control_parameters,                                                 &
            ONLY:  land_surface, message_string, urban_surface
@@ -2750 +2761 @@
           CALL message( 'netcdf_data_input_mod', 'PA0999', 2, 2, 0, 6, 0 )
        ENDIF
+!
+!--    If 3D buildings are read, check if building information is consistent 
+!--    to numeric grid. 
+       IF ( buildings_f%from_file )  THEN
+          IF ( buildings_f%lod == 2 )  THEN
+             IF ( buildings_f%nz > SIZE( zu ) )  THEN
+                message_string = 'Reading 3D building data - too much ' //     &
+                                 'data points along the vertical coordinate.' 
+                CALL message( 'netcdf_data_input_mod', 'PA0999', 2, 2, 0, 6, 0 )
+             ENDIF
+
+             IF ( ANY( buildings_f%z(0:buildings_f%nz-1) /=                    &
+                       zu(0:buildings_f%nz-1) ) )  THEN
+                message_string = 'Reading 3D building data - vertical ' //     &
+                                 'coordinate do not match numeric grid.'
+                CALL message( 'netcdf_data_input_mod', 'PA0999', 2, 2, 0, 6, 0 )
+             ENDIF
+          ENDIF
+       ENDIF
+
 !
 !--    Skip further checks concerning buildings and natural surface properties