[cig-commits] [commit] master: Added warning for pressure in geotherm (b344b96)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Fri Dec 12 16:47:11 PST 2014
Repository : https://github.com/geodynamics/burnman
On branch : master
Link : https://github.com/geodynamics/burnman/compare/498861c40268a2c440d99f9ba62ccc5ab2a79dfa...1d2396e3c9d632ea6c4c7f8cbc61788c9e87ff80
>---------------------------------------------------------------
commit b344b96f4cbe53369dc7c7975f280466e4f1b392
Author: I <kmanunterborn at gmail.com>
Date: Fri Dec 12 13:56:23 2014 -0800
Added warning for pressure in geotherm
Pressures outside of PREM should not be evaluated for temperatures
>---------------------------------------------------------------
b344b96f4cbe53369dc7c7975f280466e4f1b392
burnman/seismic.py | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/burnman/seismic.py b/burnman/seismic.py
index 552ac27..97e03ac 100644
--- a/burnman/seismic.py
+++ b/burnman/seismic.py
@@ -228,14 +228,15 @@ class SeismicRadiusTable(Seismic1DModel):
return self._lookup(depth, self.table_density)
def depth(self, pressure):
- radius = np.interp(pressure, self.table_pressure[::-1], self.table_radius[::-1] )
+ if pressure > max(self.table_pressure):
+ raise ValueError, "Pressure outside range of PREM"
+ radius = np.interp(pressure, self.table_pressure[::-1], self.table_radius[::-1] )
return self.earth_radius - radius
def _lookup(self, depth, value_table):
radius = self.earth_radius - depth
return np.interp(radius, self.table_radius, value_table)
-
class PREM(SeismicRadiusTable):
"""
Reads PREM (1s) (input_seismic/prem_table.txt, Dziewonski & Anderson 1981).
More information about the CIG-COMMITS
mailing list