gravity_anomaly
Calculates the gravity anomaly at a given latitude and longitude using different methods
Calling Sequence
from geoid_toolkit.gravity_anomaly import gravity_anomaly
ddelta_g = gravity_anomaly(lat, lon, h, 'WGS84', clm, slm, lmax, R, GM, METHOD='first')
- geoid_toolkit.gravity_anomaly(lat, lon, h, refell, clm, slm, lmax, R, GM, METHOD='first', GAUSS=0)[source]
Calculates the gravity anomaly for a given method 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
- METHOD: str
Method for calculating gravity anomalies
'first': classic first approximation method'second': classic second approximation method'molodensky': Molodensky method [9]
- GAUSS: float, default 0
Gaussian Smoothing Radius in km
- Returns:
- ddelta_g: float
gravity anomaly for a given ellipsoid in meters