anuga 1.3.5__zip

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/smf.py

@@ -0,0 +1,467 @@
+#
+# slide_tsunami function
+#
+
+"""This function returns a callable object representing an initial water
+   displacement generated by a submarine sediment slide.
+
+Using input parameters:
+
+Required
+ length  downslope slide length
+ depth   water depth to slide centre of mass
+ slope   bathymetric slope
+
+Optional
+ x0      x origin (0)
+ y0      y origin (0)
+ alpha   angular orientation of slide in xy plane (0)
+ w       slide width (0.25*length)
+ T       slide thickness (0.01*length)
+ g       acceleration due to gravity (9.8)
+ gamma   specific density of sediments (1.85)
+ Cm      added mass coefficient (1)
+ Cd      drag coefficient (1)
+ Cn      friction coefficient (0)
+ psi     (0)
+ dx      offset of second Gaussian (0.2*width of first Gaussian)
+ kappa   multiplier for sech^2 function (3.0)
+ kappad  multiplier for second Gaussian function (0.8)
+ zsmall  an amount near to zero (0.01)
+
+The following parameters are calculated:
+
+ a0      initial acceleration
+ ut      theoretical terminal velocity
+ s0      charactistic distance of motion
+ t0      characteristic time of motion
+ w       initial wavelength of tsunami
+ a2D     2D initial amplitude of tsunami
+ a3D     3D initial amplitude of tsunami
+
+The returned object is a callable double Gaussian function that represents
+the initial water displacement generated by a submarine sediment slide.
+
+Adrian Hitchman
+Geoscience Australia, June 2005
+"""
+
+import numpy as num
+
+import anuga.utilities.log as lg
+
+
+def find_min(x0, wa, kappad, dx):
+    """Determine eta_min to scale eta(x,y)
+
+       eta(x,y) = n03d/nmin*[-f(y)]*g(x)
+
+       nmin = min (-f(y)*g(x) )
+            = -f(ystar)*g(xstar)
+
+       ystar = min (-f(y) ), i.e. diff(-f(y))=0
+       xstar = min ( g(x) ), i.e. diff(g(x))=0
+
+       ystar = y0 and -f(ystar)=1.0
+    """
+    from math import exp, cosh
+
+    step = 0.05
+    x = x0+50.
+    deriv = 10.0
+    count_max = 1000000
+    c = 0
+    deriv = 10.
+    f_ystar = 1.
+   
+    while c < count_max and deriv > 0:
+        deriv = (x-x0)*exp(-((x-x0)/wa)**2.0) - \
+                kappad*(x-dx-x0)*exp(-((x-dx-x0)/wa)**2.0)
+        
+        if deriv <= 0: xstar = x
+        x -= step
+        c += 1
+    
+    g_xstar = exp(-((xstar-x0)/wa)**2)-kappad*exp(-((xstar-dx-x0)/wa)**2)
+
+    etastar = g_xstar*f_ystar
+
+    return etastar
+
+def slide_tsunami(length, depth, slope, width=None, thickness=None, \
+                  x0=0.0, y0=0.0, alpha=0.0, \
+                  gravity=9.8, gamma=1.85, \
+                  massco=1, dragco=1, frictionco=0, psi=0, \
+                  dx=None, kappa=3.0, kappad=0.8, zsmall=0.01, \
+                  scale=None,
+                  domain=None, verbose=False):
+
+    from math import sin, tan, radians, pi, sqrt, exp
+    
+    if domain is not None:
+        xllcorner = domain.geo_reference.get_xllcorner()
+        yllcorner = domain.geo_reference.get_yllcorner()
+        x0 = x0 - xllcorner  # slump origin (relative)
+        y0 = y0 - yllcorner
+        
+    #if width not provided, set to typical value
+    if width is None:
+        width = 0.25 * length
+
+    #if thickness not provided, set to typical value
+    if thickness is None:
+        thickness = 0.01 * length
+
+    #calculate some parameters of the slide
+
+    sint = sin(radians(slope))
+    tant = tan(radians(slope))
+    tanp = tan(radians(psi))
+
+    a0 = gravity * sint * ((gamma-1)/(gamma+massco)) * (1-(tanp/tant))
+    ut = sqrt((gravity*depth) * (length*sint/depth) \
+                    * (pi*(gamma-1)/(2*dragco)) * (1-(tanp/tant)))
+    s0 = ut**2 / a0
+    t0 = ut / a0
+
+    #calculate some parameters of the water displacement produced by the slide
+
+    w = t0 * sqrt(gravity*depth)
+    a2D = s0 * (0.0574 - (0.0431*sint)) \
+             * (thickness/length) \
+             * ((length*sint/depth)**1.25) \
+             * (1 - exp(-2.2*(gamma-1)))
+    a3D = a2D / (1 + (15.5*sqrt(depth/(length*sint))))
+
+    from math import sqrt, log, e
+    dx = 2.0 * (w * sqrt(-log((zsmall/a3D),e))) / 5.0
+        
+    # determine nmin for scaling of eta(x,y)
+    nmin = find_min(x0,w,kappad,dx)  
+    
+    if scale is None:
+        scale = a3D/nmin
+        
+    #a few temporary print statements
+    if verbose is True:
+        lg.critical('\nThe slide ...')
+        lg.critical('\tLength: %s' % str(length))
+        lg.critical('\tDepth: %s' % str(depth))
+        lg.critical('\tSlope: %s' % str(slope))
+        lg.critical('\tWidth: %s' % str(width))
+        lg.critical('\tThickness: %s' % str(thickness))
+        lg.critical('\tx0: %s' % str(x0))
+        lg.critical('\ty0: %s' % str(y0))
+        lg.critical('\tAlpha: %s' % str(alpha))
+        lg.critical('\tAcceleration: %s' % str(a0))
+        lg.critical('\tTerminal velocity: %s' % str(ut))
+        lg.critical('\tChar time: %s' % str(t0))
+        lg.critical('\tChar distance: %s' % str(s0))
+        lg.critical('\nThe tsunami ...')
+        lg.critical('\tWavelength: %s' % str(w))
+        lg.critical('\t2D amplitude: %s' % str(a2D))
+        lg.critical('\t3D amplitude: %s' % str(a3D))
+        lg.critical('\tscale for eta(x,y): %s' % str(scale))
+
+    #keep an eye on some of the assumptions built into the maths
+
+    if ((slope < 5) or (slope > 30)):
+        if verbose is True:
+            lg.critical('WARNING: slope out of range (5 - 30 degrees) %s'
+                         % str(slope))
+    if ((depth/length < 0.06) or (depth/length > 1.5)):
+        if verbose is True:
+            lg.critical('WARNING: d/b out of range (0.06 - 1.5) %s'
+                         % str(depth/length))
+    if ((thickness/length < 0.008) or (thickness/length > 0.2)):
+        if verbose is True:
+            lg.critical('WARNING: T/b out of range (0.008 - 0.2) %s'
+                         % str(thickness/length))
+    if ((gamma < 1.46) or (gamma > 2.93)):
+        if verbose is True:
+            lg.critical('WARNING: gamma out of range (1.46 - 2.93) %s'
+                         % str(gamma))
+
+    return Double_gaussian(a3D, w, width, x0, y0, alpha, kappa, kappad, zsmall, dx, scale)
+
+#
+# slump_tsunami function
+#
+
+"""This function returns a callable object representing an initial water
+   displacement generated by a submarine sediment slump.
+
+Using input parameters:
+
+Required
+ length  downslope slump length
+ depth   water depth to slump centre of mass
+ slope   bathymetric slope
+
+Optional
+ x0      x origin (0)
+ y0      y origin (0)
+ alpha   angular orientation of slide in xy plane (0)
+ w       slump width (1.0*length)
+ T       slump thickness (0.1*length)
+ R       slump radius of curvature (b^2/(8*T))
+ del_phi slump angular displacement (0.48)
+ g       acceleration due to gravity (9.8)
+ gamma   specific density of sediments (1.85)
+ Cm      added mass coefficient (1)
+ Cd      drag coefficient (1)
+ Cn      friction coefficient (0)
+ dx      offset of second Gaussian (0.2*width of first Gaussian)
+ kappa   multiplier for sech^2 function (3.0)
+ kappad  multiplier for second Gaussian function (0.8)
+ zsmall  an amount near to zero (0.01)
+
+The following parameters are calculated:
+
+ a0      initial acceleration
+ um      maximum velocity
+ s0      charactistic distance of motion
+ t0      characteristic time of motion
+ w       initial wavelength of tsunami
+ a2D     2D initial amplitude of tsunami
+ a3D     3D initial amplitude of tsunami
+
+The returned object is a callable double Gaussian function that represents
+the initial water displacement generated by a submarine sediment slump.
+
+Adrian Hitchman
+Geoscience Australia, June 2005
+"""
+
+def slump_tsunami(length, depth, slope, width=None, thickness=None, \
+                  radius=None, dphi=0.48, x0=0.0, y0=0.0, alpha=0.0, \
+                  gravity=9.8, gamma=1.85, \
+                  massco=1, dragco=1, frictionco=0, \
+                  dx=None, kappa=3.0, kappad=1.0, zsmall=0.01, scale=None, \
+                  domain=None,
+                  verbose=False):
+
+    from math import sin, radians, sqrt
+
+    if domain is not None:
+        xllcorner = domain.geo_reference.get_xllcorner()
+        yllcorner = domain.geo_reference.get_yllcorner()
+        x0 = x0 - xllcorner  # slump origin (relative)
+        y0 = y0 - yllcorner
+
+    #if width not provided, set to typical value
+    if width is None:
+        width = length
+
+    #if thickness not provided, set to typical value
+    if thickness is None:
+        thickness = 0.1 * length
+
+    #if radius not provided, set to typical value
+    if radius is None:
+        radius = length**2 / (8.0 * thickness)
+
+    #calculate some parameters of the slump
+
+    sint = sin(radians(slope))
+
+    s0 = radius * dphi / 2
+    t0 = sqrt((radius*(gamma+massco)) / (gravity*(gamma-1)))
+    a0 = s0 / t0**2
+    um = s0 / t0
+
+    #calculate some parameters of the water displacement produced by the slump
+
+    w = t0 * sqrt(gravity*depth)
+    a2D = s0 * (0.131/sint) \
+             * (thickness/length) \
+             * (length*sint/depth)**1.25 \
+             * (length/radius)**0.63 * dphi**0.39 \
+             * (1.47 - (0.35*(gamma-1))) * (gamma-1)
+    a3D = a2D / (1 + (2.06*sqrt(depth/length)))
+
+    from math import sqrt, log, e
+    dx = 2.0 * (w * sqrt(-log((zsmall/a3D),e))) / 5.0
+        
+    # determine nmin for scaling of eta(x,y)
+    nmin = find_min(x0,w,kappad,dx)  
+    
+    if scale is None:
+        scale = a3D/nmin
+        
+    #a few temporary print statements
+    if verbose is True:
+        lg.critical('\nThe slump ...')
+        lg.critical('\tLength: %s' % str(length))
+        lg.critical('\tDepth: %s' % str(depth))
+        lg.critical('\tSlope: %s' % str(slope))
+        lg.critical('\tWidth: %s' % str(width))
+        lg.critical('\tThickness: %s' % str(thickness))
+        lg.critical('\tRadius: %s' % str(radius))
+        lg.critical('\tDphi: %s' % str(dphi))
+        lg.critical('\tx0: %s' % str(x0))
+        lg.critical('\ty0: %s' % str(y0))
+        lg.critical('\tAlpha: %s' % str(alpha))
+        lg.critical('\tAcceleration: %s' % str(a0))
+        lg.critical('\tMaximum velocity: %s' % str(um))
+        lg.critical('\tChar time: %s' % str(t0))
+        lg.critical('\tChar distance: %s' % str(s0))
+        lg.critical('\nThe tsunami ...')
+        lg.critical('\tWavelength: %s' % str(w))
+        lg.critical('\t2D amplitude: %s' % str(a2D))
+        lg.critical('\t3D amplitude: %s' % str(a3D))
+        lg.critical('\tDelta x %s' % str(dx))
+        lg.critical('\tsmall %s' % str(zsmall))
+        lg.critical('\tKappa d  %s' % str(kappad))
+        lg.critical('\tscale for eta(x,y): %s' % str(scale))
+
+    #keep an eye on some of the assumptions built into the maths
+
+    if ((slope < 10) or (slope > 30)):        
+        if verbose is True:
+            lg.critical('WARNING: slope out of range (10 - 30 degrees) %s'
+                         % str(slope))
+    if ((depth/length < 0.34) or (depth/length > 0.5)):     
+        if verbose is True:
+            lg.critical('WARNING: d/b out of range (0.34 - 0.5) %s'
+                         % str(depth/length))
+    if ((thickness/length < 0.10) or (thickness/length > 0.15)):     
+        if verbose is True:
+            lg.critical('WARNING: T/b out of range (0.10 - 0.15) %s'
+                         % str(thickness/length))
+    if ((radius/length < 1.0) or (radius/length > 2.0)):     
+        if verbose is True:
+            lg.critical('WARNING: R/b out of range (1 - 2) %s'
+                         % str(radius/length))
+    if ((dphi < 0.10) or (dphi > 0.52)):     
+        if verbose is True:
+            lg.critical('WARNING: del_phi out of range (0.10 - 0.52) %s'
+                         % str(dphi))
+    if ((gamma < 1.46) or (gamma > 2.93)):     
+        if verbose is True:
+            lg.critical('WARNING: gamma out of range (1.46 - 2.93) %s'
+                         % str(gamma))
+
+    return Double_gaussian(a3D, w, width, x0, y0, alpha, kappa, kappad, zsmall, dx, scale)
+
+#
+# Double_gaussian class
+#
+
+"""This is a callable class representing the initial water displacment 
+   generated by a sediment slide or slump.
+
+Using input parameters:
+
+Required
+ w       initial wavelength of tsunami
+ a3D     3D initial amplitude of tsunami
+ width   width of smf
+
+Optional
+ x0      x origin of smf
+ y0      y origin of smf
+ alpha   angular orientation of smf in xy plane (0)
+ dx      offset of second Gaussian (0.2*width of first Gaussian)
+ kappa   multiplier for sech^2 function (3.0)
+ kappad  multiplier for second Gaussian function (0.8)
+ zsmall  an amount near to zero (0.01)
+
+Adrian Hitchman
+Geoscience Australia, June 2005
+"""
+
+class Double_gaussian:
+
+    def __init__(self, a3D, wavelength, width, x0, y0, alpha, \
+                 kappa, kappad, zsmall, dx, scale):
+        self.a3D = a3D
+        self.wavelength = wavelength
+        self.width = width
+        self.x0 = x0
+        self.y0 = y0
+        self.alpha = alpha
+        self.kappa = kappa
+        self.kappad = kappad
+        self.scale = scale
+
+        if dx is None:
+            from math import sqrt, log, e
+            dx = 2.0 * (self.wavelength * sqrt(-log((zsmall/self.a3D),e))) / 5.0
+        self.dx = dx
+
+    def __call__(self, x, y):
+        """Make Double_gaussian a callable object.
+
+        If called as a function, this object returns z values representing
+        the initial 3D distribution of water heights at the points (x,y)
+        produced by a submarine mass failure.
+        """
+
+        from math import sin, cos, radians, exp, cosh
+
+        #ensure vectors x and y have the same length
+        N = len(x)
+        assert N == len(y)
+
+        am = self.a3D
+        am2 = 1.0
+        wa = self.wavelength
+        wi = self.width
+        x0 = self.x0
+        y0 = self.y0
+        alpha = self.alpha
+        dx = self.dx
+        kappa = self.kappa
+        kappad = self.kappad
+        scale = self.scale
+
+        #double Gaussian calculation assumes water displacement is oriented
+        #E-W, so, for displacement at some angle alpha clockwise from the E-W
+        #direction, rotate (x,y) coordinates anti-clockwise by alpha
+
+        cosa = cos(radians(alpha))
+        sina = sin(radians(alpha))
+
+        xr = ((x-x0) * cosa - (y-y0) * sina) + x0
+        yr = ((x-x0) * sina + (y-y0) * cosa) + y0
+
+        z = num.zeros(N, num.float)
+        maxz = 0.0
+        minz = 0.0
+        for i in range(N):
+            try:
+                z[i] =  -scale / ((cosh(kappa*(yr[i]-y0)/(wi+wa)))**2) \
+                            * (exp(-((xr[i]-x0)/wa)**2) - \
+                                kappad*exp(-((xr[i]-dx-x0)/wa)**2))
+                if z[i] > maxz: maxz = z[i]
+                if z[i] < minz: minz = z[i]
+                
+            except OverflowError:
+                pass
+                
+        return z
+
+    def determineDX(self, zsmall):
+        """Determine a suitable offset for the second Gaussian function.
+
+        A suitable offset for the second Gaussian function is taken to
+        be some fraction of the 'width' of the first Gaussian function.
+
+        The 'width' of the first Gaussian is obtained from the range of
+        the x coordinates over which the function takes values from
+        'near zero', through 1, and back to 'near zero'.
+
+        The parameter zsmall passed to this function specifies how much
+        'near zero' is.
+        """
+
+        from math import sqrt, log, e
+        
+        a = self.a3D
+        c = self.wavelength
+
+        self.dx = 2.0 * (c * sqrt(-log((zsmall/a),e))) / 5.0
+
+        return self.dx
+

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_MS.lic

