GGm1=-1.0*sigmaR(i,j,k-1)
IF (GGm1.LE.SaltPlumeCriterion) THEN
tmpFac = (SaltPlumeCriterion - GGm1)
& / (GG - GGm1)
_RL GG, GGm1