mktable3.m

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

function mktable3 (file,tsat,eph,elev,pdop,cutoff)
%MKTABLE: Create a table with results
%
% The function create a table showing the satellite constellation,
% as a function of time. Each record in the table represents a
% change in the constellation.
%
% Syntax:
%    mktable3 (file,tsat,eph,elev,pdop,cutoff)
%
% Input arguments:
%    file   - Filename of file to be created
%             empty string ==> write to screen
%    tsat   - Times on which positions are computed
%    eph    - Ephemerides
%    elev   - Elevation for each satellite
%    pdop   - Computed pdops
%    cutoff - Select cutoff-elevation
%
% Output arguments:
%    none

% ----------------------------------------------------------------------
% File.....: mktable3.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 ---
% ---------------------------------

prevind = [];

for i = 1:length(tsat);

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

  if ~isequal (prevind,ind);
    
    pline   = char(32*ones(1,80));
    prevind = ind;
    
    for j = 1:length(ind);

      tmp  = num2str(eph(ind(j),1));
      pline ((j-1)*3+4-length(tmp):j*3) = tmp;
      
    end;

    tmp           = gpst2str(tsat(i));
    pline (40:47) = tmp(13:20);
    pline (50:54) = sprintf ('%5.2f',pdop(i));
    fprintf (fid,'%s\n',pline);
    
  end;
  
end;

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

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

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