SkyMapMod 0.1.7__tar.gz → 0.1.9__tar.gz

{skymapmod-0.1.7 → skymapmod-0.1.9}/PKG-INFO

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

{skymapmod-0.1.7 → skymapmod-0.1.9}/SkyMapMod/zodiac.py

@@ -176,6 +176,13 @@ 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 
@@ -223,4 +230,24 @@ def N_S10_2(Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_
     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)
+
+#предположение о том, что спектр солнечный у меня все-таки лежит в Вт / (м^2 ср нм)
+def N_S10_3(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) #Вт/ (м^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) # Вт / (м^2 ср)
+    mean_Sun = np.sum(Sun_V_fx) # Вт / (м^2 ср)
+    N_S10_to_Sun = 100 * mean_A0 / mean_Sun
     return (N_S10_to_Sun)

{skymapmod-0.1.7 → skymapmod-0.1.9}/pyproject.toml

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