[cig-commits] [commit] inversion, master, validate_MT_params: Added warning for pressure in geotherm (b344b96)

[cig_noreply at geodynamics.org]

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

On branches: inversion,master,validate_MT_params
Link       : https://github.com/geodynamics/burnman/compare/80c2a295c42dfdb38f83f6c1334bf7d8f97a8463...409647ff05dfad6a686198cac1481bd46b5e2e62

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

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).