Changeset 4329 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
- Timestamp:
- Dec 10, 2019 3:46:36 PM (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r4314 r4329 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Renamed wall_flags_0 to wall_flags_static_0 30 ! 31 ! 4314 2019-11-29 10:29:20Z suehring 29 32 ! - Bugfix, plant canopy was still considered at building edges on for the u- 30 33 ! and v-component. … … 121 124 USE indices, & 122 125 ONLY: nbgp, nxl, nxlg, nxlu, nxr, nxrg, nyn, nyng, nys, nysg, nysv, & 123 nz, nzb, nzt, topo_top_ind, wall_flags_ 0126 nz, nzb, nzt, topo_top_ind, wall_flags_static_0 124 127 125 128 USE kinds … … 1125 1128 !-- an informative message is given. 1126 1129 IF ( ANY( lad_s(:,j,i) /= 0.0_wp ) .AND. & 1127 ANY( BTEST( wall_flags_ 0(:,j,i), 6 ) ) .AND. &1128 ANY( BTEST( wall_flags_ 0(:,j,i), 4 ) ) ) THEN1130 ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) .AND. & 1131 ANY( BTEST( wall_flags_static_0(:,j,i), 4 ) ) ) THEN 1129 1132 lad_s(:,j,i) = 0.0_wp 1130 1133 WRITE( message_string, * ) & … … 1958 1961 !-- of the potential rooftop grid point, while building_edge_e is 1959 1962 !-- is set from the perspective of the non-building grid point. 1960 building_edge_w = ANY( BTEST( wall_flags_ 0(:,j,i), 6 ) ) .AND. &1961 .NOT. ANY( BTEST( wall_flags_ 0(:,j,i-1), 6 ) )1962 building_edge_e = ANY( BTEST( wall_flags_ 0(:,j,i-1), 6 ) ) .AND. &1963 .NOT. ANY( BTEST( wall_flags_ 0(:,j,i), 6 ) )1963 building_edge_w = ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) .AND. & 1964 .NOT. ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) ) 1965 building_edge_e = ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) ) .AND. & 1966 .NOT. ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) 1964 1967 ! 1965 1968 !-- Determine topography-top index on u-grid … … 2038 2041 !-- of the potential rooftop grid point, while building_edge_n is 2039 2042 !-- is set from the perspective of the non-building grid point. 2040 building_edge_s = ANY( BTEST( wall_flags_ 0(:,j,i), 6 ) ) .AND. &2041 .NOT. ANY( BTEST( wall_flags_ 0(:,j-1,i), 6 ) )2042 building_edge_n = ANY( BTEST( wall_flags_ 0(:,j-1,i), 6 ) ) .AND. &2043 .NOT. ANY( BTEST( wall_flags_ 0(:,j,i), 6 ) )2043 building_edge_s = ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) .AND. & 2044 .NOT. ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) ) 2045 building_edge_n = ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) ) .AND. & 2046 .NOT. ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) 2044 2047 ! 2045 2048 !-- Determine topography-top index on v-grid
Note: See TracChangeset
for help on using the changeset viewer.