Skip to content

Recognize TypeB #105

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Jiefangsong wants to merge 4 commits into
base: master
Choose a base branch
from
Open

Recognize TypeB #105

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
a0448f1
Update AUTHORS.md
Jiefangsong Jun 7, 2022
3e6258a
Update heart.py
Jiefangsong Jun 7, 2022
bee96ec
Update AUTHORS.md
Jiefangsong Jun 7, 2022
254617f
Update config.py
Jiefangsong Jun 7, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 9 additions & 0 deletions beat/config.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -308,6 +308,15 @@ class NonlinearGFConfig(GFConfig):
optional=True,
help='Custom Earthmodel, in case crust2 and standard model not'
' wanted. Needs to be a :py::class:cake.LayeredModel')
receiver_depth_min = Float.T(
default=0.,
help='Minimum receiver depth [km] for GF function grid.')
receiver_depth_max = Float.T(
default=0.,
help='Maximum receiver depth [km] for GF function grid.')
receiver_depth_delta = Float.T(
default=1.,
help='Depth receiver spacing [km] for GF function grid.')
source_depth_min = Float.T(
default=0.,
help='Minimum depth [km] for GF function grid.')
Expand Down
9 changes: 6 additions & 3 deletions beat/heart.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1299,6 +1299,7 @@ def init_seismic_targets(
'%i' % crust_ind, channel), # n, s, l, c
lat=station.lat,
lon=station.lon,
depth=station.depth,
azimuth=cha.azimuth,
dip=cha.dip,
interpolation=interpolation,
Expand Down Expand Up @@ -1735,15 +1736,17 @@ def get_fomosto_baseconfig(
id='slowest',
definition='0.8'))

return gf.ConfigTypeA(
return gf.ConfigTypeB(
id='%s_%s_%.3fHz_%s' % (
station.station,
get_earth_model_prefix(sf.earth_model_name),
sf.sample_rate,
crust_ind),
ncomponents=10,
sample_rate=sf.sample_rate,
receiver_depth=0. * km,
receiver_depth_min=sf.receiver_depth_min * km,
receiver_depth_max=sf.receiver_depth_max * km,
receiver_depth_delta=sf.receiver_depth_delta * km,
source_depth_min=sf.source_depth_min * km,
source_depth_max=sf.source_depth_max * km,
source_depth_delta=sf.source_depth_spacing * km,
Expand Down Expand Up @@ -2244,7 +2247,7 @@ def get_phase_arrival_time(engine, source, target, wavename=None, snap=True):
wavename, cake.m2d * dist))

try:
atime = store.t(wavename, (source.depth, dist)) + source.time
atime = store.t(wavename, (source.depth, target.depth, dist)) + source.time
except TypeError:
raise RayPathError(
'No wave-arrival for wavename "%s" distance %f [deg]! '
Expand Down