#plot quad efield and gradient computed in track1_preprocess 
set term pdf fontscale 0.75 size 6in,3in
#set term x11
outfile = 'Efield_vs_time_vs_quad.pdf'
set output outfile

set timestamp offset graph 1.18,0.7 rotate font 'Verdana,6'
set label 4 'rfquad/documentation/Efield_vs_time_vs_quad.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
!pwd > directory
!sed -n -e 's/^/name="/' -e 's/$/"/p' directory > name
load 'name'
print 'directory='.name
set label 1 at graph 0.0, 1.02 name font 'Verdana,6' noenhance
set xlabel 'time [{/Symbol m}s]'
set ylabel 'Horizontal electric field  [MV/m]'
quad1=1
quad2=2
quad3=3
quad4=4
icol=11
t1='Q1 long'
t2='Q2 long'
t3='Q3 long'
t4='Q4 long'
set pointsize 0.5
plot 'fort.91' u ($3*1.e6):($2 == quad1 ? column(icol):1/0) pt 4 t t1, ''u ($3*1.e6):($2 == quad2 ? column(icol):1/0) pt 5 t t2,
''u ($3*1.e6):($2 == quad3 ? column(icol):1/0) pt 6 t t3, ''u ($3*1.e6):($2 == quad4 ? column(icol):1/0) w pt 7 t t4
print outfile*