d1 <- read.table("tmp_1.out")
d2 <- read.table("tmp_2.out")

##d1 <- read.table("805-180_1.out")
##d2 <- read.table("805-180_2.out")






inds <- seq(201,500)
inds = c(inds, inds+(500*2), inds+(500*3), inds+(500*4), inds+(500*5))
#inds = c(1:3000)

sim.inds <- seq(300,300*5,300)


pdf("dist_hist.pdf")
par(mfcol=c(2,1))
#h <- hist(d1[inds,2], breaks = 300, main = "", xlab = "", xlim = c(3,15) )
#h <- hist(d2[inds,2], breaks = 300, main = "", xlab = "", xlim = c(3,15) )

plot(d1[inds,2], type='l')
lines(d2[inds,2], type='l', col="red")
abline(v=sim.inds, lty=2, col=4)

dev.off()