// 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] } // 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 // (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"