_RL flxFrac
flxFrac = iceFrac
IF (iceStart.EQ.1) flxFrac = 1.
dEnerg = dEnerg + thSIce_deltaT*flxFrac*(flxAtm-qleft)
dWater = dWater - thSIce_deltaT*flxFrac*(ffresh+frwAtm)
dSalt = dSalt - thSIce_deltaT*flxFrac*fsalt