Diskontinuitäten haben jetzt auch wirklich einen Sprung bei konstanter Tiefe!

AdamsWilliamsonCMB.py

@@ -78,15 +78,15 @@ def calc_density(rho_0: int, rho_CMB: int, radius: list, vp: list, vs: list) ->
             rho_ex[i] = rho_CMB
             rho_im[i] = 9500
         elif depth[i] == 5153.5 and depth[i-1] == 5153.5: # ICB
-            rho_ex[i] = 12814.285714285714
+            rho_ex[i] = 12700
 
         # more discontinuities
-        elif depth[i] == 410 and depth[i-1] == 410:   # 410km discontinuity
-           rho_ex[i] = 3921.4285714285716
-        elif depth[i] == 660 and depth[i-1] == 660:   # 660km discontinuity
-           rho_ex[i] = 4238.7
-        elif depth[i] == 2740 and depth[i-1] == 2740:   # D" discontinuity
-           rho_ex[i] = 5693.4
+        # elif depth[i] == 410 and depth[i-1] == 410:   # 410km discontinuity
+        #    rho_ex[i] = 3921.4285714285716
+        # elif depth[i] == 660 and depth[i-1] == 660:   # 660km discontinuity
+        #    rho_ex[i] = 4238.7
+        # elif depth[i] == 2740 and depth[i-1] == 2740:   # D" discontinuity
+        #    rho_ex[i] = 5693.4 
         else:
             rho_ex = density_expansion_explicit(rho_ex, vp, vs, radius, i)
             # rho_im = density_expansion_implicit(rho_im, vp, vs, r, i)
@@ -135,7 +135,7 @@ def best_inertia_model(file, n=31, rho_0_low=3500, rho_0_high=3700, rho_CMB_low=
     plt.legend()
     plt.show()
 
-def best_mass_model(file, n=31, rho_0_low=3700, rho_0_high=3900, rho_CMB_low=9900, rho_CMB_high=10050):
+def best_mass_model(file, n=31, rho_0_low=3700, rho_0_high=3900, rho_CMB_low=9800, rho_CMB_high=10050):
     depth, density, vp, vs = read_data(file)
 
     M_0 = 5.972e24  # Masse in 18km Tiefe
@@ -185,9 +185,9 @@ file = os.getcwd() + r"\VpVs-valuesak135.csv"
 
 print('Calculating optimal density model in regard to total mass. Depending on the grid size of the search domain this may take a while. \n')
 
-best_mass_model(file,n=20,rho_0_low=3750,rho_0_high=3800)
+best_mass_model(file,n=20,rho_0_low=3600,rho_0_high=3800)
 
 print('\nCalculating optimal density model in regard to inertia. Depending on the grid size of the search domain this may take a while. \n')
 
-# best_inertia_model(file,n=15)
+best_inertia_model(file,n=15)