Variational Wasserstein gradient flow