Adds the diffusion constants in the output of plotND.py

also includes volume_err and improves the output format overall