GriddedTopography2AtmospherePotentialCoefficients

Estimate interior and exterior potential coefficients for atmosphere above digital terrain models. Coefficients for interior $(1/r)^{n+1}$ and exterior ($r^n$) are computed. The density of the atmosphere is assumed to be (Sjöberg, 1998) \[ \rho_0\left(\frac{R}{R+h}\right)^\nu, \]where $R$ is the radial distance of the ellipsoid at each point, $h$ the radial height above the ellipsoid, $\rho_0$ is densitySeaLevel and nu $\nu$ is a constant factor. The density is integrated from radialLowerBound and upperAtmosphericBoundary above the ellipsoid. The radialLowerBound is typically the topography and can be computed as expression at every point from inputfileGriddedData.

NameTypeAnnotation
outputfilePotentialCoefficientsExterior
filename
outputfilePotentialCoefficientsInterior
filename
inputfileGriddedData
filenameDigital Terrain Model
densitySeaLevel
double[kg/m**3]
ny
doubleConstant for Atmosphere
radialLowerBound
expressionexpression (variables 'L', 'B', 'height', 'data', and 'area' are taken from the gridded data
upperAtmosphericBoundary
doubleconstant upper bound [m]
factor
doublethe result is multiplied by this factor, set -1 to subtract the field
minDegree
uint
maxDegree
uint
GM
doubleGeocentric gravitational constant
R
doublereference radius
This program is parallelized.