out=1

reset
clear
#if (out==1) 
set size 0.6,0.6

D=0.2 #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/V]"
set key bottom
set mxtics 2
unset mytics

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

#set label 1 "scissored GGA" at 45,-4.65

p "Sun_FeAs2_B.dat" u 1:(-$2) w l lt 1 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 2 lw lws ti "Theory : Im{/Symbol S}=cnst"

#rep "0.2newnew/Si0p/d0/inp.efshift_list" u 1:(-10.**6.*($2+0.1)/$1) w l lt 1 lw lws ti "(-1)xS, Im{/Symbol S}=cnst"

#rep "0.2newnew/Si0p/d0-0.005/response.dat_1" u 1:(-kB*10.**3.*$4/$3) w l lt 2 lw lws ti "Im{/Symbol S}=cnst, {/Symbol de}= 5meV"

#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 "Sun_FeAs2_B.dat" u 1:(-$2) w l lt 1  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 ""

#rep -a/x*10.**3. w l lt 3 lw 4 \
#ti  "phenomenological 1/T fit" 

rep (-D/(x)*10.**3.) w l lt 3 lw 2 ti "electronic limit : {/Symbol D}/T"  #mV/K !!!


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