rUpwind = ABS(rCenter)*upwindFac
& - (rCenter+rUpwind)
& + (rCenter+rUpwind)
& - (rCenter-rUpwind)
& + (rCenter-rUpwind)
_RL rCenter, rUpwind, upwindFac
rUpwind = ABS(rCenter)*upwindFac(i,j)
& - (rCenter+rUpwind)*deltaTarg(k)
& - (rCenter-rUpwind)*deltaTarg(k)
& + (rCenter+rUpwind)*deltaTarg(k-1)
& + (rCenter-rUpwind)*deltaTarg(k-1)
_RL rCenter, rUpwind
rUpwind= rU1k+rU1k
& - (rCenter + rC4km + rUpwind + rU3km)
& - (rCenter + rC4kp - rUpwind - rU3kp)
& + (rCenter + rC4km + rUpwind + rU3km)
& + (rCenter + rC4kp - rUpwind - rU3kp)
rUpwind= 0. _d 0