_RL fresh
fresh = (mwater0 - (rhos*hs + rhoi*hi))/dt
fresh = fresh + snowPr
fresh = fresh + (rhos*hs + rhoi*hi)/dt
& 'ThSI_CALC_TH: -2 : esurp,fresh=', esurp, fresh
fresh = fresh + chi*rhoi/dt
fresh = fresh - evap
& (mwater0-(rhos*hs+rhoi*hi))/dt,evap,fresh,fsalt(i,j)
fresh = fresh + (prcAtm(i,j)-snowPr)
fresh=fresh+extend/rqh*freshe
fresh=fresh+freshe
& iceFrac, qleft, fsalt(i,j), fresh
I qleft, fresh, fsalt,
frw2oc(i,j) = fresh