Redundant Zdensity output?

I noticed that the zdensity data is printed in OUTPUT and in ZDNDAT, with the addition on n(r) in OUTPUT.

Can we combine this into ZDNDAT only? It would help with #2081

# In OUTPUT
z density profiles:
  calculated using        1 configurations
  rho(r): H_
         r      rho         n(r)
    -49.1818  0.000000E+00  0.000000E+00
    -48.9627  0.000000E+00  0.000000E+00
    -48.7436  0.000000E+00  0.000000E+00
    -48.5245  0.000000E+00  0.000000E+00
    -48.3055  0.000000E+00  0.000000E+00
    -48.0864  0.000000E+00  0.000000E+00
# ZDNDAT
Chain polymer with   6 units and Carboxy Head Group
         5       450
H_
 -4.918176E+01  0.000000E+00
 -4.896269E+01  0.000000E+00
 -4.874361E+01  0.000000E+00
 -4.852454E+01  0.000000E+00
 -4.830547E+01  0.000000E+00
 -4.808640E+01  0.000000E+00