# Read first transmission spectra from fil z-a-z-6-6.hdf5 and z-a-z-20-6.hdf5
transmission_spectrum6 = nlread("z-a-z-6-6.hdf5",TransmissionSpectrum)[0]
transmission_spectrum20 = nlread("z-a-z-20-6.hdf5",TransmissionSpectrum)[0]
#make list of temperatures
temperature_list=numpy.linspace(0,2000,21)*Kelvin
#make list for holding conductance
conductance_list6=numpy.zeros(len(temperature_list))
conductance_list20=numpy.zeros(len(temperature_list))
for i in range(len(temperature_list)):
    conductance_list6[i]=transmission_spectrum6.conductance(
        electrode_temperatures=(temperature_list[i],temperature_list[i]) )
    conductance_list20[i]=transmission_spectrum20.conductance(
        electrode_temperatures=(temperature_list[i],temperature_list[i]) )
#plot the conductance as function of temperature
import pylab
pylab.figure()
pylab.semilogy(temperature_list,conductance_list6,label="short junction")
pylab.semilogy(temperature_list,conductance_list20,label="long junction")
pylab.xlabel("Temperature (K)")
pylab.ylabel("Conductance (S)")
pylab.legend(loc="lower left")
pylab.grid(True)
pylab.show()