D=0.22 #exp. gap

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

#set size 0.6,0.6

lws=3
lws2=4
set border lw 2

set xlabel "T [K]"
set ylabel "S [mK/eV]"
set key bottom

set xrange[0:140]
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 2 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.02 #0.02 #-0.02 #-0.3
dl2=0.6
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 "model : {/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.02, best fit : {/Symbol dl}=0.59"


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