O F_io_net,F_ia_net,F_ia, IcePenetSWFlux,
F_ia (I,J) = 0.0 _d 0
F_ia(I,J) = F_lwd(I,J) + F_swi(I,J) + F_lwu(I,J) +
tsurfLoc(I,J) = tsurfLoc(I,J) - F_ia(I,J) / dFiDTs1
& F_io_net(I,J), F_ia_net(I,J), F_ia(I,J)
_RL F_ia (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
& i,j,F_ia(i,j),
F_ia (I,J) = 0.0
O F_io_net,F_ia_net,F_ia, QSWI,
print *,'sig: I,J,F_ia,F_ia_net',I,J,F_ia(I,J),