saltFreeze = (- bqe - SQRT(discrim))*recip_aqe
IF ( saltFreeze .LT. 0. _d 0 )
& saltFreeze = (- bqe + SQRT(discrim))*recip_aqe
thetaFreeze = saltFreeze * a0 + eps4
& * ( saltFreeze - sloc )/saltFreeze
& * saltFreeze * convertEmP2rUnit
_RL thetaFreeze, saltFreeze