The diluent gas of Bose-Einstein condensed atoms in a non-rotationally axisymmetric harmonic trap is modeled by the time-dependent Gross-Pitaevsky equation. If the angular momentum carried by the condensate does not go to zero, the minimum energy state represents a vortex (or anti-vortex) propagating around the trap center. The number of (anti-)vortices increases with angular momentum and repel each other to form an Abrikosov lattice. In addition to eddies and anti-eddies, there are also stagnation points where superflow disappears. To our knowledge, the stagnation point has not been previously analyzed in the context of the Gross-Pitaevsky equation. The Poincaré exponent formula shows that the difference between the number of eddies and stagnation points never changes. If the number of stagnation points is low, they tend to aggregate into degenerate propagating structures. However, when the number is large enough, the stagnation points tend to pair with the eddy core and propagate around the trap center in a regular lattice arrangement. There is an analogy with the geometry of the Kosterlitz-Thouless transition, where the angular momentum of the condensate is used as the external control parameter instead of temperature.
