phxFac = pfFacMom
phyFac = pfFacMom
hxFac = pfFacMom
hyFac = pfFacMom
hx3dFac = pfFacMom
hy3dFac = pfFacMom
psFac = pfFacMom*(1. _d 0 - implicSurfPress)
pfFacMom = 1.D0
pfFacMom = 0.D0
& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
_RL pfFacMom
phxFac = pfFacMom*1.
phyFac = pfFacMom*1.