MSIS2000 atmospheric model
[res1, res2, ...] = CL_mod_atmMSIS00(cjd, lon, lat, alt, f107, f107a, ap [[, res, opts, tt_tref]])
Computes various quantities using NRLMSISE-00 atmospheric model.
The quantities that can be computed are:
- dens: total density [kg/m^3]
- temp: temperature at current position [K]
- temp_exo: exospheric temperature [K]
- dens_part: partial densities of various species [kg/m3] in the order: He, O, N2, O2, Ar, H, N, O1 (O1=anomalous oxygen).
The necessary solar and geomagnetic data are:
- f107: daily F10.7 flux for previous day (sfu)
- f107a: 81 day centered average of F10.7 flux (sfu)
- ap: either daily Ap index, or an array with 7 rows containing:
> 1 : daily Ap
> 2 : 3 hr Ap index for current time
> 3 : 3 hr Ap index for 3 hrs before current time
> 4 : 3 hr Ap index for 6 hrs before current time
> 5 : 3 hr Ap index for 9 hrs before current time
> 6 : average of eight 3 hr Ap indices from 12 to 33 hrs prior to current time
> 7 : average of eight 3 hr Ap indices from 36 to 57 hrs prior to current time
By default, the computation is in "daily Ap mode" if only one value (one row) is passed. To force the computation to only use 7 identical Ap values, use the "opts" argument.
The "opts" argument is a structure that may contain the optional fields:
- nap: Number of ap values to be used: 1 (only the daily ap) or 7. Default is the number of ap values passed (= number of rows in ap).
- lst: Type of local time passed to the internal function: "mhl": mean local time or "tlh": true local time. Default is "tlh" (true local time).
- cv_utc: %t if the input date is internally converted to UTC, %f otherwise. Default is %f.
Notes:
- By default, the TREF time scale is supposed to be identical to UTC (or close enough). If necessary, the input date can be explicitly converted to UTC by setting the "cv_utc" field in the "opts" structure to true. The argument tt_tref is used only if cv_utc is true.
- The geodetic coordinates are given in a terrestrial frame (tied to the Earth).
Modified (1950.0) julian day (Time scale: TREF) (1x1 or 1xN)
Geodetic longitude [rad] (1x1 or 1xN)
Geodetic latitude [rad] (1x1 or 1xN)
Geodetic altitude [m] (1x1 or 1xN)
Solar flux at 10.7 cm for previous day [sfu] (1x1 or 1xN)
solar flux at 10.7 cm - average value [sfu] (1x1 or 1xN)
Geomagnetic indices [-] (1x1 or 7xN) - see above for details
(optional) Name(s) of the quantities to be computed. Possible names are: 'dens','dens_part','temp','temp_exo' or "all" (1 x Noutput). Default is "all"
(optional) Computation options - see above
(optional) TT-TREF [seconds]. Default is %CL_TT_TREF (1x1 or 1xN)
Computed quantities (Noutput variables)
CNES - DCT/SB
1) NRLMSISE-00 model, M. Picone, A.E. Hedin, D. Drob, Naval Research Laboratory
2) NRLMSISE-00, C language implementation, Dominik Brodowski, http://www.brodo.de/english/pub/nrlmsise/index.html
cjd = CL_dat_cal2cjd(2000,6,21); // TREF lon = 0; lat = 0; alt = 700.e3; f107 = 150; f107a = 150; ap = 15; // // Compute all quantities in a structure res = CL_mod_atmMSIS00(cjd, lon, lat, alt, f107, f107a, ap) // // Compute total density and temperature only [dens, temp] = CL_mod_atmMSIS00(cjd, lon, lat, alt, f107, f107a, ap, res=["dens", "temp"]) // // Compute total density for 2 longitudes dens = CL_mod_atmMSIS00(cjd, [lon, lon], lat, alt, f107, f107a, ap, res="dens") // // Compute total density using 7 ap values dens = CL_mod_atmMSIS00(cjd, lon, lat, alt, f107, f107a, ap, res="dens", opts=struct("nap", 7)) | ![]() | ![]() |