Changeset 3103 for palm/trunk/SOURCE/init_grid.f90
- Timestamp:
- Jul 4, 2018 5:30:52 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r3068 r3103 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Reference lowest terrain height to zero level 28 ! 29 ! 3068 2018-06-12 14:49:41Z Giersch 27 30 ! New warning message concerning grid stretching has been introduced 28 31 ! … … 1276 1279 1277 1280 USE control_parameters, & 1278 ONLY: bc_lr_cyc, bc_ns_cyc, land_surface, ocean, urban_surface 1281 ONLY: bc_lr_cyc, bc_ns_cyc, land_surface, message_string, ocean, & 1282 urban_surface 1279 1283 1280 1284 USE indices, & … … 1313 1317 1314 1318 REAL(wp) :: ocean_offset !< offset to consider inverse vertical coordinate at topography definition 1319 REAL(wp) :: oro_min = 0.0_wp !< minimum terrain height in entire model domain, used to reference terrain to zero 1315 1320 REAL(wp), DIMENSION(:), ALLOCATABLE :: oro_max !< maximum terrain height occupied by an building with certain id 1316 1321 REAL(wp), DIMENSION(:), ALLOCATABLE :: oro_max_l !< maximum terrain height occupied by an building with certain id, on local subdomain 1317 1322 1323 1324 ! 1325 !-- Reference lowest terrain height to zero. In case the minimum terrain height 1326 !-- is non-zero, all grid points of the lower vertical grid levels might be 1327 !-- entirely below the surface, meaning a waste of computational resources. 1328 !-- In order to avoid this, remove the lowest terrain height. Please note, 1329 !-- in case of a nested run, the global minimum from all parent and childs 1330 !-- need to be remove to avoid steep edges at the child-domain boundaries. 1331 IF ( input_pids_static ) THEN 1332 1333 CALL MPI_ALLREDUCE( MINVAL( terrain_height_f%var ), oro_min, 1, & 1334 MPI_REAL, MPI_MIN, MPI_COMM_WORLD, ierr ) 1335 1336 terrain_height_f%var = terrain_height_f%var - oro_min 1337 ! 1338 !-- Give an informative message that terrain height is referenced to zero 1339 IF ( oro_min > 0.0_wp ) THEN 1340 WRITE( message_string, * ) 'Terrain height was internally shifted '//& 1341 'downwards by ', oro_min, 'meter(s) to save ' // & 1342 'computational resources.' 1343 CALL message( 'init_grid', 'PA0505', 0, 0, 0, 6, 0 ) 1344 ENDIF 1345 ENDIF 1346 1318 1347 ! 1319 1348 !-- In the following, buildings and orography are further preprocessed
Note: See TracChangeset
for help on using the changeset viewer.