I would like to compute my diagnostics using this "full" EPflux not just using the ordinary quasi-geostrophic one. The results seem quite different so far.
Addtional to the paper above, you can have a feel about possible errors in estimatint the full EP flux and EP flux divergence based on reanalysis data sets from the followng paper:
Hua Lu, Thomas J. Bracegirdle, Tony Phillips, and John Turner, 2015: A Comparative Study of Wave Forcing Derived from the ERA-40 and ERA-Interim Reanalysis Datasets. J. Climate, 28, 2291–2311. doi: http://dx.doi.org/10.1175/JCLI-D-14-00356.1