t = nlread('analysis.nc', TransmissionSpectrum)[0]
t1 = nlread('analysis1.nc', TransmissionSpectrum)[0]
import pylab as P
P.plot(t.evaluate(),t.energies(), 'r-' )
P.plot(t1.evaluate(),t1.energies(), 'b-' )
P.xlabel("Transmission T(E)")
P.ylabel("Energy / eV")
P.savefig('transmission2.png', dpi=120)