88 |
C-- Add centered & upwind contributions |
C-- Add centered & upwind contributions |
89 |
DO j=jMin,jMax |
DO j=jMin,jMax |
90 |
DO i=iMin,iMax |
DO i=iMin,iMax |
91 |
rCenter= 0.5 _d 0 *rTrans(i,j)*recip_rA(i,j,bi,bj)*rkFac |
rCenter= 0.5 _d 0 *rTrans(i,j)*recip_rA(i,j,bi,bj)*rkSign |
92 |
mskM = maskC(i,j,km2,bi,bj)*maskM2 |
mskM = maskC(i,j,km2,bi,bj)*maskM2 |
93 |
mskP = maskC(i,j,kp1,bi,bj)*maskP1 |
mskP = maskC(i,j,kp1,bi,bj)*maskP1 |
94 |
rC4km = oneSixth*rCenter*mskM |
rC4km = oneSixth*rCenter*mskM |
98 |
rU3km = 0. _d 0 |
rU3km = 0. _d 0 |
99 |
rU3kp = 0. _d 0 |
rU3kp = 0. _d 0 |
100 |
ELSE |
ELSE |
101 |
rU1k = oneSixth*abs(rCenter) |
rU1k = oneSixth*ABS(rCenter) |
102 |
rUpwind= rU1k+rU1k |
rUpwind= rU1k+rU1k |
103 |
rU3km = rU1k*mskM |
rU3km = rU1k*mskM |
104 |
rU3kp = rU1k*mskP |
rU3kp = rU1k*mskP |
105 |
ENDIF |
ENDIF |
106 |
a5d(i,j,k) = a5d(i,j,k) |
a5d(i,j,k) = a5d(i,j,k) |
107 |
& - (rC4km - rU3km) |
& + (rC4km + rU3km) |
108 |
& *deltaTarg(k) |
& *deltaTarg(k) |
109 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
110 |
b5d(i,j,k) = b5d(i,j,k) |
b5d(i,j,k) = b5d(i,j,k) |
111 |
& + (rCenter + rC4km - rUpwind - rU3km) |
& - (rCenter + rC4km + rUpwind + rU3km) |
112 |
& *deltaTarg(k) |
& *deltaTarg(k) |
113 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
114 |
c5d(i,j,k) = c5d(i,j,k) |
c5d(i,j,k) = c5d(i,j,k) |
115 |
& + (rCenter + rC4kp + rUpwind + rU3kp) |
& - (rCenter + rC4kp - rUpwind - rU3kp) |
116 |
& *deltaTarg(k) |
& *deltaTarg(k) |
117 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
118 |
d5d(i,j,k) = d5d(i,j,k) |
d5d(i,j,k) = d5d(i,j,k) |
119 |
& - (rC4kp + rU3kp) |
& + (rC4kp - rU3kp) |
120 |
& *deltaTarg(k) |
& *deltaTarg(k) |
121 |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
& *recip_hFacC(i,j,k,bi,bj)*recip_drF(k) |
122 |
b5d(i,j,k-1) = b5d(i,j,k-1) |
b5d(i,j,k-1) = b5d(i,j,k-1) |
123 |
& + (rC4km - rU3km) |
& - (rC4km + rU3km) |
124 |
& *deltaTarg(k-1) |
& *deltaTarg(k-1) |
125 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
126 |
c5d(i,j,k-1) = c5d(i,j,k-1) |
c5d(i,j,k-1) = c5d(i,j,k-1) |
127 |
& - (rCenter + rC4km - rUpwind - rU3km) |
& + (rCenter + rC4km + rUpwind + rU3km) |
128 |
& *deltaTarg(k-1) |
& *deltaTarg(k-1) |
129 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
130 |
d5d(i,j,k-1) = d5d(i,j,k-1) |
d5d(i,j,k-1) = d5d(i,j,k-1) |
131 |
& - (rCenter + rC4kp + rUpwind + rU3kp) |
& + (rCenter + rC4kp - rUpwind - rU3kp) |
132 |
& *deltaTarg(k-1) |
& *deltaTarg(k-1) |
133 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
134 |
e5d(i,j,k-1) = e5d(i,j,k-1) |
e5d(i,j,k-1) = e5d(i,j,k-1) |
135 |
& + (rC4kp + rU3kp) |
& - (rC4kp - rU3kp) |
136 |
& *deltaTarg(k-1) |
& *deltaTarg(k-1) |
137 |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
& *recip_hFacC(i,j,k-1,bi,bj)*recip_drF(k-1) |
138 |
ENDDO |
ENDDO |