To 1D acoustic-advection and Boussinesq, add a DIRK method of order 5 to standard_integrators. Coefficients can be found here: http://link.springer.com/article/10.1007%2FBF01990538#page-1