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=7
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"



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 3.5 \
ti  "phenomenological 1/T fit" 

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

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