@@ -0,0 +1,23 @@
+<?xml version='1.0' encoding='iso-8859-1'?>
+<ga_license_file>
+  <metadata>
+    <author>Herve Damlamian</author>
+  </metadata>
+  
+  <datafile>
+    <filename>fullokada_MS.txt</filename>
+    <checksum>1591397974</checksum>
+    <publishable>Yes</publishable>
+    <accountable>Jane Sexton</accountable>
+    <source>Output from fortran code (okada.f) supplied through the publication:
+    Computation of deformation induced by earthquakes in a multi-layered elastic crust - 
+    FORTRAN programs EDGRN/EDCMP by Rongjiang Wang, Francisco Lorenzo Martin and Frank Roth, 
+    in Computers and Geosciences, volume 29, 195-207, 2003. The fortran code is available 
+    from http://www.iamg.org/CGEditor/index.htm. Code run with parameters in 
+    test_tsunami_okada.py in the multiple source mode. Hence file naming of _MS.</source>
+    <IP_owner>SOPAC and Geoscience Australia.</IP_owner>
+    <IP_info>Based on freely available fortran code available at 
+    http://www.iamg.org/CGEditor/index.htm.</IP_info>
+  </datafile>
+  
+  </ga_license_file>

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_MS.txt

@@ -0,0 +1,36 @@
+-.19170511    -0.15687813    -0.25413750
+-.19499435    -0.23341475    -0.38796115
+-.14990029    -0.31190623    -0.53039482
+-.05117426    -0.36155257    -0.62515593
+0.07004436    -0.35600619    -0.62020153
+0.16342513    -0.29764163    -0.51694933
+-.21481727    -0.12919853    -0.29298271
+-.22793125    -0.20245203    -0.47511847
+-.17981488    -0.28082927    -0.67947090
+-.05933581    -0.33021106    -0.81717521
+0.09041813    -0.32123460    -0.80362204
+0.19899882    -0.25897659    -0.64599505
+-.19604864    -0.07868877    -0.27029146
+-.21197020    -0.12799217    -0.45724775
+-.16713973    -0.18179022    -0.67089015
+-.04998927    -0.21443959    -0.81021771
+0.09290397    -0.20494846    -0.78376904
+0.19004325    -0.15959199    -0.60949022
+-.12749360    -0.03840539    -0.16795117
+-.13443714    -0.06340126    -0.28812204
+-.10108254    -0.09149020    -0.42087346
+-.02361551    -0.10892580    -0.49882011
+0.06509034    -0.10412453    -0.46979073
+0.12179823    -0.08046168    -0.35499996
+-.03545205    -0.04119929    -0.02699778
+-.03056141    -0.06610209    -0.04770727
+-.01651872    -0.09493258    -0.06474178
+0.00300779    -0.11484027    -0.06708284
+0.01978183    -0.11349546    -0.05288846
+0.02871849    -0.09198804    -0.03208618
+0.03950435    -0.08214826     0.08389296
+0.04779941    -0.12461724     0.12914981
+0.04176548    -0.17027240     0.18056085
+0.01715442    -0.19991123     0.21583972
+-.01593428    -0.19704504     0.21428797
+-.03986000    -0.16381659     0.17744573

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_SP.lic

