SkyMapMod 0.1.58__tar.gz → 0.2.0__tar.gz

{skymapmod-0.1.58 → skymapmod-0.2.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: SkyMapMod
-Version: 0.1.58
+Version: 0.2.0
 Summary: for modeling
 Project-URL: Homepage, https://istina.msu.ru/profile/aminevann/
 Author-email: Anna Amineva <amineva.2219@physics.msu.ru>

{skymapmod-0.1.58 → skymapmod-0.2.0}/SkyMapMod/__init__.py

@@ -7,5 +7,4 @@ from .solar_radio_flux import *
 from .star_catalogues import *
 from .albedo_of_planets import *
 from .planets import *
-from .total_background import *
-from data.load_data import *
+from .total_background import *

{skymapmod-0.1.58 → skymapmod-0.2.0}/SkyMapMod/airglow.py

@@ -5,6 +5,7 @@ import matplotlib.pyplot as plt
 from .solar_radio_flux import fluxdate
 from .solar_radio_flux import fluxtime
 from .solar_radio_flux import fluxobsflux
+
 from .zodiac import cos
 
 from .data.load_data import load_wavelenght_kp, load_intensity_kp
@@ -27,7 +28,7 @@ def radioflux(date, time): #date -- в формате строки 'дд.мм.г
 wavelenght_kp = load_wavelenght_kp()
 intensity_kp = load_intensity_kp() 
 
-def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere_kp = wavelenght_modtran_kp, transparency_atmosphere_kp = trancparency_modtran_kp, zenith_angle = 0):
+def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere_kp = wavelenght_modtran_kp, transparency_atmosphere_kp = trancparency_modtran_kp, z = 0):
     """
     Возвращает длины волн и спектр собственного свечения за атмосферой.
 
@@ -53,6 +54,7 @@ def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = int
         - Для совпадающих длин волн значение интенсивности делится на соответствующую прозрачность.
         - Пересчёт из "Рэлей/Ангстрем" в "фотоны/(с·м²·бин·нм)" производится по формуле:
           intensity = intensity * 10^11 / (4 * pi) * pi^2 / (1800 * 1800)
+        - Учитывает толщину слоя путем деления на cos зенитного угла
     """
     wavelenght = []
     intensity = []
@@ -61,7 +63,7 @@ def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = int
             if wavelenght_airglow[i] == wavelenght_atmosphere_kp[j]:
                 wavelenght.append(wavelenght_airglow[i])
                 intensity.append(intensity_airglow[i] / transparency_atmosphere_kp[j])
-    intensity = np.array(intensity) / cos(zenith_angle)
+    intensity = np.array(intensity) / (cos(z) + 10**(-30))
     wavelenght = np.array(wavelenght)
     return(wavelenght, intensity)   #возвращает длину волны в нм и поток в фот / (сек м^2 ср нм) + от стерадиан к сетке
 

skymapmod-0.2.0/SkyMapMod/data/__init__.py

@@ -0,0 +1,3 @@
+from .load_data import load_star_temperatures
+
+__all__ = ['load_star_temperatures']

{skymapmod-0.1.58 → skymapmod-0.2.0}/SkyMapMod/planets.py

@@ -110,7 +110,7 @@ def mag_to_phot(mag): #в полосе V
     """
     F0 = 10**2
     F1 = F0 * 10 ** (0.4 * (10 - mag)) # фот / (см^2 сек)
-    F1 = F1 * 10**4 * 10**(-2) #/ перешли к м^2 и от бина как бы
+    F1 = F1 * 10**4 * 10**(-2) #/ на метр квадратный от бина
     return F1
 
 

{skymapmod-0.1.58 → skymapmod-0.2.0}/SkyMapMod/star_catalogues.py

@@ -34,7 +34,6 @@ def star_spectrum(l, b, lmbd):
             - spectrum = A * black_body_phot(lmbd, T),
             где A — амплитуда, T — температура, lmbd — длина волны.
         - Функция `black_body_phot` определена отдельно.
-        - Учитывается порядок амплитуды с помощью множителя 1e-30
 
 
     Зависимости:
@@ -50,12 +49,11 @@ def star_spectrum(l, b, lmbd):
         j =  int(1800 - round((l+10)*10))
      
     i = int(900 + round(b*10))
-    print('Техническая строчка: i = ', i, ', j = ', j, sep='')
     star_brightness = load_star_brightness()
     star_temperatures = load_star_temperatures()    
-    A = star_brightness[i][j] #учитывается порядок амплитуды, которая хранится в больших числах
+    A = star_brightness[i][j]
     T = star_temperatures[i][j]
     print('Информация по звезде в этом пикселе: T = ', T, ', A = ', A, sep='')
     print('Техническая строчка: i = ', i, ', j = ', j, sep='')
-    spectrum = A * black_body_phot(lmbd, T) * 10**(-30) 
+    spectrum = A * black_body_phot(lmbd, T)
     return spectrum

{skymapmod-0.1.58 → skymapmod-0.2.0}/SkyMapMod/total_background.py

@@ -230,7 +230,7 @@ def total_background(l, b, date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_s
 #     return(zodiac_wl, zodiac_spec)
 
     #собственное свечение достаем
-    airglow_wl, airglow_spec = airglow_spectrum(wavelenght_airglow, intensity_airglow, wavelenght_atmosphere_kp, transparency_atmosphere_kp, zenith_angle=0)
+    airglow_wl, airglow_spec = airglow_spectrum(wavelenght_airglow, intensity_airglow, wavelenght_atmosphere_kp, transparency_atmosphere_kp)
     print('ФОТОНОВ СОБСТВЕННОГО:', integral(airglow_wl, airglow_spec))
 
 #     return(airglow_wl, airglow_spec)

{skymapmod-0.1.58 → skymapmod-0.2.0}/pyproject.toml

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