We carry out extensive pseudospectral direct numerical simulations (DNSs) to
show that our minimal model yields pulsar glitches with properties that are
akin to those seen in observations. In particular, the time series of the
angular momentum \(J_c\)  of the crust shows the hallmarks of self-organized
criticality
(SOC) [1]}, [2]}, [3]}, [4]}, [5]}, [6]}, [7]}, [8]}.
To obtain these results, we develop a sophisticated algorithm to find the
ground state of the GPPE, with rotation: this alogorithm uses an ancilliary
advective real Ginzburg-Landau-Poisson equation (ARGLPE), an imaginary-time
version of the GPPE. The resulting ground states contain vortices and
yield uniformly rotating solutions of the GPPE. When we include the crust and its
dynamics, we find a transfer of the angular momentum from the star to the
crust, where it is dissipated by friction; if this transfer is large enough,
vortices move outwards and glitches are observed.
