tmpSurf = 1.
tmpSurf = deltaTMom*deltaTfreesurf
tmpSurf = tmpSurf / (1. _d 0 + tmpSurf )
& *tmpSurf
_RL tmpSurf
dPhiNH(i,j,bi,bj) = ( tmpSurf*dPhiNH(i,j,bi,bj)
& ) )/(1. _d 0 + tmpSurf )
tmpSurf = ( tmpSurf*( etaN(i,j,bi,bj)-etaH(i,j,bi,bj) )
& )/(1. _d 0 + tmpSurf )
tmpSurf = etaN(i,j,bi,bj)-etaH(i,j,bi,bj)
& +freeSurfFac*tmpSurf
_RL tmpFac, tmpSurf