#plot fields and closed orbit
field_file = 'fort.98'
closed_orbit_file = 'fort.97'
set term pdf
outfile = 'Bfields_closed_orbit.pdf'
set output outfile

set timestamp 
set label 1 '$d9/g-2/examples/example10/edm/plotting_scripts/fields_closed_orbit.gnu' at graph 0.5 ,-0.10 font 'Verdana,6' noenhance
!pwd > directory
!sed -n -e 's/^/name="/' -e 's/$/"/p' directory > name
load 'name'
print 'directory='.name
#set label 2 at graph 0.0,1.02 name font 'Verdana,10' noenhance
set label 2 at graph 0.3, -0.12 name font 'Verdana,8' noenhance


set xlabel 's [m]'
set ylabel '{/Symbol D}B/B [ppm]'
plot field_file u 3:($4*1.e6) pt 5 ps 0.5 t  'B_{radial}', '' u 3:($6*1.e6) pt 6 ps .5  t 'B_{longitudinal}'

set ylabel 'vertical displacement[mm]'
plot closed_orbit_file u 1:($4*1000) w l t 'vertical closed orbit'
print outfile