& rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF, rhoConstFresh,
_RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)
& *recip_deepFac2F(k)*recip_rhoFacF(k)
& *recip_rhoFacF(k)
recip_rhoFacF(k) = 1. _d 0/rhoFacF(k)
recip_rhoFacF(k) = 1. _d 0
& * wUnit2rVel(k)*wUnit2rVel(k)*recip_rhoFacF(k)