def dummy_homogeneous_store(self):
mod = cake.LayeredModel.from_scanlines(
cake.read_nd_model_str('''
0 6 3.46 3.0 1000 500
20 6 3.46 3.0 1000 500
'''.lstrip()))
if self._dummy_homogeneous_store is None:
conf = gf.ConfigTypeA(id='empty_homogeneous',
source_depth_min=0.,
source_depth_max=20 * km,
source_depth_delta=10 * km,
distance_min=1000 * km,
distance_max=2000 * km,
distance_delta=10 * km,
sample_rate=2.0,
ncomponents=10,
earthmodel_1d=mod)
store_dir = mkdtemp(prefix='gfstore')
self.tempdirs.append(store_dir)
gf.Store.create(store_dir, config=conf)
self._dummy_homogeneous_store = gf.Store(store_dir)
return self._dummy_homogeneous_store
def _create_regional_ttt_store(self):
conf = gf.ConfigTypeA(
id='empty_regional',
source_depth_min=0.,
source_depth_max=20 * km,
source_depth_delta=10 * km,
distance_min=1000 * km,
distance_max=2000 * km,
distance_delta=10 * km,
sample_rate=2.0,
ncomponents=10,
earthmodel_1d=cake.load_model(),
tabulated_phases=[
gf.TPDef(id=id, definition=defi)
for (id, defi) in [('depthp',
'p'), ('pS', 'pS'), ('P', 'P'), ('S', 'S')]
])
store_dir = mkdtemp(prefix='gfstore')
self.tempdirs.append(store_dir)
gf.Store.create(store_dir, config=conf)
store = gf.Store(store_dir)
store.make_ttt()
store.close()
return store_dir
def _create_benchmark_store(self):
conf = gf.ConfigTypeA(id='benchmark_store',
source_depth_min=0.,
source_depth_max=2.,
source_depth_delta=1.,
distance_min=1.0,
distance_max=5001.0,
distance_delta=5.0,
sample_rate=2.0,
ncomponents=5)
deltat = 1.0 / conf.sample_rate
store_dir = mkdtemp(prefix='gfstore')
self.tempdirs.append(store_dir)
gf.Store.create(store_dir, config=conf)
store = gf.Store(store_dir, 'w')
for args in conf.iter_nodes():
nsamples = int(round(args[1]))
data = num.ones(nsamples)
itmin = int(round(args[1]))
tr = gf.GFTrace(data=data, itmin=itmin, deltat=deltat)
store.put(args, tr)
store.close()
return store_dir
def dummy_store(self):
if self._dummy_store is None:
conf = gf.ConfigTypeA(id='empty_regional',
source_depth_min=0.,
source_depth_max=20 * km,
source_depth_delta=1 * km,
distance_min=1 * km,
distance_max=2000 * km,
distance_delta=1 * km,
sample_rate=2.0,
ncomponents=10)
store_dir = mkdtemp(prefix='gfstore')
self.tempdirs.append(store_dir)
gf.Store.create(store_dir, config=conf)
self._dummy_store = gf.Store(store_dir)
return self._dummy_store
def dummy_store(self):
if self._dummy_store is None:
conf = gf.ConfigTypeA(
id='empty_regional',
source_depth_min=0.,
source_depth_max=20 * km,
source_depth_delta=10 * km,
distance_min=1000 * km,
distance_max=2000 * km,
distance_delta=10 * km,
sample_rate=2.0,
ncomponents=10,
earthmodel_1d=cake.load_model(crust2_profile=(50., 10.)))
store_dir = mkdtemp(prefix='gfstore')
self.tempdirs.append(store_dir)
gf.Store.create(store_dir, config=conf)
self._dummy_store = gf.Store(store_dir)
return self._dummy_store
