#!/usr/bin/gnuplot

set terminal postscript eps color enhanced solid font 'Helvetica, 16pt' size 8.6cm, 6cm  dl .5

# set encoding utf8

# mean_offset = 
mean_std = 25.2
mean_offset_0 = 1
mean_offset_1 = -22.2
mean_offset_2 = 17.3

sd_offset_0 = 16.6
sd_offset_1 = 19.2
sd_offset_2 = 21.8

# point_sz = 0.7
# data_f = 'offset_rtt_18secperpt.dat'
data_0 = 'offset_10km.dat'
data_1 = 'offset_10km_1m.dat'
data_2 = 'offset_10km_10m.dat'

rtt_0 = 'RTT_10km.dat'
rtt_1 = 'RTT_10km_1m.dat'
rtt_2 = 'RTT_10km_10m.dat'

point_sz = 0.5

set output 'offsets.eps'

set multiplot

set lmargin at screen 0.13
set tmargin at screen 0.97
set bmargin at screen 0.4
set yrange [-100:120]
set xrange [-5:190]

# labels
set ylabel '{/Symbol d} - ~{/Symbol d&{.}}{.6/Symbol \55 } (ps)' offset 1.1, 0

set label 1 '(a)' at screen 0.15, 0.91
set label 2 '(b)' at screen 0.15, 0.32

# set xrange [-1:29.9]
set xtics auto
set xtics format ' '
set ytics 50
set ytics format '% 4g'

# set arrow 1 from mean_offset_0

plot \
0 notitle lc 'blue' lw 3, \
(mean_std) notitle lc 'blue' lw 2 dashtype 2, \
(- mean_std) notitle lc 'blue' lw 2 dashtype 2, \
data_0 using ($1/60):($2):($3) notitle with errorbars pt 6 ps point_sz lc rgb "black" , \
data_1 using ($1/60):($2):($3) notitle with errorbars pt 6 ps point_sz lc rgb "black" , \
data_2 using ($1/60):($2):($3) notitle with errorbars pt 6 ps point_sz lc rgb "black" 
# data_f using ($0):(($8-mean_offset)*1000) notitle with points pt 6 ps point_sz lc rgb "black"


set tmargin at screen 0.38
set bmargin at screen 0.14
set xtics 50
set xlabel "Elapsed Time (mins)" offset 0, 0.4
set ylabel '{/Symbol D}T - {/Symbol D}T_0 (ns)' offset 1.1, 0
unset ytics
set ytics 50
set yrange [0:120]
set xtics format '%g'

rtt_origin = 103300

# plot data_f using ($0):($3) notitle with points pt 6 ps point_sz lc rgb "red"
plot \
rtt_0 using ($1/60):($2-rtt_origin) notitle with points pt 6 ps point_sz lc rgb "red" , \
rtt_1 using ($1/60):($2-rtt_origin) notitle with points pt 6 ps point_sz lc rgb "red" , \
rtt_2 using ($1/60):($2-rtt_origin) notitle with points pt 6 ps point_sz lc rgb "red"