discrim = bqe*bqe - 4. _d 0*aqe*cqe
saltFreeze = (- bqe - SQRT(discrim))*recip_aqe
& saltFreeze = (- bqe + SQRT(discrim))*recip_aqe
_RL aqe, bqe, cqe, discrim, recip_aqe
IF ( discrim .LT. 0. _d 0 ) THEN
print *, 'ml-shelfice: discrim = ', discrim,aqe,bqe,cqe