mm_log10pi = -aa1 / t1 + aa2
mm_pi = 10**(mm_log10pi)
_RL mm_pi,mm_log10pi,dqhice_dTice
mm_pi = exp(mm_log10pi*lnTEN)