_RL k12, k32
k12 = 4. _d 0*kIce*kSnow / (kSnow*hi + 4. _d 0*kIce*hs)
& 'ThSI_SOLVE4T: flx0,df0dT,k12,D=', flx0,df0dT,k12,k12-df0dT
a1 = a10 - k12*df0dT / (k12-df0dT)
b1 = b10 - k12*(flx0-df0dT*Tsf) / (k12-df0dT)
dTsf = (flx0 + k12*(Tice(1)-Tsf)) / (k12-df0dT)
a1 = a10 + k12
fct = k12*(Tsf-Tice(1))
_RL k12
k12 = k1local*k2local
b = x2 + k1local*x + k12
& 2.0*diclocal*k12/b +
& 2.0*diclocal*k12*db/b2 -