SkyMapMod 0.1.6__py3-none-any.whl → 0.1.8__py3-none-any.whl

SkyMapMod/airglow.py

@@ -10,6 +10,7 @@ from .solar_radio_flux import fluxobsflux
 from .airglow_spectrum import wavelenght_kp, intensity_kp
 from .modtran_default_kp_transparency import wavelenght_modtran_kp, trancparency_modtran_kp
 
+#следующая функция для поправки на радиопоток, пока не используется
 def radioflux(date, time): #date -- в формате строки 'дд.мм.гггг', time -- в формате строки 'чч:мм:сс'
     day, month, year = date.split('.')
     date = year + month + day
@@ -31,7 +32,7 @@ def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = int
             if wavelenght_airglow[i] == wavelenght_atmosphere[j]:
                 wavelenght.append(wavelenght_airglow[i])
                 intensity.append(intensity_airglow[i] / transparency_atmosphere[j])
-    return(np.array(wavelenght), np.array(intensity))
+    return(np.array(wavelenght), np.array(intensity) * 10**10 / (4 * np.pi)) #возвращает длину волны в нм и поток в фот / (сек м^2 ср нм)
 #функция написана для спектра взятого от китт пик, прозрачность атмосферы взята из модтрана, параметры -- см. доклад, там прямо скрин. Прозрачность взята для 45 градусов, пересчитана в 30 градусов
 
 from .solar_radio_flux import fluxdate

SkyMapMod/star_catalogues.py

@@ -2,14 +2,7 @@ import numpy as np
 
 from .data.load_data import load_star_temperatures, load_star_brightness
 
-#спектр АЧТ в фотонах
-def black_body_phot(lmbd, T):
-    h = 6.63e-34 #Дж * сек
-    c = 3e8 # м / сек
-    k = 1.38e-23 # Дж / К
-    lmbd = lmbd * 1e-9 # м
-    u = 2 * c / (lmbd**4) / (np.exp(h * c / (lmbd) / k / T) - 1)
-    return u #уже в фотонах/сек (поделено на hc/lambda)
+from.zodiac import black_body_phot
 
 #на вход подаем галактические координаты (l -- долгота, b -- широта) и массив длин волн, по которым строим спектр
 def star_spectrum(l, b, lmbd):

SkyMapMod/zodiac.py

@@ -155,6 +155,16 @@ def black_body(lmbd, T):
     lmbd = lmbd * 1e-9 # м
     u = 2 * h * c**2 / (lmbd**5) / (np.exp(h * c / (lmbd) / k / T) - 1)
     return u 
+
+#спектр АЧТ в фотонах
+def black_body_phot(lmbd, T):
+    h = 6.63e-34 #Дж * сек
+    c = 3e8 # м / сек
+    k = 1.38e-23 # Дж / К
+    lmbd = lmbd * 1e-9 # м
+    u = 2 * c / (lmbd**4) / (np.exp(h * c / (lmbd) / k / T) - 1)
+    return u #уже в фотонах/сек (поделено на hc/lambda)
+
 #универсальная функция для свертки двух спектров
 def convolution(array_1, meaning_1, array_2, meaning_2):
     array = []
@@ -166,16 +176,25 @@ def convolution(array_1, meaning_1, array_2, meaning_2):
                 meaning.append(meaning_1[i] * meaning_2[j])
     return np.array(array), np.array(meaning)
 
+def integral(wl, spec):
+    shape = wl.shape[0]
+    result = 0
+    for i in range(shape - 1):
+        result += (spec[i] + spec[i+1]) / 2 * (wl[i+1] - wl[i])
+    return result
+
 def N_S10(Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
     h = 6.63e-34
     c = 3e8 
     A0_sp_wl = Sun_sp_wl
     A0_sp_fx = black_body(A0_sp_wl, T_A0)
+    
     A0_V_wl, A0_V_fx = convolution(A0_sp_wl, A0_sp_fx, V_wl, V_tr)
     A0_V_fx = A0_V_fx / (np.sum(A0_V_fx * A0_V_wl)) #надо домножать на A0_V_wl при нормировке.
 
     Sun_V_wl, Sun_V_fx = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr)
     Sun_V_fx = Sun_V_fx / (np.sum(Sun_V_fx * Sun_V_wl))
+    
     mean_A0 = np.sum(A0_V_fx * h * c / A0_V_wl * 1e9)
     mean_Sun = np.sum(Sun_V_fx * h * c / Sun_V_wl * 1e9)
     N_S10_to_Sun = 100 * mean_A0 / mean_Sun
@@ -194,4 +213,21 @@ def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003,
     return(Sun_sp_wl, counts)
     
     
