[cig-commits] r8367 - cs/spatialdata-0.1/trunk/libsrc/spatialdb

[brad at geodynamics.org]

Author: brad
Date: 2007-12-03 13:20:26 -0800 (Mon, 03 Dec 2007)
New Revision: 8367

Modified:
   cs/spatialdata-0.1/trunk/libsrc/spatialdb/SCECCVMH.cc
Log:
Fixed bug in SCECCVMH where we didn't check to see if query point was deeper than extent of domain (want to use background model in such a case).

Modified: cs/spatialdata-0.1/trunk/libsrc/spatialdb/SCECCVMH.cc
===================================================================
--- cs/spatialdata-0.1/trunk/libsrc/spatialdb/SCECCVMH.cc	2007-12-03 20:40:28 UTC (rev 8366)
+++ cs/spatialdata-0.1/trunk/libsrc/spatialdb/SCECCVMH.cc	2007-12-03 21:20:26 UTC (rev 8367)
@@ -380,7 +380,10 @@
   assert(0 != _topoElev);
   double elev = 0.0;
   const int outsideVoxet = _topoElev->query(&elev, _xyzUTM);
-  if (outsideVoxet) {
+  if (outsideVoxet || z < -10000.0) {
+    // if outside horizontal extent of domain or deeper; assume domain
+    // is at leat 10km deep (z < -10000.0)
+
     if (z < -35000.0)
       vp = 7800.0;
     else if (z < -15000.0)