The toolkit for high-order neutrino radiation hydrodynamics (thornado) aims to provide numerical solution methods with applications to relativistic astrophysical systems such as core-collapse supernovae and binary neutron star mergers.