@@ -0,0 +1,23 @@
+<?xml version='1.0' encoding='iso-8859-1'?>
+<ga_license_file>
+  <metadata>
+    <author>Herve Damlamian</author>
+  </metadata>
+  
+  <datafile>
+    <filename>fullokada_SP.txt</filename>
+    <checksum>3395079577</checksum>
+    <publishable>Yes</publishable>
+    <accountable>Jane Sexton</accountable>
+    <source>OOutput from fortran code (okada.f) supplied through the publication:
+    Computation of deformation induced by earthquakes in a multi-layered elastic crust - 
+    FORTRAN programs EDGRN/EDCMP by Rongjiang Wang, Francisco Lorenzo Martin and Frank Roth, 
+    in Computers and Geosciences, volume 29, 195-207, 2003. The fortran code is available 
+    from http://www.iamg.org/CGEditor/index.htm. Code run with parameters in 
+    test_tsunami_okada.py in the single point mode. Hence file naming of _SP.</source>
+    <IP_owner>SOPAC and Geoscience Australia.  </IP_owner>
+    <IP_info>Based on freely available fortran code available at 
+    http://www.iamg.org/CGEditor/index.htm.</IP_info>
+  </datafile>
+  
+  </ga_license_file>

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_SP.txt

@@ -0,0 +1,36 @@
+-.00234481    -0.00165609    -0.00391997
+-.00213709    -0.00259559    -0.00628137
+-.00090816    -0.00336007    -0.00824223
+0.00090816    -0.00336007    -0.00824223
+0.00213709    -0.00259559    -0.00628137
+0.00234481    -0.00165609    -0.00391997
+-.00217612    -0.00091269    -0.00361610
+-.00208233    -0.00151934    -0.00617730
+-.00091648    -0.00204983    -0.00846255
+0.00091648    -0.00204983    -0.00846255
+0.00208233    -0.00151934    -0.00617730
+0.00217612    -0.00091269    -0.00361610
+-.00127284    -0.00035749    -0.00195464
+-.00120174    -0.00060274    -0.00339635
+-.00052565    -0.00082338    -0.00469174
+0.00052565    -0.00082338    -0.00469174
+0.00120174    -0.00060274    -0.00339635
+0.00127284    -0.00035749    -0.00195465
+-.00012261    -0.00053631     0.00011798
+-.00005774    -0.00085618     0.00018346
+-.00001092    -0.00113253     0.00026295
+0.00001092    -0.00113253     0.00026295
+0.00005774    -0.00085618     0.00018346
+0.00012261    -0.00053631     0.00011798
+0.00064093    -0.00117261     0.00140100
+0.00061642    -0.00174294     0.00214396
+0.00026790    -0.00219196     0.00275077
+-.00026790    -0.00219196     0.00275077
+-.00061642    -0.00174294     0.00214396
+-.00064093    -0.00117261     0.00140100
+0.00088998    -0.00169357     0.00171918
+0.00076592    -0.00233657     0.00242693
+0.00031156    -0.00279843     0.00294873
+-.00031156    -0.00279843     0.00294873
+-.00076592    -0.00233657     0.00242693
+-.00088998    -0.00169357     0.00171918

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_SS.lic

