[cig-commits] [commit] inversion: further cleaning (eb459ac)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/burnman

On branch  : inversion
Link       : https://github.com/geodynamics/burnman/compare/1098a3bf1b7cf8d21e6a06c1cca7639f1b95764a...eb459ac5beb7a02a772b89cf937922150abf01af

>---------------------------------------------------------------

commit eb459ac5beb7a02a772b89cf937922150abf01af
Author: sannecottaar <sanne.cottaar at gmail.com>
Date:   Sat Dec 13 14:30:35 2014 -0800

    further cleaning


>---------------------------------------------------------------

eb459ac5beb7a02a772b89cf937922150abf01af
 .../setup_isochemical_newmisfit_moduli_realdata.py  | 21 ++-------------------
 1 file changed, 2 insertions(+), 19 deletions(-)

diff --git a/inversion/setup_isochemical_newmisfit_moduli_realdata.py b/inversion/setup_isochemical_newmisfit_moduli_realdata.py
index cb9056f..16f2971 100644
--- a/inversion/setup_isochemical_newmisfit_moduli_realdata.py
+++ b/inversion/setup_isochemical_newmisfit_moduli_realdata.py
@@ -20,31 +20,14 @@ from burnman import minerals
 
 
 
-number_of_points = 10 #set on how many depth slices the computations should be done
+number_of_points = 5 #set on how many depth slices the computations should be done
 
 
 
 # velocity constraints from seismology
 seismic_model = burnman.seismic.PREM() # pick from .prem() .slow() .fast() (see code/seismic.py)
 depths = np.linspace(1000e3,2500e3, number_of_points)
-seis_p, seis_rho, seis_vp, seis_vs, seis_vphi = seismic_model.evaluate_all_at(depths)
-
-
-seismic_model2 = burnman.seismic.PREM() # pick from .prem() .slow() .fast() (see code/seismic.py)
-depths = np.linspace(1000e3,2500e3, number_of_points)
-seis_p2, seis_rho2, seis_vp2, seis_vs2, seis_vphi2 = seismic_model2.evaluate_all_at(depths)
-
-'''
-print seis_vs-seis_vs2
-plt.plot(depths,seis_vs,'r')
-plt.plot(depths,seis_vs2,'k')
-plt.plot(depths,seis_vp,'r')
-plt.plot(depths,seis_vp2,'k')
-#plt.plot(depths,seis_rho)
-#plt.plot(depths,seis_rho2)
-plt.show()
-stophere
-'''
+seis_p, seis_rho, seis_vp, seis_vs, seis_vphi = seismic_model.evaluate_all_at(depths
 
 seis_G= seis_vs**2.*seis_rho
 seis_K= seis_vphi**2*seis_rho