# SM file to plot the average delpk spectra of the 1D DM "spectra".
#
erase
window 1 1 1 1
ctype default
expand 1.5
lweight 2
ltype 0
ptype 4 3
location 4500 31000 3500 24000
notation -4 4 -4 4
define TeX_strings 1
#
limits -1.09 2.0 -2.0 2.00
ticksize -1 10 0 0
box
xlabel k (h/Mpc)
ylabel log_{10}[k P(k)/\pi]
#
data OS_z3_Kf025.dat
read {k 1 pk 5} SET lk=LG(k) SET pk=LG(pk)
connect lk pk
#
ctype cyan
relocate -0.9 1.5 putlabel 6 30Mpc/h
ctype red
relocate -0.9 1.2 putlabel 6 60Mpc/h
#
ctype cyan
data "OS30A/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 1 points lk d2
ctype default relocate -0.9 -1.00 dot putlabel 6  128^3/ 256^3
#
ctype cyan
data "OS30B/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
ctype default relocate -0.9 -1.25 dot putlabel 6  256^3/ 512^3
#
ctype cyan
data "OS30C/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 3 0 points lk d2
ctype default relocate -0.9 -1.50 dot putlabel 6  512^3/1024^3
#
ctype red
data "OS60A/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 1 points lk d2
#
ctype red
data "OS60B/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype red
data "OS60C/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 3 0 points lk d2
#