@@ -0,0 +1,23 @@
+<?xml version='1.0' encoding='iso-8859-1'?>
+<ga_license_file>
+  <metadata>
+    <author>Herve Damlamian</author>
+  </metadata>
+  
+  <datafile>
+    <filename>fullokada_SS.txt</filename>
+    <checksum>3315294267</checksum>
+    <publishable>Yes</publishable>
+    <accountable>Jane Sexton</accountable>
+    <source>Output from fortran code (okada.f) supplied through the publication:
+    Computation of deformation induced by earthquakes in a multi-layered elastic crust - 
+    FORTRAN programs EDGRN/EDCMP by Rongjiang Wang, Francisco Lorenzo Martin and Frank Roth, 
+    in Computers and Geosciences, volume 29, 195-207, 2003. The fortran code is available 
+    from http://www.iamg.org/CGEditor/index.htm. Code run with parameters in 
+    test_tsunami_okada.py in the single source mode. Hence file naming of _SS.</source>
+    <IP_owner>SOPAC and Geoscience Australia.  </IP_owner>
+    <IP_info>Based on freely available fortran code available at 
+    http://www.iamg.org/CGEditor/index.htm.</IP_info>
+  </datafile>
+  
+  </ga_license_file>

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/tsunami_source/tests/data/fullokada_SS.txt

