The output of the code 'script_used.py' is 'density.dat' with the coordinates in the units of the box length.
The density parameter 'rs' must be used to get the box length for the number of electrons in the main box.  
The density is in units of Angstrom^{-3}.