Vacuum level only printed if VH or VT is written
This is a little unfortunate since I think it would be nice to have them printed out for all final dhscf calculations
- add flag to DHSCF to signal a final run
- always write out vacuum level in case we have a chain/slab calculation, does it make sense for molecules?
vacuum_levelto be passed the rhoatm array this will remove two forall loops that adds/removes
consider always having the spin factor for
rhopcc, nearly all places a factor of
rnsd == 1/spinoris used, these are superfluous multiplications and could easily be changed