@@ -0,0 +1,36 @@
+-.09896401    -0.06266776    -0.11488581
+-.11014517    -0.09974463    -0.18732074
+-.09764119    -0.14380126    -0.27540737
+-.05117426    -0.17945972    -0.34791887
+0.01778525    -0.18790121    -0.36521408
+0.07857595    -0.16397151    -0.31630892
+-.10198684    -0.04392203    -0.11887774
+-.11901769    -0.07401557    -0.20685036
+-.11013730    -0.11228523    -0.32191411
+-.05933581    -0.14482059    -0.42171589
+0.02074054    -0.15269056    -0.44606525
+0.09008526    -0.13054013    -0.37772694
+-.08209710    -0.02150989    -0.09234827
+-.09747238    -0.03748287    -0.16744398
+-.09175030    -0.05869899    -0.26925331
+-.04998927    -0.07729888    -0.35980535
+0.01751454    -0.08185723    -0.38213220
+0.07554544    -0.06908268    -0.31968644
+-.04271989    -0.01189426    -0.03980895
+-.04862055    -0.02033263    -0.07372370
+-.04423248    -0.03174617    -0.11865463
+-.02361551    -0.04188265    -0.15795134
+0.00824028    -0.04438050    -0.16757187
+0.03598163    -0.03739305    -0.14060162
+-.00083686    -0.02226588     0.01429015
+0.00232566    -0.03596514     0.02151249
+0.00435985    -0.05347275     0.03187149
+0.00300778    -0.06841496     0.04135995
+-.00109674    -0.07203563     0.04372482
+-.00416858    -0.06185109     0.03713359
+0.02803219    -0.04347391     0.04853847
+0.03399120    -0.06635202     0.07537680
+0.03182526    -0.09328642     0.10832243
+0.01715442    -0.11494618     0.13554835
+-.00599406    -0.12005907     0.14204955
+-.02605180    -0.10555138     0.12367272