  We numerically investigate the quench expansion dynamics of an initially
confined state in a two-dimensional Gross-Pitaevskii lattice in the presence of
external disorder. The expansion dynamics is conveniently described in the
control parameter space of the energy and norm densities. The expansion can
slow down substantially if the expected final state is a non-ergodic non-Gibbs
one, regardless of the disorder strength. Likewise stronger disorder delays
expansion. We compare our results with recent studies for quantum many body
quench experiments.
