[cig-commits] [commit] inversion: rename file (635dc1a)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Sat Dec 13 15:14:26 PST 2014
Repository : https://github.com/geodynamics/burnman
On branch : inversion
Link : https://github.com/geodynamics/burnman/compare/eb459ac5beb7a02a772b89cf937922150abf01af...635dc1a53d1fabfd5d1b2a13b45421af951ad5e7
>---------------------------------------------------------------
commit 635dc1a53d1fabfd5d1b2a13b45421af951ad5e7
Author: sannecottaar <sanne.cottaar at gmail.com>
Date: Sat Dec 13 15:13:46 2014 -0800
rename file
>---------------------------------------------------------------
635dc1a53d1fabfd5d1b2a13b45421af951ad5e7
...chemical_newmisfit_moduli_realdata.py => basic_inversion_setup.py} | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/inversion/setup_isochemical_newmisfit_moduli_realdata.py b/inversion/basic_inversion_setup.py
similarity index 99%
copy from inversion/setup_isochemical_newmisfit_moduli_realdata.py
copy to inversion/basic_inversion_setup.py
index 16f2971..cfe8c4f 100644
--- a/inversion/setup_isochemical_newmisfit_moduli_realdata.py
+++ b/inversion/basic_inversion_setup.py
@@ -27,11 +27,13 @@ number_of_points = 5 #set on how many depth slices the computations should be do
# 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
+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
+
+
temperature = burnman.geotherm.brown_shankland(seis_p)
'''
More information about the CIG-COMMITS
mailing list