Variational time discretization for free surface flows