bqe = - eps1*eps6 -sLoc*a0*eps3 + eps3*eps7 + eps2
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
bqe = eps1*eps6
print *, 'ml-shelfice: discrim = ', discrim,aqe,bqe,cqe