suzerain: #3127 output Dy1, Dy2, etc.

Tested with

import h5py
globals().update(h5py.File('foo.h5','r'))
Dy1[()].dot(bar_u.attrs["mu"])
Dy1[()].dot(bar_u.attrs["mu"]) - bar_u__y.attrs["mu"]
Dy2[()].dot(bar_u.attrs["mu"]) - bar_u__yy.attrs["mu"]

which shows O(1e-14) discrepancies for Dy1 and O(1e-10) for Dy2
relative to non-dense computation of bar_u__y, bar_u__yy.
I'll take it.

git-svn-id: https://svn.ices.utexas.edu/repos/pecos/suzerain/trunk@45397 0bd85f0b-c151-43b9-a1e9-27231104551f