# Initial and final density profiles
#

location 4350 31000 3500 31000
erase
ex 1.4
lw 5
lt 0
mycolors
notation -2 2 -2 2

window 1 -2 1 1

ticksize -1 10 -1 10
lim 0 1.4771 -.69897 2

# M=1.e10, z=2
data Infall4.out
lines 22 172
read RoRv 5
read delta 7
set RoRv=LG(RoRv)
set delta=LG(1.+delta)
blue lt 0
connect RoRv delta

#-1sig
data Infall5.out
lines 22 172
read RoRv 5
read delta 7
set RoRv=LG(RoRv)
set delta=LG(1.+delta)
blue lt 0
connect RoRv delta

#+1sig
data Infall6.out
lines 22 172
read RoRv 5
read delta 7
set RoRv=LG(RoRv)
set delta=LG(1.+delta)
blue lt 0
connect RoRv delta

rel -.04576 0
green lt 1
draw 2 0

lw 4 blue
xlabel r/r_M or R/R_{vir}
ylabel "\rho(R)/\bar{\rho}"
lw 5
black
box 0 0 0 0

rel -.04 0
putl 4 1
rel -.04 1
putl 4 10
rel -.04 1.69897
putl 4 50

rel 0 -.95
putl 5 1
rel 1 -.95
putl 5 10
rel 1.4771 -.95
putl 5 30
rel .4771 -.95
putl 5 3

window 1 -2 1 2

ticksize -1 10 0 0
lim 0 1.4771 -.4 1.1

# M=1.e10, z=2
data Infall1.out
lines 20 155
read RoRv 1
read PSr 2
# ePSk plus ansatz:
set RoRv=LG(RoRv)
read ePSr 3
set ePSr = ePSr
blue lt 0
connect RoRv ePSr

data Infall2.out
lines 20 155
read RoRv 1
read PSr 2
# ePSk plus ansatz:
set RoRv=LG(RoRv)
read ePSr 3
set ePSr = ePSr
blue lt 0
connect RoRv ePSr

data Infall3.out
lines 20 155
read RoRv 1
read PSr 2
# ePSk plus ansatz:
set RoRv=LG(RoRv)
read ePSr 3
set ePSr = ePSr
blue lt 0
connect RoRv ePSr

rel -.04576 0
green lt 1
draw 1.4772 0

lw 4 blue
ylabel "\langle\delta(r)\rangle/\nu"
lw 5
black
box 0 2 0 0

location 3500 31000 3500 31000
window 1 1 1 1