hvasbath / beat

Bayesian Earthquake Analysis Tool
GNU General Public License v3.0
132 stars 42 forks source link

pyrocko.gf.meta.OutOfBounds: out of bounds #35

Closed braunfuss closed 5 years ago

braunfuss commented 5 years ago

sergintav1 opened the following issue: Hi I use local model to produce the green function and the codes gives the following error:

beat check bodrum --what=traces config - INFO All parameter-priors ok! config - INFO All hyper-parameters ok! config - INFO All hierarchical-parameters ok! models - INFO ... Initialising Geometry Optimizer ...

models - INFO Analysing problem ... models - INFO ---------------------

heart - INFO Consistent number of datasets and targets in any_P_0 wavemap! heart - INFO Number of seismic datasets for wavemap: any_P_0: 22 seismic - INFO The waveform defined in "slowest 1" config is not included in the optimization! seismic - INFO The waveform defined in "any_S 2" config is not included in the optimization! Traceback (most recent call last): File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/gf/meta.py", line 482, in evaluate times = [phase(args) for phase in phases] File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/gf/meta.py", line 482, in times = [phase(args) for phase in phases] File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/spit.py", line 354, in call return self.interpolate(x) File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/spit.py", line 349, in interpolate raise OutOfBounds() pyrocko.spit.OutOfBounds

During handling of the above exception, another exception occurred:

Traceback (most recent call last): File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/bin/beat", line 11, in load_entry_point('beat==1.0rc1', 'console_scripts', 'beat')() File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x8664.egg/beat/apps/beat.py", line 1600, in main globals()['command' + command](args) File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/apps/beat.py", line 1382, in command_check outmode='stacked_traces') File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/heart.py", line 2196, in prepare_data target=target, wavename=self.name) File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/heart.py", line 1974, in get_phase_arrival_time atime = store.t(wavename, (source.depth, dist)) + source.time File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/gf/store.py", line 1623, in t return timing.evaluate(self.get_phase, args) File "/Users/sergintav/anaconda3/envs/GFZ_py3.7/lib/python3.6/site-packages/pyrocko/gf/meta.py", line 496, in evaluate raise OutOfBounds(args) pyrocko.gf.meta.OutOfBounds: out of bounds: (4300,163488)

Also, I attached my config file:

config_geometry.yaml.txt

How can I solve this problem?

braunfuss commented 5 years ago

Hi sergintav1, it seems like a very simple issue. You calculated the the greens function in the limits: source_depth_min: 5.0 source_depth_max: 20.0 meaning that you calculated it only for source points between 5 and 20km. However as prior bounds for depth you give a source depth between 0 and 5 km: depth: !beat.heart.Parameter name: depth form: Uniform lower: [0.0] upper: [5.0] testvalue: [2.5]

So there is simply no greens functions available at the depths you want so sample.

sergintav1 commented 5 years ago

Thanks a lot for the fast answer.

hvasbath commented 5 years ago

Closing here, as issue seems resolved ...