_RL k12
k12 = k1local*k2local
b = x2 + k1local*x + k12
& 2.0*diclocal*k12/b +
& 2.0*diclocal*k12*db/b2 -
_RL k12, k32 ! thermal conductivity terms
k12 = 4. _d 0*kice*ksnow / (ksnow*hi + 4. _d 0*kice*hs)
& 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))