& afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
_RL pfFacMom
psFac = pfFacMom*implicSurfPress
nhFac = pfFacMom*implicitNHPress
pfFacMom = 1. _d 0
pfFacMom = 0. _d 0
nhFac = pfFacMom*(1. _d 0 - implicitNHPress)
psFac = pfFacMom*(1. _d 0 - implicSurfPress)
phxFac = pfFacMom
phyFac = pfFacMom
tmpFac = pfFacMom*(1. _d 0 - implicitNHPress)