set print "fit_spin.dat" append 
file =  'tbt_0._bl_bz_0.0_br_bz_0._0' ; omega = 1.43934
harmonic = '0'
set title  'B_l/B_y = 0'
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega_ref = omega
print " Bl/Bz "," harmonic ","    omega","          (omega-omega_ref)/omega_ref "
print " 0       ","   0            ", omega_ref, (omega_ref-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 0;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega
file =  'tbt_0._bl_bz_0.01_br_bz_0._0'; omega=1.50100
harmonic = '0'
set title 'B_l/B_y = 0.01'
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega0 = omega
print " 0.01    ","   0            ", omega0, (omega0-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 1;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega
file =  'tbt_0._bl_bz_0.01_br_bz_0._1';harmonic = '1';omega=1.43931
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega1 = omega
print " 0.01    ","   1            ", omega1, (omega1-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 2;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega
file =  'tbt_0._bl_bz_0.01_br_bz_0._2';harmonic = '2'
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega2 = omega
print " 0.01    ","   2            ", omega2, (omega2-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 3;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega
file =  'tbt_0._bl_bz_0.01_br_bz_0._3';harmonic = '3'
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega3 = omega
print " 0.01    ","   3            ", omega3, (omega3-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 4;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega

file =  'tbt_0._bl_bz_0.01_br_bz_0._4';harmonic = '4'
fit f(x) file u ($14*1.e6):13 via a,omega,c,d
omega4 = omega
print " 0.01    ","   4            ", omega4, (omega4-omega_ref)/omega_ref
word_omega=sprintf("%8.7f",omega)
set term x11 5;plot file u ($14*1.e6):13 w l t 'data', f(x) t 'harmonic = '.harmonic.' {/Symbol w}_a = '.word_omega