& thSIce_deltaT, ocean_deltaT, tauRelax_MxL,
_RL thSIce_deltaT, ocean_deltaT, tauRelax_MxL
ocean_deltaT = dTtracerLev(1)
WRITE(iUnit,*) 'ThSI: ocean_deltaT =',ocean_deltaT
frzmltMxL = (TFrzOce-oceTs)*cphm/ocean_deltaT
dtFac = ocean_deltaT/rhosw
fwFac = ocean_deltaT*sMxL_default*rhofw/rhosw
heatFac = ocean_deltaT/(cpwater*rhosw)
dtFac = ocean_deltaT/tauRelax_MxL