& +phyFac*dPhiHydY(i,j) )
_RL phxFac, phyFac
phyFac = pfFacMom
phyFac = 0.
& - phyFac*dPhiHydY(i,j)
_RL phxFac,phyFac, psFac
_RL phyFac
phyFac = pfFacMom*1.
gV(i,j,k,bi,bj) = vCf(i,j) - phyFac*dPhiHydY(i,j)
gV(i,j,k,bi,bj) = -phyFac*dPhiHydY(i,j)