- Timestamp:
- Jan 13, 2016 8:19:40 AM (9 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_u.f90
r1692 r1740 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! unnecessary calculations of kmzm and kmzp in wall bounded parts removed 22 22 ! 23 23 ! Former revisions: … … 225 225 kmzp = 0.25_wp * & 226 226 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) ) 227 kmzm = 0.25_wp * &228 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )229 227 230 228 tend(k,j,i) = tend(k,j,i) & … … 243 241 ! 244 242 !-- Interpolate eddy diffusivities on staggered gridpoints 245 kmzp = 0.25_wp * &246 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )247 243 kmzm = 0.25_wp * & 248 244 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) ) … … 410 406 kmzp = 0.25_wp * & 411 407 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) ) 412 kmzm = 0.25_wp * &413 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )414 408 415 409 tend(k,j,i) = tend(k,j,i) & … … 436 430 ! 437 431 !-- Interpolate eddy diffusivities on staggered gridpoints 438 kmzp = 0.25_wp * &439 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )440 432 kmzm = 0.25_wp * & 441 433 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) ) … … 574 566 !-- Interpolate eddy diffusivities on staggered gridpoints 575 567 kmzp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) ) 576 kmzm = 0.25_wp * ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) )577 568 578 569 tend(k,j,i) = tend(k,j,i) & … … 591 582 ! 592 583 !-- Interpolate eddy diffusivities on staggered gridpoints 593 kmzp = 0.25_wp * &594 ( km(k,j,i)+km(k+1,j,i)+km(k,j,i-1)+km(k+1,j,i-1) )595 584 kmzm = 0.25_wp * & 596 585 ( km(k,j,i)+km(k-1,j,i)+km(k,j,i-1)+km(k-1,j,i-1) ) -
palm/trunk/SOURCE/diffusion_v.f90
r1683 r1740 19 19 ! Current revisions: 20 20 ! ----------------- 21 ! 21 ! unnecessary calculations of kmzm and kmzp in wall bounded parts removed 22 22 ! 23 23 ! Former revisions: … … 220 220 kmzp = 0.25_wp * & 221 221 ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) ) 222 kmzm = 0.25_wp * &223 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )224 222 225 223 tend(k,j,i) = tend(k,j,i) & … … 238 236 ! 239 237 !-- Interpolate eddy diffusivities on staggered gridpoints 240 kmzp = 0.25_wp * &241 ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )242 238 kmzm = 0.25_wp * & 243 239 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) ) … … 405 401 kmzp = 0.25_wp * & 406 402 ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) ) 407 kmzm = 0.25_wp * &408 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )409 403 410 404 tend(k,j,i) = tend(k,j,i) & … … 431 425 ! 432 426 !-- Interpolate eddy diffusivities on staggered gridpoints 433 kmzp = 0.25_wp * &434 ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )435 427 kmzm = 0.25_wp * & 436 428 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) ) … … 571 563 !-- Interpolate eddy diffusivities on staggered gridpoints 572 564 kmzp = 0.25_wp * ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) ) 573 kmzm = 0.25_wp * ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )574 565 575 566 tend(k,j,i) = tend(k,j,i) & … … 588 579 ! 589 580 !-- Interpolate eddy diffusivities on staggered gridpoints 590 kmzp = 0.25_wp * &591 ( km(k,j,i)+km(k+1,j,i)+km(k,j-1,i)+km(k+1,j-1,i) )592 581 kmzm = 0.25_wp * & 593 582 ( km(k,j,i)+km(k-1,j,i)+km(k,j-1,i)+km(k-1,j-1,i) )
Note: See TracChangeset
for help on using the changeset viewer.