SkyMapMod 0.0.1__tar.gz → 0.0.3__tar.gz

{skymapmod-0.0.1 → skymapmod-0.0.3}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: SkyMapMod
-Version: 0.0.1
+Version: 0.0.3
 Project-URL: Homepage, https://istina.msu.ru/profile/aminevann/
 Author-email: Anna Amineva <amineva.2219@physics.msu.ru>
 License-File: LICENSE

{skymapmod-0.0.1 → skymapmod-0.0.3}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "SkyMapMod"
-version = "0.0.1"
+version = "0.0.3"
 authors = [
  { name="Anna Amineva", email="amineva.2219@physics.msu.ru" },
 ]

skymapmod-0.0.3/src/SkyMapMod/__init__.py

@@ -0,0 +1 @@
+from .zodiac import *

skymapmod-0.0.3/src/SkyMapMod/zodiac.py

@@ -0,0 +1,139 @@
+import numpy as np
+import math
+import matplotlib.pyplot as plt
+
+def gret():
+    print('Hello')
+
+#переопределяю тригонометрические функции, чтобы работать с градусами
+def cos(x):
+    return math.cos(x * math.pi / 180)
+
+def sin(x):
+    return math.sin(x * math.pi / 180)
+
+def tan(x):
+    return math.tan(x * math.pi / 180)
+
+def acos(x):
+    return math.acos(x) * 180 / math.pi
+
+def asin(x):
+    return math.asin(x) * 180 / math.pi
+
+def atan(x):
+    return math.atan(x) * 180 / math.pi
+
+#переопределяю прочие функции для красоты записи
+def exp(x):
+    return math.exp(x)
+
+def sqrt(x):
+    return math.sqrt(x)
+
+#функция Хевисайда (в контексте статьи определяется именно так)
+def u(x):
+    if x <= 0:
+        return 0
+    elif x > 0:
+        return 1
+    
+#функции для перехода из одной СК в другую
+#сферические в декартовы
+def spheral_to_decart(r, lmbd, beta):
+    x = r * cos(lmbd) * cos(beta)
+    y = r * sin(lmbd) * cos(beta)
+    z = r * sin(beta)
+    return (x, y, z)
+
+#декартовы в сферические
+def decart_to_spheral(x, y, z):
+    r = sqrt(x**2 + y**2 + z**2)
+    lmbd = atan(y/(x + 1e-30))
+    beta = asin(z/r)
+    return(r, lmbd, beta)
+
+#функция для поворота эклиптических координат
+def rotate_ekl(lmbd, beta, lmbd0, beta0):
+    x, y, z = spheral_to_decart(1, lmbd-lmbd0, beta)
+    x1 = x * cos(beta0) + z * sin(beta0) 
+    z1 = - x * sin(beta0) + z * cos(beta0)
+    y1 = y
+    r, lmbd2, beta2 = decart_to_spheral(x1, y1, z1)
+    return(lmbd2, beta2)
+
+#Далее идут функции, которые входят в общую формулу для расчета зод. света
+#Компонента S
+def S(lmbd, beta, lmbd_sun, Omega):
+    c = cos(lmbd - lmbd_sun) * cos(beta)
+    eps = acos(c)
+    first = 6 * abs(sin(lmbd_sun - Omega))
+    second = ((1-u(eps-90)) * (sin(eps) + 1e-30)**(-2.3) + u(eps - 90) * sin(eps))
+    #third = (1 - u(Omega - lmbd_sun) / 4 + 2 * (1-u(eps-90)) * c)
+    third = (1 - (1-u(lmbd_sun - Omega))/4 + 2*(1-u(eps-90)) * c)
+    #equap = ((u(Omega - lmbd_sun) - u(lmbd_sun - Omega)) * beta + 5) / 10
+    equap = ((1 - u(lmbd_sun - Omega) - u(lmbd_sun - Omega)) * beta + 5) / 10
+    fourth = max(0, min(equap, 1))
+    S = first * second * third * fourth
+    return S
+
+#Компонента D -- "эмпирический вклад в форме гантели"
+def D(lmbd, beta, lmbd_sun):
+    
+    d = abs(lmbd - lmbd_sun) / 6.5 - abs(beta) + 15 + 5 * u(beta)
+    eps = acos(cos(lmbd - lmbd_sun) * cos(beta))
+    if beta > 0:
+        lmbd1, beta1 = rotate_ekl(lmbd, beta, lmbd_sun, 21)
+    elif beta <= 0:
+        lmbd1, beta1 = rotate_ekl(lmbd, beta, lmbd_sun, -15)
+    gamma1 = atan(sin(lmbd1)/( 1e-30 + tan(beta1)))
+    eps1 = acos(cos(lmbd1) * cos(beta1))
+    h = gamma1 * (1 - u(beta)) + (gamma1) * u(beta)  
+    A = (1 + 0.5 * (1 - u(beta))) * (1-u(eps-90)) * ((1-u(eps-75)) + u(eps - 75) * exp(-(eps - 75)**2 / 120))
+#     f = 60 * exp(- eps1 / 16) * (0.15 + 0.85 * exp(-6 * cos(0.6 * h)**4))
+    f = 60 * exp(- eps1 / 16) * (0.85 + 0.15 * exp(-6 * cos(0.6 * h)**4))
+    g = 25 * exp(- acos(cos(2 * lmbd1) * cos(0.7 * beta1)) / 10)
+    D = A * ((f * math.exp(-d / 10) + g * math.exp(-d / 8)) * u(d) + (f + g) * (1 - u(d)))
+    return D
+
+#Компонента E -- остаточная, не учитывалась
+def E(lmbd, beta):
+    return 0
+
+#Компонента F -- остаточная, не учитывалась
+def F(lmbd, beta, lmbd_sun):
+    return 0
+
+#Компонента G -- антисолнечная точка
+def G(lmbd, beta, lmbd_sun):
+    lmbd_tild = lmbd - lmbd_sun - 180
+    if lmbd_tild > 180:
+        lmbd_tild = lmbd_tild - 360
+    elif lmbd_tild < -180:
+        lmbd_tild = lmbd_tild + 360        
+    eps_tild = 180 - acos(cos(lmbd_tild + 180) * cos(beta))    
+    G_b         = 7.5 * exp(-eps_tild / 4) + 39.5 * exp(-eps_tild / 25)
+    G_lmbd_tild = 7.5 * exp(-eps_tild / 4) + 39.5 * exp(-eps_tild / 35)  
+    first = 1 - 0.02 * beta * lmbd_tild**2 / (eps_tild**3 + 1e-300)
+    second = (beta**2 * G_b + lmbd_tild**2 * G_lmbd_tild) / (eps_tild**2 + 1e-300)
+    third = (1 - u(eps_tild - 60) * (1 - exp(-(eps_tild - 60)**2 / 300)))
+    G = first * second * third
+    return G
+
+#Общая формула для зод. света
+def zodiacal_light(lmbd, beta, lmbd_sun):
+    b = 1.5 * (sqrt(1 + (beta / 1.5)**2) - 1)
+    c = cos(lmbd - lmbd_sun) * cos(beta)
+    eps = acos(c)
+    eps_tild = 180 - eps
+    lmbd_tild = lmbd - lmbd_sun - 180
+    abs_gm = abs(atan(sin(lmbd - lmbd_sun) / (tan(beta) + 1e-30)))
+    Omega = 78.25
+    first = 7 + 8 * (1 - cos(b)) + 6 * exp(- beta**2 / 512)
+    second_1 = (1-u(eps-90)) * (65 + 120 * c - 185 * c**2 + 65 * c**3) * (sin(eps) + 1e-30)**(-2.3)
+    second_2 = u(eps - 90) * (65 + 120 * c + 154 * c**2 + 88 * c**3)
+    second_mult = 10 ** (- sin(b) / (0.009 * (eps + 40)))  
+    third = (1-u(eps-90)) * (30 * ((sin(eps) + 1e-30)**(-2.3) - 1) * cos(b))
+    fourth = (8800 * exp((1 - sqrt(1 + ((abs_gm - 90)/3)**2))/10) - 1200) * exp(-eps/10)
+    zod = first  + (second_1 + second_2) * second_mult + third + fourth + S(lmbd, beta, lmbd_sun, Omega) + D(lmbd, beta, lmbd_sun) + E(lmbd, beta) + F(lmbd, beta, lmbd_sun) + G(lmbd, beta, lmbd_sun)
+    return zod/0.415 #деление на 0.415 осуществляет переход из ед. ADU в ед. S10