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.59115 
#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 100
set mytics 2


# 0.3/13.6=+0.0221


plot  "<awk '$1<0.59115+0.00005+0.0221 && $1>0.59115-0.00005+0.0221' FeSb2.condtens" u ($2):(10**6*$13) w p ps 3 lw 4 ti "{/Symbol dm}= +0.3eV - Sx"\
,"<awk '$1<0.59115+0.00005+0.0221 && $1>0.59115-0.00005+0.0221' FeSb2.condtens" u ($2):(10**6*$17) w p ps 2 lw 3 ti "{/Symbol dm}= +0.3eV - Sy"\
,"<awk '$1<0.59115+0.00005+0.0221 && $1>0.59115-0.00005+0.0221' FeSb2.condtens" u ($2):(10**6*$21) w p ps 1 lw 2 ti "{/Symbol dm}= +0.3eV - Sz"\

#,"<awk '$3<0.0008 && $3>-0.0008' FeSb2.condtens" u ($2):(10**6*$17) w p lw 2 ti "{/Symbol d}n_e= 0.00 - Sy"\
#,"<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"\

#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"\




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