%% Function makes array of H_proj 
%  5-min interval in the year
%  direction 0 - 355 deg with 5 deg step 
%%

function Hproj_arr = fi_arr(year);
  
  % load dist file
  load(['dist_' num2str(year) '.mat']);
  eval (['dist = dist_' num2str(year)]);
  
  %load time file
  load('yeartime.mat');
  
  
  % calculate fi array
  
    for fi=0:5:355
       Hproj_arr((fi+5)/5,1:105120) = hproj(yeartime(1,:),yeartime(2,:),fi, dist.x(1,1:105120), dist.y(1,1:105120));
    end % for fi   
 
