pip install qe_rho

# import superposition of atomic density
from qe_rho import SAD

# initialize SAD density
rho = SAD('tests/pwscf.in')

# save charge-density.dat to folder pwscf.save
rho.saverhog('pwscf.save')

# output real-space charge density in a 3D numpy array
rhor = rho.rho_g2r()