110 |
rUpwind = ABS(rCenter)*upwindFac(i,j) |
rUpwind = ABS(rCenter)*upwindFac(i,j) |
111 |
a3d(i,j,k) = a3d(i,j,k) |
a3d(i,j,k) = a3d(i,j,k) |
112 |
& - (rCenter+rUpwind)*deltaTarg(k) |
& - (rCenter+rUpwind)*deltaTarg(k) |
113 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
114 |
b3d(i,j,k) = b3d(i,j,k) |
b3d(i,j,k) = b3d(i,j,k) |
115 |
& - (rCenter-rUpwind)*deltaTarg(k) |
& - (rCenter-rUpwind)*deltaTarg(k) |
116 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *_recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
117 |
b3d(i,j,k-1) = b3d(i,j,k-1) |
b3d(i,j,k-1) = b3d(i,j,k-1) |
118 |
& + (rCenter+rUpwind)*deltaTarg(k-1) |
& + (rCenter+rUpwind)*deltaTarg(k-1) |
119 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
120 |
c3d(i,j,k-1) = c3d(i,j,k-1) |
c3d(i,j,k-1) = c3d(i,j,k-1) |
121 |
& + (rCenter-rUpwind)*deltaTarg(k-1) |
& + (rCenter-rUpwind)*deltaTarg(k-1) |
122 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *_recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
123 |
ENDDO |
ENDDO |
124 |
ENDDO |
ENDDO |
125 |
|
|