Compares atmosphere models

See also AtmNRLMSISE

Contents

%--------------------------------------------------------------------------
%   Copyright (c) 2021 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------
%   Since 2021.1
%--------------------------------------------------------------------------

NRL Model

input(1) = struct('year',0,'doy',0,'sec',0,'alt',0,'lst',0,'g_lat',0,'g_long',0,...
  'f107A',0,'f107',0,'ap',0,'ap_a',aph);

z     = zeros(1,24);
flags = struct('switches',z,'sw',z,'swc',z);

flags.switches(1) = 0;
flags.switches(10)=1; % Does not use daily ap

for i = 2:24
	flags.switches(i)=1;
end

alt = linspace(0,2000);
rho = zeros(1,length(alt));
input(1).ap_a = aph;
for k = 1:length(alt)
  input(1).alt = alt(k);
  flags.switches(1) = 0;
	output = AtmNRLMSISE( 'gtd7', [], input(1), flags );
  rho(k) = output.d(6);
end

Plot2D(alt,rho,'Altitude (km)','\rho (gm/cm^3)','Atmosphere','ylog')


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

Published with MATLAB® R2021a