I qleft, ffresh, fsalt, flxAtm, frwAtm,
_RL qleft, fsalt, ffresh
& -ffresh-frwAtm,-fsalt
dWater = dWater - thSIce_deltaT*flxFrac*(ffresh+frwAtm)