-    
+def N_S10_2(Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
+    h = 6.63e-34
+    c = 3e8 
+    A0_sp_wl = Sun_sp_wl #нм
+    A0_sp_fx = black_body_phot(A0_sp_wl, T_A0) #фот/ (м^2 ср сек м)
+    A0_V_wl, A0_V_fx = convolution(A0_sp_wl, A0_sp_fx, V_wl, V_tr) #нм, фот / (м^2 сек ср м), стерадианы уйдут при нормировке
+    
+    A0_V_fx = A0_V_fx / (integral(A0_V_wl * 10**(-9), A0_V_fx)) #беру спектр и нормирую его на общее количество фотонов под спектром
+    #на выходе: фот / (м^2 сек ср)
+
+    Sun_V_wl, Sun_V_fx = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr) #нм, фот / (м^2  сек ср нм)
+    
+    Sun_V_fx = Sun_V_fx / (integral(Sun_V_fx, Sun_V_wl)) #на выходе: фот / (м^2 сек ср)
+    
+    mean_A0 = np.sum(A0_V_fx * h * c / A0_V_wl * 1e9) # Вт / (м^2 ср)
+    mean_Sun = np.sum(Sun_V_fx * h * c / Sun_V_wl * 1e9) # Вт / (м^2 ср)
+    N_S10_to_Sun = 100 * mean_A0 / mean_Sun
+    return (N_S10_to_Sun)

{skymapmod-0.1.6.dist-info → skymapmod-0.1.8.dist-info}/METADATA

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

{skymapmod-0.1.6.dist-info → skymapmod-0.1.8.dist-info}/RECORD

@@ -1,18 +1,18 @@
 SkyMapMod/__init__.py,sha256=Q0vaTOPF5a5ZVCmJpg6VYOnHAjGpoheLKM_6SewIuQM,271
-SkyMapMod/airglow.py,sha256=UvUeXz_9682CEZWsbmh36bvRNpjQGddll2tsIR9tX1M,2043
+SkyMapMod/airglow.py,sha256=4Uo4e2sk7Ux2Hc4MwgaZ7UrZwl2qe_ms0crCLOgXHKM,2293
 SkyMapMod/airglow_spectrum.py,sha256=4pY4QSGjIDvjs_jvxrowsZ9WLtIUwR0DesoCK6moWjs,127314
 SkyMapMod/band_V_data.py,sha256=fOj9H2IdZuqErAxnoLq4A02UXfYMgLNBiLjwrzjfWTs,15223
 SkyMapMod/modtran_default_kp_transparency.py,sha256=QeLYQA9FNUsKp-eRtOMQFxfq0pdwCnqVIl0X6Vvf7GU,78306
 SkyMapMod/solar_radio_flux.py,sha256=I02O6Mltv3lfqxif3BHUNAzeElQ5AZax_2L09A3cXb8,842
 SkyMapMod/solar_spectrum.py,sha256=LjBPSTI6Yq86pwF7haUja7RjRYH4C8HZL8SjyI7zp6g,24433
-SkyMapMod/star_catalogues.py,sha256=O1CpcGLZ1wDoFliMgetzjp68hp4eVB_nWRne0ognJ00,1034
+SkyMapMod/star_catalogues.py,sha256=6J1SzpCq93xG13mmw8rtABRgGHX3xwCk7TT75BQFs9M,740
 SkyMapMod/transparency.py,sha256=HiSbbVX80neh3n284f64sWnrbxpRzPGKvFyid9sVo7I,376
-SkyMapMod/zodiac.py,sha256=qbyYaciyGyQ_AVR5bGR7H0d1w00tze5_4jQi1siF1HU,8574
+SkyMapMod/zodiac.py,sha256=YEzm7RKab1azf-SOGlkYRgXIF7Vc3NoC4Swb32uJ-gs,10236
 SkyMapMod/data/__init__.py,sha256=_292qgM72hPa2vxDNSaI9QvPunvi-aSXUsZki99q63o,83
 SkyMapMod/data/load_data.py,sha256=p-Nq5St713HTTDunSWu9ANrF9V0Slaak17IhD0n_bSY,931
 SkyMapMod/data/star_brightness.npy,sha256=waKIZ0fIBM-92MgjGfgutQpXQrzTZkEMJXy8Kluvopg,51840128
 SkyMapMod/data/star_temperatures.npy,sha256=Yjp252bzqmKQz00K0ADTkstjzGFulRMgyW1b8kaDOPE,51840128
-skymapmod-0.1.6.dist-info/METADATA,sha256=X6WOCWj2DSvPY48MSVUwqZNxCePMoRmj3UC-eZRVKh4,468
-skymapmod-0.1.6.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-skymapmod-0.1.6.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
-skymapmod-0.1.6.dist-info/RECORD,,
+skymapmod-0.1.8.dist-info/METADATA,sha256=uM7bodtyGQn-8Gjq8EkAc7SDLrmL7VQyQvBRfXy9aXw,468
+skymapmod-0.1.8.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+skymapmod-0.1.8.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
+skymapmod-0.1.8.dist-info/RECORD,,