set term pdf fontscale 0.5 size 4in,2.6in
outfile = 'plot_fit.pdf'
set output outfile
load 'define_fit.gnu'
#set term x11 0
#set format xy "$%g$"
set timestamp offset graph 1.172,0.7 rotate font ',6'
set label 'fixed_pitch/plot_fit.gnu' at graph 1.02,0.02 rotate left font ',6' noenhance
g=29.3
f(x) = -0.5*g**2/(x**2*(g**2-1)**2-g**2)*1.e6
h001= 0.0420528422585867 
  h0005=0.0106766052315869
  set key bottom right Left
  set xtics 1,1
  set xrange [0.75:4.75]
  set xlabel 'harmonic'
  set ylabel '{/Symbol Dw}(n)/{/Symbol Dw}(0) [ppm]'     # C_z(h)/C_z(0) [ppm]'
#  set title 'C_z(n,B_z), normalized by {C_z(0,B_z)'
# set label at 2.5, -270 'C_l (n=0,B_l) = (B_z)^2{/Symbol g}^2/2'
#set label at 2.5,-350 'C_l (n=0, B_z/B_y=0.005) = 0.0107'
#  set label at 2.5,-310 'C_l (n=0, B_z/B_y= 0.01) = 0.0421'
#set term x11 1;  
set timestamp offset graph 1.18,0.7 rotate font ',6'
set output 'long_dw_dn_d0.pdf'; plot 'fit_spin.dat'  u 2:($4*1.e6/h001) i 1:1 pt 6 ps 1.5 t '(b_l/B_z)=0.01', '' u 2:($4*1.e6/h0005) i 3:3 pt 7 t '(b_l/B_z)=0.005', f(x) t 'theory'
print 'long_dw_dn_d0.pdf'
set ylabel '{/Symbol Dw}(0)/{/Symbol w}_a [ppm]'
#set title 
set key bottom right Left
set autoscale y
#set term x11 2;  
set output 'long_dw_dn.pdf'; plot 'fit_spin.dat'  u 2:($4*1.e6) i 1:1 pt 6 ps 1.5 t '(b_l/B_z)=0.01', '' u 2:($4*1.e6) i 3:3 pt 7 t '(b_l/B_z)=0.005 '
print 'long_dw_dn.pdf'
set xlabel 'b_l/B_z [ppm]'
set xrange [0:100]
set xtics 0,10
unset timestamp
set timestamp offset graph 1.172,0.7 rotate font ',6'
set output 'long_dw_bl.pdf'; plot hfenergy(x)/1.e6 not
print 'long_dw_bl.pdf'