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

SkyMapMod/zodiac.py

@@ -4,7 +4,7 @@ import random
 
 #коэффициент пересчета из единиц S10 в количество фотонов от Солнца
 #вычислен на основе спектра АЧТ от звезды класса A0, солнечного спектра и факта, что 1 ед. S10 соответствует 100 фотонов от звезды класса A0
-N_S10 = 100.63910099562007
+
 
 #переопределяю тригонометрические функции, чтобы работать с градусами
 def cos(x):
@@ -182,18 +182,25 @@ def integral(wl, spec):
         result += (spec[i] + spec[i+1]) / 2 * (wl[i+1] - wl[i])
     return result
 
-
+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):
+    #коэффициент пересчета из единиц S10 в количество фотонов от Солнца
+    #вычислен на основе спектра АЧТ от звезды класса A0, солнечного спектра и факта, что 1 ед. S10 соответствует 100 фотонов от звезды класса A0
+    #вот столько фотонов/(см^2 сек) от Солнца несет 1 ед S10:
+    N_S10 = 100.6 
+    N_S10 = N_S10 * 10**4 / (4 * np.pi) # фот / (м^2 сек ср) -- НАДО ЛИ ДЕЛИТЬ НА 4 ПИ?
+    S10 = zodiacal_light(lmbd, beta, lmbd_sun) 
+    Sun_V_wl, Sun_V_sp = convolution(Sun_sp_wl, Sun_sp_fx, V_wl, V_tr) #свертка Солнечного спектра с полосой V
+    integr = integral(Sun_V_wl, Sun_V_sp)
+    A = S10 * N_S10 / integr
     
-
-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
-    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])
-    for i in range(N_phot):
-        spectrum = random.choices(Sun_sp_wl, weights=probabilities)
-        counts[np.where(Sun_sp_wl == spectrum)[0][0]] += 1
-    return(Sun_sp_wl, counts)
+    
+    #############################################
+    a, b = smm.convolution(Sun_sp_wl, A * Sun_sp_fx, V_wl, V_tr)
+    c = integr = smm.integral(a, b)
+    print('Техническая проверка: S10*N_S10 =', S10 * N_S10)
+    print('А при нормировке и интекрированию получается:', c)
+    #############################################
+    
+    return Sun_sp_wl, A * Sun_sp_fx #нм, фот / (м^2 сек нм ср)
     
     

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

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

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

@@ -9,12 +9,12 @@ SkyMapMod/solar_radio_flux.py,sha256=I02O6Mltv3lfqxif3BHUNAzeElQ5AZax_2L09A3cXb8
 SkyMapMod/solar_spectrum.py,sha256=LjBPSTI6Yq86pwF7haUja7RjRYH4C8HZL8SjyI7zp6g,24433
 SkyMapMod/star_catalogues.py,sha256=hl8-EMMCfJGz1MGBqdwqIY_d1-9G8j8cZGyWmE3OSBA,1166
 SkyMapMod/transparency.py,sha256=HiSbbVX80neh3n284f64sWnrbxpRzPGKvFyid9sVo7I,376
-SkyMapMod/zodiac.py,sha256=pgOvKoMoz6P0JFyFDihhxTzIgbLFlM-Ey6D08RnzxYU,8280
+SkyMapMod/zodiac.py,sha256=TKrJleZg7AWVFKNQD9RUS-5OkCvORq90f1S_I0VJexw,9103
 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.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,,
+skymapmod-0.1.19.dist-info/METADATA,sha256=6jGzhwGjxSGeyOi6qvxrdoQdE6YHz3wRs0nD-xrgi9w,469
+skymapmod-0.1.19.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+skymapmod-0.1.19.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
+skymapmod-0.1.19.dist-info/RECORD,,