SkyMapMod 0.1.17__py3-none-any.whl → 0.1.18__py3-none-any.whl

SkyMapMod/star_catalogues.py

@@ -5,6 +5,7 @@ from .data.load_data import load_star_temperatures, load_star_brightness
 from.zodiac import black_body_phot
 
 #на вход подаем галактические координаты (l -- долгота, b -- широта) и массив длин волн, по которым строим спектр
+#РАЗОБРАТЬСЯ С КОЭФФИЦИЕНТАМИ А -- НАДО ЛИ НОРМИРОВАТЬ НА C (КАЖЕТСЯ НАДО) А ЕЩЕ ВРОДЕ Я ЕГО ОТНОРМИРОВАЛА И ПОДГРУЗИЛА -- ПРОВЕРИТЬ
 def star_spectrum(l, b, lmbd):
     j = int((1800 - round(l*10)) + 3600) #пересчет из галактических к позиции в массиве
     i = int(900 + round(b*10))

SkyMapMod/zodiac.py

@@ -4,7 +4,7 @@ import random
 
 #коэффициент пересчета из единиц S10 в количество фотонов от Солнца
 #вычислен на основе спектра АЧТ от звезды класса A0, солнечного спектра и факта, что 1 ед. S10 соответствует 100 фотонов от звезды класса A0
-k = 100.63910099562007
+N_S10 = 100.63910099562007
 
 #переопределяю тригонометрические функции, чтобы работать с градусами
 def cos(x):
@@ -182,27 +182,12 @@ def integral(wl, spec):
         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
-    return (N_S10_to_Sun) #ЭТО НА САНТИМЕТР КВАДРАТНЫЙ В СЕКУНДУ!!!! ВАЖНО!!!!!
     
 
-def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800):
+def zodiacal_spectrum(lmbd, beta, lmbd_sun, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, T_A0 = 9800, N_S10 = N_S10):
     S10 = zodiacal_light(lmbd, beta, lmbd_sun)
-    N_S10_to_Sun = N_S10(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr, T_A0)
+    N_S10_to_Sun = N_S10
     N_phot = round(S10 * N_S10_to_Sun)
     probabilities = Sun_sp_fx / np.sum(Sun_sp_fx) #или все-таки нормируем на интеграл?
     counts = np.zeros(Sun_sp_wl.shape[0])
@@ -212,30 +197,3 @@ 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 
-    Sun_V_wl, Sun_V_fx = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr) #сворачиваю с полосой V
-    
-    Sun_norm = integral(Sun_V_wl, Sun_V_fx)#интегрирую спектр, photons/(m^2 sec sr)
-    Sun_V_fx_norm = Sun_V_fx / Sun_norm #спектр отнормирован на 1 фотон
-
-    Sun_V_fx_norm_Vt = Sun_V_fx_norm * h * c / (Sun_V_wl * 10**(-9)) #перевела в мощности
-    Sun = np.sum(Sun_V_fx_norm_Vt)
-    
-    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)
-    norm_S10 = integral(A0_V_wl * 10**(-9), A0_V_fx)
-
-    A = 100/norm_S10
-
-    A0_V_fx = A * A0_V_fx
-    A0_norm = integral(A0_V_wl, A0_V_fx)
-    A0_V_fx_norm = A0_V_fx / A0_norm #photons/(m^2 sec sr)
-
-    A0_V_fx_norm_Vt = A0_V_fx_norm * h * c / (A0_V_wl * 10**(-9)) #Vt/(m^2 sec sr)
-    A0 = np.sum(A0_V_fx_norm_Vt)
-    N_S10_to_Sun = 100 * A0 / Sun
-    return (N_S10_to_Sun) #ЭТО НА САНТИМЕТР КВАДРАТНЫЙ В СЕКУНДУ

{skymapmod-0.1.17.dist-info → skymapmod-0.1.18.dist-info}/METADATA

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

{skymapmod-0.1.17.dist-info → skymapmod-0.1.18.dist-info}/RECORD

@@ -7,14 +7,14 @@ SkyMapMod/modtran_default_kp_transparency.py,sha256=QeLYQA9FNUsKp-eRtOMQFxfq0pdw
 SkyMapMod/planets.py,sha256=7xzna0cybGgDPOaNzx97hQ08jWYN9Ogbtk89nkId_RE,6702
 SkyMapMod/solar_radio_flux.py,sha256=I02O6Mltv3lfqxif3BHUNAzeElQ5AZax_2L09A3cXb8,842
 SkyMapMod/solar_spectrum.py,sha256=LjBPSTI6Yq86pwF7haUja7RjRYH4C8HZL8SjyI7zp6g,24433
-SkyMapMod/star_catalogues.py,sha256=QnSvH0dkRTQD87kGoG1pu7J5l6JP5Afmk61GTtij9K8,930
+SkyMapMod/star_catalogues.py,sha256=hl8-EMMCfJGz1MGBqdwqIY_d1-9G8j8cZGyWmE3OSBA,1166
 SkyMapMod/transparency.py,sha256=HiSbbVX80neh3n284f64sWnrbxpRzPGKvFyid9sVo7I,376
-SkyMapMod/zodiac.py,sha256=VcSOCsl1NL3MDqdfrl39EfBZwymFTQFDWylmBqRH69A,10406
+SkyMapMod/zodiac.py,sha256=pgOvKoMoz6P0JFyFDihhxTzIgbLFlM-Ey6D08RnzxYU,8280
 SkyMapMod/data/__init__.py,sha256=_292qgM72hPa2vxDNSaI9QvPunvi-aSXUsZki99q63o,83
 SkyMapMod/data/load_data.py,sha256=p-Nq5St713HTTDunSWu9ANrF9V0Slaak17IhD0n_bSY,931
 SkyMapMod/data/star_brightness.npy,sha256=OVrYGTl3I2wbWAMmFnDkunkRyPLL2-g-nD_JX_0FBdQ,51840128
 SkyMapMod/data/star_temperatures.npy,sha256=Yjp252bzqmKQz00K0ADTkstjzGFulRMgyW1b8kaDOPE,51840128
-skymapmod-0.1.17.dist-info/METADATA,sha256=57wHb2KX6HWXHoUtJJiGv1iL9CZnBecu2uCGD0-u2Do,469
-skymapmod-0.1.17.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-skymapmod-0.1.17.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
-skymapmod-0.1.17.dist-info/RECORD,,
+skymapmod-0.1.18.dist-info/METADATA,sha256=gsBcBZwbqwerABc92sEPE09CbeBIn8NUQRJTYdLo6Z4,469
+skymapmod-0.1.18.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+skymapmod-0.1.18.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
+skymapmod-0.1.18.dist-info/RECORD,,