_RL const1,const2
const2=SEAICE_lhFusion*SEAICE_rhoIce
& - (-1. _d 0 *epsilonHEFF*const2)
& const2 ) - (epsilonTemp*const1)