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

Timestamp:

Apr 9, 2018 3:14:01 PM (6 years ago)

Author:

suehring

Message:

Revise topography filter; extend checks for consistent setting of building ID and type; add cpu measures in netcdf-data input

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Improve topography filter routine and add ghost-point exchange for building 
+! ID and building type. 
+! 
+! 2927 2018-03-23 15:13:00Z suehring
 ! Bugfix, setting boundary conditions for topography index array.
 ! 
@@ -1186 +1190 @@
     INTEGER(iwp) ::  num_wall   !< number of surrounding vertical walls for a single grid point
 
-    INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE ::  topo_tmp      !< temporary 3D-topography used to fill holes
-    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  topo_3d       !< 3D-topography array merging buildings and orography
+    INTEGER(iwp), DIMENSION(nysg:nyng,nxlg:nxrg)           ::  var_exchange_int  !< dummy array for exchanging ghost-points 
+    INTEGER(iwp), DIMENSION(:,:,:), ALLOCATABLE            ::  topo_tmp          !< temporary 3D-topography used to fill holes
+    INTEGER(iwp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  topo_3d           !< 3D-topography array merging buildings and orography
 !
 !-- Before checking for holes, set lateral boundary conditions for 
@@ -1201 +1206 @@
        num_hole = 0    
        CALL exchange_horiz_int( topo_3d, nys, nyn, nxl, nxr, nzt, nbgp )
+!
+!--    Exchange also building ID and type. Note, building_type is an one-byte 
+!--    variable.
+       IF ( building_id_f%from_file )                                          &
+          CALL exchange_horiz_2d_int( building_id_f%var, nys, nyn, nxl, nxr, nbgp )
+       IF ( building_type_f%from_file )  THEN
+          var_exchange_int = INT( building_type_f%var, KIND = 4 )
+          CALL exchange_horiz_2d_int( var_exchange_int, nys, nyn, nxl, nxr, nbgp )
+          building_type_f%var = INT( var_exchange_int, KIND = 1 )
+       ENDIF
 
        topo_tmp = topo_3d
@@ -1343 +1358 @@
        IF ( nxl == 0  )  topo_3d(:,:,-1)   = topo_3d(:,:,0)
        IF ( nxr == nx )  topo_3d(:,:,nx+1) = topo_3d(:,:,nx)          
+    ENDIF
+!
+!-- Exchange building ID and type. Note, building_type is an one-byte variable.
+    IF ( building_id_f%from_file )                                             &
+       CALL exchange_horiz_2d_int( building_id_f%var, nys, nyn, nxl, nxr, nbgp )
+    IF ( building_type_f%from_file )  THEN
+       var_exchange_int = INT( building_type_f%var, KIND = 4 )
+       CALL exchange_horiz_2d_int( var_exchange_int, nys, nyn, nxl, nxr, nbgp )
+       building_type_f%var = INT( var_exchange_int, KIND = 1 )
     ENDIF