StandardAtmosphere.frink

View or download 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:
// http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19770009539_1977009539.pdf
// See: http://www.atmosculator.com/The%20Standard%20Atmosphere.html
//
//   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:
   //   [geopotentialHeight, seaLevelTemperature, seaLevelPressure]
   //  The latter two arguments have defaults at the defaults for the
   // Standard Atmposphere 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, 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] :=
   {
      return secantInterpolate[{|alt| getTemperatureAndPressure[alt]@1}, pressure, .01 m]
   }

   // Calculate the equivalent standard altitude that corresponds to a given
   // density. 
   class getDensityAltitude[density] :=
   {
      return secantInterpolate[{|alt|
                                  [t,p] = getTemperatureAndPressure[alt]
                                  getDensity[alt, t, p]
                               }, density, .01 m]
   }

   // Secant method solver.  Given a function f[x], returns the value of x
   // which yields the desired value.  maxDelta is a value that indicates the
   // largest error in x.
   class secantInterpolate[f, desired, maxDelta] :=
   {
      x1 = 0 km
      x2 = 86 km
      f1 = f[x1] - desired
      f2 = f[x2] - desired

      while (true)
      {
         correction = (f1 * (x1 - x2)) / (f1-f2)
         x = x1 - correction
//         println["correction is $correction"]
         if (x < -5000 m)
            x = -5000 m
         if (x > 1000 km)
            x = 1000 km

//         println["x1 = $x1"]
//         println["x2 = $x2"]

         if abs[correction] < maxDelta
            return x

         x2 = x1
         x1 = x
         f2 = f1

         f1 = f[x] - desired
      }

      return f1
   }

   // 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"


View or download 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 17645 days, 8 hours, 34 minutes ago.