Restart from ground state and fix potential:
calc = GPAW('Cd_gs.gpw',
 nbands=16,
 fixdensity=True,
    symmetry='off',
 kpts={'path': 'GXWKL', 'npoints': 60},
 convergence={'bands': 8})
calc.get_potential_energy()