Download or view StandardAtmosphere.frink in plain text format
// Implementation of U.S. Standard Atmosphere, 1976 in Frink.
//
// Alan Eliasen, eliasen@mindspring.com
//
// The official book of the Standard Atmosphere is available in PDF format
// at:
// https://ntrs.nasa.gov/citations/19770009539
//
// See: http://www.atmosculator.com/
//
// This is the same as the ICAO Standard Atmosphere as far as the ICAO
// atmosphere is defined, but the U.S. Standard Atmosphere extends much higher.
// The ICAO Standard Atmosphere extends up to 86 km, while the U.S. Standard
// Atmosphere extends up to 1000 km. This model defines up to 1000 km.
//
// This model defines the atmosphere from 0 to 86 km to be in 7 layers with
// different (constant) temperature lapse rates in each layer.
//
// The class contains all class-level methods and variables, so you do *not*
// need to create an instance of this class.
//
// You will generally want to start with a call to
// [t,p] = StandardAtmosphere.getTemperatureAndPressure[alt]
// and then call other functions with that temperature and pressure if
// necessary.
//
// See the file StandardAtmosphereTest.frink for more samples.
//
// Section numbers and equations refer to the book
// _U.S. Standard Atmosphere, 1976_ (available at the first URL above.)
// A fundamental basic of the Standard Atmosphere is that the atmosphere
// below 86 km is made up of 7 different layers, each of them having a
// different (constant) temperature lapse rate, that is, the rate at which
// temperature varies with (geopotential) altitude.
class AtmosphereLayer
{
var h_base // The base geopotential altitude of the layer
var lambda // The temperature lapse rate in this layer
new[base, lapse] :=
{
h_base = base
lambda = lapse
}
}
// This is the class that contains the methods of the standard atmosphere
// model.
class StandardAtmosphere
{
// An array of AtmosphereLayer objects.
class var layers = StandardAtmosphere.initializeLayers[]
// Standard temperature at sea level.
class var standardTemp = Celsius[15]
// Standard pressure at sea level
class var standardPressure = 760 mmHg // Model says 1.013250e5 N/m^2
// Mean molecular weight of air, M0
class var M0 = 0.0289644 kg/mol
// The radius of the earth as defined in the Standard Atmosphere 1976
class var r0 = 6356766 m
// The Standard Atmosphere 1976 actually defines a different gas constant
// than is currently known to be good. It is also inconsistent with the
// values they use for Boltzmann's constant or Avogadro's number, which
// are also inconsistent with current best-known values. Sigh. For
// verifiability, and adherence to the model, this file uses the Standard
// Atmosphere's definition of the constants. (The best-known gas constant
// is called "gasconstant" or "R" in Frink.)
class var Rstar = 8.31432e+3 N*m/(kmol*K)
// This is the mean effective collision diameter sigma of molecules. This
// is no longer valid above 85 km
class var meanSigma = 3.65e-10 m
// This is the definition of Avogadro's number as defined in the Standard
// Atmosphere 1976. Note that it is not quite the same as the current
// best-known value for Avogadro's number, called "avogadro" in Frink.
class var N_avogadro = 6.022169e23 / mol
// The temperature in the layer starting at 86 km (Z7) to 91 km. See eq. 25
class var T7 = 186.8673 K
// The column where the geometric altitude resides.
class var geometricAltCol = 0
// The column where the N2 number density resides.
class var N2col = 2
// The column where the O number density resides.
class var Ocol = 3
// The column where the O2 number density resides.
class var O2col = 4
// The column where the Ar number density resides.
class var Arcol = 5
// The column where the He number density resides.
class var Hecol = 6
// The column where the H number density resides.
class var Hcol = 7
// The column where the sum of numeric densities resides.
class var sumcol = 8
class var N2mass = 28.0134 g/mol
class var O2mass = 31.9988 g/mol
class var Omass = 1/2 O2mass
class var Armass = 39.948 g/mol
class var Hemass = 4.0026 g/mol
class var Hmass = 1.00794 g/mol
// Initialize the layers.
class initializeLayers[] :=
{
layers = new array[8]
layers@0 = new AtmosphereLayer[ 0 km, -6.5 K/km]
layers@1 = new AtmosphereLayer[11 km, 0 K/km]
layers@2 = new AtmosphereLayer[20 km, +1.0 K/km]
layers@3 = new AtmosphereLayer[32 km, +2.8 K/km]
layers@4 = new AtmosphereLayer[47 km, 0 K/km]
layers@5 = new AtmosphereLayer[51 km, -2.8 K/km]
layers@6 = new AtmosphereLayer[71 km, -2.0 K/km]
layers@7 = new AtmosphereLayer[84.8520 km, undef]
return layers
}
// Get the atmospheric conditions for the specified geometric altitude.
//
// parameters:
// [geometricHeight, seaLevelTemperature, seaLevelPressure]
// The latter two arguments have defaults at the defaults for the
// Standard Atmosphere model.
//
// returns:
// [temperature, pressure]
class getTemperatureAndPressure[geometricHeight,
T0 = StandardAtmosphere.standardTemp,
P0 = StandardAtmosphere.standardPressure] :=
{
if (geometricHeight > 86 km)
return getHighAltitudeTemperatureAndPressure[geometricHeight]
geopotentialHeight = geometricToGeopotentialHeight[geometricHeight]
return getTemperatureAndPressureGeopotential[geopotentialHeight, T0, P0]
}
/** Convenience method to get the temperature, pressure *and* density at a given
geometric height. If you need to call further methods, you will need to save
the temperature and pressure.
parameters:
[geometricHeight, seaLevelTemperature, seaLevelPressure]
The latter two arguments have defaults at the defaults for the
Standard Atmosphere model.
returns:
[temperature, pressure, density]
*/
class getTPD[geometricHeight,
T0 = StandardAtmosphere.standardTemp,
P0 = StandardAtmosphere.standardPressure] :=
{
[temp, pressure] = getTemperatureAndPressure[geometricHeight, T0, P0]
density = getDensity[geometricHeight, temp, pressure]
return [temp, pressure, density]
}
// Get the atmospheric conditions for the specified *geopotential* altitude.
// Note that the geopotential altitude is a somewhat fictitious quantity
// and is *not* the true geometrical altitude as would be measured by
// a very long tape measure. Users should generally use the
// getTemperatureAndPressure[] function instead.
//
// parameters:
// [geopotentialHeight, temperature, pressure]
// The latter two arguments have defaults at the defaults for the
// Standard Atmposphere model.
//
// returns:
// [temperature, pressure]
class getTemperatureAndPressureGeopotential[geopotentialHeight, T0 = StandardAtmosphere.standardTemp, P0 = StandardAtmosphere.standardPressure] :=
{
if (geopotentialHeight > 84.8520 km)
return getHighAltitudeTemperatureAndPressure[geopotentialToGeometricHeight[geopotentialHeight]]
layerNum = 0
lambda = 1
theta = 1
delta = 1
T = T0
do
{
thetaPrev = theta
lambdaLayer = layers@layerNum.lambda
base = layers@layerNum.h_base
nextBase = layers@(layerNum+1).h_base
if (geopotentialHeight > nextBase) // In next layer?
h = nextBase
else
h = geopotentialHeight
theta = theta + (h-base) lambdaLayer/T0
// println["Theta: $theta"]
T = T0 * theta
// println["T: $T"]
// Isothermic layer? If so, use a different equation
if (lambdaLayer == 0 K/km)
delta = delta e^(-(h-base) gravity M0/Rstar/T)
else
delta = delta (theta/thetaPrev)^(-gravity M0/lambdaLayer/Rstar)
// println["Delta: $delta"]
layerNum = layerNum + 1
} while (geopotentialHeight > nextBase)
p = P0 * delta
// rho = p M0 / (Rstar T)
return[T, p]
}
// Returns the high altitude temperature and pressure at the specified
// geometric height.
// Returns: [temperature, pressure]
class getHighAltitudeTemperatureAndPressure[geometricHeight] :=
{
// Z7=86 km to Z8=91 km
if geometricHeight < 91 km
{
// TODO: Adjust for different base temperature?
temp = T7 // Constant (eq. 25)
return [temp, getHighAltitudePressure[geometricHeight, temp]]
}
if geometricHeight < 110 km
{
// Eq. 27
Tc = 263.1905 K
A = -76.3232 K
a = -19.9429 km
temp = Tc + A * sqrt[1 - ((geometricHeight - 91 km)/a)^2]
return [temp, getHighAltitudePressure[geometricHeight, temp]]
}
// Layer Z9=110 km to Z10=120 km
if geometricHeight < 120 km
{
// Eq. 29
T9 = 240.0 K
LK9 = 12.0 K/km
temp = T9 + LK9 (geometricHeight - 110 km)
return [temp, getHighAltitudePressure[geometricHeight, temp]]
}
// Layer Z10 = 120 km to Z12=1000 km
if geometricHeight <= 1000 km
{
// Eq. 31
Tinf = 1000 K
T10 = 360.0 K
lambda = 0.01875 km^-1
Z10 = 120 km
xi = (geometricHeight - Z10) * (r0 + Z10) / (r0 + geometricHeight)
temp = Tinf - (Tinf - T10) e^(-lambda xi)
return [temp, getHighAltitudePressure[geometricHeight, temp]]
}
println["Error: Standard Atmosphere 1976 is only defined to altitude of 1000 km, requested altitude of " + (geometricHeight -> "km")]
return [undef, undef]
}
// Convert a geometric height to a geopotential height.
// The geometric height is the normal height that you'd measure, and
// the geopotential height is a somewhat fictitious height which is used
// in the Standard Atmosphere.
// This is equation 18 in the Standard Atmosphere.
class geometricToGeopotentialHeight[Z] := (r0 * Z) / (r0 + Z)
// Convert a geopotential height to a geometric height.
// The geometric height is the normal height that you'd measure, and
// the geopotential height is a somewhat fictitious height which is used
// in the Standard Atmosphere.
// This is equation 19 in the Standard Atmosphere.
class geopotentialToGeometricHeight[H] := (r0 * H) / (r0 - H)
// The model usually assumes a standard temperature of 288.15 K (59 F) at
// sea level, but obviously that's not always very applicable. This
// function calculates the new base temperature T0 to feed to the model
// given the temperature at a specified *geometric* (true) altitude,
// assuming that the lapse rate remains the same.
class getBaseTemperature[geometricAltitude, temp] :=
{
// Temperature in standard model at that altitude
[standardT] = getTemperatureAndPressure[geometricAltitude]
return standardTemp + temp - standardT
}
// The model usually assumes a standard temperature of 288.15 K (59 F) at
// sea level, but obviously that's not always very applicable. This
// function calculates the conditions at a particular geometric elevation,
// specified by the argument targetAlt, given the temperature at another
// point, whose geometric altitude and temperature are given by givenAlt
// and givenTemp respectively, assuming that the lapse rate remains the
// same as in the standard model.
//
// returns:
// [temperature, pressure]
class getTemperatureAndPressureGiven[targetAlt, givenAlt, givenTemp] :=
{
baseT = getBaseTemperature[givenAlt, givenTemp]
return getTemperatureAndPressure[targetAlt, baseT]
}
// Calculates the mass density. This requires calculating partial pressures
// of each gas, which varies above 86 km, so you have to pass it the
// geometric height as well.
// See Section 1.3.4 and equation 42.
class getDensity[geometricHeight, temp, pressure] :=
{
if (geometricHeight < 86 km)
return pressure * M0 / (Rstar * temp)
else
{
// Do a binary search on the first column in numberData, which is the
// geometric height.
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
if (index == 0)
return pressure * M0 / (Rstar * temp)
else
return getMolarMass[geometricHeight, index] / N_avogadro
}
}
// Returns the mass of a mole of gas at the specified height.
class getMolarMass[geometricHeight] :=
{
if geometricHeight <= 86 km
return M0
else
{
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
return getMolarMass[geometricHeight, index] / logarithmicInterpolation[geometricHeight, index, sumcol]
}
}
// Get molecular weight. This is the sum of all gas masses multiplied by
// their respective number densities.
class getMolarMass[geometricHeight, index] :=
{
sum = logarithmicInterpolation[geometricHeight, index, N2col] * N2mass
sum = sum + logarithmicInterpolation[geometricHeight, index, Ocol] * Omass
sum = sum + logarithmicInterpolation[geometricHeight, index, O2col] * O2mass
sum = sum + logarithmicInterpolation[geometricHeight, index, Arcol] * Armass
sum = sum + logarithmicInterpolation[geometricHeight, index, Hecol] * Hemass
if geometricHeight > 150 km
sum = sum + logarithmicInterpolation[geometricHeight, index, Hcol] * Hmass
return sum
}
// Calculates the mole volume from temperature and pressure.
// See section 1.3.5
class getMoleVolume[geometricHeight, temp, pressure] :=
{
if geometricHeight <= 86 km
return Rstar * temp / pressure
else
{
// Do a binary search on the first column in numberData, which is the
// geometric height.
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
return N_avogadro / getNumberDensity[geometricHeight, temp, pressure]
}
}
// Calculates the number density from temperature and pressure.
// The number density is the number of neutral atmospheric gas particles
// per volume of the atmosphere.
// See Section 1.3.3
class getNumberDensity[geometricHeight, temp, pressure] :=
{
if geometricHeight <= 86 km
return N_avogadro * pressure / (Rstar * temp)
else
{
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
return logarithmicInterpolation[geometricHeight, index, sumcol]
}
}
// Calculates the mean air-particle speed from temperature and pressure.
// Since this requires knowledge of each gas species above 86 km, it
// needs to be passed the geometricHeight.
// See section 1.3.7
class getMeanAirParticleSpeed[geometricHeight, temp, pressure] :=
{
if geometricHeight <= 86 km
return sqrt[8 Rstar temp / (pi M0)]
else
{
// Do a binary search on the first column in numberData, which is the
// geometric height.
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
return sqrt[8 Rstar temp logarithmicInterpolation[geometricHeight, index, sumcol] / (pi getMolarMass[geometricHeight, index])]
}
}
// Calculates the mean free path from temperature and pressure.
// See section 1.3.8
class getMeanFreePath[geometricHeight, temp, pressure] :=
{
if geometricHeight <= 86 km
return sqrt[2] Rstar temp / (2 pi N_avogadro meanSigma^2 pressure)
else
return sqrt[2]/(2 pi meanSigma^2 getNumberDensity[geometricHeight, temp, pressure])
}
// Calculate the mean collision frequency from temperature and pressure.
// See section 1.3.9
class getMeanCollisionFrequency[geometricHeight, temp, pressure] :=
{
if geometricHeight <= 86 km
return 4 N_avogadro * meanSigma^2 * sqrt[pi * pressure^2 / (Rstar M0 temp)]
else
{
// Do a binary search on the first column in numberData, which is the
// geometric height.
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
return 4 N_avogadro * meanSigma^2 * sqrt[pi * pressure^2 * logarithmicInterpolation[geometricHeight, index, sumcol] / (Rstar temp * getMolarMass[geometricHeight, index])]
}
}
// Calculate the speed of sound from temperature and pressure.
// See section 1.3.10
class getSpeedOfSound[geometricHeight, temp, pressure] :=
{
if (geometricHeight > 86 km)
{
println["Speed of sound cannot be calculated over 86 km."]
return undef
}
return sqrt[1.40 Rstar temp / M0]
}
// Calculate the dynamic viscosity from temperature and pressure.
// This should not be used above 86 km.
// See section 1.3.11
class getDynamicViscosity[geometricHeight, temp, pressure] :=
{
if (geometricHeight > 86 km)
{
println["Dynamic viscosity cannot be calculated over 86 km."]
return undef
}
// Frink can now handle fractional exponents.
beta = 1.458e-6 kg/(s m K^(1/2))
// Sutherland's constant S
S = 110.4 K
return beta * temp^(3/2) / (temp + S)
}
// Calculate the kinematic viscosity from temperature and pressure.
// This should not be used above 86 km.
// See section 1.3.12
class getKinematicViscosity[geometricHeight, temp, pressure] :=
{
if (geometricHeight > 86 km)
{
println["Kinematic viscosity cannot be calculated over 86 km."]
return undef
}
getDynamicViscosity[geometricHeight, temp, pressure] / getDensity[geometricHeight, temp, pressure]
}
// Gets the high altitude (> 86 km) pressure.
class getHighAltitudePressure[geometricHeight, temp] :=
{
if (geometricHeight < 86 km || geometricHeight > 1000 km)
{
println["Error: StandardAtmosphere.getHighAltitudePressure called with an inappropriate altitude of " + (geometricHeight + "km") ". Valid values are between 86 km and 1000 km."]
return undef
}
// Do a binary search on the first column in numberData, which is the
// geometric height.
index = binarySearch[numberData, geometricHeight, {|a,b| a <=> b@0}]
if (index == 0)
interpval = numberData@index@sumcol
else
interpval = logarithmicInterpolation[geometricHeight, index, sumcol]
return interpval * k * temp // Eq. 33c
}
// This performs a logarithmic interpolation from values in the upper
// atmosphere table, searching on the geometricHeight field.
class logarithmicInterpolation[geometricHeight, index, columnNum] :=
{
// Do a logarithmic interpolation.
lowerSum = ln[numberData@(index-1)@columnNum / m^-3]
upperSum = ln[numberData@index@columnNum / m^-3]
lowerAlt = numberData@(index-1)@geometricAltCol
upperAlt = numberData@index@geometricAltCol
frac = (geometricHeight - lowerAlt) / (upperAlt - lowerAlt)
interp = lowerSum + frac(upperSum - lowerSum)
interpval = exp[interp] m^-3
return interpval
}
// Calculates the number density of atomic Nitrogen (N2) at heights above
// 86 km (geometric height)
class getHighAltitudeN2NumberDensity[geometricHeight, temp] :=
{
// From table 9, number density of N2 at the 86 km level.
nN2_7 = 1.129794e+20 / m^3
if geometricHeight < 86 km
{
println["Error: StandardAtmosphere.getHighAltitudeN2NumberDensity called with a too-low altitude of $geometricHeight"]
return undef
}
// Temperature is constant from 86 km to 91 km (eq. 25)
// This simplifies the integral in eq. 38.
if geometricHeight < 91 km
{
integral = gravity M0 r0^2 (geometricHeight - 86 km) / (Rstar T7) / ((r0+geometricHeight)(r0+86 km))
// println["Integral is $integral"]
return nN2_7 * e^-integral
}
if geometricHeight < 100 km
{
Tc = 263.1905 K
A = -76.3232 K // Ooh, negative Kelvin!
a = -19.9429 km
Z8 = 91 km // Geometric
Z = geometricHeight
satc = sqrt[-A^2 + Tc^2]
return (E^((gravity*M0*(a*(satc*arcsin[(5 km)/a] - satc*arcsin[(Z8 - Z)/a] +
Tc*(arctan[(5 km*A)/(a*satc)] -a^2*arctan[(5 km*Tc)/(sqrt[-25 km^2 + a^2]*satc)] +
arctan[(A*(Z - Z8))/(a*satc)] km^2 -
a^2*arctan[(Tc*(Z - Z8))/(satc*sqrt[a^2 - (Z - Z8)^2])])))/(A*satc))/Rstar)*nN2_7*T7)/ (Tc + A*sqrt[1 - (Z - Z8)^2/a^2])
}
}
// This is an algorithm for doing binary search of a list in Frink.
//
// arguments:
// list: an ordered list of items
//
// item: the item to search for in the list
//
// orderFunction: a two-argument function of the form {|a,b| ... } where
// the function returns -1 if a<b, 0 if a==b, and 1 if a > b
// (such as the <=> operator returns. This is the default
// comparison if no order function is passed in.)
//
// This returns the index of the item to insert *before* (for example,
// using the method list.insert[pos, item] ) to put the item in the right
// order. If the item is found anywhere in the list, this returns the index
// of the matching item in the list. If the item is not found in the
// list, this still returns the index before which it should be inserted.
class binarySearch[list, item, orderFunction] :=
{
start = 0
end = length[list] - 1
var current
var comp
while (start <= end)
{
current = floor[(end + start) / 2]
comp = orderFunction[item, list@current]
if comp == 0
return current
if (comp == -1)
end = current - 1
else
start = current + 1
}
return start
}
// Calculate the equivalent standard altitude that corresponds to a given
// pressure.
class getPressureAltitude[pressure, T0 = StandardAtmosphere.standardTemp, P0 = StandardAtmosphere.standardPressure] :=
{
return secantInvert[{|alt, data|
[t0, p0] = data
getTemperatureAndPressure[alt, t0, p0]@1
}, pressure, 0 km, 86 km, 1e-8, [T0, P0]]
}
// Calculate the equivalent standard altitude that corresponds to a given
// density.
class getDensityAltitude[density, T0 = StandardAtmosphere.standardTemp, P0 = StandardAtmosphere.standardPressure] :=
{
return secantInvert[{|alt, data|
[t0, p0] = data
[t,p] = getTemperatureAndPressure[alt, t0, p0]
getDensity[alt, t, p]
}, density, 0 km, 86 km, 1e-8, [T0, P0]]
}
/** This calculates corrections to "station" barometric pressures to give
"true" barometric pressures or vice-versa.
When most weather sites give you an atmospheric pressure reading, this
pressure is a "corrected" pressure that represents what the pressure would
be if taken at sea level at a specific temperature. This is rather useless
for local physical calculations.
The following functions convert between the "true" "station" pressure as
would be measured by an actual device at that altitude and the "corrected"
sea-level pressure that meteorologists like to publish.
*/
/** Calculate the "station" pressure given the geometric elevation and,
optionally, the base pressure and temperature to correct from. By
default, the temperature and pressure are the standard temperature and
pressure of the Standard Atmosphere.
While some references will give you a simple equation based on a single
assumption, (that is, a fixed pressure lapse rate and maybe a
temperature lapse rate or no temperature lapse rate, and some hidden
base temperature,) those assumptions
are totally hidden and probably invalid over more than one layer of
the Standard Atmosphere model.
This uses the full direct calculations of the Standard Atmosphere model
to find the pressure in any atmospheric layer.
arguments:
elevation: The station elevation above sea level
P0: The pressure at sea level (default is standard pressure)
T0: The temperature at sea level (default is standard temperature)
*/
class stationPressure[elevation is length, P0 is pressure = StandardAtmosphere.standardPressure, T0 is temperature = StandardAtmosphere.standardTemp] :=
{
[t,p] = getTemperatureAndPressure[elevation, T0, P0]
return p
}
/** Calculate the "corrected" sea level pressure given an elevation and a
"true" "station" pressure reading as read at the station.
While some references will give you a simple equation based on a single
assumption, (that is, a fixed pressure lapse rate and maybe a
temperature lapse rate or no temperature lapse rate, and some hidden
base temperature,) those assumptions
are totally hidden and probably invalid over more than one layer of
the Standard Atmosphere model.
This solution is very powerful in that it inverts all of the
calculations of the full Standard Atmosphere model (using the secant
method to invert a function) so it works over multiple atmospheric
layers and the results are completely compatible with Standard
Atmosphere calculations.
Arguments:
[elevation, stationPressure]
where
elevation is the geometric height above sea level
stationPressure is the true pressure as read at that location
returns:
pressure at sea level, assuming standardTemp
*/
class seaLevelPressure[elevation is length, stationPressure is pressure] :=
{
return secantInvert[{|p0, data|
[alt, press] = data
StandardAtmosphere.getTemperatureAndPressure[alt, StandardAtmosphere.standardTemp, p0]@1
}, stationPressure, 0 bar, 1.777 bar, 1e-8, [elevation, stationPressure]]
}
// This uses the secant method to invert the function y = f[x].
// This will essentially find an inverse function for f[x] and return a value
// of x for which f[x] = y.
// other parameters:
// x1,x2: initial guesses that hopefully bound the desired result.
// maxDelta: maximum error in y
// TODO: Use interval techniques to make this more rigorous and powerful?
// TODO: Automatically make guesses for x1 and x2? Somehow?
class secantInvert[f, y, xmin, xmax, maxDelta, data=undef] :=
{
x1 = xmin
x2 = xmax
y1 = f[x1,data]
y2 = f[x2,data]
xnew = (x2-x1)/2 + x1
while true
{
ydiff = y2 - y1
if ydiff == 0 y // Degenerate case to avoid dividing by zero.
return xnew // This may not be always a correct solution?
invSlope = (x2-x1) / ydiff
xnew = x1 + (y - y1) invSlope
if xnew < xmin
xnew = xmin
if xnew > xmax
xnew = xmax
ynew = f[xnew,data]
// println["xnew=$xnew\tynew=$ynew"]
if y == 0 y // Degenerate case to avoid dividing by zero.
return xnew // This may not be always a correct solution?
if abs[(ynew - y) / y] < maxDelta
return xnew
y2 = y1
y1 = ynew
x2 = x1
x1 = xnew
}
}
// The following is the tabulation of upper-atmosphere number data for
// each major element.
// Columns are:
// Z (geometric height), H (geopotential height), N2, O, O2, Ar, He, H, sum
class var numberData=[[86000 m,84852 m,1.130e20 m^-3,8.600e16 m^-3,3.031e19 m^-3,1.351e18 m^-3,7.582e14 m^-3,0 m^-3,1.4475e+20 m^-3],
[86500 m,85339 m,1.034e20 m^-3,9.939e16 m^-3,2.772e19 m^-3,1.236e18 m^-3,6.976e14 m^-3,0 m^-3,1.3246e+20 m^-3],
[87000 m,85825 m,9.456e19 m^-3,1.147e17 m^-3,2.535e19 m^-3,1.130e18 m^-3,6.422e14 m^-3,0 m^-3,1.2116e+20 m^-3],
[87500 m,86312 m,8.651e19 m^-3,1.320e17 m^-3,2.319e19 m^-3,1.033e18 m^-3,5.915e14 m^-3,0 m^-3,1.1086e+20 m^-3],
[88000 m,86798 m,7.915e19 m^-3,1.513e17 m^-3,2.120e19 m^-3,9.437e17 m^-3,5.453e14 m^-3,0 m^-3,1.0144e+20 m^-3],
[88500 m,87285 m,7.242e19 m^-3,1.724e17 m^-3,1.938e19 m^-3,8.624e17 m^-3,5.031e14 m^-3,0 m^-3,9.2835e+19 m^-3],
[89000 m,87771 m,6.626e19 m^-3,1.952e17 m^-3,1.772e19 m^-3,7.880e17 m^-3,4.647e14 m^-3,0 m^-3,8.4963e+19 m^-3],
[89500 m,88257 m,6.062e19 m^-3,2.193e17 m^-3,1.619e19 m^-3,7.198e17 m^-3,4.296e14 m^-3,0 m^-3,7.7749e+19 m^-3],
[90000 m,88744 m,5.547e19 m^-3,2.443e17 m^-3,1.479e19 m^-3,6.574e17 m^-3,3.976e14 m^-3,0 m^-3,7.1161e+19 m^-3],
[90500 m,89230 m,5.075e19 m^-3,2.699e17 m^-3,1.351e19 m^-3,6.002e17 m^-3,3.685e14 m^-3,0 m^-3,6.5130e+19 m^-3],
[91000 m,89716 m,4.643e19 m^-3,2.953e17 m^-3,1.234e19 m^-3,5.478e17 m^-3,3.419e14 m^-3,0 m^-3,5.9613e+19 m^-3],
[91500 m,90202 m,4.248e19 m^-3,3.200e17 m^-3,1.126e19 m^-3,4.998e17 m^-3,3.177e14 m^-3,0 m^-3,5.4560e+19 m^-3],
[92000 m,90688 m,3.886e19 m^-3,3.434e17 m^-3,1.027e19 m^-3,4.557e17 m^-3,2.956e14 m^-3,0 m^-3,4.9929e+19 m^-3],
[92500 m,91173 m,3.553e19 m^-3,3.651e17 m^-3,9.361e18 m^-3,4.152e17 m^-3,2.753e14 m^-3,0 m^-3,4.5671e+19 m^-3],
[93000 m,91659 m,3.249e19 m^-3,3.846e17 m^-3,8.527e18 m^-3,3.781e17 m^-3,2.568e14 m^-3,0 m^-3,4.1780e+19 m^-3],
[93500 m,92145 m,2.970e19 m^-3,4.016e17 m^-3,7.761e18 m^-3,3.441e17 m^-3,2.399e14 m^-3,0 m^-3,3.8207e+19 m^-3],
[94000 m,92630 m,2.715e19 m^-3,4.159e17 m^-3,7.060e18 m^-3,3.129e17 m^-3,2.244e14 m^-3,0 m^-3,3.4939e+19 m^-3],
[94500 m,93116 m,2.481e19 m^-3,4.275e17 m^-3,6.418e18 m^-3,2.844e17 m^-3,2.103e14 m^-3,0 m^-3,3.1940e+19 m^-3],
[95000 m,93601 m,2.268e19 m^-3,4.365e17 m^-3,5.830e18 m^-3,2.583e17 m^-3,1.973e14 m^-3,0 m^-3,2.9205e+19 m^-3],
[95500 m,94087 m,2.072e19 m^-3,4.429e17 m^-3,5.293e18 m^-3,2.345e17 m^-3,1.854e14 m^-3,0 m^-3,2.6691e+19 m^-3],
[96000 m,94572 m,1.894e19 m^-3,4.471e17 m^-3,4.801e18 m^-3,2.127e17 m^-3,1.745e14 m^-3,0 m^-3,2.4401e+19 m^-3],
[96500 m,95057 m,1.730e19 m^-3,4.493e17 m^-3,4.353e18 m^-3,1.928e17 m^-3,1.645e14 m^-3,0 m^-3,2.2295e+19 m^-3],
[97000 m,95542 m,1.581e19 m^-3,4.500e17 m^-3,3.943e18 m^-3,1.746e17 m^-3,1.553e14 m^-3,0 m^-3,2.0378e+19 m^-3],
[97500 m,96027 m,1.445e19 m^-3,4.494e17 m^-3,3.570e18 m^-3,1.581e17 m^-3,1.468e14 m^-3,0 m^-3,1.8627e+19 m^-3],
[98000 m,96512 m,1.320e19 m^-3,4.476e17 m^-3,3.230e18 m^-3,1.430e17 m^-3,1.390e14 m^-3,0 m^-3,1.7021e+19 m^-3],
[98500 m,96997 m,1.206e19 m^-3,4.447e17 m^-3,2.920e18 m^-3,1.292e17 m^-3,1.317e14 m^-3,0 m^-3,1.5554e+19 m^-3],
[99000 m,97482 m,1.102e19 m^-3,4.408e17 m^-3,2.639e18 m^-3,1.167e17 m^-3,1.251e14 m^-3,0 m^-3,1.4217e+19 m^-3],
[99500 m,97967 m,1.008e19 m^-3,4.358e17 m^-3,2.383e18 m^-3,1.053e17 m^-3,1.190e14 m^-3,0 m^-3,1.3004e+19 m^-3],
[100000 m,98451 m,9.210e18 m^-3,4.298e17 m^-3,2.151e18 m^-3,9.501e16 m^-3,1.133e14 m^-3,0 m^-3,1.1886e+19 m^-3],
[101000 m,99420 m,7.740e18 m^-3,4.168e17 m^-3,1.756e18 m^-3,7.735e16 m^-3,1.034e14 m^-3,0 m^-3,9.9903e+18 m^-3],
[102000 m,100389 m,6.508e18 m^-3,4.007e17 m^-3,1.430e18 m^-3,6.279e16 m^-3,9.497e13 m^-3,0 m^-3,8.4016e+18 m^-3],
[103000 m,101358 m,5.475e18 m^-3,3.821e17 m^-3,1.163e18 m^-3,5.082e16 m^-3,8.776e13 m^-3,0 m^-3,7.0710e+18 m^-3],
[104000 m,102326 m,4.609e18 m^-3,3.619e17 m^-3,9.434e17 m^-3,4.101e16 m^-3,8.160e13 m^-3,0 m^-3,5.9554e+18 m^-3],
[105000 m,103294 m,3.883e18 m^-3,3.406e17 m^-3,7.645e17 m^-3,3.299e16 m^-3,7.633e13 m^-3,0 m^-3,5.0212e+18 m^-3],
[106000 m,104261 m,3.273e18 m^-3,3.188e17 m^-3,6.189e17 m^-3,2.645e16 m^-3,7.181e13 m^-3,0 m^-3,4.2373e+18 m^-3],
[107000 m,105229 m,2.760e18 m^-3,2.968e17 m^-3,5.005e17 m^-3,2.113e16 m^-3,6.789e13 m^-3,0 m^-3,3.5785e+18 m^-3],
[108000 m,106196 m,2.327e18 m^-3,2.748e17 m^-3,4.045e17 m^-3,1.681e16 m^-3,6.443e13 m^-3,0 m^-3,3.0232e+18 m^-3],
[109000 m,107162 m,1.959e18 m^-3,2.528e17 m^-3,3.263e17 m^-3,1.331e16 m^-3,6.128e13 m^-3,0 m^-3,2.5515e+18 m^-3],
[110000 m,108129 m,1.641e18 m^-3,2.303e17 m^-3,2.621e17 m^-3,1.046e16 m^-3,5.821e13 m^-3,0 m^-3,2.1440e+18 m^-3],
[111000 m,109095 m,1.373e18 m^-3,2.083e17 m^-3,2.104e17 m^-3,8.200e15 m^-3,5.526e13 m^-3,0 m^-3,1.8000e+18 m^-3],
[112000 m,110061 m,1.158e18 m^-3,1.889e17 m^-3,1.706e17 m^-3,6.481e15 m^-3,5.271e13 m^-3,0 m^-3,1.5241e+18 m^-3],
[113000 m,111026 m,9.841e17 m^-3,1.718e17 m^-3,1.398e17 m^-3,5.169e15 m^-3,5.044e13 m^-3,0 m^-3,1.3010e+18 m^-3],
[114000 m,111992 m,8.422e17 m^-3,1.565e17 m^-3,1.156e17 m^-3,4.163e15 m^-3,4.838e13 m^-3,0 m^-3,1.1185e+18 m^-3],
[115000 m,112957 m,7.254e17 m^-3,1.428e17 m^-3,9.646e16 m^-3,3.386e15 m^-3,4.648e13 m^-3,0 m^-3,9.6810e+17 m^-3],
[116000 m,113921 m,6.285e17 m^-3,1.305e17 m^-3,8.120e16 m^-3,2.779e15 m^-3,4.473e13 m^-3,0 m^-3,8.4302e+17 m^-3],
[117000 m,114885 m,5.475e17 m^-3,1.194e17 m^-3,6.891e16 m^-3,2.301e15 m^-3,4.310e13 m^-3,0 m^-3,7.3815e+17 m^-3],
[118000 m,115849 m,4.794e17 m^-3,1.096e17 m^-3,5.892e16 m^-3,1.920e15 m^-3,4.160e13 m^-3,0 m^-3,6.4988e+17 m^-3],
[119000 m,116813 m,4.217e17 m^-3,1.007e17 m^-3,5.072e16 m^-3,1.614e15 m^-3,4.019e13 m^-3,0 m^-3,5.7477e+17 m^-3],
[120000 m,117777 m,3.726e17 m^-3,9.275e16 m^-3,4.395e16 m^-3,1.366e15 m^-3,3.888e13 m^-3,0 m^-3,5.1071e+17 m^-3],
[121000 m,118740 m,3.306e17 m^-3,8.562e16 m^-3,3.832e16 m^-3,1.164e15 m^-3,3.766e13 m^-3,0 m^-3,4.5574e+17 m^-3],
[122000 m,119703 m,2.947e17 m^-3,7.925e16 m^-3,3.360e16 m^-3,9.979e14 m^-3,3.652e13 m^-3,0 m^-3,4.0859e+17 m^-3],
[123000 m,120665 m,2.637e17 m^-3,7.354e16 m^-3,2.963e16 m^-3,8.606e14 m^-3,3.547e13 m^-3,0 m^-3,3.6777e+17 m^-3],
[124000 m,121627 m,2.368e17 m^-3,6.840e16 m^-3,2.625e16 m^-3,7.460e14 m^-3,3.448e13 m^-3,0 m^-3,3.3223e+17 m^-3],
[125000 m,122589 m,2.135e17 m^-3,6.376e16 m^-3,2.336e16 m^-3,6.498e14 m^-3,3.356e13 m^-3,0 m^-3,3.0130e+17 m^-3],
[126000 m,123551 m,1.930e17 m^-3,5.956e16 m^-3,2.087e16 m^-3,5.685e14 m^-3,3.270e13 m^-3,0 m^-3,2.7403e+17 m^-3],
[127000 m,124512 m,1.750e17 m^-3,5.576e16 m^-3,1.871e16 m^-3,4.994e14 m^-3,3.189e13 m^-3,0 m^-3,2.5000e+17 m^-3],
[128000 m,125473 m,1.592e17 m^-3,5.229e16 m^-3,1.683e16 m^-3,4.403e14 m^-3,3.112e13 m^-3,0 m^-3,2.2879e+17 m^-3],
[129000 m,126434 m,1.451e17 m^-3,4.914e16 m^-3,1.519e16 m^-3,3.896e14 m^-3,3.040e13 m^-3,0 m^-3,2.0985e+17 m^-3],
[130000 m,127395 m,1.326e17 m^-3,4.625e16 m^-3,1.375e16 m^-3,3.458e14 m^-3,2.972e13 m^-3,0 m^-3,1.9298e+17 m^-3],
[131000 m,128355 m,1.215e17 m^-3,4.361e16 m^-3,1.247e16 m^-3,3.078e14 m^-3,2.907e13 m^-3,0 m^-3,1.7792e+17 m^-3],
[132000 m,129315 m,1.116e17 m^-3,4.118e16 m^-3,1.134e16 m^-3,2.748e14 m^-3,2.846e13 m^-3,0 m^-3,1.6442e+17 m^-3],
[133000 m,130274 m,1.026e17 m^-3,3.894e16 m^-3,1.034e16 m^-3,2.460e14 m^-3,2.787e13 m^-3,0 m^-3,1.5216e+17 m^-3],
[134000 m,131234 m,9.460e16 m^-3,3.688e16 m^-3,9.444e15 m^-3,2.207e14 m^-3,2.732e13 m^-3,0 m^-3,1.4117e+17 m^-3],
[135000 m,132193 m,8.735e16 m^-3,3.497e16 m^-3,8.645e15 m^-3,1.985e14 m^-3,2.679e13 m^-3,0 m^-3,1.3120e+17 m^-3],
[136000 m,133151 m,8.080e16 m^-3,3.320e16 m^-3,7.927e15 m^-3,1.789e14 m^-3,2.629e13 m^-3,0 m^-3,1.2214e+17 m^-3],
[137000 m,134110 m,7.487e16 m^-3,3.156e16 m^-3,7.283e15 m^-3,1.616e14 m^-3,2.581e13 m^-3,0 m^-3,1.1390e+17 m^-3],
[138000 m,135068 m,6.947e16 m^-3,3.004e16 m^-3,6.702e15 m^-3,1.463e14 m^-3,2.535e13 m^-3,0 m^-3,1.0639e+17 m^-3],
[139000 m,136026 m,6.456e16 m^-3,2.862e16 m^-3,6.177e15 m^-3,1.326e14 m^-3,2.491e13 m^-3,0 m^-3,9.9515e+16 m^-3],
[140000 m,136983 m,6.009e16 m^-3,2.729e16 m^-3,5.702e15 m^-3,1.205e14 m^-3,2.449e13 m^-3,0 m^-3,9.3227e+16 m^-3],
[141000 m,137940 m,5.600e16 m^-3,2.605e16 m^-3,5.272e15 m^-3,1.096e14 m^-3,2.408e13 m^-3,0 m^-3,8.7456e+16 m^-3],
[142000 m,138897 m,5.225e16 m^-3,2.489e16 m^-3,4.881e15 m^-3,9.989e13 m^-3,2.369e13 m^-3,0 m^-3,8.2145e+16 m^-3],
[143000 m,139854 m,4.881e16 m^-3,2.380e16 m^-3,4.524e15 m^-3,9.118e13 m^-3,2.332e13 m^-3,0 m^-3,7.7248e+16 m^-3],
[144000 m,140810 m,4.565e16 m^-3,2.278e16 m^-3,4.199e15 m^-3,8.335e13 m^-3,2.296e13 m^-3,0 m^-3,7.2735e+16 m^-3],
[145000 m,141766 m,4.275e16 m^-3,2.183e16 m^-3,3.903e15 m^-3,7.630e13 m^-3,2.261e13 m^-3,0 m^-3,6.8582e+16 m^-3],
[146000 m,142722 m,4.007e16 m^-3,2.092e16 m^-3,3.631e15 m^-3,6.994e13 m^-3,2.228e13 m^-3,0 m^-3,6.4713e+16 m^-3],
[147000 m,143677 m,3.760e16 m^-3,2.007e16 m^-3,3.382e15 m^-3,6.420e13 m^-3,2.196e13 m^-3,0 m^-3,6.1138e+16 m^-3],
[148000 m,144633 m,3.531e16 m^-3,1.927e16 m^-3,3.153e15 m^-3,5.900e13 m^-3,2.165e13 m^-3,0 m^-3,5.7814e+16 m^-3],
[149000 m,145587 m,3.320e16 m^-3,1.852e16 m^-3,2.943e15 m^-3,5.428e13 m^-3,2.135e13 m^-3,0 m^-3,5.4738e+16 m^-3],
[150000 m,146542 m,3.124e16 m^-3,1.780e16 m^-3,2.750e15 m^-3,5.000e13 m^-3,2.106e13 m^-3,3.767e11 m^-3,5.1861e+16 m^-3],
[151000 m,147496 m,2.942e16 m^-3,1.712e16 m^-3,2.572e15 m^-3,4.611e13 m^-3,2.078e13 m^-3,3.659e11 m^-3,4.9179e+16 m^-3],
[152000 m,148450 m,2.773e16 m^-3,1.648e16 m^-3,2.407e15 m^-3,4.256e13 m^-3,2.051e13 m^-3,3.557e11 m^-3,4.6681e+16 m^-3],
[153000 m,149404 m,2.616e16 m^-3,1.587e16 m^-3,2.255e15 m^-3,3.933e13 m^-3,2.024e13 m^-3,3.461e11 m^-3,4.4344e+16 m^-3],
[154000 m,150357 m,2.469e16 m^-3,1.530e16 m^-3,2.114e15 m^-3,3.638e13 m^-3,1.999e13 m^-3,3.369e11 m^-3,4.2160e+16 m^-3],
[155000 m,151311 m,2.333e16 m^-3,1.475e16 m^-3,1.984e15 m^-3,3.368e13 m^-3,1.974e13 m^-3,3.283e11 m^-3,4.0118e+16 m^-3],
[156000 m,152263 m,2.206e16 m^-3,1.423e16 m^-3,1.863e15 m^-3,3.121e13 m^-3,1.950e13 m^-3,3.201e11 m^-3,3.8204e+16 m^-3],
[157000 m,153216 m,2.087e16 m^-3,1.373e16 m^-3,1.751e15 m^-3,2.895e13 m^-3,1.927e13 m^-3,3.123e11 m^-3,3.6399e+16 m^-3],
[158000 m,154168 m,1.975e16 m^-3,1.326e16 m^-3,1.647e15 m^-3,2.687e13 m^-3,1.905e13 m^-3,3.049e11 m^-3,3.4703e+16 m^-3],
[159000 m,155120 m,1.871e16 m^-3,1.281e16 m^-3,1.550e15 m^-3,2.496e13 m^-3,1.883e13 m^-3,2.978e11 m^-3,3.3114e+16 m^-3],
[160000 m,156072 m,1.774e16 m^-3,1.238e16 m^-3,1.460e15 m^-3,2.321e13 m^-3,1.861e13 m^-3,2.911e11 m^-3,3.1622e+16 m^-3],
[161000 m,157023 m,1.682e16 m^-3,1.197e16 m^-3,1.376e15 m^-3,2.159e13 m^-3,1.841e13 m^-3,2.847e11 m^-3,3.0206e+16 m^-3],
[162000 m,157974 m,1.596e16 m^-3,1.158e16 m^-3,1.297e15 m^-3,2.011e13 m^-3,1.820e13 m^-3,2.786e11 m^-3,2.8875e+16 m^-3],
[163000 m,158925 m,1.516e16 m^-3,1.120e16 m^-3,1.224e15 m^-3,1.874e13 m^-3,1.801e13 m^-3,2.728e11 m^-3,2.7621e+16 m^-3],
[164000 m,159875 m,1.440e16 m^-3,1.085e16 m^-3,1.156e15 m^-3,1.747e13 m^-3,1.782e13 m^-3,2.672e11 m^-3,2.6441e+16 m^-3],
[165000 m,160826 m,1.369e16 m^-3,1.050e16 m^-3,1.092e15 m^-3,1.630e13 m^-3,1.763e13 m^-3,2.619e11 m^-3,2.5316e+16 m^-3],
[166000 m,161775 m,1.302e16 m^-3,1.018e16 m^-3,1.032e15 m^-3,1.522e13 m^-3,1.745e13 m^-3,2.568e11 m^-3,2.4264e+16 m^-3],
[167000 m,162725 m,1.239e16 m^-3,9.863e15 m^-3,9.757e14 m^-3,1.422e13 m^-3,1.727e13 m^-3,2.520e11 m^-3,2.3260e+16 m^-3],
[168000 m,163674 m,1.179e16 m^-3,9.562e15 m^-3,9.232e14 m^-3,1.329e13 m^-3,1.710e13 m^-3,2.473e11 m^-3,2.2305e+16 m^-3],
[169000 m,164623 m,1.123e16 m^-3,9.273e15 m^-3,8.739e14 m^-3,1.243e13 m^-3,1.693e13 m^-3,2.429e11 m^-3,2.1406e+16 m^-3],
[170000 m,165572 m,1.070e16 m^-3,8.996e15 m^-3,8.277e14 m^-3,1.163e13 m^-3,1.676e13 m^-3,2.386e11 m^-3,2.0553e+16 m^-3],
[171000 m,166521 m,1.020e16 m^-3,8.730e15 m^-3,7.843e14 m^-3,1.089e13 m^-3,1.660e13 m^-3,2.345e11 m^-3,1.9742e+16 m^-3],
[172000 m,167469 m,9.724e15 m^-3,8.474e15 m^-3,7.435e14 m^-3,1.020e13 m^-3,1.644e13 m^-3,2.306e11 m^-3,1.8968e+16 m^-3],
[173000 m,168417 m,9.277e15 m^-3,8.228e15 m^-3,7.051e14 m^-3,9.565e12 m^-3,1.629e13 m^-3,2.268e11 m^-3,1.8236e+16 m^-3],
[174000 m,169364 m,8.853e15 m^-3,7.992e15 m^-3,6.690e14 m^-3,8.970e12 m^-3,1.614e13 m^-3,2.232e11 m^-3,1.7539e+16 m^-3],
[175000 m,170311 m,8.452e15 m^-3,7.765e15 m^-3,6.350e14 m^-3,8.417e12 m^-3,1.599e13 m^-3,2.197e11 m^-3,1.6876e+16 m^-3],
[176000 m,171258 m,8.072e15 m^-3,7.546e15 m^-3,6.030e14 m^-3,7.901e12 m^-3,1.585e13 m^-3,2.163e11 m^-3,1.6245e+16 m^-3],
[177000 m,172205 m,7.712e15 m^-3,7.335e15 m^-3,5.728e14 m^-3,7.420e12 m^-3,1.571e13 m^-3,2.131e11 m^-3,1.5643e+16 m^-3],
[178000 m,173151 m,7.371e15 m^-3,7.132e15 m^-3,5.443e14 m^-3,6.971e12 m^-3,1.557e13 m^-3,2.100e11 m^-3,1.5070e+16 m^-3],
[179000 m,174098 m,7.047e15 m^-3,6.936e15 m^-3,5.174e14 m^-3,6.553e12 m^-3,1.543e13 m^-3,2.070e11 m^-3,1.4522e+16 m^-3],
[180000 m,175043 m,6.740e15 m^-3,6.747e15 m^-3,4.921e14 m^-3,6.162e12 m^-3,1.530e13 m^-3,2.041e11 m^-3,1.4000e+16 m^-3],
[181000 m,175989 m,6.448e15 m^-3,6.565e15 m^-3,4.681e14 m^-3,5.797e12 m^-3,1.517e13 m^-3,2.013e11 m^-3,1.3502e+16 m^-3],
[182000 m,176934 m,6.170e15 m^-3,6.389e15 m^-3,4.455e14 m^-3,5.456e12 m^-3,1.504e13 m^-3,1.987e11 m^-3,1.3025e+16 m^-3],
[183000 m,177879 m,5.907e15 m^-3,6.220e15 m^-3,4.241e14 m^-3,5.136e12 m^-3,1.492e13 m^-3,1.961e11 m^-3,1.2571e+16 m^-3],
[184000 m,178824 m,5.656e15 m^-3,6.056e15 m^-3,4.039e14 m^-3,4.838e12 m^-3,1.479e13 m^-3,1.936e11 m^-3,1.2136e+16 m^-3],
[185000 m,179768 m,5.417e15 m^-3,5.897e15 m^-3,3.847e14 m^-3,4.558e12 m^-3,1.467e13 m^-3,1.911e11 m^-3,1.1719e+16 m^-3],
[186000 m,180712 m,5.190e15 m^-3,5.744e15 m^-3,3.666e14 m^-3,4.296e12 m^-3,1.456e13 m^-3,1.888e11 m^-3,1.1320e+16 m^-3],
[187000 m,181656 m,4.974e15 m^-3,5.596e15 m^-3,3.494e14 m^-3,4.050e12 m^-3,1.444e13 m^-3,1.866e11 m^-3,1.0937e+16 m^-3],
[188000 m,182600 m,4.768e15 m^-3,5.453e15 m^-3,3.331e14 m^-3,3.820e12 m^-3,1.433e13 m^-3,1.844e11 m^-3,1.0572e+16 m^-3],
[189000 m,183543 m,4.572e15 m^-3,5.315e15 m^-3,3.177e14 m^-3,3.604e12 m^-3,1.421e13 m^-3,1.823e11 m^-3,1.0223e+16 m^-3],
[190000 m,184486 m,4.385e15 m^-3,5.181e15 m^-3,3.031e14 m^-3,3.401e12 m^-3,1.410e13 m^-3,1.802e11 m^-3,9.8868e+15 m^-3],
[191000 m,185428 m,4.207e15 m^-3,5.051e15 m^-3,2.892e14 m^-3,3.211e12 m^-3,1.400e13 m^-3,1.782e11 m^-3,9.5646e+15 m^-3],
[192000 m,186371 m,4.037e15 m^-3,4.926e15 m^-3,2.760e14 m^-3,3.033e12 m^-3,1.389e13 m^-3,1.763e11 m^-3,9.2561e+15 m^-3],
[193000 m,187313 m,3.875e15 m^-3,4.804e15 m^-3,2.635e14 m^-3,2.865e12 m^-3,1.379e13 m^-3,1.745e11 m^-3,8.9594e+15 m^-3],
[194000 m,188255 m,3.720e15 m^-3,4.686e15 m^-3,2.517e14 m^-3,2.707e12 m^-3,1.368e13 m^-3,1.727e11 m^-3,8.6743e+15 m^-3],
[195000 m,189196 m,3.572e15 m^-3,4.572e15 m^-3,2.404e14 m^-3,2.558e12 m^-3,1.358e13 m^-3,1.709e11 m^-3,8.4008e+15 m^-3],
[196000 m,190137 m,3.430e15 m^-3,4.461e15 m^-3,2.297e14 m^-3,2.419e12 m^-3,1.348e13 m^-3,1.692e11 m^-3,8.1368e+15 m^-3],
[197000 m,191078 m,3.295e15 m^-3,4.354e15 m^-3,2.195e14 m^-3,2.288e12 m^-3,1.339e13 m^-3,1.676e11 m^-3,7.8844e+15 m^-3],
[198000 m,192019 m,3.166e15 m^-3,4.249e15 m^-3,2.098e14 m^-3,2.164e12 m^-3,1.329e13 m^-3,1.660e11 m^-3,7.6405e+15 m^-3],
[199000 m,192959 m,3.043e15 m^-3,4.148e15 m^-3,2.006e14 m^-3,2.047e12 m^-3,1.319e13 m^-3,1.645e11 m^-3,7.4070e+15 m^-3],
[200000 m,193899 m,2.925e15 m^-3,4.050e15 m^-3,1.918e14 m^-3,1.938e12 m^-3,1.310e13 m^-3,1.630e11 m^-3,7.1820e+15 m^-3],
[201000 m,194839 m,2.812e15 m^-3,3.955e15 m^-3,1.834e14 m^-3,1.834e12 m^-3,1.301e13 m^-3,1.615e11 m^-3,6.9654e+15 m^-3],
[202000 m,195779 m,2.704e15 m^-3,3.862e15 m^-3,1.755e14 m^-3,1.737e12 m^-3,1.292e13 m^-3,1.601e11 m^-3,6.7563e+15 m^-3],
[203000 m,196718 m,2.601e15 m^-3,3.773e15 m^-3,1.679e14 m^-3,1.645e12 m^-3,1.283e13 m^-3,1.587e11 m^-3,6.5565e+15 m^-3],
[204000 m,197657 m,2.502e15 m^-3,3.685e15 m^-3,1.607e14 m^-3,1.558e12 m^-3,1.274e13 m^-3,1.574e11 m^-3,6.3622e+15 m^-3],
[205000 m,198595 m,2.407e15 m^-3,3.600e15 m^-3,1.538e14 m^-3,1.477e12 m^-3,1.266e13 m^-3,1.561e11 m^-3,6.1752e+15 m^-3],
[206000 m,199534 m,2.316e15 m^-3,3.518e15 m^-3,1.473e14 m^-3,1.399e12 m^-3,1.257e13 m^-3,1.548e11 m^-3,5.9955e+15 m^-3],
[207000 m,200472 m,2.229e15 m^-3,3.438e15 m^-3,1.410e14 m^-3,1.327e12 m^-3,1.249e13 m^-3,1.536e11 m^-3,5.8220e+15 m^-3],
[208000 m,201410 m,2.146e15 m^-3,3.360e15 m^-3,1.351e14 m^-3,1.258e12 m^-3,1.240e13 m^-3,1.524e11 m^-3,5.6550e+15 m^-3],
[209000 m,202347 m,2.066e15 m^-3,3.284e15 m^-3,1.294e14 m^-3,1.193e12 m^-3,1.232e13 m^-3,1.512e11 m^-3,5.4931e+15 m^-3],
[210000 m,203284 m,1.989e15 m^-3,3.211e15 m^-3,1.239e14 m^-3,1.131e12 m^-3,1.224e13 m^-3,1.501e11 m^-3,5.3374e+15 m^-3],
[211000 m,204221 m,1.915e15 m^-3,3.139e15 m^-3,1.188e14 m^-3,1.073e12 m^-3,1.216e13 m^-3,1.490e11 m^-3,5.1862e+15 m^-3],
[212000 m,205158 m,1.845e15 m^-3,3.069e15 m^-3,1.138e14 m^-3,1.019e12 m^-3,1.208e13 m^-3,1.479e11 m^-3,5.0410e+15 m^-3],
[213000 m,206094 m,1.777e15 m^-3,3.001e15 m^-3,1.091e14 m^-3,9.666e11 m^-3,1.201e13 m^-3,1.468e11 m^-3,4.9002e+15 m^-3],
[214000 m,207030 m,1.712e15 m^-3,2.935e15 m^-3,1.046e14 m^-3,9.176e11 m^-3,1.193e13 m^-3,1.458e11 m^-3,4.7645e+15 m^-3],
[215000 m,207966 m,1.650e15 m^-3,2.871e15 m^-3,1.003e14 m^-3,8.711e11 m^-3,1.185e13 m^-3,1.448e11 m^-3,4.6342e+15 m^-3],
[216000 m,208902 m,1.590e15 m^-3,2.808e15 m^-3,9.617e13 m^-3,8.272e11 m^-3,1.178e13 m^-3,1.439e11 m^-3,4.5069e+15 m^-3],
[217000 m,209837 m,1.533e15 m^-3,2.747e15 m^-3,9.224e13 m^-3,7.856e11 m^-3,1.171e13 m^-3,1.429e11 m^-3,4.3848e+15 m^-3],
[218000 m,210772 m,1.477e15 m^-3,2.688e15 m^-3,8.848e13 m^-3,7.463e11 m^-3,1.163e13 m^-3,1.420e11 m^-3,4.2659e+15 m^-3],
[219000 m,211706 m,1.424e15 m^-3,2.630e15 m^-3,8.489e13 m^-3,7.090e11 m^-3,1.156e13 m^-3,1.411e11 m^-3,4.1513e+15 m^-3],
[220000 m,212641 m,1.373e15 m^-3,2.573e15 m^-3,8.145e13 m^-3,6.737e11 m^-3,1.149e13 m^-3,1.402e11 m^-3,4.0398e+15 m^-3],
[221000 m,213575 m,1.324e15 m^-3,2.518e15 m^-3,7.816e13 m^-3,6.402e11 m^-3,1.142e13 m^-3,1.393e11 m^-3,3.9323e+15 m^-3],
[222000 m,214509 m,1.277e15 m^-3,2.465e15 m^-3,7.502e13 m^-3,6.085e11 m^-3,1.135e13 m^-3,1.385e11 m^-3,3.8291e+15 m^-3],
[223000 m,215442 m,1.232e15 m^-3,2.412e15 m^-3,7.201e13 m^-3,5.785e11 m^-3,1.128e13 m^-3,1.377e11 m^-3,3.7280e+15 m^-3],
[224000 m,216375 m,1.188e15 m^-3,2.361e15 m^-3,6.913e13 m^-3,5.500e11 m^-3,1.122e13 m^-3,1.369e11 m^-3,3.6300e+15 m^-3],
[225000 m,217308 m,1.147e15 m^-3,2.312e15 m^-3,6.637e13 m^-3,5.230e11 m^-3,1.115e13 m^-3,1.361e11 m^-3,3.5372e+15 m^-3],
[226000 m,218241 m,1.106e15 m^-3,2.263e15 m^-3,6.373e13 m^-3,4.974e11 m^-3,1.108e13 m^-3,1.353e11 m^-3,3.4444e+15 m^-3],
[227000 m,219173 m,1.068e15 m^-3,2.216e15 m^-3,6.121e13 m^-3,4.731e11 m^-3,1.102e13 m^-3,1.345e11 m^-3,3.3568e+15 m^-3],
[228000 m,220105 m,1.030e15 m^-3,2.170e15 m^-3,5.879e13 m^-3,4.501e11 m^-3,1.095e13 m^-3,1.338e11 m^-3,3.2704e+15 m^-3],
[229000 m,221037 m,9.945e14 m^-3,2.125e15 m^-3,5.647e13 m^-3,4.282e11 m^-3,1.089e13 m^-3,1.331e11 m^-3,3.1874e+15 m^-3],
[230000 m,221969 m,9.600e14 m^-3,2.081e15 m^-3,5.425e13 m^-3,4.075e11 m^-3,1.083e13 m^-3,1.324e11 m^-3,3.1066e+15 m^-3],
[231000 m,222900 m,9.268e14 m^-3,2.038e15 m^-3,5.212e13 m^-3,3.878e11 m^-3,1.076e13 m^-3,1.317e11 m^-3,3.0282e+15 m^-3],
[232000 m,223831 m,8.948e14 m^-3,1.996e15 m^-3,5.009e13 m^-3,3.691e11 m^-3,1.070e13 m^-3,1.310e11 m^-3,2.9521e+15 m^-3],
[233000 m,224762 m,8.640e14 m^-3,1.955e15 m^-3,4.813e13 m^-3,3.514e11 m^-3,1.064e13 m^-3,1.304e11 m^-3,2.8782e+15 m^-3],
[234000 m,225692 m,8.343e14 m^-3,1.915e15 m^-3,4.626e13 m^-3,3.345e11 m^-3,1.058e13 m^-3,1.297e11 m^-3,2.8066e+15 m^-3],
[235000 m,226622 m,8.058e14 m^-3,1.876e15 m^-3,4.446e13 m^-3,3.185e11 m^-3,1.052e13 m^-3,1.291e11 m^-3,2.7372e+15 m^-3],
[236000 m,227552 m,7.782e14 m^-3,1.838e15 m^-3,4.274e13 m^-3,3.033e11 m^-3,1.046e13 m^-3,1.285e11 m^-3,2.6698e+15 m^-3],
[237000 m,228481 m,7.517e14 m^-3,1.801e15 m^-3,4.109e13 m^-3,2.888e11 m^-3,1.040e13 m^-3,1.279e11 m^-3,2.6046e+15 m^-3],
[238000 m,229411 m,7.262e14 m^-3,1.765e15 m^-3,3.951e13 m^-3,2.751e11 m^-3,1.034e13 m^-3,1.273e11 m^-3,2.5414e+15 m^-3],
[239000 m,230340 m,7.016e14 m^-3,1.729e15 m^-3,3.799e13 m^-3,2.621e11 m^-3,1.029e13 m^-3,1.267e11 m^-3,2.4793e+15 m^-3],
[240000 m,231268 m,6.778e14 m^-3,1.695e15 m^-3,3.653e13 m^-3,2.497e11 m^-3,1.023e13 m^-3,1.261e11 m^-3,2.4198e+15 m^-3],
[241000 m,232197 m,6.550e14 m^-3,1.661e15 m^-3,3.513e13 m^-3,2.379e11 m^-3,1.017e13 m^-3,1.256e11 m^-3,2.3616e+15 m^-3],
[242000 m,233125 m,6.329e14 m^-3,1.628e15 m^-3,3.379e13 m^-3,2.267e11 m^-3,1.012e13 m^-3,1.250e11 m^-3,2.3051e+15 m^-3],
[243000 m,234053 m,6.117e14 m^-3,1.595e15 m^-3,3.251e13 m^-3,2.160e11 m^-3,1.006e13 m^-3,1.245e11 m^-3,2.2496e+15 m^-3],
[244000 m,234980 m,5.912e14 m^-3,1.564e15 m^-3,3.127e13 m^-3,2.059e11 m^-3,1.001e13 m^-3,1.240e11 m^-3,2.1968e+15 m^-3],
[245000 m,235908 m,5.714e14 m^-3,1.533e15 m^-3,3.008e13 m^-3,1.962e11 m^-3,9.953e12 m^-3,1.234e11 m^-3,2.1448e+15 m^-3],
[246000 m,236835 m,5.523e14 m^-3,1.503e15 m^-3,2.895e13 m^-3,1.871e11 m^-3,9.899e12 m^-3,1.229e11 m^-3,2.0945e+15 m^-3],
[247000 m,237761 m,5.339e14 m^-3,1.473e15 m^-3,2.785e13 m^-3,1.783e11 m^-3,9.846e12 m^-3,1.224e11 m^-3,2.0449e+15 m^-3],
[248000 m,238688 m,5.162e14 m^-3,1.444e15 m^-3,2.680e13 m^-3,1.700e11 m^-3,9.794e12 m^-3,1.219e11 m^-3,1.9971e+15 m^-3],
[249000 m,239614 m,4.991e14 m^-3,1.416e15 m^-3,2.579e13 m^-3,1.621e11 m^-3,9.741e12 m^-3,1.215e11 m^-3,1.9509e+15 m^-3],
[250000 m,240540 m,4.826e14 m^-3,1.388e15 m^-3,2.482e13 m^-3,1.546e11 m^-3,9.690e12 m^-3,1.210e11 m^-3,1.9054e+15 m^-3],
[251000 m,241466 m,4.666e14 m^-3,1.361e15 m^-3,2.389e13 m^-3,1.474e11 m^-3,9.638e12 m^-3,1.205e11 m^-3,1.8613e+15 m^-3],
[252000 m,242391 m,4.512e14 m^-3,1.335e15 m^-3,2.300e13 m^-3,1.406e11 m^-3,9.587e12 m^-3,1.201e11 m^-3,1.8190e+15 m^-3],
[253000 m,243316 m,4.364e14 m^-3,1.309e15 m^-3,2.214e13 m^-3,1.341e11 m^-3,9.537e12 m^-3,1.196e11 m^-3,1.7772e+15 m^-3],
[254000 m,244241 m,4.221e14 m^-3,1.284e15 m^-3,2.132e13 m^-3,1.280e11 m^-3,9.487e12 m^-3,1.192e11 m^-3,1.7371e+15 m^-3],
[255000 m,245165 m,4.082e14 m^-3,1.259e15 m^-3,2.052e13 m^-3,1.221e11 m^-3,9.438e12 m^-3,1.188e11 m^-3,1.6973e+15 m^-3],
[256000 m,246089 m,3.949e14 m^-3,1.235e15 m^-3,1.976e13 m^-3,1.165e11 m^-3,9.389e12 m^-3,1.183e11 m^-3,1.6593e+15 m^-3],
[257000 m,247013 m,3.820e14 m^-3,1.211e15 m^-3,1.903e13 m^-3,1.111e11 m^-3,9.340e12 m^-3,1.179e11 m^-3,1.6215e+15 m^-3],
[258000 m,247937 m,3.695e14 m^-3,1.188e15 m^-3,1.832e13 m^-3,1.060e11 m^-3,9.292e12 m^-3,1.175e11 m^-3,1.5853e+15 m^-3],
[259000 m,248860 m,3.575e14 m^-3,1.165e15 m^-3,1.765e13 m^-3,1.012e11 m^-3,9.244e12 m^-3,1.171e11 m^-3,1.5496e+15 m^-3],
[260000 m,249784 m,3.459e14 m^-3,1.143e15 m^-3,1.700e13 m^-3,9.658e10 m^-3,9.196e12 m^-3,1.167e11 m^-3,1.5153e+15 m^-3],
[261000 m,250706 m,3.347e14 m^-3,1.121e15 m^-3,1.637e13 m^-3,9.218e10 m^-3,9.149e12 m^-3,1.163e11 m^-3,1.4814e+15 m^-3],
[262000 m,251629 m,3.238e14 m^-3,1.100e15 m^-3,1.577e13 m^-3,8.799e10 m^-3,9.103e12 m^-3,1.159e11 m^-3,1.4489e+15 m^-3],
[263000 m,252551 m,3.134e14 m^-3,1.079e15 m^-3,1.519e13 m^-3,8.399e10 m^-3,9.056e12 m^-3,1.156e11 m^-3,1.4169e+15 m^-3],
[264000 m,253473 m,3.033e14 m^-3,1.059e15 m^-3,1.463e13 m^-3,8.019e10 m^-3,9.010e12 m^-3,1.152e11 m^-3,1.3861e+15 m^-3],
[265000 m,254395 m,2.935e14 m^-3,1.039e15 m^-3,1.410e13 m^-3,7.655e10 m^-3,8.965e12 m^-3,1.148e11 m^-3,1.3558e+15 m^-3],
[266000 m,255316 m,2.841e14 m^-3,1.019e15 m^-3,1.358e13 m^-3,7.309e10 m^-3,8.920e12 m^-3,1.145e11 m^-3,1.3258e+15 m^-3],
[267000 m,256237 m,2.749e14 m^-3,9.998e14 m^-3,1.309e13 m^-3,6.979e10 m^-3,8.875e12 m^-3,1.141e11 m^-3,1.2969e+15 m^-3],
[268000 m,257158 m,2.661e14 m^-3,9.811e14 m^-3,1.261e13 m^-3,6.665e10 m^-3,8.830e12 m^-3,1.138e11 m^-3,1.2688e+15 m^-3],
[269000 m,258079 m,2.576e14 m^-3,9.627e14 m^-3,1.215e13 m^-3,6.365e10 m^-3,8.786e12 m^-3,1.134e11 m^-3,1.2415e+15 m^-3],
[270000 m,258999 m,2.494e14 m^-3,9.447e14 m^-3,1.171e13 m^-3,6.078e10 m^-3,8.743e12 m^-3,1.131e11 m^-3,1.2147e+15 m^-3],
[271000 m,259919 m,2.414e14 m^-3,9.270e14 m^-3,1.128e13 m^-3,5.805e10 m^-3,8.699e12 m^-3,1.127e11 m^-3,1.1886e+15 m^-3],
[272000 m,260839 m,2.337e14 m^-3,9.097e14 m^-3,1.088e13 m^-3,5.545e10 m^-3,8.656e12 m^-3,1.124e11 m^-3,1.1632e+15 m^-3],
[273000 m,261758 m,2.263e14 m^-3,8.928e14 m^-3,1.048e13 m^-3,5.297e10 m^-3,8.613e12 m^-3,1.121e11 m^-3,1.1384e+15 m^-3],
[274000 m,262678 m,2.191e14 m^-3,8.762e14 m^-3,1.010e13 m^-3,5.060e10 m^-3,8.571e12 m^-3,1.118e11 m^-3,1.1142e+15 m^-3],
[275000 m,263597 m,2.121e14 m^-3,8.599e14 m^-3,9.739e12 m^-3,4.834e10 m^-3,8.529e12 m^-3,1.115e11 m^-3,1.0903e+15 m^-3],
[276000 m,264515 m,2.054e14 m^-3,8.440e14 m^-3,9.388e12 m^-3,4.618e10 m^-3,8.487e12 m^-3,1.112e11 m^-3,1.0674e+15 m^-3],
[277000 m,265434 m,1.989e14 m^-3,8.284e14 m^-3,9.050e12 m^-3,4.412e10 m^-3,8.445e12 m^-3,1.109e11 m^-3,1.0449e+15 m^-3],
[278000 m,266352 m,1.926e14 m^-3,8.131e14 m^-3,8.725e12 m^-3,4.216e10 m^-3,8.404e12 m^-3,1.106e11 m^-3,1.0229e+15 m^-3],
[279000 m,267269 m,1.865e14 m^-3,7.981e14 m^-3,8.412e12 m^-3,4.029e10 m^-3,8.363e12 m^-3,1.103e11 m^-3,1.0015e+15 m^-3],
[280000 m,268187 m,1.806e14 m^-3,7.834e14 m^-3,8.110e12 m^-3,3.850e10 m^-3,8.322e12 m^-3,1.100e11 m^-3,9.8058e+14 m^-3],
[281000 m,269104 m,1.750e14 m^-3,7.691e14 m^-3,7.820e12 m^-3,3.679e10 m^-3,8.282e12 m^-3,1.097e11 m^-3,9.6035e+14 m^-3],
[282000 m,270021 m,1.695e14 m^-3,7.549e14 m^-3,7.540e12 m^-3,3.516e10 m^-3,8.242e12 m^-3,1.094e11 m^-3,9.4033e+14 m^-3],
[283000 m,270938 m,1.641e14 m^-3,7.411e14 m^-3,7.271e12 m^-3,3.360e10 m^-3,8.202e12 m^-3,1.091e11 m^-3,9.2081e+14 m^-3],
[284000 m,271854 m,1.590e14 m^-3,7.276e14 m^-3,7.011e12 m^-3,3.212e10 m^-3,8.163e12 m^-3,1.088e11 m^-3,9.0191e+14 m^-3],
[285000 m,272771 m,1.540e14 m^-3,7.143e14 m^-3,6.761e12 m^-3,3.070e10 m^-3,8.124e12 m^-3,1.086e11 m^-3,8.8332e+14 m^-3],
[286000 m,273686 m,1.492e14 m^-3,7.012e14 m^-3,6.521e12 m^-3,2.935e10 m^-3,8.085e12 m^-3,1.083e11 m^-3,8.6515e+14 m^-3],
[287000 m,274602 m,1.445e14 m^-3,6.885e14 m^-3,6.289e12 m^-3,2.805e10 m^-3,8.046e12 m^-3,1.080e11 m^-3,8.4748e+14 m^-3],
[288000 m,275517 m,1.400e14 m^-3,6.759e14 m^-3,6.065e12 m^-3,2.682e10 m^-3,8.008e12 m^-3,1.078e11 m^-3,8.3012e+14 m^-3],
[289000 m,276432 m,1.356e14 m^-3,6.637e14 m^-3,5.850e12 m^-3,2.564e10 m^-3,7.969e12 m^-3,1.075e11 m^-3,8.1326e+14 m^-3],
[290000 m,277347 m,1.314e14 m^-3,6.516e14 m^-3,5.643e12 m^-3,2.451e10 m^-3,7.931e12 m^-3,1.073e11 m^-3,7.9670e+14 m^-3],
[291000 m,278262 m,1.273e14 m^-3,6.398e14 m^-3,5.443e12 m^-3,2.344e10 m^-3,7.894e12 m^-3,1.070e11 m^-3,7.8056e+14 m^-3],
[292000 m,279176 m,1.234e14 m^-3,6.282e14 m^-3,5.251e12 m^-3,2.241e10 m^-3,7.856e12 m^-3,1.067e11 m^-3,7.6484e+14 m^-3],
[293000 m,280090 m,1.195e14 m^-3,6.169e14 m^-3,5.065e12 m^-3,2.143e10 m^-3,7.819e12 m^-3,1.065e11 m^-3,7.4942e+14 m^-3],
[294000 m,281004 m,1.158e14 m^-3,6.058e14 m^-3,4.886e12 m^-3,2.049e10 m^-3,7.782e12 m^-3,1.063e11 m^-3,7.3440e+14 m^-3],
[295000 m,281917 m,1.122e14 m^-3,5.948e14 m^-3,4.714e12 m^-3,1.960e10 m^-3,7.746e12 m^-3,1.060e11 m^-3,7.1959e+14 m^-3],
[296000 m,282830 m,1.088e14 m^-3,5.841e14 m^-3,4.548e12 m^-3,1.874e10 m^-3,7.709e12 m^-3,1.058e11 m^-3,7.0529e+14 m^-3],
[297000 m,283743 m,1.054e14 m^-3,5.736e14 m^-3,4.388e12 m^-3,1.792e10 m^-3,7.673e12 m^-3,1.055e11 m^-3,6.9119e+14 m^-3],
[298000 m,284656 m,1.021e14 m^-3,5.633e14 m^-3,4.234e12 m^-3,1.714e10 m^-3,7.637e12 m^-3,1.053e11 m^-3,6.7740e+14 m^-3],
[299000 m,285568 m,9.898e13 m^-3,5.532e14 m^-3,4.085e12 m^-3,1.639e10 m^-3,7.602e12 m^-3,1.051e11 m^-3,6.6400e+14 m^-3],
[300000 m,286480 m,9.593e13 m^-3,5.433e14 m^-3,3.942e12 m^-3,1.568e10 m^-3,7.566e12 m^-3,1.049e11 m^-3,6.5086e+14 m^-3],
[302000 m,288303 m,9.011e13 m^-3,5.241e14 m^-3,3.670e12 m^-3,1.435e10 m^-3,7.496e12 m^-3,1.044e11 m^-3,6.2549e+14 m^-3],
[304000 m,290125 m,8.466e13 m^-3,5.055e14 m^-3,3.418e12 m^-3,1.313e10 m^-3,7.427e12 m^-3,1.040e11 m^-3,6.0112e+14 m^-3],
[306000 m,291946 m,7.954e13 m^-3,4.877e14 m^-3,3.184e12 m^-3,1.202e10 m^-3,7.358e12 m^-3,1.035e11 m^-3,5.7789e+14 m^-3],
[308000 m,293766 m,7.474e13 m^-3,4.705e14 m^-3,2.966e12 m^-3,1.100e10 m^-3,7.290e12 m^-3,1.031e11 m^-3,5.5561e+14 m^-3],
[310000 m,295585 m,7.024e13 m^-3,4.540e14 m^-3,2.763e12 m^-3,1.007e10 m^-3,7.224e12 m^-3,1.027e11 m^-3,5.3433e+14 m^-3],
[312000 m,297403 m,6.602e13 m^-3,4.380e14 m^-3,2.574e12 m^-3,9.223e9 m^-3,7.157e12 m^-3,1.023e11 m^-3,5.1386e+14 m^-3],
[314000 m,299220 m,6.206e13 m^-3,4.227e14 m^-3,2.399e12 m^-3,8.447e9 m^-3,7.092e12 m^-3,1.019e11 m^-3,4.9436e+14 m^-3],
[316000 m,301035 m,5.834e13 m^-3,4.079e14 m^-3,2.236e12 m^-3,7.737e9 m^-3,7.028e12 m^-3,1.015e11 m^-3,4.7562e+14 m^-3],
[318000 m,302850 m,5.485e13 m^-3,3.937e14 m^-3,2.084e12 m^-3,7.087e9 m^-3,6.964e12 m^-3,1.012e11 m^-3,4.5770e+14 m^-3],
[320000 m,304663 m,5.158e13 m^-3,3.800e14 m^-3,1.942e12 m^-3,6.493e9 m^-3,6.901e12 m^-3,1.008e11 m^-3,4.4053e+14 m^-3],
[322000 m,306470 m,4.850e13 m^-3,3.668e14 m^-3,1.811e12 m^-3,5.950e9 m^-3,6.839e12 m^-3,1.004e11 m^-3,4.2406e+14 m^-3],
[324000 m,308287 m,4.561e13 m^-3,3.541e14 m^-3,1.688e12 m^-3,5.452e9 m^-3,6.777e12 m^-3,1.001e11 m^-3,4.0829e+14 m^-3],
[326000 m,310097 m,4.290e13 m^-3,3.418e14 m^-3,1.574e12 m^-3,4.997e9 m^-3,6.717e12 m^-3,9.971e10 m^-3,3.9309e+14 m^-3],
[328000 m,311906 m,4.035e13 m^-3,3.300e14 m^-3,1.468e12 m^-3,4.580e9 m^-3,6.657e12 m^-3,9.937e10 m^-3,3.7858e+14 m^-3],
[330000 m,313714 m,3.796e13 m^-3,3.186e14 m^-3,1.369e12 m^-3,4.199e9 m^-3,6.597e12 m^-3,9.903e10 m^-3,3.6463e+14 m^-3],
[332000 m,315521 m,3.571e13 m^-3,3.076e14 m^-3,1.277e12 m^-3,3.850e9 m^-3,6.538e12 m^-3,9.869e10 m^-3,3.5123e+14 m^-3],
[334000 m,317327 m,3.360e13 m^-3,2.970e14 m^-3,1.191e12 m^-3,3.530e9 m^-3,6.480e12 m^-3,9.836e10 m^-3,3.3837e+14 m^-3],
[336000 m,319132 m,3.162e13 m^-3,2.868e14 m^-3,1.111e12 m^-3,3.237e9 m^-3,6.423e12 m^-3,9.804e10 m^-3,3.2605e+14 m^-3],
[338000 m,320935 m,2.975e13 m^-3,2.770e14 m^-3,1.037e12 m^-3,2.969e9 m^-3,6.366e12 m^-3,9.772e10 m^-3,3.1426e+14 m^-3],
[340000 m,322738 m,2.800e13 m^-3,2.675e14 m^-3,9.674e11 m^-3,2.723e9 m^-3,6.310e12 m^-3,9.741e10 m^-3,3.0288e+14 m^-3],
[342000 m,324539 m,2.635e13 m^-3,2.583e14 m^-3,9.027e11 m^-3,2.498e9 m^-3,6.254e12 m^-3,9.710e10 m^-3,2.9190e+14 m^-3],
[344000 m,326340 m,2.480e13 m^-3,2.495e14 m^-3,8.424e11 m^-3,2.292e9 m^-3,6.199e12 m^-3,9.680e10 m^-3,2.8144e+14 m^-3],
[346000 m,328139 m,2.335e13 m^-3,2.410e14 m^-3,7.862e11 m^-3,2.103e9 m^-3,6.145e12 m^-3,9.650e10 m^-3,2.7139e+14 m^-3],
[348000 m,329938 m,2.198e13 m^-3,2.328e14 m^-3,7.338e11 m^-3,1.929e9 m^-3,6.091e12 m^-3,9.620e10 m^-3,2.6170e+14 m^-3],
[350000 m,331735 m,2.069e13 m^-3,2.249e14 m^-3,6.850e11 m^-3,1.771e9 m^-3,6.038e12 m^-3,9.591e10 m^-3,2.5242e+14 m^-3],
[352000 m,333531 m,1.948e13 m^-3,2.172e14 m^-3,6.394e11 m^-3,1.625e9 m^-3,5.985e12 m^-3,9.562e10 m^-3,2.4341e+14 m^-3],
[354000 m,335326 m,1.834e13 m^-3,2.099e14 m^-3,5.969e11 m^-3,1.491e9 m^-3,5.933e12 m^-3,9.534e10 m^-3,2.3487e+14 m^-3],
[356000 m,337120 m,1.727e13 m^-3,2.027e14 m^-3,5.573e11 m^-3,1.369e9 m^-3,5.881e12 m^-3,9.505e10 m^-3,2.2651e+14 m^-3],
[358000 m,338913 m,1.627e13 m^-3,1.959e14 m^-3,5.204e11 m^-3,1.257e9 m^-3,5.830e12 m^-3,9.478e10 m^-3,2.1861e+14 m^-3],
[360000 m,340705 m,1.532e13 m^-3,1.893e14 m^-3,4.859e11 m^-3,1.154e9 m^-3,5.779e12 m^-3,9.450e10 m^-3,2.1098e+14 m^-3],
[362000 m,342496 m,1.443e13 m^-3,1.829e14 m^-3,4.538e11 m^-3,1.059e9 m^-3,5.729e12 m^-3,9.423e10 m^-3,2.0360e+14 m^-3],
[364000 m,344286 m,1.359e13 m^-3,1.767e14 m^-3,4.238e11 m^-3,9.728e8 m^-3,5.680e12 m^-3,9.397e10 m^-3,1.9648e+14 m^-3],
[366000 m,346074 m,1.280e13 m^-3,1.707e14 m^-3,3.958e11 m^-3,8.934e8 m^-3,5.631e12 m^-3,9.370e10 m^-3,1.8962e+14 m^-3],
[368000 m,347862 m,1.206e13 m^-3,1.650e14 m^-3,3.697e11 m^-3,8.205e8 m^-3,5.582e12 m^-3,9.344e10 m^-3,1.8310e+14 m^-3],
[370000 m,349648 m,1.136e13 m^-3,1.594e14 m^-3,3.454e11 m^-3,7.536e8 m^-3,5.534e12 m^-3,9.318e10 m^-3,1.7673e+14 m^-3],
[372000 m,351434 m,1.070e13 m^-3,1.541e14 m^-3,3.226e11 m^-3,6.922e8 m^-3,5.487e12 m^-3,9.293e10 m^-3,1.7070e+14 m^-3],
[374000 m,353218 m,1.008e13 m^-3,1.489e14 m^-3,3.014e11 m^-3,6.359e8 m^-3,5.439e12 m^-3,9.268e10 m^-3,1.6481e+14 m^-3],
[376000 m,355002 m,9.498e12 m^-3,1.439e14 m^-3,2.816e11 m^-3,5.842e8 m^-3,5.393e12 m^-3,9.243e10 m^-3,1.5916e+14 m^-3],
[378000 m,356784 m,8.950e12 m^-3,1.391e14 m^-3,2.631e11 m^-3,5.367e8 m^-3,5.347e12 m^-3,9.218e10 m^-3,1.5375e+14 m^-3],
[380000 m,358565 m,8.434e12 m^-3,1.344e14 m^-3,2.459e11 m^-3,4.932e8 m^-3,5.301e12 m^-3,9.193e10 m^-3,1.4847e+14 m^-3],
[382000 m,360346 m,7.948e12 m^-3,1.300e14 m^-3,2.297e11 m^-3,4.532e8 m^-3,5.256e12 m^-3,9.169e10 m^-3,1.4353e+14 m^-3],
[384000 m,362125 m,7.490e12 m^-3,1.256e14 m^-3,2.147e11 m^-3,4.165e8 m^-3,5.211e12 m^-3,9.145e10 m^-3,1.3860e+14 m^-3],
[386000 m,363903 m,7.059e12 m^-3,1.214e14 m^-3,2.006e11 m^-3,3.827e8 m^-3,5.167e12 m^-3,9.121e10 m^-3,1.3392e+14 m^-3],
[388000 m,365680 m,6.653e12 m^-3,1.174e14 m^-3,1.875e11 m^-3,3.518e8 m^-3,5.123e12 m^-3,9.098e10 m^-3,1.2945e+14 m^-3],
[390000 m,367456 m,6.271e12 m^-3,1.135e14 m^-3,1.753e11 m^-3,3.234e8 m^-3,5.079e12 m^-3,9.074e10 m^-3,1.2512e+14 m^-3],
[392000 m,369231 m,5.911e12 m^-3,1.097e14 m^-3,1.638e11 m^-3,2.972e8 m^-3,5.036e12 m^-3,9.051e10 m^-3,1.2090e+14 m^-3],
[394000 m,371005 m,5.572e12 m^-3,1.061e14 m^-3,1.532e11 m^-3,2.733e8 m^-3,4.993e12 m^-3,9.028e10 m^-3,1.1690e+14 m^-3],
[396000 m,372778 m,5.253e12 m^-3,1.025e14 m^-3,1.432e11 m^-3,2.512e8 m^-3,4.951e12 m^-3,9.005e10 m^-3,1.1293e+14 m^-3],
[398000 m,374549 m,4.952e12 m^-3,9.913e13 m^-3,1.339e11 m^-3,2.310e8 m^-3,4.909e12 m^-3,8.983e10 m^-3,1.0921e+14 m^-3],
[400000 m,376320 m,4.669e12 m^-3,9.584e13 m^-3,1.252e11 m^-3,2.124e8 m^-3,4.868e12 m^-3,8.960e10 m^-3,1.0560e+14 m^-3],
[402000 m,378090 m,4.402e12 m^-3,9.267e13 m^-3,1.170e11 m^-3,1.953e8 m^-3,4.827e12 m^-3,8.938e10 m^-3,1.0211e+14 m^-3],
[404000 m,379858 m,4.151e12 m^-3,8.960e13 m^-3,1.094e11 m^-3,1.796e8 m^-3,4.786e12 m^-3,8.916e10 m^-3,9.8735e+13 m^-3],
[406000 m,381626 m,3.914e12 m^-3,8.664e13 m^-3,1.023e11 m^-3,1.652e8 m^-3,4.746e12 m^-3,8.894e10 m^-3,9.5491e+13 m^-3],
[408000 m,383392 m,3.691e12 m^-3,8.378e13 m^-3,9.568e10 m^-3,1.519e8 m^-3,4.706e12 m^-3,8.872e10 m^-3,9.2362e+13 m^-3],
[410000 m,385158 m,3.480e12 m^-3,8.101e13 m^-3,8.948e10 m^-3,1.397e8 m^-3,4.666e12 m^-3,8.851e10 m^-3,8.9334e+13 m^-3],
[412000 m,386922 m,3.282e12 m^-3,7.834e13 m^-3,8.369e10 m^-3,1.285e8 m^-3,4.627e12 m^-3,8.829e10 m^-3,8.6421e+13 m^-3],
[414000 m,388686 m,3.095e12 m^-3,7.576e13 m^-3,7.827e10 m^-3,1.182e8 m^-3,4.588e12 m^-3,8.808e10 m^-3,8.3609e+13 m^-3],
[416000 m,390448 m,2.919e12 m^-3,7.327e13 m^-3,7.321e10 m^-3,1.088e8 m^-3,4.550e12 m^-3,8.787e10 m^-3,8.0900e+13 m^-3],
[418000 m,392210 m,2.754e12 m^-3,7.086e13 m^-3,6.848e10 m^-3,1.001e8 m^-3,4.512e12 m^-3,8.766e10 m^-3,7.8282e+13 m^-3],
[420000 m,393970 m,2.597e12 m^-3,6.853e13 m^-3,6.406e10 m^-3,9.207e7 m^-3,4.474e12 m^-3,8.745e10 m^-3,7.5752e+13 m^-3],
[422000 m,395729 m,2.450e12 m^-3,6.628e13 m^-3,5.993e10 m^-3,8.472e7 m^-3,4.437e12 m^-3,8.725e10 m^-3,7.3314e+13 m^-3],
[424000 m,397487 m,2.311e12 m^-3,6.410e13 m^-3,5.606e10 m^-3,7.796e7 m^-3,4.399e12 m^-3,8.704e10 m^-3,7.0953e+13 m^-3],
[426000 m,399245 m,2.180e12 m^-3,6.200e13 m^-3,5.245e10 m^-3,7.174e7 m^-3,4.363e12 m^-3,8.684e10 m^-3,6.8682e+13 m^-3],
[428000 m,401001 m,2.057e12 m^-3,5.997e13 m^-3,4.907e10 m^-3,6.602e7 m^-3,4.326e12 m^-3,8.663e10 m^-3,6.6489e+13 m^-3],
[430000 m,402756 m,1.940e12 m^-3,5.800e13 m^-3,4.592e10 m^-3,6.076e7 m^-3,4.290e12 m^-3,8.643e10 m^-3,6.4362e+13 m^-3],
[432000 m,404510 m,1.831e12 m^-3,5.611e13 m^-3,4.297e10 m^-3,5.593e7 m^-3,4.255e12 m^-3,8.623e10 m^-3,6.2325e+13 m^-3],
[434000 m,406263 m,1.727e12 m^-3,5.427e13 m^-3,4.020e10 m^-3,5.148e7 m^-3,4.219e12 m^-3,8.603e10 m^-3,6.0342e+13 m^-3],
[436000 m,408015 m,1.630e12 m^-3,5.250e13 m^-3,3.762e10 m^-3,4.739e7 m^-3,4.184e12 m^-3,8.583e10 m^-3,5.8438e+13 m^-3],
[438000 m,409766 m,1.538e12 m^-3,5.079e13 m^-3,3.521e10 m^-3,4.362e7 m^-3,4.150e12 m^-3,8.564e10 m^-3,5.6599e+13 m^-3],
[440000 m,411516 m,1.451e12 m^-3,4.913e13 m^-3,3.295e10 m^-3,4.016e7 m^-3,4.115e12 m^-3,8.544e10 m^-3,5.4814e+13 m^-3],
[442000 m,413265 m,1.369e12 m^-3,4.753e13 m^-3,3.084e10 m^-3,3.698e7 m^-3,4.081e12 m^-3,8.525e10 m^-3,5.3096e+13 m^-3],
[444000 m,415013 m,1.292e12 m^-3,4.598e13 m^-3,2.887e10 m^-3,3.404e7 m^-3,4.047e12 m^-3,8.505e10 m^-3,5.1433e+13 m^-3],
[446000 m,416760 m,1.220e12 m^-3,4.448e13 m^-3,2.702e10 m^-3,3.135e7 m^-3,4.014e12 m^-3,8.486e10 m^-3,4.9826e+13 m^-3],
[448000 m,418505 m,1.151e12 m^-3,4.303e13 m^-3,2.529e10 m^-3,2.887e7 m^-3,3.981e12 m^-3,8.467e10 m^-3,4.8272e+13 m^-3],
[450000 m,420250 m,1.086e12 m^-3,4.164e13 m^-3,2.368e10 m^-3,2.658e7 m^-3,3.948e12 m^-3,8.448e10 m^-3,4.6782e+13 m^-3],
[452000 m,421994 m,1.025e12 m^-3,4.028e13 m^-3,2.216e10 m^-3,2.448e7 m^-3,3.915e12 m^-3,8.429e10 m^-3,4.5326e+13 m^-3],
[454000 m,423737 m,9.679e11 m^-3,3.898e13 m^-3,2.075e10 m^-3,2.255e7 m^-3,3.883e12 m^-3,8.410e10 m^-3,4.3936e+13 m^-3],
[456000 m,425478 m,9.136e11 m^-3,3.771e13 m^-3,1.943e10 m^-3,2.077e7 m^-3,3.851e12 m^-3,8.391e10 m^-3,4.2578e+13 m^-3],
[458000 m,427219 m,8.625e11 m^-3,3.649e13 m^-3,1.819e10 m^-3,1.913e7 m^-3,3.819e12 m^-3,8.373e10 m^-3,4.1274e+13 m^-3],
[460000 m,428959 m,8.142e11 m^-3,3.531e13 m^-3,1.703e10 m^-3,1.762e7 m^-3,3.788e12 m^-3,8.354e10 m^-3,4.0013e+13 m^-3],
[462000 m,430698 m,7.686e11 m^-3,3.416e13 m^-3,1.595e10 m^-3,1.623e7 m^-3,3.757e12 m^-3,8.336e10 m^-3,3.8785e+13 m^-3],
[464000 m,432435 m,7.256e11 m^-3,3.306e13 m^-3,1.493e10 m^-3,1.495e7 m^-3,3.726e12 m^-3,8.317e10 m^-3,3.7610e+13 m^-3],
[466000 m,434172 m,6.851e11 m^-3,3.199e13 m^-3,1.398e10 m^-3,1.377e7 m^-3,3.695e12 m^-3,8.299e10 m^-3,3.6467e+13 m^-3],
[468000 m,435907 m,6.468e11 m^-3,3.096e13 m^-3,1.309e10 m^-3,1.269e7 m^-3,3.665e12 m^-3,8.281e10 m^-3,3.5368e+13 m^-3],
[470000 m,437642 m,6.107e11 m^-3,2.996e13 m^-3,1.226e10 m^-3,1.169e7 m^-3,3.635e12 m^-3,8.263e10 m^-3,3.4301e+13 m^-3],
[472000 m,439376 m,5.766e11 m^-3,2.899e13 m^-3,1.148e10 m^-3,1.077e7 m^-3,3.605e12 m^-3,8.245e10 m^-3,3.3265e+13 m^-3],
[474000 m,441108 m,5.445e11 m^-3,2.806e13 m^-3,1.076e10 m^-3,9.929e6 m^-3,3.576e12 m^-3,8.227e10 m^-3,3.2274e+13 m^-3],
[476000 m,442840 m,5.142e11 m^-3,2.715e13 m^-3,1.007e10 m^-3,9.149e6 m^-3,3.547e12 m^-3,8.209e10 m^-3,3.1303e+13 m^-3],
[478000 m,444570 m,4.855e11 m^-3,2.628e13 m^-3,9.436e9 m^-3,8.432e6 m^-3,3.518e12 m^-3,8.191e10 m^-3,3.0375e+13 m^-3],
[480000 m,446300 m,4.585e11 m^-3,2.543e13 m^-3,8.839e9 m^-3,7.771e6 m^-3,3.489e12 m^-3,8.173e10 m^-3,2.9469e+13 m^-3],
[482000 m,448028 m,4.330e11 m^-3,2.461e13 m^-3,8.280e9 m^-3,7.162e6 m^-3,3.461e12 m^-3,8.155e10 m^-3,2.8594e+13 m^-3],
[484000 m,449756 m,4.090e11 m^-3,2.382e13 m^-3,7.757e9 m^-3,6.602e6 m^-3,3.432e12 m^-3,8.138e10 m^-3,2.7750e+13 m^-3],
[486000 m,451482 m,3.863e11 m^-3,2.306e13 m^-3,7.267e9 m^-3,6.085e6 m^-3,3.404e12 m^-3,8.120e10 m^-3,2.6938e+13 m^-3],
[488000 m,453208 m,3.648e11 m^-3,2.232e13 m^-3,6.808e9 m^-3,5.609e6 m^-3,3.377e12 m^-3,8.103e10 m^-3,2.6150e+13 m^-3],
[490000 m,454932 m,3.446e11 m^-3,2.160e13 m^-3,6.378e9 m^-3,5.171e6 m^-3,3.349e12 m^-3,8.085e10 m^-3,2.5381e+13 m^-3],
[492000 m,456656 m,3.255e11 m^-3,2.091e13 m^-3,5.976e9 m^-3,4.767e6 m^-3,3.322e12 m^-3,8.068e10 m^-3,2.4645e+13 m^-3],
[494000 m,458378 m,3.075e11 m^-3,2.024e13 m^-3,5.599e9 m^-3,4.395e6 m^-3,3.295e12 m^-3,8.051e10 m^-3,2.3930e+13 m^-3],
[496000 m,460100 m,2.904e11 m^-3,1.959e13 m^-3,5.247e9 m^-3,4.052e6 m^-3,3.268e12 m^-3,8.034e10 m^-3,2.3233e+13 m^-3],
[498000 m,461820 m,2.744e11 m^-3,1.896e13 m^-3,4.917e9 m^-3,3.737e6 m^-3,3.242e12 m^-3,8.017e10 m^-3,2.2561e+13 m^-3],
[500000 m,463540 m,2.592e11 m^-3,1.836e13 m^-3,4.607e9 m^-3,3.445e6 m^-3,3.215e12 m^-3,8.000e10 m^-3,2.1919e+13 m^-3],
[505000 m,467834 m,2.249e11 m^-3,1.693e13 m^-3,3.917e9 m^-3,2.814e6 m^-3,3.151e12 m^-3,7.959e10 m^-3,2.0390e+13 m^-3],
[510000 m,472122 m,1.951e11 m^-3,1.561e13 m^-3,3.331e9 m^-3,2.299e6 m^-3,3.087e12 m^-3,7.918e10 m^-3,1.8974e+13 m^-3],
[515000 m,476404 m,1.694e11 m^-3,1.440e13 m^-3,2.834e9 m^-3,1.878e6 m^-3,3.026e12 m^-3,7.878e10 m^-3,1.7677e+13 m^-3],
[520000 m,480679 m,1.470e11 m^-3,1.328e13 m^-3,2.411e9 m^-3,1.535e6 m^-3,2.965e12 m^-3,7.838e10 m^-3,1.6472e+13 m^-3],
[525000 m,484949 m,1.277e11 m^-3,1.225e13 m^-3,2.052e9 m^-3,1.255e6 m^-3,2.906e12 m^-3,7.798e10 m^-3,1.5364e+13 m^-3],
[530000 m,489212 m,1.109e11 m^-3,1.130e13 m^-3,1.747e9 m^-3,1.027e6 m^-3,2.848e12 m^-3,7.758e10 m^-3,1.4339e+13 m^-3],
[535000 m,493469 m,9.633e10 m^-3,1.043e13 m^-3,1.487e9 m^-3,8.400e5 m^-3,2.791e12 m^-3,7.719e10 m^-3,1.3395e+13 m^-3],
[540000 m,497719 m,8.370e10 m^-3,9.624e12 m^-3,1.267e9 m^-3,6.875e5 m^-3,2.735e12 m^-3,7.680e10 m^-3,1.2521e+13 m^-3],
[545000 m,501964 m,7.274e10 m^-3,8.883e12 m^-3,1.079e9 m^-3,5.628e5 m^-3,2.681e12 m^-3,7.641e10 m^-3,1.1714e+13 m^-3],
[550000 m,506202 m,6.323e10 m^-3,8.200e12 m^-3,9.196e8 m^-3,4.609e5 m^-3,2.628e12 m^-3,7.602e10 m^-3,1.0968e+13 m^-3],
[555000 m,510435 m,5.497e10 m^-3,7.570e12 m^-3,7.838e8 m^-3,3.775e5 m^-3,2.576e12 m^-3,7.564e10 m^-3,1.0278e+13 m^-3],
[560000 m,514661 m,4.781e10 m^-3,6.989e12 m^-3,6.682e8 m^-3,3.093e5 m^-3,2.525e12 m^-3,7.526e10 m^-3,9.6378e+12 m^-3],
[565000 m,518881 m,4.158e10 m^-3,6.454e12 m^-3,5.697e8 m^-3,2.535e5 m^-3,2.475e12 m^-3,7.488e10 m^-3,9.0461e+12 m^-3],
[570000 m,523095 m,3.617e10 m^-3,5.960e12 m^-3,4.859e8 m^-3,2.079e5 m^-3,2.426e12 m^-3,7.451e10 m^-3,8.4972e+12 m^-3],
[575000 m,527303 m,3.148e10 m^-3,5.505e12 m^-3,4.146e8 m^-3,1.705e5 m^-3,2.379e12 m^-3,7.413e10 m^-3,7.9900e+12 m^-3],
[580000 m,531505 m,2.740e10 m^-3,5.085e12 m^-3,3.537e8 m^-3,1.398e5 m^-3,2.332e12 m^-3,7.376e10 m^-3,7.5186e+12 m^-3],
[585000 m,535701 m,2.385e10 m^-3,4.698e12 m^-3,3.019e8 m^-3,1.147e5 m^-3,2.286e12 m^-3,7.339e10 m^-3,7.0816e+12 m^-3],
[590000 m,539890 m,2.076e10 m^-3,4.341e12 m^-3,2.578e8 m^-3,9.419e4 m^-3,2.241e12 m^-3,7.303e10 m^-3,6.6761e+12 m^-3],
[595000 m,544074 m,1.808e10 m^-3,4.011e12 m^-3,2.201e8 m^-3,7.733e4 m^-3,2.197e12 m^-3,7.267e10 m^-3,6.2990e+12 m^-3],
[600000 m,548252 m,1.575e10 m^-3,3.707e12 m^-3,1.880e8 m^-3,6.351e4 m^-3,2.154e12 m^-3,7.231e10 m^-3,5.9493e+12 m^-3],
[605000 m,552424 m,1.372e10 m^-3,3.426e12 m^-3,1.606e8 m^-3,5.217e4 m^-3,2.112e12 m^-3,7.195e10 m^-3,5.6239e+12 m^-3],
[610000 m,556589 m,1.196e10 m^-3,3.167e12 m^-3,1.372e8 m^-3,4.287e4 m^-3,2.071e12 m^-3,7.159e10 m^-3,5.3217e+12 m^-3],
[615000 m,560749 m,1.042e10 m^-3,2.928e12 m^-3,1.173e8 m^-3,3.524e4 m^-3,2.031e12 m^-3,7.124e10 m^-3,5.0407e+12 m^-3],
[620000 m,564903 m,9.085e9 m^-3,2.707e12 m^-3,1.003e8 m^-3,2.896e4 m^-3,1.991e12 m^-3,7.089e10 m^-3,4.7781e+12 m^-3],
[625000 m,569051 m,7.921e9 m^-3,2.503e12 m^-3,8.573e7 m^-3,2.383e4 m^-3,1.953e12 m^-3,7.054e10 m^-3,4.5345e+12 m^-3],
[630000 m,573193 m,6.908e9 m^-3,2.315e12 m^-3,7.332e7 m^-3,1.961e4 m^-3,1.915e12 m^-3,7.019e10 m^-3,4.3072e+12 m^-3],
[635000 m,577329 m,6.025e9 m^-3,2.141e12 m^-3,6.272e7 m^-3,1.613e4 m^-3,1.878e12 m^-3,6.985e10 m^-3,4.0950e+12 m^-3],
[640000 m,581459 m,5.257e9 m^-3,1.981e12 m^-3,5.367e7 m^-3,1.328e4 m^-3,1.842e12 m^-3,6.950e10 m^-3,3.8979e+12 m^-3],
[645000 m,585583 m,4.587e9 m^-3,1.832e12 m^-3,4.593e7 m^-3,1.094e4 m^-3,1.806e12 m^-3,6.916e10 m^-3,3.7118e+12 m^-3],
[650000 m,589701 m,4.003e9 m^-3,1.695e12 m^-3,3.932e7 m^-3,9.006e3 m^-3,1.771e12 m^-3,6.883e10 m^-3,3.5388e+12 m^-3],
[655000 m,593614 m,3.495e9 m^-3,1.569e12 m^-3,3.367e7 m^-3,7.420e3 m^-3,1.737e12 m^-3,6.849e10 m^-3,3.3780e+12 m^-3],
[660000 m,597920 m,3.051e9 m^-3,1.452e12 m^-3,2.883e7 m^-3,6.114e3 m^-3,1.704e12 m^-3,6.816e10 m^-3,3.2273e+12 m^-3],
[665000 m,602021 m,2.665e9 m^-3,1.344e12 m^-3,2.470e7 m^-3,5.040e3 m^-3,1.671e12 m^-3,6.782e10 m^-3,3.0855e+12 m^-3],
[670000 m,606116 m,2.327e9 m^-3,1.244e12 m^-3,2.116e7 m^-3,4.155e3 m^-3,1.639e12 m^-3,6.749e10 m^-3,2.9528e+12 m^-3],
[675000 m,610205 m,2.033e9 m^-3,1.151e12 m^-3,1.813e7 m^-3,3.427e3 m^-3,1.608e12 m^-3,6.717e10 m^-3,2.8282e+12 m^-3],
[680000 m,614288 m,1.777e9 m^-3,1.066e12 m^-3,1.554e7 m^-3,2.827e3 m^-3,1.577e12 m^-3,6.684e10 m^-3,2.7116e+12 m^-3],
[685000 m,618365 m,1.553e9 m^-3,9.870e11 m^-3,1.333e7 m^-3,2.333e3 m^-3,1.547e12 m^-3,6.652e10 m^-3,2.6021e+12 m^-3],
[690000 m,622437 m,1.357e9 m^-3,9.140e11 m^-3,1.143e7 m^-3,1.926e3 m^-3,1.518e12 m^-3,6.620e10 m^-3,2.4996e+12 m^-3],
[695000 m,626503 m,1.187e9 m^-3,8.465e11 m^-3,9.802e6 m^-3,1.590e3 m^-3,1.489e12 m^-3,6.588e10 m^-3,2.4026e+12 m^-3],
[700000 m,630563 m,1.038e9 m^-3,7.840e11 m^-3,8.410e6 m^-3,1.313e3 m^-3,1.461e12 m^-3,6.556e10 m^-3,2.3117e+12 m^-3],
[705000 m,634617 m,9.075e8 m^-3,7.263e11 m^-3,7.216e6 m^-3,1.085e3 m^-3,1.433e12 m^-3,6.524e10 m^-3,2.2254e+12 m^-3],
[710000 m,638666 m,7.939e8 m^-3,6.728e11 m^-3,6.194e6 m^-3,8.964e2 m^-3,1.406e12 m^-3,6.493e10 m^-3,2.1445e+12 m^-3],
[715000 m,642709 m,6.946e8 m^-3,6.234e11 m^-3,5.317e6 m^-3,7.409e2 m^-3,1.379e12 m^-3,6.462e10 m^-3,2.0677e+12 m^-3],
[720000 m,646746 m,6.078e8 m^-3,5.777e11 m^-3,4.566e6 m^-3,6.126e2 m^-3,1.353e12 m^-3,6.431e10 m^-3,1.9956e+12 m^-3],
[725000 m,650778 m,5.320e8 m^-3,5.354e11 m^-3,3.921e6 m^-3,5.066e2 m^-3,1.328e12 m^-3,6.400e10 m^-3,1.9279e+12 m^-3],
[730000 m,654803 m,4.658e8 m^-3,4.962e11 m^-3,3.368e6 m^-3,4.191e2 m^-3,1.303e12 m^-3,6.370e10 m^-3,1.8634e+12 m^-3],
[735000 m,658824 m,4.078e8 m^-3,4.599e11 m^-3,2.894e6 m^-3,3.467e2 m^-3,1.278e12 m^-3,6.339e10 m^-3,1.8017e+12 m^-3],
[740000 m,662838 m,3.572e8 m^-3,4.264e11 m^-3,2.487e6 m^-3,2.870e2 m^-3,1.254e12 m^-3,6.309e10 m^-3,1.7439e+12 m^-3],
[745000 m,666847 m,3.129e8 m^-3,3.953e11 m^-3,2.138e6 m^-3,2.376e2 m^-3,1.231e12 m^-3,6.279e10 m^-3,1.6894e+12 m^-3],
[750000 m,670850 m,2.741e8 m^-3,3.666e11 m^-3,1.838e6 m^-3,1.967e2 m^-3,1.208e12 m^-3,6.249e10 m^-3,1.6374e+12 m^-3],
[755000 m,674848 m,2.402e8 m^-3,3.399e11 m^-3,1.581e6 m^-3,1.630e2 m^-3,1.185e12 m^-3,6.220e10 m^-3,1.5873e+12 m^-3],
[760000 m,678840 m,2.105e8 m^-3,3.153e11 m^-3,1.360e6 m^-3,1.350e2 m^-3,1.163e12 m^-3,6.190e10 m^-3,1.5404e+12 m^-3],
[765000 m,682826 m,1.845e8 m^-3,2.924e11 m^-3,1.170e6 m^-3,1.119e2 m^-3,1.141e12 m^-3,6.161e10 m^-3,1.4952e+12 m^-3],
[770000 m,686807 m,1.618e8 m^-3,2.712e11 m^-3,1.007e6 m^-3,9.276e1 m^-3,1.120e12 m^-3,6.132e10 m^-3,1.4527e+12 m^-3],
[775000 m,690782 m,1.419e8 m^-3,2.516e11 m^-3,8.664e5 m^-3,7.692e1 m^-3,1.099e12 m^-3,6.103e10 m^-3,1.4117e+12 m^-3],
[780000 m,694751 m,1.244e8 m^-3,2.335e11 m^-3,7.458e5 m^-3,6.380e1 m^-3,1.079e12 m^-3,6.074e10 m^-3,1.3733e+12 m^-3],
[785000 m,698715 m,1.092e8 m^-3,2.166e11 m^-3,6.422e5 m^-3,5.293e1 m^-3,1.059e12 m^-3,6.046e10 m^-3,1.3362e+12 m^-3],
[790000 m,702674 m,9.577e7 m^-3,2.011e11 m^-3,5.531e5 m^-3,4.392e1 m^-3,1.039e12 m^-3,6.017e10 m^-3,1.3004e+12 m^-3],
[795000 m,706627 m,8.404e7 m^-3,1.866e11 m^-3,4.764e5 m^-3,3.646e1 m^-3,1.020e12 m^-3,5.989e10 m^-3,1.2666e+12 m^-3],
[800000 m,710574 m,7.377e7 m^-3,1.732e11 m^-3,4.105e5 m^-3,3.027e1 m^-3,1.001e12 m^-3,5.961e10 m^-3,1.2339e+12 m^-3],
[805000 m,714516 m,6.476e7 m^-3,1.608e11 m^-3,3.537e5 m^-3,2.514e1 m^-3,9.826e11 m^-3,5.933e10 m^-3,1.2028e+12 m^-3],
[810000 m,718452 m,5.686e7 m^-3,1.493e11 m^-3,3.049e5 m^-3,2.088e1 m^-3,9.645e11 m^-3,5.905e10 m^-3,1.1730e+12 m^-3],
[815000 m,722383 m,4.993e7 m^-3,1.386e11 m^-3,2.628e5 m^-3,1.735e1 m^-3,9.468e11 m^-3,5.878e10 m^-3,1.1443e+12 m^-3],
[820000 m,726309 m,4.386e7 m^-3,1.287e11 m^-3,2.267e5 m^-3,1.442e1 m^-3,9.294e11 m^-3,5.851e10 m^-3,1.1166e+12 m^-3],
[825000 m,730229 m,3.853e7 m^-3,1.195e11 m^-3,1.955e5 m^-3,1.199e1 m^-3,9.124e11 m^-3,5.823e10 m^-3,1.0901e+12 m^-3],
[830000 m,734143 m,3.386e7 m^-3,1.110e11 m^-3,1.686e5 m^-3,9.970e0 m^-3,8.957e11 m^-3,5.796e10 m^-3,1.0647e+12 m^-3],
[835000 m,738052 m,2.975e7 m^-3,1.031e11 m^-3,1.455e5 m^-3,8.293e0 m^-3,8.793e11 m^-3,5.769e10 m^-3,1.0401e+12 m^-3],
[840000 m,741956 m,2.615e7 m^-3,9.580e10 m^-3,1.256e5 m^-3,6.900e0 m^-3,8.632e11 m^-3,5.743e10 m^-3,1.0165e+12 m^-3],
[845000 m,745854 m,2.299e7 m^-3,8.901e10 m^-3,1.084e5 m^-3,5.742e0 m^-3,8.475e11 m^-3,5.716e10 m^-3,9.9369e+11 m^-3],
[850000 m,749747 m,2.022e7 m^-3,8.270e10 m^-3,9.358e4 m^-3,4.780e0 m^-3,8.320e11 m^-3,5.690e10 m^-3,9.7162e+11 m^-3],
[855000 m,753634 m,1.778e7 m^-3,7.685e10 m^-3,8.081e4 m^-3,3.980e0 m^-3,8.169e11 m^-3,5.664e10 m^-3,9.5041e+11 m^-3],
[860000 m,757516 m,1.564e7 m^-3,7.142e10 m^-3,6.979e4 m^-3,3.314e0 m^-3,8.021e11 m^-3,5.637e10 m^-3,9.2991e+11 m^-3],
[865000 m,761393 m,1.376e7 m^-3,6.638e10 m^-3,6.029e4 m^-3,2.761e0 m^-3,7.875e11 m^-3,5.612e10 m^-3,9.1001e+11 m^-3],
[870000 m,765264 m,1.211e7 m^-3,6.171e10 m^-3,5.210e4 m^-3,2.301e0 m^-3,7.733e11 m^-3,5.586e10 m^-3,8.9088e+11 m^-3],
[875000 m,769130 m,1.066e7 m^-3,5.737e10 m^-3,4.503e4 m^-3,1.918e0 m^-3,7.593e11 m^-3,5.560e10 m^-3,8.7228e+11 m^-3],
[880000 m,772991 m,9.380e6 m^-3,5.334e10 m^-3,3.892e4 m^-3,1.599e0 m^-3,7.456e11 m^-3,5.535e10 m^-3,8.5430e+11 m^-3],
[885000 m,776846 m,8.258e6 m^-3,4.959e10 m^-3,3.365e4 m^-3,1.333e0 m^-3,7.321e11 m^-3,5.509e10 m^-3,8.3679e+11 m^-3],
[890000 m,780696 m,7.271e6 m^-3,4.612e10 m^-3,2.910e4 m^-3,1.112e0 m^-3,7.189e11 m^-3,5.484e10 m^-3,8.1987e+11 m^-3],
[895000 m,784541 m,6.404e6 m^-3,4.289e10 m^-3,2.517e4 m^-3,9.277e-1 m^-3,7.060e11 m^-3,5.459e10 m^-3,8.0349e+11 m^-3],
[900000 m,788380 m,5.641e6 m^-3,3.989e10 m^-3,2.177e4 m^-3,7.742e-1 m^-3,6.933e11 m^-3,5.434e10 m^-3,7.8754e+11 m^-3],
[905000 m,792214 m,4.970e6 m^-3,3.711e10 m^-3,1.884e4 m^-3,6.462e-1 m^-3,6.809e11 m^-3,5.410e10 m^-3,7.7212e+11 m^-3],
[910000 m,796043 m,4.379e6 m^-3,3.452e10 m^-3,1.631e4 m^-3,5.396e-1 m^-3,6.687e11 m^-3,5.385e10 m^-3,7.5707e+11 m^-3],
[915000 m,799866 m,3.859e6 m^-3,3.212e10 m^-3,1.411e4 m^-3,4.506e-1 m^-3,6.567e11 m^-3,5.361e10 m^-3,7.4243e+11 m^-3],
[920000 m,803685 m,3.402e6 m^-3,2.989e10 m^-3,1.222e4 m^-3,3.764e-1 m^-3,6.450e11 m^-3,5.336e10 m^-3,7.2825e+11 m^-3],
[925000 m,807498 m,2.999e6 m^-3,2.781e10 m^-3,1.058e4 m^-3,3.145e-1 m^-3,6.335e11 m^-3,5.312e10 m^-3,7.1443e+11 m^-3],
[930000 m,811305 m,2.645e6 m^-3,2.588e10 m^-3,9.165e3 m^-3,2.629e-1 m^-3,6.222e11 m^-3,5.288e10 m^-3,7.0096e+11 m^-3],
[935000 m,815108 m,2.332e6 m^-3,2.409e10 m^-3,7.940e3 m^-3,2.197e-1 m^-3,6.111e11 m^-3,5.264e10 m^-3,6.8783e+11 m^-3],
[940000 m,818905 m,2.057e6 m^-3,2.242e10 m^-3,6.880e3 m^-3,1.837e-1 m^-3,6.003e11 m^-3,5.241e10 m^-3,6.7513e+11 m^-3],
[945000 m,822697 m,1.815e6 m^-3,2.088e10 m^-3,5.962e3 m^-3,1.537e-1 m^-3,5.896e11 m^-3,5.217e10 m^-3,6.6265e+11 m^-3],
[950000 m,826484 m,1.602e6 m^-3,1.944e10 m^-3,5.168e3 m^-3,1.286e-1 m^-3,5.792e11 m^-3,5.194e10 m^-3,6.5058e+11 m^-3],
[955000 m,830266 m,1.414e6 m^-3,1.810e10 m^-3,4.481e3 m^-3,1.076e-1 m^-3,5.689e11 m^-3,5.170e10 m^-3,6.3870e+11 m^-3],
[960000 m,834043 m,1.248e6 m^-3,1.685e10 m^-3,3.886e3 m^-3,9.004e-2 m^-3,5.589e11 m^-3,5.147e10 m^-3,6.2722e+11 m^-3],
[965000 m,837814 m,1.102e6 m^-3,1.569e10 m^-3,3.370e3 m^-3,7.538e-2 m^-3,5.490e11 m^-3,5.124e10 m^-3,6.1593e+11 m^-3],
[970000 m,841580 m,9.726e5 m^-3,1.462e10 m^-3,2.924e3 m^-3,6.312e-2 m^-3,5.393e11 m^-3,5.101e10 m^-3,6.0493e+11 m^-3],
[975000 m,845342 m,8.590e5 m^-3,1.362e10 m^-3,2.537e3 m^-3,5.287e-2 m^-3,5.298e11 m^-3,5.078e10 m^-3,5.9420e+11 m^-3],
[980000 m,849098 m,7.587e5 m^-3,1.268e10 m^-3,2.201e3 m^-3,4.430e-2 m^-3,5.205e11 m^-3,5.056e10 m^-3,5.8374e+11 m^-3],
[985000 m,852849 m,6.703e5 m^-3,1.182e10 m^-3,1.911e3 m^-3,3.712e-2 m^-3,5.114e11 m^-3,5.033e10 m^-3,5.7355e+11 m^-3],
[990000 m,856594 m,5.922e5 m^-3,1.101e10 m^-3,1.659e3 m^-3,3.111e-2 m^-3,5.024e11 m^-3,5.011e10 m^-3,5.6352e+11 m^-3],
[995000 m,860335 m,5.234e5 m^-3,1.026e10 m^-3,1.440e3 m^-3,2.609e-2 m^-3,4.936e11 m^-3,4.989e10 m^-3,5.5375e+11 m^-3],
[1000000 m,864071 m,4.626e5 m^-3,9.562e9 m^-3,1.251e3 m^-3,2.188e-2 m^-3,4.850e11 m^-3,4.967e10 m^-3,5.4423e+11 m^-3]]
}
"StandardAtmosphere.frink included"
Download or view StandardAtmosphere.frink in plain text format
This is a program written in the programming language Frink.
For more information, view the Frink
Documentation or see More Sample Frink Programs.
Alan Eliasen was born 20136 days, 5 hours, 8 minutes ago.