mktable1.m

Sky lab calculate position of satellite in GNSS systems and calculate elevation and azimuth them and

function mktable1 (file,eph,tsat,azim,elev,xsat,ysat,zsat,cutoff)
%MKTABLE: Create a table with results
%
% The function create a table showing numerical values for azimuth 
% and elevation of the satellites, as well as the positions of the
% satellites in earth-fixed WGS'84 coordinates.
%
% Syntax:
%    mktable1 (file,eph,tsat,azim,elev,xsat,ysat,zsat,cutoff)
%
% Input arguments:
%    file   - Filename of file to be created
%             empty string ==> write to screen
%    eph    - Ephemerides
%    tsat   - Times on which positions are computed
%    azim   - Azimuth of satellites
%    elev   - Elevation for each satellite
%    xsat   - X-coordinates in earth-fixed WGS'84 coordinates
%    ysat   - Y-coordinates in earth-fixed WGS'84 coordinates
%    zsat   - Z-coordinates in earth-fixed WGS'84 coordinates
%    cutoff - Cutoff-elevation
%
% Output arguments:
%    none

% ----------------------------------------------------------------------
% File.....: mktable1.m
% Date.....: 12-OCT-2000
% Version..: 2.0
% Author...: Peter Joosten
%            Mathematical Geodesy and Positioning
%            Delft University of Technology
% ----------------------------------------------------------------------

% -----------------
% --- Open file ---
% -----------------

if length(file) ~= 0;
  fid = fopen (file,'wt'); 
else;
  fid = 1;
end;

% ---------------------------------
% --- Write satellite positions ---
% ---------------------------------

for i = 1:size(eph,1);

  ind = find (elev(i,:) >= cutoff);

  if length(ind) > 0;
    
    fprintf (fid,'-----------------------------\n');
    fprintf (fid,'--- Satellite number: %3d ---\n',eph(i,1));
    fprintf (fid,'-----------------------------\n\n');
    fprintf (fid,'|--------------------------|-----------|-----------|\n');
    fprintf (fid,'|          Date/time [UTC] |   Azimuth | Elevation |\n');
    fprintf (fid,'|--------------------------|-----------|-----------|\n');

    for j = ind;
      
      str = gpst2str (tsat(j));
      fprintf (fid,'| %24s | %9.2f | %9.2f |\n',str,azim(i,j),elev(i,j));
    
    end;

    fprintf (fid,'|--------------------------|-----------|-----------|\n\n\n');

  end;

end;

% ------------------
% --- Close file ---
% ------------------

if fid ~= 1;
  fclose (fid);
end;

% ----------------------
% --- End of routine ---
% ----------------------