This is the QG model as was used for the Asselin et al. (2016ab,2018) and my PhD thesis. The transform method is used to compute horizontal derivatives with spectral accuracy. Time and vertical derivatives are approximated by second-order centered finite differences.