# SM file to plot the average power spectra of the fake spectra.
#
erase
window 1 1 1 1
ctype default
expand 1.1
lweight 2
ltype 0
ptype 4 3
location 4500 31000 3500 31000
notation -4 4 -5 -1
define TeX_strings 1
#
window 1 -2 1 1
limits -2.99 -.5 -4.5 -1.0
ticksize -1 10 0 0
box
xlabel k (s/km)
ylabel log_{10}[k P(k)/\pi]
expand 0.9
#
ctype default
relocate -2.8 -2.50 putlabel 6 128^3/ 256^3, 30Mpc/h
data "OS30A/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype blue
relocate -2.8 -2.75 putlabel 6 256^3/ 512^3, 30Mpc/h
data "OS30B/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype cyan
relocate -2.8 -3.00 putlabel 6 512^3/1024^3, 30Mpc/h
data "OS30C/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
#ctype magenta
#relocate -2.8 -3.25 putlabel 6 512^3/1024^3, 20Mpc/h
#data "OS20A/power_3.dat"
#read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
#ptype 4 0 points lk d2
#
#ctype red
#relocate -2.8 -3.50 putlabel 6 256^3/ 512^3, 15Mpc/h
#data "OS15B/power_3.dat"
#read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
#ptype 4 0 points lk d2
#
ctype default
expand 1.1
window 1 -2 1 2
limits -2.99 -.5 0.00001 0.05
ticksize -1 10 0 0
box 0 2 0 0
ylabel k P(k)/\pi
expand 0.9
#
relocate -2.8 0.045 putlabel 6  O&S, z=5
relocate -2.8 0.040 putlabel 6  Subrandom ICs
#
ctype default
data "OS30A/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype blue
data "OS30B/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype cyan
data "OS30C/power_nj_3.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
#ctype magenta
#data "OS20A/power_3.dat"
#read {k 1 d2 2} SET lk=LG(k)
#ptype 4 0 points lk d2
#
#ctype red
#data "OS15B/power_3.dat"
#read {k 1 d2 2} SET lk=LG(k)
#ptype 4 0 points lk d2
#