# 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.60172
set xrange [-0.2:0.2]
 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 d}n_e"

set ylabel "S [{/Symbol m}V/K]"
set xtics 0.1
set mxtics 2
set ytics 20
set mytics 2



#plot 'FeSb2.trace' u (13.6*($1-ef)):(10**6*$5) every 10::0  w lp lw 2 ti '20K'\

#plot "<awk '$2==20' FeSb2.trace" u (13.6*($1-ef)):(10**6*$5) w l lw 2 ti "20K - <S>"\
#, "<awk '$2==200' FeSb2.trace" u (13.6*($1-ef)):(10**6*$5) w l lw 2 ti "200K - <S>"

plot  "<awk '$2==1200' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "1200K - <S>"\
,"<awk '$2==600' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "600K - <S>"\
,"<awk '$2==300' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "300K - <S>"\
,"<awk '$2==120' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "120K - <S>"\
,"<awk '$2==60' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "60K - <S>"\
,"<awk '$2==20' FeSb2.trace" u (-($3+0.01128)):(10**6*$5) w l lw 2 ti "20K - <S>"\



set terminal postscript enhanced color
set output "FeSb2.S_dn_T.eps"
rep