Am I missing something here? You have available the natural (or some other) orbitals and occupations, right? Surely it's a simple matter to write about ten lines of Fortran to read these in and calculate the total electron density at any point?