Aerospace Toolbox

atmosnrlmsise00

Implement mathematical representation of 2001 United States Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere

Syntax

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec) [T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, lst) [T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, f107a, f107, aph) [T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, flags) [T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, otype) [T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, action)

Description

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec) and implement the mathematical representation of the 2001 United States Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere (NRLMSISE-00). NRLMSISE-00 calculates the neutral atmosphere empirical model from the surface to lower exosphere (0 to 1,000,000 meters) with the option of including contributions from anomalous oxygen which can affect satellite drag above 500,000 meters.

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, lst) lets you specify an array of m local apparent solar time (hours).

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, f107a, f107, aph) lets you specify an array of m 81 day average of F10.7 flux (centered on doy), an array of m daily F10.7 flux for previous day, and an array of m-by-7 of magnetic index information.

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, flags) lets you specify an array of 23 to enable or disable particular variations for the outputs.

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, otype) lets you specify a string for total mass density output.

[T rho] = atmosnrlmsise00(h, lat, lon, year, doy, sec, action) lets you specify out-of-range input action.

Inputs for atmosnrlmsise00 are:

`h`	An array of `m` altitudes, in meters.
`lat`	An array of `m` geodetic latitudes, in meters.
`long`	An array of `m` longitudes, in degrees.
`year`	Depending on the value of `ctype`, this argument specifies one of the following arrays: `m` pressures in pascal `m` geopotential height in meters
`doy`	An array `m` day of year.
`sec`	An array of `m` seconds in day in universal time (UT)
`lst`	An array of `m` local apparent solar time (hours). To obtain a physically realistic value, lst is set to (`sec/3600 + lon/15`) by default. See Limitations for more information.
`F107a`	An array of `m` 81 day average of F10.7 flux (centered on day of year (`doy`)). If `F107A` is input, `F107` and `aph` must also be input. The effects of `F107A` are neither large nor well established below 80,000 meters; therefore, the default value is set to 150. See Limitations for more information.
`F107`	An array of `m` daily F10.7 flux for previous day. If `F107` is input, `F107A` and `aph` must also be input. The effects of F107 are neither large nor well established below 80,000 meters; therefore, the default value is set to 150. See Limitations for more information.
`aph`	An array of `m`-by-7 of magnetic index information. If `aph` is input, `F107a` and `F107` must also be input. This information consists of daily magnetic index (AP), 3 hour AP for current time, 3 hour AP for 3 hours before current time, 3 hour AP for 6 hours before current time, 3 hour AP for 9 hours before current time, average of eight 3 hour AP indices from 12 to 33 hours prior to current time, and average of eight 3 hour AP indices from 36 to 57 hours prior to current time. The effects of daily magnetic index are neither large nor well established below 80,000 meters. As a result, the default value is set to 4. See Limitations for more information.
`flags`	An array of 23 to enable or disable particular variations for the outputs. See following table.
`otype`	A string for total mass density output. `‘Oxygen'` Total mass density outputs include anomalous oxygen number density. `‘NoOxygen'` Total mass density outputs do not include anomalous oxygen number density.
`action`	A string to determine action for out-of-range input. Specify if out-of-range input invokes a `'Warning'`, `'Error'`, or no action (`'None'`). The default is `'Warning'`.

These flags, associated with the eleventh input, enable or disable particular variations for the outputs.

Field	Description
`Flags(1)`	F10.7 effect on mean
`Flags(2)`	Independent of time
`Flags(3)`	Symmetrical annual
`Flags(4)`	Symmetrical semiannual
`Flags(5)`	Asymmetrical annual
`Flags(6)`	Asymmetrical semiannual
`Flags(7)`	Diurnal
`Flags(8)`	Semidiurnal
`Flags(9)`	Daily AP. If you set this field to -1, the block uses the entire matrix of magnetic index information (APH) instead of `APH(:,1)`
`Flags(10)`	All UT, longitudinal effects
`Flags(11)`	Longitudinal
`Flags(12)`	UT and mixed UT, longitudinal
`Flags(13)`	Mixed AP, UT, longitudinal
`Flags(14)`	Terdiurnal
`Flags(15)`	Departures from diffusive equilibrium
`Flags(16)`	All exospheric temperature variations
`Flags(17)`	All variations from 120,000 meter temperature (TLB)
`Flags(18)`	All lower thermosphere (TN1) temperature variations
`Flags(19)`	All 120,000 meter gradient (S) variations
`Flags(20)`	All upper stratosphere (TN2) temperature variations
`Flags(21)`	All variations from 120,000 meter values (ZLB)
`Flags(22)`	All lower mesosphere temperature (TN3) variations
`Flags(23)`	Turbopause scale height variations

