_RL rec_pLaw
rec_pLaw = powerLaw
rec_pLaw = 1. _d 0 / rec_pLaw
c1Mlt = fracEnMelt**rec_pLaw
c2Mlt = (1. _d 0 - fracEnMelt)**rec_pLaw
c1Frz = fracEnFreez**rec_pLaw
c2Frz = (1. _d 0 -fracEnFreez)**rec_pLaw