Comment (by Roland Haas):

@Vaishak Prasad Turns out that you can also use:

iohdf5::out_var = "QuasiLocalMeasures::qlm_newman_penrose"

which handles output of 2D arrays fine (it just cannot create a 2D slice through a 3D grid function).

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2451/reg-hdf5-output-of-quasilocalmeasures