_RL saltice
& saltice, S_winton, mu_Tf,
msalt0 = rhoi*hi*saltice
fsalt = (msalt0 - rhoi*hi*saltice)/dt
salte=(rhoi*hi*saltice)/dt
dSalt = rhoi*iceHeight(i,j,bi,bj)*saltice
& - compact* rhoi*hIce*saltice
saltice = 4. _d 0
WRITE(iUnit,*) 'ThSI: saltice =',saltice
fsalt=-(rhoi*iceThick*saltice)*newIceFrac/deltaTice