453 |
& +maskC( i ,j,kp1,bi,bj)* |
& +maskC( i ,j,kp1,bi,bj)* |
454 |
& (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,bi,bj)) |
& (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,bi,bj)) |
455 |
& ) ) * 0.25 _d 0 |
& ) ) * 0.25 _d 0 |
456 |
tmp1k(i,j) = dyG(i,j,bi,bj)*drF(k)*hFacW(i,j,k,bi,bj) |
tmp1k(i,j) = dyG(i,j,bi,bj)*drF(k) |
457 |
|
& * _hFacW(i,j,k,bi,bj) |
458 |
& * tmp1k(i,j) * dTdz |
& * tmp1k(i,j) * dTdz |
459 |
ENDDO |
ENDDO |
460 |
ENDDO |
ENDDO |
572 |
& +maskC(i, j ,kp1,bi,bj)* |
& +maskC(i, j ,kp1,bi,bj)* |
573 |
& (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj)) |
& (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj)) |
574 |
& ) ) * 0.25 _d 0 |
& ) ) * 0.25 _d 0 |
575 |
tmp1k(i,j) = dxG(i,j,bi,bj)*drF(k)*hFacS(i,j,k,bi,bj) |
tmp1k(i,j) = dxG(i,j,bi,bj)*drF(k) |
576 |
|
& * _hFacS(i,j,k,bi,bj) |
577 |
& * tmp1k(i,j) * dTdz |
& * tmp1k(i,j) * dTdz |
578 |
ENDDO |
ENDDO |
579 |
ENDDO |
ENDDO |