Changeset 153 for palm/trunk/SOURCE/plant_canopy_model.f90
- Timestamp:
- Mar 19, 2008 9:41:30 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model.f90
r142 r153 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! Consideration of heat sources within the forest canopy which represent the 7 ! rate of heat input into the air from the forest leaves 8 ! 9 ! Evaluation of sinks and sources for scalar concentration due to canopy 10 ! elements 11 ! 6 12 ! Bugfix: remove IF statement in plant_canopy_model_ij 7 13 ! … … 15 21 ! Description: 16 22 ! ------------ 17 ! Evaluation of the drag due to vegetation 23 ! Evaluation of sinks and sources of momentum, heat and scalar concentration 24 ! due to canopy elements 18 25 !------------------------------------------------------------------------------! 19 26 … … 59 66 v(k,j,i) + & 60 67 v(k,j+1,i) + & 61 v(k,j+1,i +1) ) &68 v(k,j+1,i-1) ) & 62 69 / 4.0 )**2 + & 63 70 ( ( w(k-1,j,i-1) + & … … 120 127 121 128 ! 129 !-- potential temperature 130 CASE ( 4 ) 131 DO i = nxl, nxr 132 DO j = nys, nyn 133 DO k = nzb_s_inner(j,i)+1, pch_index 134 tend(k,j,i) = tend(k,j,i) + & 135 ( canopy_heat_flux(k,j,i) - & 136 canopy_heat_flux(k-1,j,i) ) / & 137 dzw(k) 138 ENDDO 139 ENDDO 140 ENDDO 141 142 ! 143 !-- scalar concentration 144 CASE ( 5 ) 145 DO i = nxl, nxr 146 DO j = nys, nyn 147 DO k = nzb_s_inner(j,i)+1, pch_index 148 tend(k,j,i) = tend(k,j,i) - & 149 sec(k,j,i) * lad_s(k,j,i) * & 150 SQRT( ( ( u(k,j,i) + & 151 u(k,j,i+1) ) & 152 / 2.0 )**2 + & 153 ( ( v(k,j,i) + & 154 v(k,j+1,i) ) & 155 / 2.0 )**2 + & 156 ( ( w(k-1,j,i) + & 157 w(k,j,i) ) & 158 / 2.0 )**2 ) * & 159 ( q(k,j,i) - sls(k,j,i) ) 160 ENDDO 161 ENDDO 162 ENDDO 163 164 ! 122 165 !-- sgs-tke 123 CASE ( 4)166 CASE ( 6 ) 124 167 DO i = nxl, nxr 125 168 DO j = nys, nyn … … 180 223 v(k,j,i) + & 181 224 v(k,j+1,i) + & 182 v(k,j+1,i +1) ) &225 v(k,j+1,i-1) ) & 183 226 / 4.0 )**2 + & 184 227 ( ( w(k-1,j,i-1) + & … … 232 275 233 276 ! 277 !-- potential temperature 278 CASE ( 4 ) 279 DO k = nzb_s_inner(j,i)+1, pch_index 280 tend(k,j,i) = tend(k,j,i) + & 281 ( canopy_heat_flux(k,j,i) - & 282 canopy_heat_flux(k-1,j,i) ) / & 283 dzw(k) 284 ENDDO 285 286 287 ! 288 !-- scalar concentration 289 CASE ( 5 ) 290 DO k = nzb_s_inner(j,i)+1, pch_index 291 tend(k,j,i) = tend(k,j,i) - & 292 sec(k,j,i) * lad_s(k,j,i) * & 293 SQRT( ( ( u(k,j,i) + & 294 u(k,j,i+1) ) & 295 / 2.0 )**2 + & 296 ( ( v(k,j,i) + & 297 v(k,j+1,i) ) & 298 / 2.0 )**2 + & 299 ( ( w(k-1,j,i) + & 300 w(k,j,i) ) & 301 / 2.0 )**2 ) * & 302 ( q(k,j,i) - sls(k,j,i) ) 303 ENDDO 304 305 ! 234 306 !-- sgs-tke 235 CASE ( 4)307 CASE ( 6 ) 236 308 DO k = nzb_s_inner(j,i)+1, pch_index 237 309 tend(k,j,i) = tend(k,j,i) - &
Note: See TracChangeset
for help on using the changeset viewer.