set term pdf fontscale 0.5 size 6in,4in

file1='tbt_9.1584e-06_bl_bz_0._br_bz_2._0_-7000000._eta_0._radial_frac_2.'
file1='tbt_0._bl_bz_100.e-6_br_bz_0.0_1_-7000000._eta_0.0_radial_frac_0.0'
outfile = 'spin_vs_fields_long.pdf'
set output outfile

set timestamp 
set label 1 '$d9/g-2/examples/example10/edm/plotting_scripts/spin_vs_fields.gnu' at graph 0.6 ,-0.12 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 title ' '
set xrange [0:10]
set xlabel 'time [us]'
set ylabel 's_{long}'
set y2label 's_{vert}'
set yrange [-1:1]
set ytics -1,.2
set  y2range [-1.:1.]
set y2tics -1.,0.2
set label 3 at graph 0.05,0.95 '{/Symbol h}=10^{-6}'
set label 5 at graph 0.05,0.9 'd = 4.67 X 10^{-20} e-cm'
factory2=1000.
print log10(factory2)
#set label 1 at graph 1, 0 sprintf("× 10^{%d}",log10(factorX)) offset 0.5,0.3
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./10000.)) offset 0,0.7
plot 'tbt_0.0_bl_bz_0._br_bz_0._0_-7000000._eta_1.e-6' u ($14*1.e6):($11*factory2) w l axes x1y2 t 's_{vert}', '' u ($14*1.e6):12 w l dt 2 t 's_{long}'

factory2=10.
set label 3 at graph 0.05,0.95 '{/Symbol h}=1.158 X 10^{-4}'
set label 5 at graph 0.05,0.9 'd = 5.4 X 10^{-18} e-cm'
set label 4 at graph 1.04, .99 sprintf("× 10^{%i}",log10(1./factory2*.99999999)) offset 0,0.7
plot 'tbt_0.0_bl_bz_0._br_bz_0._0_-7000000._eta_1.15848e-4_radial_frac_1.'  u ($14*1.e6):($11*factory2) w l axes x1y2 t 'sy', '' u ($14*1.e6):12 w l dt 2 t 'sz'


print outfile