Hello,
I noticed a small issue when plotting molecular orbitals using the new cube format (a very useful addition, by the way!). If a system contains atoms that use an ECP, dscf and ridft do not take the ECP electrons into account when writing the cube header. For example, gold (79 electrons, 60-electron ECP) is written as "19" in the header. Due to this, visualization programs interpret the gold atoms as potassium atoms (19 electrons). This is not a big issue, as it can be easily remedied with a simple sed command wrapped in a shell script. Maybe it can be corrected in a future release...
