from gpaw.transport.analysor import Transport_Plotter
import numpy as np
from pylab import *
import sys
plotter = Transport_Plotter()
plotter.plot_setup()
if len(sys.argv) <= 2:
if len(sys.argv[1]) <= 2:
nt = plotter.get_info('nty', int(sys.argv[1]), 0)
else:
tmp = sys.argv[1].split('-')
sam = int(tmp[0])
ref = int(tmp[1])
nt = plotter.get_info('nty', sam, 0) - plotter.get_info('nty', ref, 0)
else:
nt = plotter.get_info('nty', int(sys.argv[1]), int(sys.argv[2]))
matshow((nt[0]))
colorbar()
xlabel('Transport Direction')
ylabel('Pseudo Density')
show()
from gpaw.transport.analysor import Transport_Plotter
from pylab import *
plotter = Transport_Plotter('bias', 'Pt_H2_nsc.dat')
plotter.set_ele_steps(None, 0)
plotter.plot_ele_step_info('ham', [0], 0, 0)
