Tsf = tSrf(i,j)
I iceMask, hIce, hSnow, tSrf, ageSnw,
_RL tSrf (siLo:siHi,sjLo:sjHi)
_RL tSrf (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
U icFrac, hIce, hSnow, tSrf, qIc1, qIc2,
tSrf(i,j) = 0. _d 0
U tSrf, tIc1, tIc2, qIc1, qIc2,
tSrf(i,j) = tFrz(i,j)
U flxSW, tSrf, qIc1, qIc2,
& 'ThSI_SOLVE4T: k, Ts, Tice=',0,tSrf(i,j),tIc1(i,j),tIc2(i,j)
Tsf(i,j) = tSrf(i,j)
tSrf(i,j) = Tsf(i,j)