10 lines
311 B
Python
10 lines
311 B
Python
|
import numpy as np
|
||
|
|
|
||
|
|
|
||
|
|
def grav(lat, h):
|
||
|
|
g0 = 9.7803253359 * ((1 + 0.00193185265241 * np.sin(np.deg2rad(lat)) ** 2) / (
|
||
|
|
np.sqrt(1 - 0.00669437999013 * np.sin(np.deg2rad(lat)) ** 2)))
|
||
|
|
|
||
|
|
res = g0 - (3.0877 * 1e-6 - 4.3 * 1e-9 * np.sin(np.deg2rad(lat)) ** 2) * h + 7.2 * 1e-13 * h ** 2
|
||
|
|
|
||
|
|
return res
|