Changeset 8943

Timestamp:

Jul 3, 2013, 11:56:10 PM (12 years ago)

Author:

davies

Message:

Added data check in Okada tsunami

File:

• trunk/anuga_core/source/anuga/shallow_water/okada_tsunami.py (modified) (1 diff)

trunk/anuga_core/source/anuga/shallow_water/okada_tsunami.py

@@ -209 +209 @@
         raise Exception, 'ERROR: Dip should be between 0 and 90 degrees'
 
+    depth=source[:,2]
+    width=source[:,6]
+    topdepth=depth-width/2.0*numpy.sin(dip/180.*numpy.pi)
+    if(topdepth.min() < 0.0):
+        print '  easting   northing   depth strike    dip length  width  disl1   disl2  disl3 topdepth'
+        for i in range(len(source[:,0])):
+            if(topdepth[i]<0.0):
+                print source[i,:], topdepth[i]
+
+        raise Exception, 'ERROR: The tops of some sub-faults are above the ground (printed above)'
+
     #A few print statements
     if verbose is True: