sstdegC = tsurfLoc(I,J) - TMELT
_RL QS1, D1, D1W, D3, TMELT
TMELT = celsius2K
_RL TMELT, XKI, XKS, HCUT, recip_HCUT, XIO
SurfMeltTemp = TMELT + SEAICE_wetAlbTemp
& 'ibi iter, TMELT ',myIter,TMELT
tsurfLoc(I,J) = MIN( tsurfLoc(I,J), TMELT )