Thanks Edoapra. So the numbers currently being output in the cube file are simply the square roots of the density at each point (times a factor of +/- 1)?
No, the correspondence is not so direct.
Here is the formula for the charge density ρ(r) in term of the orbitals ψa(r)