Outputs calculated for the neutral atmosphere empirical model are:

Field Description

T

An array of N-by-2 values of temperature, in kelvin. The first column is exospheric temperature in kelvin, the second column is temperature at altitude, in kelvin.

rho

An array of N-by-9 values of densities (kg/m³ or 1/m³) in selected density units. The column order is:

Density of He, in 1/m³

Density of O, in 1/m³

Density of N2, in 1/m³

Density of O2, in 1/m³

Density of Ar, in 1/m³

Total mass density, in 1/m³

Density of H, in 1/m³

Density of N, in 1/m³

Anomalous oxygen number density, in 1/m³

rho(6), total mass density, is defined as the sum of the mass densities of He, O, N2, O2, Ar, H, and N. Optionally, rho(6) can include the mass density of anomalous oxygen making rho(6), the effective total mass density for drag.

Examples

Calculate the temperatures, densities not including anomalous oxygen using the NRLMSISE-00 model at 10,000 meters, 45 degrees latitude, -50 degrees longitude, on January 4, 2007 at 0 UT using default values for flux, magnetic index data, and local solar time with out of range actions generating warnings:

[T rho] = atmosnrlmsise00( 10000, 45, -50, 2007, 4, 0)

Calculate the temperatures, densities not including anomalous oxygen using the NRLMSISE-00 model at 10,000 meters, 45 degrees latitude, -50 degrees longitude, and at 25,000 meters, 47 degrees latitude, -55 degrees longitude on January 4, 2007 at 0 UT using default values for flux, magnetic index data, and local solar time with out of range actions generating warnings:

[T rho] = atmosnrlmsise00( [10000; 25000], [45; 47], 
[-50; -55], [2007; 2007], [4; 4], [0; 0])

Calculate the temperatures, densities including anomalous oxygen using the NRLMSISE-00 model at 10,000 meters, 45 degrees latitude, -50 degrees longitude, on January 4, 2007 at 0 UT using default values for flux, magnetic index data, and local solar time with out of range actions generating errors:

[T rho] = atmosnrlmsise00( 10000, 45, -50, 2007, 
4, 0, 'Oxygen', 'Error')

Calculate the temperatures, densities including anomalous oxygen using the NRLMSISE-00 model at 100,000 meters, 45 degrees latitude, -50 degrees longitude, on January 4, 2007 at 0 UT using defined values for flux, and magnetic index data, and default local solar time with out of range actions generating no message:

aph = [17.375 15 20 15 27 (32+22+15+22+9+18+12+15)/8 (39+27+9+32+39+9+7+12)/8]
f107 = 87.7
nov_6days  = [ 78.6 78.2 82.4 85.5 85.0 84.1]
dec_31daymean = 84.5
jan_31daymean = 83.5
feb_13days = [ 89.9 90.3 87.3 83.7 83.0 81.9 82.0 78.4 76.7 75.9 74.7 73.6 72.7]
f107a = (sum(nov_6days) + sum(feb_13days) + (dec_31daymean + jan_31daymean)*31)/81
flags = ones(1,23)
flags(9) = -1
[T rho] = atmosnrlmsise00( 100000, 45, -50, 2007, 4, 0, f107a, f107,   
aph, flags, 'Oxygen', 'None')

Limitations

If flags array length, m, is 23 and all available inputs are not specified, this function assumes that flags is set.

This function has the limitations of the NRLMSISE-00 model. For more information, see the NRLMSISE-00 model documentation.

The NRLMSISE-00 model uses sec, lst, and lon independently. These arguments are not of equal importance for every situation. For the most physically realistic calculation, choose these three variables to be consistent by default:

 lst = sec/3600 + lon/15

Departures from this equation for lst can be included if available but are of minor importance.

The F107 and F107A values that are used to generate the model correspond to the 10.7 cm radio flux at the actual distance of the Earth from the Sun rather than the radio flux at 1 AU. The following site provides both classes of values: ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/FLUX/

References