set terminal x11

# 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:1200]
 ymin=-0.00012*10**6
 ymax= 0.00008*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 "T [K]"

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





plot  "<awk '$3<0.0108 && $3>0.0092' FeSb2.condtens" u ($2):(10**6*$21) w p lw 2 ti "{/Symbol d}n_e=-0.01 - Sz"\
,"<awk '$3<0.0008 && $3>-0.0008' FeSb2.condtens" u ($2):(10**6*$21) w p lw 2 ti "{/Symbol d}n_e= 0.00 - Sz"\
,"<awk '$3<-0.0092 && $3>-0.0108' FeSb2.condtens" u ($2):(10**6*$21) w p lw 2 ti "{/Symbol d}n_e=+0.01 - Sz"\




#,"<awk '$2==600' FeSb2.condtens" u ($5):(10**6*$21) w l lw 2 ti "600K - Sz"\
#,"<awk '$2==300' FeSb2.condtens" u ($5):(10**6*$21) w l lw 2 ti "300K - Sz"\
#,"<awk '$2==120' FeSb2.condtens" u ($5):(10**6*$21) w l lw 2 ti "120K - Sz"\
#,"<awk '$2==60' FeSb2.condtens" u ($5):(10**6*$21) w l lw 2 ti "60K - Sz"\
#,"<awk '$2==20' FeSb2.condtens" u ($5):(10**6*$21) w l lw 2 ti "20K - Sz"\



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