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