calc = GPAW('Al_gs.gpw',
            nbands=16,
            fixdensity=True,
            symmetry='off',
            kpts={'path': 'GXWKL', 'npoints': 60},
            convergence={'bands': 8})
calc.get_potential_energy()
bs = calc.band_structure()
bs.plot(filename='bandstructure.png', show=True, emax=10.0)