Changeset 3855 for palm/trunk/SOURCE/init_grid.f90
- Timestamp:
- Apr 3, 2019 10:00:59 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_grid.f90
r3763 r3855 25 25 ! ----------------- 26 26 ! $Id$ 27 ! In projection of non-building 3D objects onto numerical grid remove 28 ! dependency on building_type 29 ! 30 ! 3763 2019-02-25 17:33:49Z suehring 27 31 ! Replace work-around for ghost point exchange of 1-byte arrays with specific 28 32 ! routine as already done in other routines … … 1606 1610 !-- exclusively by building type 7, terrain shape should be 1607 1611 !-- maintained. 1608 IF ( building_type_f%var(j,i) /= 7 ) THEN 1609 DO k = topo_top_index + 1, nzt + 1 1610 IF ( zw(k) - ocean_offset <= oro_max(nr) ) THEN 1611 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1612 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) 1613 ENDIF 1614 ENDDO 1615 ! 1616 !-- After surface irregularities are smoothen, determine lower 1617 !-- start index where building starts. 1618 DO k = nzb, nzt 1619 IF ( zw(k) - ocean_offset <= oro_max(nr) ) & 1620 topo_top_index = k 1621 ENDDO 1612 IF ( building_type_f%from_file ) THEN 1613 IF ( building_type_f%var(j,i) /= 7 ) THEN 1614 DO k = topo_top_index + 1, nzt + 1 1615 IF ( zw(k) - ocean_offset <= oro_max(nr) ) THEN 1616 topo_3d(k,j,i) = IBCLR( topo_3d(k,j,i), 0 ) 1617 topo_3d(k,j,i) = IBSET( topo_3d(k,j,i), 2 ) 1618 ENDIF 1619 ENDDO 1620 ! 1621 !-- After surface irregularities are smoothen, determine 1622 !-- lower start index where building starts. 1623 DO k = nzb, nzt 1624 IF ( zw(k) - ocean_offset <= oro_max(nr) ) & 1625 topo_top_index = k 1626 ENDDO 1627 ENDIF 1622 1628 ENDIF 1623 1629 !
Note: See TracChangeset
for help on using the changeset viewer.