Variational structure of the Vlasov equation in multidimensional systems