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 17683 days, 14 hours, 52 minutes ago.