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