& +phxFac*dPhiHydX(i,j) )
_RL phxFac, phyFac
phxFac = pfFacMom
phxFac = 0.
& - phxFac*dPhiHydX(i,j)
_RL phxFac,phyFac, psFac
_RL phxFac
phxFac = pfFacMom*1.
gU(i,j,k,bi,bj) = uCf(i,j) - phxFac*dPhiHydX(i,j)
gU(i,j,k,bi,bj) = -phxFac*dPhiHydX(i,j)