#set size 0.6,0.6

D=0.22 #exp. gap

kB=8.9*10.**(-5.)  #eV/K
kBJ=1.38*10.**(-23.) #J/K
e=1.602*10.**(-19.)

lws=3
lws2=5
set border lw 2

set xlabel "T [K]"
set ylabel "S [mK/eV]"
set key bottom
set mxtics 2
unset mytics

set xrange[0:100]
set yrange[-8:0]

set label 1 "FeAs_2" at 12,-2

p "Sun_FeAs2_B.dat" u 1:(-$2) w l lw lws2 ti "experiment (Sun et al)"


efT=0.62883 #Ry ?? 0.6388
def=0.0  #Ry


#rep "0.2/d0/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lw 1
#rep "0.2/d.0005/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lw 2

rep "0.2new/d.0002/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lw lws ti "theory (scissored GGA)"

###rep "0.2TEST/response.dat_1" u 1:(kB*10.**3.*$4/$3) w l lw 2
#rep "0.2/d.0010/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lw 1



#ziman type of S

de2=0.03 #0.02 #-0.02 #-0.3
dl2=0.5
#fit [14.5:50.] ((-(kBJ/e+(D/2)/x)*dl2 -(1./x)*de2 ) *10.**3) 'Sun_FeAs2_B.dat' u 1:(-$2) via dl2


rep (-(kBJ/e+(D/2)/x)*dl2 -(1./x)*de2 )   *10.**3. lw lws \
ti  "fit for 2-band model, Eq. (1)" 
#: {/Symbol de}=0.03, best fit : {/Symbol dl}=0.5"
#: {/Symbol de}=0.02, best fit : {/Symbol dl}=0.59"

de=0.035
dl=0.45
#rep (-(kBJ/e+(D/2)/x)*dl -(1./x)*de )   *10.**3. lw lws \
#ti "model : {/Symbol de}=0.035, best fit : {/Symbol dl}=0.45"


de=-0.1
dL=-0.1	

#fit [14.5:50.]  (-((5./2.*kBJ/e+(D/2)/x)*dL  -de/x))   'Sun_FeAs2_B.dat' u 1:(-$2) via de

rep -((5./2.*kBJ/e+(D/2)/x)*dL  -de/x) *10.**3. 



#set terminal postscript enhanced color
#set output "S_FeAs2.eps"
#rep