def getRepresentedSource(sandbox):
length = 5000.0
return OkadaSource(length=length,
width=15.0 * length**0.66,
strike=45.0,
rake=0,
slip=2)
def getRepresentedSource(sandbox):
length = 5000.
src = OkadaSource(
easting=num.mean(sandbox.frame.E),
northing=num.mean(sandbox.frame.N),
depth=4000,
length=length,
width=15. * length**.66,
strike=45.,
rake=0,
slip=2,
)
return src
def testOkadaSource(self):
nsources = 2
def r(lo, hi):
return num.random.randint(lo, high=hi, size=1).astype(num.float)
for s in xrange(nsources):
length = r(5000, 15000)
self.ms.addSource(
OkadaSource(
easting=r(0., self.ms.frame.E.max()), # ok
northing=r(0., self.ms.frame.N.max()), # ok
depth=r(0, 8000), # ok
strike=r(0, 360), # ok
dip=r(0, 170),
slip=r(1, 5), # ok
rake=r(0, 180),
length=length,
width=15. * length**.66,
))
self.plotDisplacement(self.ms)
from kite.sources import OkadaSource
km = 1e3
sandbox = SandboxScene()
nframes = 20
depths = num.linspace(4 * km, 4 * km, nframes)
strikes = num.linspace(0, 180, nframes)
okada = OkadaSource(
northing=40 * km,
easting=40 * km,
depth=4 * km,
length=8 * km,
width=4 * km,
strike=63.0,
rake=0,
dip=0.0,
slip=3.0,
opening=0,
)
sandbox.addSource(okada)
sandbox.processSources()
fig, axis = plt.subplots(nrows=2, ncols=2)
# Plot the resulting surface displacements
import matplotlib.pyplot as plt
import numpy as num
km = 1e3
sandbox = SandboxScene()
# Set the LOS incident angles, remember :class:`kite.Scene`
sandbox.phi.fill(num.rad2deg(100.))
sandbox.theta.fill(num.rad2deg(23.))
okada = OkadaSource(
northing=40*km,
easting=40*km,
depth=4*km,
length=8*km,
width=4*km,
strike=63.,
dip=33.,
slip=3.,
opening=1)
pyrocko_rectangular = PyrockoRectangularSource(
northing=40*km,
easting=40*km,
depth=4*km,
length=8*km,
width=4*km,
strike=63.,
dip=33.,
slip=3.,
store_dir='gfstore_halfspace')