# Final velocity profiles
#

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

ticksize -1 10 0 0
lim 0 2 -1.3 .1

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

rel -.04576 0
green lt 1
draw 2 0

rel 0 0
magenta lt 3
set Vhub=-(10.**RoRv)/(3.*3.1415926536)
connect RoRv Vhub

lw 4 blue
xlabel R/R_{vir}
ylabel "V_{pec}(R)/V_c"
lw 5
black
box 

location 3500 31000 3500 31000
window 1 1 1 1