From ae08f496d7d8fcda119196033f803ac8dddc9110 Mon Sep 17 00:00:00 2001
From: Marcel Frommelt <lrt86919@stud.uni-stuttgart.de>
Date: Mon, 3 May 2021 17:44:16 +0900
Subject: [PATCH] seperate script for gravity model

---
 models/gravity.py | 10 ++++++++++
 1 file changed, 10 insertions(+)
 create mode 100644 models/gravity.py

diff --git a/models/gravity.py b/models/gravity.py
new file mode 100644
index 0000000..e8aa102
--- /dev/null
+++ 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
\ No newline at end of file