%% herat.m displays the relationship between F10.7 and the 
%% scaling factor required by the He+ at night, so that
%% GPID results agree with FLIP.
%%
%% NOTE - results for L=2.5 only so far.
%%

clear all

f107v=[166 189 219 157 223 218 217 278 251 89 93 94 82 121 116 122 128];
sfv=[1.1 1.1 1.1 1.2 1.1 1.1 1.1 1.3 1.3 1.5 1.4 1.4 1.4 1.2 1.0 1.1 1.3];

f107f=linspace(50,350);
sf=hesf(f107f,zeros(size(f107f)));

figure(1)
plot(f107v,sfv,'o',f107f,sf)
axis([50 350 1 1.7])
grid
xlabel('F10.7')
ylabel('sf')

