Demonstrate the J70 atmosphere model.

------------------------------------------------------------------------
See also NPlot, Plot2D, Date2JD, JD2DN, AtmJ70, SolarFluxPrediction,
GMSTime
------------------------------------------------------------------------

%-------------------------------------------------------------------------------
%   Copyright (c) 2000-2003 Princeton Satellite Systems, Inc.
%   All rights reserved.
%-------------------------------------------------------------------------------

longitude = linspace(0,360);
nPlot     = length(longitude);
rhoPlot   = zeros(5,nPlot);
degToMin  = 24*60/360;
altitude  = 600; % km

for j = 1:5
  year = j + 2001;
  jD   = Date2JD([year 3 16 0 0 0]);
  d    = SolarFluxPrediction( jD, 'nominal', true );
  d.z  = altitude;
  for k = 1:length(longitude)
    d.lng = longitude(k);
    rhoPlot(j,k) = AtmJ70( d );
  end
end

[~, hA] = Plot2D(longitude,rhoPlot*1000,'Longitude (deg)', 'Density (kg/m^3)',...
  sprintf('Density at %g km Altitude, Jacchia 1970',altitude),'ylog');

legend(hA.h,'2002','2003','2004','2005','2006');


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