# Read the DeviceConfiguration corresponding to different gate bias voltages
configurations = nlread('gatescan_conf.nc', DeviceConfiguration)

# Generate and save transmission spectrum
for device_configuration in configurations:
    transmission_spectrum = TransmissionSpectrum(
        configuration=device_configuration,
        energies=numpy.linspace(-0.5,0.5,401)*eV,
        kpoints=MonkhorstPackGrid(1,15),
        energy_zero_parameter=AverageFermiLevel,
        infinitesimal=1e-09*eV,
        self_energy_calculator=RecursionSelfEnergy(),
        )
    nlsave("gatescan_TS.nc", transmission_spectrum)
    nlprint(transmission_spectrum)