Changeset 3115
- Timestamp:
- Jul 10, 2018 12:49:26 PM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 2 edited
-
init_grid.f90 (modified) (3 diffs)
-
urban_surface_mod.f90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r3103 r3115 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Referencing of buildings onto top of terrain - special treatment for bridges. 28 ! 29 ! 3103 2018-07-04 17:30:52Z suehring 27 30 ! Reference lowest terrain height to zero level 28 31 ! … … 1287 1290 1288 1291 USE netcdf_data_input_mod, & 1289 ONLY: buildings_f, building_id_f, input_pids_static,&1292 ONLY: buildings_f, building_id_f, building_type_f, input_pids_static, & 1290 1293 terrain_height_f 1291 1294 … … 1577 1580 !-- terrain top is defined at upper bound of the grid box. 1578 1581 !-- Hence, check for zw in this case. 1579 DO k = topo_top_index + 1, nzt + 1 1580 IF ( zw(k) - ocean_offset <= oro_max(nr) ) THEN 1581 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1582 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) 1583 ENDIF 1584 ENDDO 1585 ! 1586 !-- After surface irregularities are smoothen, determine lower 1587 !-- start index where building starts. 1588 DO k = nzb, nzt 1589 IF ( zw(k) - ocean_offset <= oro_max(nr) ) & 1590 topo_top_index = k 1591 ENDDO 1582 !-- Note, do this only for buildings which are surface mounted, 1583 !-- i.e. building types 1-6. Below bridges, which are represented 1584 !-- exclusively by building type 7, terrain shape should be 1585 !-- maintained. 1586 IF ( building_type_f%var(j,i) /= 7 ) THEN 1587 DO k = topo_top_index + 1, nzt + 1 1588 IF ( zw(k) - ocean_offset <= oro_max(nr) ) THEN 1589 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1590 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) 1591 ENDIF 1592 ENDDO 1593 ! 1594 !-- After surface irregularities are smoothen, determine lower 1595 !-- start index where building starts. 1596 DO k = nzb, nzt 1597 IF ( zw(k) - ocean_offset <= oro_max(nr) ) & 1598 topo_top_index = k 1599 ENDDO 1600 ENDIF 1592 1601 ! 1593 1602 !-- Finally, map building on top. -
palm/trunk/SOURCE/urban_surface_mod.f90
r3091 r3115 28 28 ! ----------------- 29 29 ! $Id$ 30 ! Additional building type to represent bridges 31 ! 32 ! 3091 2018-06-28 16:20:35Z suehring 30 33 ! - Limit aerodynamic resistance at vertical walls. 31 34 ! - Add check for local roughness length not exceeding surface-layer height and … … 428 431 429 432 430 CHARACTER(37), DIMENSION(0: 6), PARAMETER :: building_type_name = (/ &433 CHARACTER(37), DIMENSION(0:7), PARAMETER :: building_type_name = (/ & 431 434 'user-defined ', & ! 0 432 435 'residential - 1950 ', & ! 1 … … 435 438 'office - 1950 ', & ! 4 436 439 'office 1951 - 2000 ', & ! 5 437 'office 2001 - ' & ! 6 440 'office 2001 - ', & ! 6 441 'bridges ' & ! 7 438 442 /) 439 443 ! … … 460 464 !-- 45 - heat capacity of the wall surface, 46 - heat conductivity 461 465 !-- Please note, only preleminary dummy values so far! 462 REAL(wp), DIMENSION(0:46,1: 6), PARAMETER :: building_pars = RESHAPE( (/ &466 REAL(wp), DIMENSION(0:46,1:7), PARAMETER :: building_pars = RESHAPE( (/ & 463 467 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, 1.0_wp, & !parameter 0-5 464 468 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 465 469 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 466 0.001_wp, 0.0001_wp, 4.0_wp, 470 0.001_wp, 0.0001_wp, 4.0_wp, & !parameter 18-20 467 471 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 468 472 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 … … 475 479 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 476 480 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 477 0.001_wp, 0.0001_wp, 4.0_wp, 481 0.001_wp, 0.0001_wp, 4.0_wp, & !parameter 18-20 478 482 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 479 483 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 … … 486 490 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 487 491 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 488 0.001_wp, 0.0001_wp, 4.0_wp, 492 0.001_wp, 0.0001_wp, 4.0_wp, & !parameter 18-20 489 493 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 490 494 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 … … 508 512 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 509 513 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 510 0.001_wp, 0.0001_wp, 4.0_wp, 514 0.001_wp, 0.0001_wp, 4.0_wp, & !parameter 18-20 511 515 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 512 516 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 … … 519 523 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 520 524 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 521 0.001_wp, 0.0001_wp, 4.0_wp, 525 0.001_wp, 0.0001_wp, 4.0_wp, & !parameter 18-20 522 526 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 523 527 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 … … 526 530 24.0_wp, 24.0_wp, 24.0_wp, & !parameter 38-40 527 531 0.0242_wp, 0.0969_wp, 0.346_wp, 1.0_wp, & !parameter 41-44 528 20000.0_wp, 10.0_wp & !parameter 45-46 - end of type 6 532 20000.0_wp, 10.0_wp, & !parameter 45-46 - end of type 6 533 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 1.0_wp, 1.0_wp, & !parameter 0-5 534 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 6-11 535 296.15_wp, 293.15_wp, 0.9_wp, 0.9_wp, 0.01_wp, 0.99_wp, & !parameter 12-17 536 0.001_wp, 0.0001_wp, 0.0_wp, & !parameter 18-20 537 1.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 3.0_wp, & !parameter 21-25 538 1000000.0_wp, 1000000.0_wp, 1000000.0_wp, & !parameter 26-28 539 0.3_wp, 0.3_wp, 0.3_wp, & !parameter 29-31 540 0.4_wp, 0.4_wp, 0.4_wp, 0.4_wp, 0.01_wp, 0.001_wp, & !parameter 32-37 541 24.0_wp, 24.0_wp, 24.0_wp, & !parameter 38-40 542 0.0242_wp, 0.0969_wp, 0.346_wp, 1.0_wp, & !parameter 41-44 543 20000.0_wp, 10.0_wp & !parameter 45-46 - end of type 7 (bridges) 529 544 /), & 530 (/47, 6/) )545 (/47, 7/) ) 531 546 532 547 !
Note: See TracChangeset
for help on using the changeset viewer.
