# SM file to plot the 3D \Delta^2(k) vs k.
#
window 1 1 1 1
expand 1.5
lweight 2
ptype 4 3
location 3500 31000 3500 24000
notation 1 1 1 1
define TeX_strings 1
#
ctype default
limits -1.0 2.0 -2.0 2.0
ticksize -1 10 -1 10
erase box
xlabel k (h/Mpc)
ylabel \Delta^2_{DM}(k)
#
relocate -0.8 1.50 putlabel 6  O&S, z=3
relocate -0.8 1.25 putlabel 6  Subrandom ICs
#
data pdnl_OS_z3.dat
read {k 1 d2 3} SET lk=LG(k) SET ld2=LG(d2)
connect lk ld2
#
ctype default
relocate 0.5 -0.25 putlabel 6 128^3/ 256^3 30Mpc/h
data "OS30A/pk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET ld2=LG(d2)
ptype 4 0 points lk ld2
#
ctype blue
relocate 0.5 -0.50 putlabel 6 256^3/ 512^3 30Mpc/h
data "OS30B/pk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET ld2=LG(d2)
ptype 4 0 points lk ld2
#
ctype cyan
relocate 0.5 -0.75 putlabel 6 512^3/1024^3 30Mpc/h
data "OS30C/pk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET ld2=LG(d2)
ptype 4 0 points lk ld2
#
ctype magenta
relocate 0.5 -1.00 putlabel 6 512^3/1024^3 20Mpc/h
data "OS20A/pk_1.dat"
read {k 1 d2 2} SET lk=LG(k) SET ld2=LG(d2)
ptype 4 0 points lk ld2
#
ctype red
relocate 0.5 -1.25 putlabel 6 256^3/ 512^3 15Mpc/h
data "OS15B/pk_1.dat"
read {k 1 d2 2} SET lk=LG(k) SET ld2=LG(d2)
ptype 4 0 points lk ld2
#