real_potential
Calculates the real potential at a given latitude and height using coefficients from a gravity model
Calling Sequence
from geoid_toolkit.real_potential import real_potential
W, dW_dr, dW_dtheta = real_potential(lat, lon, h, clm, slm, lmax, R, GM)
- geoid_toolkit.real_potential(lat, lon, h, refell, clm, slm, lmax, R, GM, GAUSS=0)[source]
Calculates the real potential using gravity model coefficients following Barthelmes [1], Hofmann-Wellenhof and Moritz [4], Moazezi and Zomorrodian [8], Molodensky [9]
- Parameters:
- lat: float
latitude in degrees
- lon: float
longitude in degrees
- h: float
ellipsoidal height in meters
- refell: str
Reference ellipsoid name
'CLK66': Clarke 1866'GRS67': Geodetic Reference System 1967'GRS80': Geodetic Reference System 1980'HGH80': Hughes 1980 Ellipsoid'WGS72': World Geodetic System 1972'WGS84': World Geodetic System 1984'ATS77': Quasi-earth centred ellipsoid for ATS77'NAD27': North American Datum 1927'NAD83': North American Datum 1983'INTER': International'KRASS': Krassovsky (USSR)'MAIRY': Modified Airy (Ireland 1965/1975)'TOPEX': TOPEX/POSEIDON ellipsoid'EGM96': EGM 1996 gravity model
- clm: float
cosine spherical harmonics for a gravity model
- slm: float
sine spherical harmonics for a gravity model
- lmax: int
maximum spherical harmonic degree
- R: float
average radius used in gravity model
- GM: float
geocentric gravitational constant used in gravity model
- GAUSS: float, default 0
Gaussian Smoothing Radius in km
- Returns:
- W: float
real potential at height h
- dW_dr: float
derivative of real potential with respect to radius
- dW_dtheta: float
derivative of real potential with respect to theta