Computes eclipses.

------------------------------------------------------------------------
See also Plot2D, GeoEclps
------------------------------------------------------------------------

%--------------------------------------------------------------------------
%   Copyright (c) 1994 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------
%   Since version 1.

% Number of increments in alpha

nAlpha     = 3600;

% Solar constant in lbf/in^2

pSolar     = 4.5e-6/(4.448*39.37^2);

degToRad   = pi/180;

% Sun out-of-plane angle (ß)

BetaSun   =  0*degToRad;

dsunAlpha = 0.12*pi/nAlpha;

cBeta     = cos(BetaSun);
sBeta     = sin(BetaSun);

nPSolar  = zeros(1,nAlpha);
eType    = zeros(1,nAlpha);
sunAlpha = zeros(1,nAlpha);

for k = 1:nAlpha;

  sunAlpha(k)           = (k-1)*dsunAlpha + 0.94*pi;

  [nPSolar(k),eType(k)] = GeoEclps(sunAlpha(k),BetaSun);

end

Plot2D(sunAlpha*180/pi,[nPSolar;eType],'Alpha (deg)',['Norm Flux   ';'Eclipse Type'])


%--------------------------------------