seperate script for gravity model

2021-05-03 17:44:16 +09:00 · 2021-05-03 17:44:16 +09:00 · ae08f496d7
parent 16e5236591
commit ae08f496d7
1 changed files with 10 additions and 0 deletions
--- a/models/gravity.py
+++ b/models/gravity.py
@ -0,0 +1,10 @@
+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