out=0

reset
clear
if (out==1) 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=5
lws2=7
set border lw 2

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

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

set label 1 "scissored GGA" at 55,-4.65

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

rep "0.2newnew/Si0p/d0/response.dat_1" u 1:(kB*10.**3.*$4/$3) w l lt 1 lw lws ti "(-1)xS, Im{/Symbol S}=cnst, {/Symbol d}n=0"

rep "0.2newnew/Si0p/d.00005/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt 2 lw lws ti "Im{/Symbol S}=cnst, {/Symbol d}n=  50ppm"

rep "0.2newnew/Si0p/d.0001/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l  lt 3 lw lws ti "Im{/Symbol S}=cnst, {/Symbol d}n=100ppm"

rep "0.2newnew/Si0p/d.0002/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt 4  lw lws ti "Im{/Symbol S}=cnst, {/Symbol d}n=200ppm"

#rep "0.2newnew/Si0p/d.0003/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt 5 lw lws ti "Im{/Symbol S}=cnst, {/Symbol d}n=300ppm"

rep 1000 w l ti ""

#rep "0.2newnew/w2-008003/d.0002/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt -1 lw 7 ti "TTT"

a=0.085
fit [14.5:50.] (-a/x *10.**3) 'Sun_FeAs2_B.dat' u 1:(-$2) via a

rep -a/x*10.**3. w l lt 0 lw 4 \
ti  "fit : 2-band model" 


rep "Sun_FeAs2_B.dat" u 1:(-$2) w l lt 7  lw lws2 ti "Experiment (Sun et al.)"

rep "0.2newnew/Si0p/d.0001/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt 3 lw lws ti ""


if (out==1) set terminal postscript enhanced color
if (out==1) set output "S_FeAs2.eps"
if (out==1) rep