%% Function calculates Hproj
%
%
%%

% =======================================================
%  hour and minutes - UT time
%
%
%
% =======================================================


function H_proj = hproj(hour,minute,fi, dist_x, dist_y);
% !!!!!!!!!!!
  lon =291; % THL
% lon =106.9; % VOS
% !!!!!!!!!!!
 
 
    UT = hour.*15 + minute.*0.25; % TIME to DEG
    y = lon + fi + UT;
    H_proj = dist_x.*sin(y.*pi./180)-dist_y.*cos(y.*pi./180);
