# Unit cell volume in a m^3 unit
#( 2 appears because it's a BCC structure : 4 for FCC : 3 for HEX (?) )
#( unit cell means unit cell in the structure file )
#V=(17.080296*0.529*10.**(-10.))**3./2.
# Charge 
 C=1.602*10.**(-19.)
#plot
ef=0.59115
set xrange [-1.:1.5]
 ymin=-0.00015*10**6
 ymax= 0.00015*10**6
 set yrange [ymin:ymax]
# set arrow nohead from 0,ymin to 0,ymax
set yzeroaxis lt 0
set xzeroaxis lt 0

#set xlabel "{/Symbol dm} [eV]"
set xlabel "{/Symbol dm} [eV]"

set ylabel "S [{/Symbol m}V/K]"
set xtics 0.2
set mxtics 2
#set ytics 25
#set mytics 2


#!!!ef=ef/13.6


plot  "<awk '$2==1200' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "1200K - Sx"\
,"<awk '$2==600' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "600K - Sx"\
,"<awk '$2==300' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "300K - Sx"\
,"<awk '$2==120' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "120K - Sx"\
,"<awk '$2==60' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "60K - Sx"\
,"<awk '$2==20' FeSb2.condtens" u (($1-ef)*13.6):(10**6*$13) w l lw 2 ti "20K - Sx"\



set terminal postscript enhanced color
set output "FeSb2.Hyb.Sx_dn_T2.eps"
rep