forig='orig_DA.clo.txt'
fbump='test_DA.clo.txt'

descript='SRmin1'

array zvals[17];
zvals[1]=0.
zvals[2]=2.0582968
zvals[3]=2.902871
zvals[4]=3.682871
zvals[5]=4.302871
zvals[6]=4.847871
zvals[7]=5.952871
zvals[8]=6.907871
zvals[9]=7.862871
zvals[10]=8.817871
zvals[11]=9.772871
zvals[12]=10.727871
zvals[13]=11.642871
zvals[14]=12.262871
zvals[15]=13.032871
zvals[16]=13.652871
zvals[17]=16.375742

ff(x)=zvals[x+1]

set terminal postscript eps enhanced color lw 1 "Helvetica" 16
set border lw 0.75
set style line 1 lw 3 lt 1 lc rgb '#ee0101'
set style line 2 lw 3 lt 1 lc rgb 'dark-green'
set style line 3 lw 3 lt 1 lc rgb 'medium-blue'
set style line 4 lw 3 lt 1 lc rgb 'brown4'
set style line 5 lw 3 lt 1 lc rgb 'plum'
set style line 6 lw 3 lt 1 lc rgb 'light-green'
set style line 7 lw 3 lt 1 lc rgb 'skyblue'
set style line 8 lw 3 dt 1 lc rgb 'orange'
set style line 9 lw 3 dt 1 lc rgb 'gold'
set style line 10 lw 3 dt 1 lc rgb 'grey80'
set style line 11 lw 3 dt 2 lc rgb '#ee0101'
set style line 12 lw 3 dt 2 lc rgb 'dark-green'
set style line 13 lw 3 dt 2 lc rgb 'medium-blue'
set style line 14 lw 3 dt 2 lc rgb 'brown4'
set style line 15 lw 3 dt 2 lc rgb 'plum'
set style line 16 lw 3 dt 2 lc rgb 'light-green'
set style line 17 lw 3 dt 2 lc rgb 'skyblue'
set style line 18 lw 3 dt 2 lc rgb 'orange'
set style line 19 lw 3 dt 2 lc rgb 'gold'
set style line 20 lw 3 dt 2 lc rgb 'grey80'

set style increment user

unset title
set key bottom left
set size 1,1
set title font ",20"

set output 'closed.eps'
set multiplot
set xrange [0:*]
set size 1,0.5
set origin 0,0.5
set key top right
set xlabel 'Z (m)'
set ylabel 'X (mm)'
set xtics auto
set ytics auto
plot \
  forig u ($1):($2*1e3) w lp ls 1 pt 1 t 'original', \
  fbump u ($1):($2*1e3) w lp ls 2 pt 2 t 'with bump', \

set origin 0,0
set key top right
set xlabel 'Z (m)'
set ylabel 'Y (mm)'
set xtics auto
set ytics auto
plot \
  forig u ($1):($4*1e3) w lp ls 1 pt 1 t 'original', \
  fbump u ($1):($4*1e3) w lp ls 2 pt 2 t 'with bump', \

unset multiplot

unset title
unset origin
set size 1,1

set output 'closed_zoom.eps'
set multiplot
set xrange [0:*]
set size 1,0.5
set origin 0,0.5
set key top right
set xlabel 'Z (m)'
set ylabel 'X (mm)'
set xtics auto
set ytics auto
plot [0:20] \
  forig u ($1):($2*1e3) w lp ls 1 pt 1 t 'original', \
  fbump u ($1):($2*1e3) w lp ls 2 pt 2 t 'with bump', \

set origin 0,0
set key top right
set xlabel 'Z (m)'
set ylabel 'Y (mm)'
set xtics auto
set ytics auto
plot [0:20] \
  forig u ($1):($4*1e3) w lp ls 1 pt 1 t 'original', \
  fbump u ($1):($4*1e3) w lp ls 2 pt 2 t 'with bump', \

unset multiplot

unset title
unset origin
set size 1,1

