effConduct=
A3(I,J)=4.0 _d +00*D3*t3 + effConduct + D1*UG(I,J)
B(I,J)=effConduct*(TB-tsurfLoc(I,J))
_RL effConduct
effConduct = XKI * XKS /
dFiDTs1 = 4.0 * D3*t3 + effConduct + D1*UG(I,J) +
F_c(I,J) = -effConduct * (TB - t1)