& thSIce_deltaT, thSIce_dtTemp, ocean_deltaT,
_RL thSIce_deltaT, thSIce_dtTemp, ocean_deltaT
ocean_deltaT = dTtracerLev(1)
WRITE(iUnit,fmtR) 'ThSI:','ocean_deltaT =',ocean_deltaT
& ocean_deltaT, tauRelax_MxL, tauRelax_MxL_salt,
dtFac = ocean_deltaT/tauRelax_MxL
dtFac = ocean_deltaT/rhosw
fwFac = ocean_deltaT*sMxL_default/rhosw
heatFac = ocean_deltaT/(cpwater*rhosw)
dtFac = ocean_deltaT/tauRelax_MxL_salt
& * cphm/ocean_deltaT
ocean_deltaT = dtcouplo