diff --git a/core.py b/core.py
index 4970129..09fba3b 100644
--- a/core.py
+++ b/core.py
@@ -274,3 +274,7 @@ def tangent_line_k(line_x, line_y, center_x, center_y, radius, init_k=10.0):
             k_angle.append(math.pi + math.atan(kk))
     # 返回相对x轴最大的角度k
     return np.array(k_candidate)[np.max(k_angle) == k_angle].tolist()[-1]
+
+
+def func_ng(td):  # 地闪密度
+    return 0.023 * (td ** 1.3)