mktable4.m

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

function mktable4 (file,eph,tsat,azim,elev,mdbs,snrx,nabx,naby,nabz,nabt,cutoff)
%MKTABLE4: Create a table with results (reliability measures)
%
% This function create a table showing numerical values for azimuth 
% and elevation of the satellites, as well as the computed reliability
% measures
%
% Syntax:
%   function mktable4 (file,eph,tsat,azim,elev, ...
%                      mdbs,snrx,nabx,naby,nabz,nabt,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
%    mdbs   - Internal reliability, minimal detectable bias
%    snrx   - External reliability, norm of the influence of an
%             observational error the size of the mdb on all unknown 
%             parameters (square root)
%    nabx   - External reliability, effect of an error the size of the
%             mdb on the X-coordinate
%    naby   - External reliability, effect of an error the size of the
%             mdb on the Y-coordinate
%    nabz   - External reliability, effect of an error the size of the
%             mdb on the Z-coordinate
%    nabt   - External reliability, effect of an error the size of the
%             mdb on the time-parameter
%    cutoff - Cutoff elevation
%
% Output arguments:
%    none

% ----------------------------------------------------------------------
% File.....: mktable4.m
% Date.....: 12-OCT-2000
% Version..: 1.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 |    MDB |   SNRX |   NABX |   NABY |   NABZ |   NABT |\n');
    fprintf (fid,'|--------------------------|-----------|-----------|--------|--------|--------|--------|--------|--------\n');

    for j = ind;
      
      str = gpst2str (tsat(j));
      fprintf (fid,'| %24s | %9.2f | %9.2f | %6.2f | %6.2f | %6.2f | %6.2f | %6.2f | %6.2f |\n', ...
	       str,azim(i,j),elev(i,j),mdbs(i,j),snrx(i,j),nabx(i,j),naby(i,j),nabz(i,j),nabt(i,j));
    
    end;

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

  end;

end;

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

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

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