Birch gave an extremely efficient algorithm to compute a certain subspace of classical modular forms using the Hecke action on classes of ternary quadratic forms. We extend this method to compute all forms of non-square level using the spinor norm, and we exhibit an implementation that is very fast in practice. This is joint work with Jeffery Hein and Gonzalo Tornaria.