unclear GF parameters, GF database

Hi, when I running "beat build_gfs Laquila_kinematic --datatypes=geodetic --execute --mode=ffo" in Scenario III, I meet an error:
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
beat - INFO Discretized fault geometry exists! Use --force to overwrite!
beat - INFO Loading existing discretized fault
beat - INFO Calculating linear Green's Functions
beat - INFO for seismic data ...
beat - INFO crust_ind 0
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
ffo - INFO Calculating GFs for starttimes: "-10.0", "-9.0", "-8.0", "-7.0", "-6.0", "-5.0", "-4.0", "-3.0", "-2.0", "-1.0", "0.0", "1.0", "2.0", "3.0", "4.0", "5.0", "6.0", "7.0", "8.0", "9.0", "10.0", "11.0", "12.0", "13.0", "14.0", "15.0", "16.0", "17.0", "18.0", "19.0", "20.0", "21.0", "22.0", "23.0", "24.0", "25.0", "26.0", "27.0", "28.0", "29.0", "30.0", "31.0", "32.0"
durations: "0.0", "1.0", "2.0", "3.0", "4.0"
ffo - INFO Using 1 workers ...
ffo - INFO For slip component: uparr
ffo - INFO Allocating GF Library
ffo - INFO Setting up Green's Function Library:
Seismic GF Library

--- !beat.SeismicGFLibraryConfig
component: uparr
event: !pf.Event
lat: 42.29
lon: 13.35
time: 2009-04-06 01:32:49.190000
name: 200904060132A
depth: 12000.0
magnitude: 6.343080192483292
region: CENTRAL ITALY
catalog: gCMT
moment_tensor: !pf.MomentTensor
mnn: 1.43e+18
mee: 1.87e+18
mdd: -3.3e+18
mne: 1.77e+18
mnd: -1.43e+18
med: 2.6900000000000003e+17
strike1: 120.23408298515041
dip1: 54.240869089580485
rake1: -112.81739742081386
strike2: 335.98575923255856
dip2: 41.58440373860804
rake2: -61.69749587601104
moment: 3.6696131948749036e+18
magnitude: 6.343080192483292
duration: 7.0
datatype: seismic
crust_ind: 0
reference_sources:

!beat.sources.RectangularSource
lat: 42.29
lon: 13.35
north_shift: 6210.765329902994
east_shift: 12215.633404840579
elevation: 0.0
depth: 371.66597868550434
time: 2009-04-06 01:32:49.190000
stf: !pf.HalfSinusoidSTF
duration: 0.0
anchor: -1.0
stf_mode: post
strike: 145.40298276603247
dip: 54.22973164156763
rake: -112.79284665995873
length: 20000.0
width: 15000.0
velocity: 3500.0
slip: 1.0
wave_config: !beat.WaveformFitConfig
include: true
preprocess_data: true
name: any_P
channels: [Z]
filterer: !beat.heart.Filter
lower_corner: 0.001
upper_corner: 0.1
order: 4
distances: [30.0, 90.0]
interpolation: multilinear
arrival_taper: !beat.heart.ArrivalTaper
a: -15.0
b: -10.0
c: 50.0
d: 55.0
starttime_sampling: 1.0
duration_sampling: 1.0
starttime_min: -10.0
duration_min: 0.0
dimensions: [35, 12, 5, 43, 120]

ntargets: 35
npatches: 12
ndurations: 5
nstarttimes: 43
nsamples: 120
size: 10836000
filesize [MB]: 82.672119
filename: seismic_uparr_any_P_0

ffo - INFO Patch Number 0
Traceback (most recent call last):
File "/beat/python35/bin/beat", line 11, in
load_entry_point('beat==1.0rc1', 'console_scripts', 'beat')()
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 1588, in main
globals()'command_' + command
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 1144, in command_build_gfs
force=options.force)
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/ffo/base.py", line 1088, in seis_construct_gf_linear
for res in p:
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/parallel.py", line 224, in paripool
yield [function(*work)]
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/ffo/base.py", line 927, in _process_patch_seismic
outmode='data')
File "/.local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/heart.py", line 2574, in seis_synthetics
targets=targets, nprocs=nprocs)
File "/beat/python35/lib/python3.5/site-packages/pyrocko/gf/seismosizer.py", line 3114, in process
self.get_store(store_id)
File "/beat/python35/lib/python3.5/site-packages/pyrocko/gf/seismosizer.py", line 2898, in get_store
store_dir = self.get_store_dir(store_id)
File "/beat/python35/lib/python3.5/site-packages/pyrocko/gf/seismosizer.py", line 2856, in get_store_dir
raise NoSuchStore(store_id, self.iter_store_dirs())
pyrocko.gf.seismosizer.NoSuchStore: no GF store with id "TSUM_ak135_2.000Hz_0" found. Searched folders:

I noticed that in "Elementray GFs" of Scenario III, the green function is calculated with "beat build_gfs Laquila --datatypes='seismic' --execute". Should it use "Laquila_kinematic" in this command? Why is the "Laquila" is used here? And I calculated the GFs with 0.5HZ and the error said there is no 2HZ GFs. Can you give me some suggestions?

BAD_REQUEST, beat sample --hypers

Hello,

I have successfully installed BEAT on another machine, and Examples 3 and 4 turned out beautifully. I am now trying to perform a RectangularSource inversion for my own data.

I have processed two InSAR scenes (one ascending, one descending) for a medium-sized earthquake in central Tibet. Because of the remote location of the earthquake, I am trying to do the inversion using only geodetic (InSAR) data.

I prepared the scenes using Kite and its quadtree approach. I initialized the project with:

beat init Xizang 2020-07-22 --datatypes=geodetic --n_sources=1 --source_type=RectangularSource --min_mag=6.2

This gave me a large EQ in Alaska, and the central Tibet event that I am interested in. In the config_geometry.yaml file I replaced the pre-loaded Alaska event with the Xizang, Tibet info. I then set the data directory and Green's Function directory and did:

beat update Xizang

It gave me the message that the covariances had not been calculated, then calculated them. Next I gave:

beat build_gfs Xizang --datatypes=geodetic
beat build_gfs Xizang --force --execute

It successfully build the Green's Functions, and displayed: "beat - INFO geodetic GF calculations successful!"
I did:

beat update Xizang --parameters='hierarchicals,hypers'

It gave the OK for the parameters (1 hyperparameter, 0 hierarchicals). So next I did:

beat sample Xizang --hypers

and it gave the following error:

(pyrocko) rzinke@mpzinke-ubuntu:~/Titan/Projects/Xizang/BeatTest1$ beat sample Xizang --hypers
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
models       - INFO     ... Initialising Geometry Optimizer ...

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

geodetic     - INFO     Number of geodetic datasets: 2
geodetic     - INFO     Number of geodetic data points: 233
geodetic     - WARNING  Covariance estimation not implemented (yet)! Using imported covariances!
models       - INFO     ... Building Hyper model ...

geodetic     - INFO     Initialized 0 hierarchical parameters (ramps).
models       - INFO     Optimization for 1 hyperparameters in total!
Traceback (most recent call last):
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/bin/beat", line 33, in <module>
    sys.exit(load_entry_point('beat==1.0', 'console_scripts', 'beat')())
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
    globals()['command_' + command](args)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 768, in command_sample
    project_dir, options.mode, options.hypers)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/problems.py", line 842, in load_model
    problem.built_hyper_model()
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/problems.py", line 280, in built_hyper_model
    composite.update_llks(point)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/geodetic.py", line 322, in update_llks
    results = self.assemble_results(point)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/geodetic.py", line 190, in assemble_results
    processed_synts = self.get_synthetics(point, outmode='stacked_arrays')
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/geodetic.py", line 479, in get_synthetics
    **kwargs)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/heart.py", line 2873, in geo_synthetics
    response = engine.process(sources, targets)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/pyrocko/gf/seismosizer.py", line 3577, in process
    nthreads=nthreads):
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/pyrocko/gf/seismosizer.py", line 3086, in process_static
    source, target, components, nthreads)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/pyrocko/gf/seismosizer.py", line 3423, in base_statics
    nthreads)
  File "/home/rzinke/Miniconda/miniconda3/envs/pyrocko/lib/python3.6/site-packages/pyrocko/gf/store.py", line 1888, in statics
    nthreads)
pyrocko.gf.store_ext.StoreExtError: BAD_REQUEST

I reran the build_gfs commands and checked that the stores are empty. I searched the issues page and did not find an answer to this issue. I'm happy to provide any information if you can help me solve it!

Thanks!

geodetic composite data loading

Needs error checking if the types configs indeed are filled with meaningful values and if defined hyper priors indeed refer to related data...

ffi: allow for number of subfault prior vector length for patch related priors

For now it is allowed, to use either a single value for ALL the patches, or patch individual priors ergo length of the vector as number of patches.

trans-D Lune source

trans-D Tape & Tape lune parametrization

How to obtain the model resolution matrix in beat

Hi again,
Does beat have the function that can show the model resolution matrix? If it does, could you please tell me the method to get it?
Thank you very much!
Best regards,

Error when running beat sample Laquila

Hi,
I meet this error when I running the Scnario 3: Rectangular source.
Before I ran beat sample Laquila, I ran: beat update Laquila --parameters="hierarchicals" , beat update Laquila --parameters="hypers" and beat sample Laquila --hypers. And then I changed n_jobs in sampler_config to 50.
The error message is:
backend - INFO Loading multitrace from /data/lxrtest/insarlab/BEAT/Laquila/geometry/trans_stage_5
smc - INFO Beta: 0.015582 Stage: 6
backend - INFO Loading multitrace from /data/lxrtest/insarlab/BEAT/Laquila/geometry/stage_6
Traceback (most recent call last):
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/bin/beat", line 11, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.0-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 1688, in main
globals()'command_' + command
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.0-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 777, in command_sample
sample(step, problem)
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.0-py3.7-linux-x86_64.egg/beat/models/base.py", line 176, in sample
rm_flag=pa.rm_flag)
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.0-py3.7-linux-x86_64.egg/beat/sampler/smc.py", line 460, in smc_sample
step.select_end_points(mtrace)
File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.0-py3.7-linux-x86_64.egg/beat/sampler/smc.py", line 210, in select_end_points
array_population[:, slc] = slc_population0/1 [00:00<?, ?it/s]
ValueError: could not broadcast input array from shape (902,1) into shape (1000,1)

And my config_geometry.yaml file is:
name: /data/lxrtest/insarlab/BEAT/Laquila
date: '2009-04-06'
event: !pf.Event
lat: 42.29
lon: 13.35
time: 2009-04-06 01:32:49.190000
depth: 12000.0
name: '200904060132A'
magnitude: 6.343080192483292
region: 'CENTRAL ITALY'
catalog: 'gCMT'
moment_tensor: !pf.MomentTensor
mnn: 1.43e+18
mee: 1.87e+18
mdd: -3.3e+18
mne: 1.77e+18
mnd: -1.43e+18
med: 2.6900000000000003e+17
strike1: 120.23408298515041
dip1: 54.240869089580485
rake1: -112.81739742081386
strike2: 335.98575923255856
dip2: 41.58440373860804
rake2: -61.69749587601104
moment: 3.6696131948749036e+18
magnitude: 6.343080192483292
duration: 7.0
project_dir: /data/lxrtest/insarlab/BEAT/Laquila
problem_config: !beat.ProblemConfig
mode: geometry
source_type: RectangularSource
stf_type: HalfSinusoid
decimation_factors:
geodetic: 4
seismic: 1
n_sources: 1
datatypes:

geodetic
hyperparameters:
h_SAR: !beat.heart.Parameter
name: h_SAR
form: Uniform
lower:
- 0.0
  upper:
- 5.0
  testvalue:
- 2.5
  priors:
  depth: !beat.heart.Parameter
  name: depth
  form: Uniform
  lower:
- 0.0
  upper:
- 5.0
  testvalue:
- 2.95
  dip: !beat.heart.Parameter
  name: dip
  form: Uniform
  lower:
- 35.0
  upper:
- 65.0
  testvalue:
- 59.46
  east_shift: !beat.heart.Parameter
  name: east_shift
  form: Uniform
  lower:
- -5.0
  upper:
- 20.0
  testvalue:
- 10.64
  length: !beat.heart.Parameter
  name: length
  form: Uniform
  lower:
- 5.0
  upper:
- 25.0
  testvalue:
- 12.4
  north_shift: !beat.heart.Parameter
  name: north_shift
  form: Uniform
  lower:
- -15.0
  upper:
- 10.0
  testvalue:
- 5.624
  nucleation_x: !beat.heart.Parameter
  name: nucleation_x
  form: Uniform
  lower:
- -1.0
  upper:
- 1.0
  testvalue:
- -0.94
  nucleation_y: !beat.heart.Parameter
  name: nucleation_y
  form: Uniform
  lower:
- -1.0
  upper:
- 1.0
  testvalue:
- -0.09508065
  rake: !beat.heart.Parameter
  name: rake
  form: Uniform
  lower:
- -180.0
  upper:
- 0.0
  testvalue:
- -125.0
  slip: !beat.heart.Parameter
  name: slip
  form: Uniform
  lower:
- 0.1
  upper:
- 1.0
  testvalue:
- 0.634
  strike: !beat.heart.Parameter
  name: strike
  form: Uniform
  lower:
- 60.0
  upper:
- 180.0
  testvalue:
- 143.0
  time: !beat.heart.Parameter
  name: time
  form: Uniform
  lower:
- -20.0
  upper:
- 10.0
  testvalue:
- -10.0
  width: !beat.heart.Parameter
  name: width
  form: Uniform
  lower:
- 5.0
  upper:
- 20.0
  testvalue:
- 10.0
  hierarchicals:
  Laquila_ascxn_offset: !beat.heart.Parameter
  name: Laquila_ascxn_offset
  form: Uniform
  lower:
- -0.05
  upper:
- 0.05
  testvalue:
- 0.0
  Laquila_ascxn_ramp: !beat.heart.Parameter
  name: Laquila_ascxn_ramp
  form: Uniform
  lower:
- -0.005
- -0.005
  upper:
- 0.005
- 0.005
  testvalue:
- 0.0
- 0.0
  Laquila_dscxn_offset: !beat.heart.Parameter
  name: Laquila_dscxn_offset
  form: Uniform
  lower:
- -0.05
  upper:
- 0.05
  testvalue:
- 0.0
  Laquila_dscxn_ramp: !beat.heart.Parameter
  name: Laquila_dscxn_ramp
  form: Uniform
  lower:
- -0.005
- -0.005
  upper:
- 0.005
- 0.005
  testvalue:
- 0.0
- 0.0
  geodetic_config: !beat.GeodeticConfig
  datadir: /home/vasyurhm/BEATS/RectangularStatic/insar
  names:
Laquila_dscxn
Laquila_ascxn
types:
SAR
calc_data_cov: true
interpolation: multilinear
fit_plane: true
dataset_specific_residual_noise_estimation: true
gf_config: !beat.GeodeticGFConfig
store_superdir: /data/lxrtest/insarlab/GF/Laquila
reference_model_idx: 0
n_variations:
- 0
- 1
  earth_model_name: ak135-f-continental.m
  nworkers: 60
  use_crust2: false
  replace_water: false
  source_depth_min: 0.0
  source_depth_max: 35.0
  source_depth_spacing: 1.0
  source_distance_radius: 100.0
  source_distance_spacing: 1.0
  error_depth: 0.1
  error_velocities: 0.1
  depth_limit_variation: 600.0
  code: psgrn
  sample_rate: 1.1574074074074073e-05
  sampling_interval: 1.0
  medium_depth_spacing: 1.0
  medium_distance_spacing: 1.0
  sampler_config: !beat.SamplerConfig
  name: SMC
  backend: csv
  progressbar: true
  buffer_size: 5000
  buffer_thinning: 1
  parameters: !beat.SMCConfig
  tune_interval: 10
  check_bnd: true
  rm_flag: false
  n_jobs: 50
  n_steps: 100
  n_chains: 1000
  coef_variation: 1.0
  stage: 0
  proposal_dist: MultivariateCauchy
  update_covariances: true
  hyper_sampler_config: !beat.SamplerConfig
  name: Metropolis
  backend: csv
  progressbar: true
  buffer_size: 5000
  buffer_thinning: 1
  parameters: !beat.MetropolisConfig
  tune_interval: 50
  proposal_dist: Normal
  check_bnd: true
  rm_flag: false
  n_jobs: 4
  n_steps: 25000
  n_chains: 12
  thin: 5
  burn: 0.5

Could you please give me some suggestions?
Thank you !

error import geometry to ffi v1.1_dev_rc

Hi @hvasbath
Sorry for disturbing you again......
The beat work well for the uniform inversion and I got a very good results.
But when I want to try the distributed slip inversion, I got an error like following (I updated the beat software before I used it in case there are new changes):

[lxrtest@dqcentos7insar BEAT]$ beat import BalochistanEQ7 --results=BalochistanEQ7 --mode='geometry' --datatypes=geodetic
beat         - INFO     Attempting to load results with mode ffi to config_geometry.yaml from directory: BalochistanEQ7
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
models       - INFO     ... Initialising Distribution Optimizer ... 

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

geodetic     - INFO     Number of geodetic datasets: 2 
geodetic     - INFO     Number of geodetic data points: 1190 
geodetic     - WARNING  Covariance estimation not implemented (yet)! Using imported covariances!
geodetic     - INFO     Initialising corrections ...
heart        - INFO     Not correcting Sen201410_201712ascMeter_sub for Euler Pole
heart        - INFO     Setting up Ramps correction for Sen201410_201712ascMeter_sub
heart        - INFO     Not correcting Sen201410_201712desMeter_sub for Euler Pole
heart        - INFO     Setting up Ramps correction for Sen201410_201712desMeter_sub
models       - INFO     Using "none" regularization ...
Traceback (most recent call last):
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/bin/beat", line 11, in <module>
    load_entry_point('beat==1.1', 'console_scripts', 'beat')()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 2041, in main
    globals()['command_' + command](args)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 499, in command_import
    load='trace', chains=[-1])
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/models/base.py", line 387, in load_results
    stage_number = self.handler.highest_sampled_stage()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/backend.py", line 906, in highest_sampled_stage
    return max(self.stage_number(s) for s in glob(self.stage_path('*')))
ValueError: max() arg is an empty sequence

Before running this command, I ran:
beat init BalochistanEQ7 --mode='ffi' --datatypes=geodetic

Could you please give me some suggestion?
Thank you very much!
Best regards

station naming and store id

Station names may start with numbers, but store IDs may NOT! leading to:

beat build_gfs CalifBajaFullMT --datatypes='seismic'. And then, I have the following error message:

heart - INFO Using crust2 profile
heart - INFO Station 035Z
heart - INFO ---------------------
Traceback (most recent call last):
File "/usr/local/bin/beat", line 11, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/usr/local/lib/python3.8/dist-packages/beat-1.0-py3.8-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
globals()'command_' + command
File "/usr/local/lib/python3.8/dist-packages/beat-1.0-py3.8-linux-x86_64.egg/beat/apps/beat.py", line 1032, in command_build_gfs
heart.seis_construct_gf(
File "/usr/local/lib/python3.8/dist-packages/beat-1.0-py3.8-linux-x86_64.egg/beat/heart.py", line 1663, in seis_construct_gf
fomosto_config = get_fomosto_baseconfig(
File "/usr/local/lib/python3.8/dist-packages/beat-1.0-py3.8-linux-x86_64.egg/beat/heart.py", line 1481, in get_fomosto_baseconfig
return gf.ConfigTypeA(
File "/usr/local/lib/python3.8/dist-packages/pyrocko/gf/meta.py", line 1606, in init
self.validate(regularize=True, depth=1)
File "/usr/local/lib/python3.8/dist-packages/pyrocko/guts.py", line 846, in validate
self.T.instance.validate(self, regularize, depth)
File "/usr/local/lib/python3.8/dist-packages/pyrocko/guts.py", line 545, in validate
val = validator.validate_children(val, regularize, depth)
File "/usr/local/lib/python3.8/dist-packages/pyrocko/guts.py", line 557, in validate_children
newpropval = prop.validate(propval, regularize, depth-1)
File "/usr/local/lib/python3.8/dist-packages/pyrocko/guts.py", line 542, in validate
validator.validate_extra(val)
File "/usr/local/lib/python3.8/dist-packages/pyrocko/guts.py", line 1306, in validate_extra
raise ValidationError('%s: "%s" does not match pattern %s' % (
pyrocko.guts.ValidationError: id: "035Z_ak135_2.000Hz_0" does not match pattern '^[A-Za-z][A-Za-z0-9._]{0,64}$'

Update hierarchicals

Hi again. I am trying to set the fit_plane option to true for two interferograms (asc/dsc) for the July 2020 Xizang earthquake. Here I use only geodetic data. I initialize the problem and compute the greens functions.

I then use beat update Xizang --parameters=hierarchicals,hypers to update the parameters, with fit_plane = true. When I run beat sample Xizang --hypers I get a message that says:

(pyrocko) [rzinke@mpzinke-arch BeatTest6]$ beat sample Xizang --hypers
Duplicate key in file '/home/rzinke/.config/matplotlib/matplotlibrc', line 44 ('backend: Qt4Cairo')
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
models       - INFO     ... Initialising Geometry Optimizer ... 

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

geodetic     - INFO     Number of geodetic datasets: 2 
geodetic     - INFO     Number of geodetic data points: 1111 
geodetic     - WARNING  Covariance estimation not implemented (yet)! Using imported covariances!
geodetic     - INFO     Fit residual ramp selected!
models       - INFO     ... Building Hyper model ...

geodetic     - INFO     Estimating ramp for each dataset...
Traceback (most recent call last):
  File "/home/rzinke/miniconda3/envs/pyrocko/bin/beat", line 33, in <module>
    sys.exit(load_entry_point('beat==1.0', 'console_scripts', 'beat')())
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1688, in main
    globals()['command_' + command](args)
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 770, in command_sample
    project_dir, options.mode, options.hypers)
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/problems.py", line 842, in load_model
    problem.built_hyper_model()
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/problems.py", line 261, in built_hyper_model
    self.init_hierarchicals()
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/problems.py", line 546, in init_hierarchicals
    composite.init_hierarchicals(self.config.problem_config)
  File "/home/rzinke/miniconda3/envs/pyrocko/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/geodetic.py", line 242, in init_hierarchicals
    'Plane corrections enabled, but they are'
beat.models.base.ConfigInconsistentError: Plane corrections enabled, but they are not defined in the problem configuration! (hierarchicals)
 Please run: "beat update <project_dir> --parameters="hierarchicals"

I re-run beat update Xizang --parameters=hierarchicals,hypers to be sure. But I still get the same error message.

I checked the config_geometry.yaml file and the hierarchicals are definitely there:

  hierarchicals:
    Tr121dsc_ct03_e002n05m001M01_kite_scene_offset: !beat.heart.Parameter
      name: Tr121dsc_ct03_e002n05m001M01_kite_scene_offset
      form: Uniform
      lower:
      - -0.05
      upper:
      - 0.05
      testvalue:
      - 0.0
    Tr121dsc_ct03_e002n05m001M01_kite_scene_ramp: !beat.heart.Parameter
      name: Tr121dsc_ct03_e002n05m001M01_kite_scene_ramp
      form: Uniform
      lower:
      - -0.005
      - -0.005
      upper:
      - 0.005
      - 0.005
      testvalue:
      - 0.0
      - 0.0
    Tr12asc_ct03_e002n05m001M01_kite_scene_offset: !beat.heart.Parameter
      name: Tr12asc_ct03_e002n05m001M01_kite_scene_offset
      form: Uniform
      lower:
      - -0.05
      upper:
      - 0.05
      testvalue:
      - 0.0
    Tr12asc_ct03_e002n05m001M01_kite_scene_ramp: !beat.heart.Parameter
      name: Tr12asc_ct03_e002n05m001M01_kite_scene_ramp
      form: Uniform
      lower:
      - -0.005
      - -0.005
      upper:
      - 0.005
      - 0.005
      testvalue:
      - 0.0
      - 0.0

Not sure what to do.

Error：running "beat sample Laquila --hypers" in Scnario I

Hi, I want to re-run the Scnario I, but when I run the "beat sample Laquila --hypers" in Scnario I, I meet an error
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: 35
geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
geodetic - INFO Fit residual ramp selected!
models - INFO ... Building Hyper model ...

geodetic - INFO Estimating ramp for each dataset...
geodetic - INFO Initialized 4 hierarchical parameters (ramps).
Traceback (most recent call last):
File "beat/python35/bin/beat", line 11, in
load_entry_point('beat==1.0rc1', 'console_scripts', 'beat')()
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 1588, in main
globals()'command_' + command
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 748, in command_sample
project_dir, options.mode, options.hypers)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/models/problems.py", line 831, in load_model
problem.built_hyper_model()
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/models/problems.py", line 260, in built_hyper_model
self.init_hyperparams()
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/models/problems.py", line 400, in init_hyperparams
raise NotImplementedError('Not fully implemented!')
NotImplementedError: Not fully implemented!

I did not change any values in configure file. Can you give me some suggestion?

pyrocko.gf.meta.OutOfBounds: out of bounds

I am trying to run beat for MT calculations. I am using a downloaded Green's function (ak135_2000km_1Hz) that was converted to beat format and calculated the phases arrivals. But I am getting the error below when I try to compute the hypers:
beat sample --hypers
config - INFO All parameter-priors ok!
models - INFO ... Initialising Geometry Optimizer ...

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

heart - INFO Consistent number of datasets and targets in any_S_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_S_0: 6
models - INFO ... Building Hyper model ...

models - INFO Optimization for 1 hyperparameters in total!
seismic - INFO Retrieving seismic data-covariances with structure "identity" for any_S_0 ...
Traceback (most recent call last):
File "/opt/anaconda2/bin/beat", line 4, in
import('pkg_resources').run_script('beat==1.0rc1', 'beat')
File "/opt/anaconda2/lib/python2.7/site-packages/pkg_resources/init.py", line 666, in run_script
self.require(requires)[0].run_script(script_name, ns)
File "/opt/anaconda2/lib/python2.7/site-packages/pkg_resources/init.py", line 1469, in run_script
exec(script_code, namespace, namespace)
File "/opt/anaconda2/lib/python2.7/site-packages/beat-1.0rc1-py2.7-macosx-10.7-x86_64.egg/EGG-INFO/scripts/beat", line 1485, in

File "/opt/anaconda2/lib/python2.7/site-packages/beat-1.0rc1-py2.7-macosx-10.7-x86_64.egg/EGG-INFO/scripts/beat", line 664, in command_sample

File "build/bdist.macosx-10.7-x86_64/egg/beat/models/problems.py", line 898, in load_model
File "build/bdist.macosx-10.7-x86_64/egg/beat/models/problems.py", line 420, in built_hyper_model
File "build/bdist.macosx-10.7-x86_64/egg/beat/models/seismic.py", line 123, in analyse_noise
File "build/bdist.macosx-10.7-x86_64/egg/beat/covariance.py", line 272, in get_data_covariances
File "build/bdist.macosx-10.7-x86_64/egg/beat/covariance.py", line 221, in do_variance_estimate
File "build/bdist.macosx-10.7-x86_64/egg/beat/heart.py", line 1949, in get_phase_arrival_time
File "/opt/anaconda2/lib/python2.7/site-packages/pyrocko/gf/store.py", line 1663, in t
return timing.evaluate(self.get_phase, args)
File "/opt/anaconda2/lib/python2.7/site-packages/pyrocko/gf/meta.py", line 528, in evaluate
raise OutOfBounds(args)
pyrocko.gf.meta.OutOfBounds: out of bounds: (15,85047)
config_geometry.txt

BUG in slip-distribution plot

But when I run "beat plot Laquila_kinematic scene_fits --mode=ffo" in Scenario III, I meet an error:
beat - INFO Loading problem ...
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
models - INFO ... Initialising Distribution Optimizer ...

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

geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
seismic - INFO Fault(s) discretized to "5.0" [km] patches.
seismic - INFO Preparing data of "any_P" for optimization
models - INFO Using "none" regularization ...
backend - INFO Loading multitrace from beat/test_beat/Laquila_kinematic/ffo/stage_-1
geodetic - INFO Initialized 0 hierarchical parameters (ramps).
plotting - INFO FFO scene fit, using reference source ...
plotting - INFO Loading full resolution kite scene: beat/test_beat/Laquila_insar_data/laquila_kite_container-master/Laquila_dscxn
plotting - INFO Plotting obs
plotting - INFO Plotting syn
plotting - INFO Plotting res
plotting - INFO Plotting reference sources
plotting - INFO Loading full resolution kite scene: beat/test_beat/Laquila_insar_data/laquila_kite_container-master/Laquila_ascxn
plotting - INFO Plotting obs
plotting - INFO Plotting syn
plotting - INFO Plotting res
plotting - INFO Plotting reference sources
plotting - INFO saving figures to beat/test_beat/Laquila_kinematic/ffo/figures/scenes_-1_max_local
[gyxu@insar test_beat]$ beat plot Laquila_kinematic waveform_fits --mode=ffo
beat - INFO Loading problem ...
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
models - INFO ... Initialising Distribution Optimizer ...

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

geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
seismic - INFO Fault(s) discretized to "5.0" [km] patches.
seismic - INFO Preparing data of "any_P" for optimization
models - INFO Using "none" regularization ...
backend - INFO Loading multitrace from beat/test_beat/Laquila_kinematic/ffo/stage_-1
plotting - INFO FFO waveform fit, using reference source ...
plotting - INFO Plotting waveforms ...
plotting - INFO saving figures to beat/test_beat/Laquila_kinematic/ffo/figures/waveforms_-1_max
[gyxu@insar test_beat]$ beat plot Laquila_kinematic moment_rate --mode=ffo
beat - INFO Loading problem ...
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
models - INFO ... Initialising Distribution Optimizer ...

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

geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
seismic - INFO Fault(s) discretized to "5.0" [km] patches.
seismic - INFO Preparing data of "any_P" for optimization
models - INFO Using "none" regularization ...
backend - INFO Loading multitrace from beat/test_beat/Laquila_kinematic/ffo/stage_-1
plotting - INFO Drawing ensemble of moment rate functions ...
plotting - INFO saving figure to beat/test_beat/Laquila_kinematic/ffo/figures/moment_rate_-1_0_max.pdf
[gyxu@insar test_beat]$ beat plot Laquila_kinematic slip_distribution --mode=ffo
beat - INFO Loading problem ...
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
models - INFO ... Initialising Distribution Optimizer ...

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

geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
seismic - INFO Fault(s) discretized to "5.0" [km] patches.
seismic - INFO Preparing data of "any_P" for optimization
models - INFO Using "none" regularization ...
backend - INFO Loading multitrace from beat/test_beat/Laquila_kinematic/ffo/stage_-1
Traceback (most recent call last):
File "beat/python35/bin/beat", line 11, in
load_entry_point('beat==1.0rc1', 'console_scripts', 'beat')()
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 1588, in main
globals()'command_' + command
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/apps/beat.py", line 1295, in command_plot
plotting.plots_catalog[plot](problem, po)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/plotting.py", line 2655, in draw_slip_dist
fault, mtrace, transform=transform, reference=reference)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/plotting.py", line 2532, in fault_slip_distribution
alpha=1., linewidth=7, zorder=-1)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/plotting.py", line 2368, in fuzzy_rupture_fronts
linewidth=linewidth)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/plotting.py", line 2819, in draw_line_on_array
check_grid_shape(ngr, naim, axis)
File ".local/lib/python3.5/site-packages/beat-1.0rc1-py3.5-linux-x86_64.egg/beat/plotting.py", line 2777, in check_grid_shape
' Expected %i got %i' % (axis, naim, ngr))
TypeError: Gridsize of given grid is inconistent for axis 0! Expected 250 got 375

Other command can run successfully except the plot of slip distribution. Can you give me some suggestion?

error in beat summarize

Hi again.
I am trying to summarize the finial results after the sampling successfully finished. When I ran the command "beat summarize BalochistanEQ7 --stage_number=-1 --mode=ffi", I meet an error: "KeyError: 'uparr' at uparr = self.vector2subfault(index, point['uparr'])". I tried to debug it and found that this line worked well for index=0 and reported error when index=1. I have 7 subfaults.
When index=0, the keys in point dict are:
{uparr, uperp, h_SAR, h_laplacian, Sen201410_201712ascMeter_sub_azimuth_ramp, Sen201410_201712ascMeter_sub_range_ramp, Sen201410_201712ascMeter_sub_offset, Sen201410_201712desMeter_sub_azimuth_ramp, Sen201410_201712desMeter_sub_range_ramp, Sen201410_201712desMeter_sub_offset, geo_like, laplacian_like, like}
When index=1, the point dict is:
{'slip': 0.8816954847133593, 'rake': 12.847378570297858}

The error message is:

beat summarize BalochistanEQ7 --stage_number=-1 --mode=ffi
beat         - INFO     Loading problem ...
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
models       - INFO     ... Initialising Distribution Optimizer ... 

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

geodetic     - INFO     Number of geodetic datasets: 2 
geodetic     - INFO     Number of geodetic data points: 1190 
geodetic     - WARNING  Covariance estimation not implemented (yet)! Using imported covariances!
geodetic     - INFO     Initialising corrections ...
heart        - INFO     Not correcting Sen201410_201712ascMeter_sub for Euler Pole
heart        - INFO     Setting up Ramps correction for Sen201410_201712ascMeter_sub
heart        - INFO     Not correcting Sen201410_201712desMeter_sub for Euler Pole
heart        - INFO     Setting up Ramps correction for Sen201410_201712desMeter_sub
models       - INFO     ... Building model ...

models       - INFO     Initialised hyperparameter h_SAR with size 1 
models       - INFO     Initialised hyperparameter h_laplacian with size 1 
models       - INFO     Optimization for 2 hyperparameters in total!
geodetic     - INFO     Loading geodetic Green's Functions
ffi          - INFO     Setting geodetic_uparr_static_0 GF Library to optimization mode.
ffi          - INFO     Setting geodetic_uperp_static_0 GF Library to optimization mode.
geodetic     - INFO     Evaluating config for Euler Pole corrections for datasets...
geodetic     - INFO     No Euler Pole correction!
geodetic     - INFO     Evaluating config for Ramps corrections for datasets...
geodetic     - INFO     Initialized 6 hierarchical parameters.
ffi.laplacian - INFO     Initialising Laplacian smoothing operator ...
models       - INFO     Model building was successful! 

beat         - INFO     Summarizing stage(s): -1
beat         - INFO     Summarizing stage under: /data/lxrtest/insarlab/BEAT/BalochistanEQ7/ffi/stage_-1
backend      - INFO     Loading multitrace from /data/lxrtest/insarlab/BEAT/BalochistanEQ7/ffi/stage_-1
backend      - INFO     Adding derived variables magnitude to trace.
  0%|                                                                                                                            | 0/1200 [00:00<?, ?it/s]
Traceback (most recent call last):
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/bin/beat", line 11, in <module>
    load_entry_point('beat==1.1', 'console_scripts', 'beat')()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 2041, in main
    globals()['command_' + command](args)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 1053, in command_summarize
    composite.point2sources(point)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/models/geodetic.py", line 821, in point2sources
    return self.fault.point2sources(tpoint, events=events)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/ffi/fault.py", line 591, in point2sources
    uparr = self.vector2subfault(index, point['uparr'])
KeyError: 'uparr'

Did I mess something or forget to set something? I have no idea what to do...
Thank you!

v1.1 features need tutorials

FFI discretization options
Laplacian configuration
hierarchical corrections: Euler Pole, Ramp
multi segment setup in Talpa
GNSS plotting (also needs improvement)
Multi-event setup

Error: running QSSP with CRUST2

Transferred from personal mail from @fangjinSGG
When I calculate seismic GF functions replacing the crust from the ak135 earth model with crust from the crust2 model based on QSSP backend, an error called "STOP Error: earth radius smaller than pre-defined! statement executed" is thrown out. I have no idea why it happened. The page (e359b98#diff-0a7529cc0f3035dcda0414a3cb4d0b06R1569) says that you have changed the qssp version to '2010beta' from '2010'. But I can't find the qssp2010beta version.

Question about beatdown

Hi, I look at the tutorial about the beatdown, I want to know that how can I download the seismic data of stations with epicentral distance 30-90 degree? Whether I should download the available data in a radius of 1000 km and then select data with epicentral distance 30-90 degree?

Config Inconsistency and data import

Error when running beat sample Scenario II, regularization, ramps

Hi,
When I running "beat sample Laquila --hypers --mode=ffo" in Scenario II, I meet a error that
beat sample Laquila --mode=ffo
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
models - INFO ... Initialising Distribution Optimizer ...
models - INFO Analysing problem ...
models - INFO ---------------------
geodetic - INFO Number of geodetic datasets: 2
geodetic - INFO Number of geodetic data points: 419
geodetic - WARNING Covariance estimation not implemented (yet)! Using imported covariances!
geodetic - INFO Fit residual ramp selected!
Traceback (most recent call last):
File "/beat/python35/lib/python3.5/site-packages/pymc3/model.py", line 179, in get_context
return cls.get_contexts()[-1]
IndexError: list index out of range

Can you help me solve this problem? Thank you!

docs: rework existing tutorials for new functionality

import function
filterers
beat check geometry vs discretization
FFI shear vs tensile

TypeError: Datatype "seismic" not included in FaultGeometry

Hi Hannes,
when I run the code "beat build_gfs Laquila_kinematic --datatypes=seismic --execute --mode=ffo
" in Scenario III, it appears the error:
fangnandeMacBook-Pro:beat_learning NAN$ beat build_gfs Laquila_kinematic --datatypes=seismic --execute --mode=ffo
config - INFO All parameter-priors ok!
config - INFO All hyper-parameters ok!
config - INFO All hierarchical-parameters ok!
beat - INFO Discretized fault geometry exists! Use --force to overwrite!
beat - INFO Loading existing discretized fault
beat - INFO Calculating linear Green's Functions
beat - INFO for seismic data ...
beat - INFO crust_ind 0
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
ffo - INFO Calculating GFs for starttimes: "-13.0", "-12.0", "-11.0", "-10.0", "-9.0", "-8.0", "-7.0", "-6.0", "-5.0", "-4.0", "-3.0", "-2.0", "-1.0", "0.0", "1.0", "2.0", "3.0", "4.0", "5.0", "6.0", "7.0", "8.0", "9.0", "10.0", "11.0", "12.0", "13.0", "14.0", "15.0"
durations: "0.0", "0.25", "0.5", "0.75", "1.0", "1.25", "1.5", "1.75", "2.0", "2.25", "2.5", "2.75", "3.0", "3.25", "3.5", "3.75", "4.0"
ffo - INFO Using 4 workers ...
ffo - INFO For slip component: uparr
Traceback (most recent call last):
File "/anaconda3/bin/beat", line 11, in
load_entry_point('beat==1.0rc1', 'console_scripts', 'beat')()
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/apps/beat.py", line 1588, in main
globals()'command_' + command
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/apps/beat.py", line 1144, in command_build_gfs
force=options.force)
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/ffo/base.py", line 1043, in seis_construct_gf_linear
datatype='seismic', component=var),
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/ffo/fault.py", line 161, in get_all_subfaults
index=i, datatype=datatype, component=component))
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/ffo/fault.py", line 146, in get_subfault
source_key = self.get_subfault_key(index, datatype, component)
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/ffo/fault.py", line 88, in get_subfault_key
self._check_datatype(datatype)
File "/anaconda3/lib/python3.6/site-packages/beat-1.0rc1-py3.6-macosx-10.7-x86_64.egg/beat/ffo/fault.py", line 75, in _check_datatype
'Datatype "%s" not included in FaultGeometry' % datatype)
TypeError: Datatype "seismic" not included in FaultGeometry

Can you give me some suggestions? Thank you very much.

How to get FULL original InSAR data in Scenario I?

Hi,
When I run "beat plot Laquila scene_fits", I can not get the figures with content. The result figure is blank. I noticed that the Scenario I said to generate the results the FULL original data needs to be downloaded and put into the specified data path. But I don't know how to get the FULL original InSAR data? Can you help me? Thank you!

Polarity / Ampspec fitting

first motions (polarity) likelihood & amplitude ratio likelihood

Error when running 'beat build_gfs FullMT --datatypes='seismic' --force --execute'

Hi, when I run 'beat build_gfs FullMT --datatypes='seismic' --force --execute' with the dataset that I download using 'beatdown', I meet an error said that:

beat build_gfs FullMT --datatypes='seismic' --force --execute
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors ok!
beat - INFO Creating Green's Function stores individually for each station!
heart - INFO Using crust2 profile
heart - INFO Water layer 1329.000000 in CRUST model! Remove and add to lower crust
heart - INFO New Lower crust layer thickness 11995.000000
heart - INFO Station SOCY
heart - INFO ---------------------
heart - INFO Creating Store at ./SOCY_ak135_2.000Hz_0
heart - INFO Using crust2 profile
heart - INFO Water layer 1792.000000 in CRUST model! Remove and add to lower crust
heart - INFO New Lower crust layer thickness 6958.000000
heart - INFO Filling store ...
pyrocko.gf.store - INFO making travel time table for phasegroup "any_P"
pyrocko.spit - INFO at level 0: 100.0% covered, 1 cell
pyrocko.gf.store - INFO making travel time table for phasegroup "any_S"
pyrocko.spit - INFO at level 0: 0.0% covered, 1 cell
pyrocko.spit - INFO at level 1: 100.0% covered, 3 cells
pyrocko.fomosto.qssp - INFO Starting step 1 / 2, block 1 / 11
pyrocko.fomosto.qssp - INFO Starting step 1 / 2, block 3 / 11
pyrocko.fomosto.qssp - INFO Starting step 1 / 2, block 2 / 11
pyrocko.fomosto.qssp - INFO Starting step 1 / 2, block 4 / 11
pyrocko.fomosto.qssp - WARNING qssp emitted something via stderr:

Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

pyrocko.fomosto.qssp - INFO Done with step 1 / 2, block 1 / 11, wallclock time: 3229 s
pyrocko.fomosto.qssp - INFO Starting step 1 / 2, block 5 / 11
pyrocko.fomosto.qssp - WARNING qssp emitted something via stderr:

Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

pyrocko.fomosto.qssp - INFO Done with step 1 / 2, block 2 / 11, wallclock time: 3274 s
pyrocko.fomosto.qssp - WARNING qssp emitted something via stderr:

Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

pyrocko.fomosto.qssp - INFO Done with step 1 / 2, block 3 / 11, wallclock time: 3278 s
pyrocko.fomosto.qssp - WARNING qssp emitted something via stderr:

Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

pyrocko.fomosto.qssp - INFO Done with step 1 / 2, block 4 / 11, wallclock time: 3279 s
Traceback (most recent call last):
File "/usr/local/bin/beat", line 9, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/usr/local/lib/python3.6/dist-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1641, in main
globals()'command_' + command
File "/usr/local/lib/python3.6/dist-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1005, in command_build_gfs
force=options.force)
File "/usr/local/lib/python3.6/dist-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/heart.py", line 1709, in seis_construct_gf
store_dir, nworkers=sf.nworkers, force=force)
File "/usr/local/lib/python3.6/dist-packages/pyrocko/fomosto/qssp.py", line 890, in build
step=step, iblock=iblock)
File "/usr/local/lib/python3.6/dist-packages/pyrocko/gf/builder.py", line 166, in build
Interrupted, store.StoreError)):
File "/usr/local/lib/python3.6/dist-packages/pyrocko/parimap.py", line 136, in parimap
raise exc
pyrocko.fomosto.qssp.QSSPError: ===== begin qssp input =====

autogenerated QSSP input by qssp.py

This is the input file of FORTRAN77 program "qssp2010" for calculating

synthetic seismograms of a self-gravitating, spherically symmetric,

isotropic and viscoelastic earth.

by

Rongjiang Wang [email protected]

Helmholtz-Centre Potsdam

GFZ German Reseach Centre for Geosciences

Telegrafenberg, D-14473 Potsdam, Germany

Last modified: Potsdam, July, 2010

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

If not specified, SI Unit System is used overall!

Coordinate systems:

spherical (r,t,p) with r = radial,

t = co-latitude,

p = east longitude.

= = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =

UNIFORM RECEIVER DEPTH

======================

1. uniform receiver depth [km]

#-------------------------------------------------------------------------------------------
0.000000e+00
#-------------------------------------------------------------------------------------------

TIME (FREQUENCY) SAMPLING

=========================

1. time window [sec], sampling interval [sec]

2. max. frequency [Hz] of Green's functions

3. max. slowness [s/km] of Green's functions

Note: if the near-field static displacement is desired, the maximum slowness should not

be smaller than the S wave slowness in the receiver layer

4. anti-aliasing factor (> 0 & < 1), if it is <= 0 or >= 1/e (~ 0.4), then

default value of 1/e is used (e.g., 0.1 = alias phases will be suppressed

to 10% of their original amplitude)

Note: The computation effort increases linearly the time window and

quadratically with the cut-off frequency.

#-------------------------------------------------------------------------------------------
7.385000e+02 5.000000e-01
1.000000e+00
1.326313e-01
1.000000e-01
#-------------------------------------------------------------------------------------------

SELF-GRAVITATING EFFECT

=======================

1. the critical frequency [Hz] and the critical harmonic degree, below which

the self-gravitating effect should be included

#-------------------------------------------------------------------------------------------
0.000000e+00 0
#-------------------------------------------------------------------------------------------

WAVE TYPES

==========

1. selection (1/0 = yes/no) of speroidal modes (P-SV waves), selection of toroidal modes

(SH waves), minimum and maximum cutoff harmonic degrees

Note: if the near-field static displacement is desired, the minimum

cutoff harmonic degree should not be smaller than, e.g., 2000.

#-------------------------------------------------------------------------------------------
1 1 0 25000
#-------------------------------------------------------------------------------------------

GREEN'S FUNCTION FILES

======================

1. number of discrete source depths, estimated radius of each source patch [km] and

directory for Green's functions

2. list of the source depths [km], the respective file names of the Green's

functions (spectra) and the switch number (0/1) (0 = do not calculate

this Green's function because it exists already, 1 = calculate or update

this Green's function. Note: update is required if any of the above

parameters is changed)

#-------------------------------------------------------------------------------------------
1 9.500000e-01 '/home/matrix/my_beat_case/Alor_seismic/AlorEQ/SOCY_ak135_2.000Hz_0/qssp_green/'
4.000000e+00 'TEMP_4km_0km' 1
#-------------------------------------------------------------------------------------------

MULTI-EVENT SOURCE PARAMETERS

=============================

1. number of discrete point sources and selection of the source data format

(1 or 2)

2. list of the multi-event sources

Format 1:

M-Unit Mrr Mtt Mpp Mrt Mrp Mtp Lat Lon Depth T_origin T_rise

[Nm] [deg] [deg] [km] [sec] [sec]

Format 2:

Moment Strike Dip Rake Lat Lon Depth T_origin T_rise

[Nm] [deg] [deg] [deg] [deg] [deg] [km] [sec] [sec]

#-------------------------------------------------------------------------------------------
1 1
1.000000e+00 1.000000e+00 1.000000e+00 1.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 8.999999e+01 0.000000e+00 1.000000e+01 0.000000e+00 5.000000e-04
#-------------------------------------------------------------------------------------------

RECEIVER PARAMETERS

===================

1. output file name and and selection of output format:

1 = cartesian: vertical(z)/north(n)/east(e);

2 = spherical: radial(r)/theta(t)/phi(p)

(Note: if output format 2 is selected, the epicenter (T_origin = 0)

2. output time window [sec] (<= Green's function time window)

3. selection of order of Butterworth bandpass filter (if <= 0, then no

filtering), lower and upper corner frequencies (smaller than the cut-off

frequency defined above)

4. lower and upper slowness cut-off [s/km] (slowness band-pass filter)

5. number of receiver

6. list of the station parameters

Format:

Lat Lon Name Time_reduction

[deg] [deg] [sec]

(Note: Time_reduction = start time of the time window)

#-------------------------------------------------------------------------------------------
'receivers' 1
7.385000e+02
0 1.000000 1.000000
0.000000e+00 1.326313e-01
1
1.664765e+01 1.800000e+02 '' 6.240000e+02
#-------------------------------------------------------------------------------------------

LAYERED EARTH MODEL (IASP91)

============================

1. number of data lines of the layered model and selection for including

the physical dispersion according Kamamori & Anderson (1977)

#-------------------------------------------------------------------------------------------
47 0
#-------------------------------------------------------------------------------------------

MULTILAYERED MODEL PARAMETERS (source site)

===========================================

no depth[km] vp[km/s] vs[km/s] ro[g/cm^3] qp qs

#-------------------------------------------------------------------------------------------
1 0.000000e+00 2.200000e+00 1.100000e+00 2.200000e+00 5.000000e+01 5.000000e+01
2 3.340000e-01 2.200000e+00 1.100000e+00 2.200000e+00 5.000000e+01 5.000000e+01
3 3.340000e-01 5.000000e+00 2.500000e+00 2.600000e+00 6.000000e+02 4.000000e+02
4 3.034000e+00 5.000000e+00 2.500000e+00 2.600000e+00 6.000000e+02 4.000000e+02
5 3.034000e+00 6.600000e+00 3.650000e+00 2.900000e+00 1.471000e+03 6.000000e+02
6 6.334000e+00 6.600000e+00 3.650000e+00 2.900000e+00 1.471000e+03 6.000000e+02
7 6.334000e+00 7.100000e+00 3.900000e+00 3.050000e+00 1.491000e+03 6.000000e+02
8 1.329000e+01 7.100000e+00 3.900000e+00 3.050000e+00 1.491000e+03 6.000000e+02
9 1.329000e+01 8.036000e+00 4.486000e+00 3.638000e+00 9.497000e+02 3.946000e+02
10 8.000000e+01 8.040000e+00 4.481000e+00 3.500000e+00 1.008000e+03 4.176000e+02
11 8.000000e+01 8.045000e+00 4.490000e+00 3.502000e+00 1.820000e+02 7.560000e+01
12 1.200000e+02 8.050000e+00 4.500000e+00 3.427000e+00 1.826000e+02 7.606000e+01
13 1.200000e+02 8.050000e+00 4.500000e+00 3.427000e+00 1.826000e+02 7.606000e+01
14 1.650000e+02 8.175000e+00 4.509000e+00 3.371000e+00 1.887000e+02 7.655000e+01
15 2.100000e+02 8.301000e+00 4.518000e+00 3.324000e+00 2.010000e+02 7.940000e+01
16 2.100000e+02 8.300000e+00 4.519000e+00 3.323000e+00 3.382000e+02 1.337000e+02
17 3.000000e+02 8.628000e+00 4.679000e+00 3.401000e+00 3.536000e+02 1.387000e+02
18 4.100000e+02 9.030000e+00 4.870000e+00 3.506000e+00 3.775000e+02 1.465000e+02
19 4.100000e+02 9.360000e+00 5.080000e+00 3.929000e+00 4.141000e+02 1.627000e+02
20 6.600000e+02 1.020000e+01 5.611000e+00 3.918000e+00 4.285000e+02 1.729000e+02
21 6.600000e+02 1.079000e+01 5.965000e+00 4.240000e+00 1.349000e+03 5.495000e+02
22 7.640000e+02 1.107000e+01 6.215000e+00 4.359000e+00 1.276000e+03 5.371000e+02
23 8.491000e+02 1.121000e+01 6.272000e+00 4.463000e+00 1.263000e+03 5.273000e+02
24 1.038000e+03 1.152000e+01 6.407000e+00 4.616000e+00 1.230000e+03 5.069000e+02
25 1.227000e+03 1.181000e+01 6.527000e+00 4.710000e+00 1.198000e+03 4.880000e+02
26 1.416000e+03 1.208000e+01 6.636000e+00 4.812000e+00 1.168000e+03 4.704000e+02
27 1.605000e+03 1.233000e+01 6.735000e+00 4.909000e+00 1.141000e+03 4.544000e+02
28 1.795000e+03 1.255000e+01 6.827000e+00 5.004000e+00 1.111000e+03 4.377000e+02
29 1.984000e+03 1.278000e+01 6.913000e+00 5.096000e+00 1.065000e+03 4.156000e+02
30 2.173000e+03 1.300000e+01 6.997000e+00 5.186000e+00 1.037000e+03 4.006000e+02
31 2.362000e+03 1.321000e+01 7.079000e+00 5.273000e+00 9.954000e+02 3.812000e+02
32 2.551000e+03 1.343000e+01 7.162000e+00 5.357000e+00 9.681000e+02 3.673000e+02
33 2.740000e+03 1.365000e+01 7.248000e+00 5.439000e+00 9.325000e+02 3.506000e+02
34 2.740000e+03 1.365000e+01 7.248000e+00 5.693000e+00 7.227000e+02 2.717000e+02
35 2.790000e+03 1.365000e+01 7.259000e+00 5.720000e+00 7.269000e+02 2.740000e+02
36 2.839000e+03 1.366000e+01 7.270000e+00 5.746000e+00 7.251000e+02 2.740000e+02
37 2.892000e+03 1.366000e+01 7.282000e+00 5.772000e+00 7.231000e+02 2.740000e+02
38 2.892000e+03 7.972000e+00 0.000000e+00 9.928000e+00 5.782000e+04 0.000000e+00
39 3.344000e+03 8.744000e+00 0.000000e+00 1.059000e+01 5.782000e+04 0.000000e+00
40 3.797000e+03 9.338000e+00 0.000000e+00 1.112000e+01 5.782000e+04 0.000000e+00
41 4.249000e+03 9.759000e+00 0.000000e+00 1.155000e+01 5.782000e+04 0.000000e+00
42 4.702000e+03 1.009000e+01 0.000000e+00 1.188000e+01 5.782000e+04 0.000000e+00
43 5.154000e+03 1.032000e+01 0.000000e+00 1.214000e+01 5.782000e+04 0.000000e+00
44 5.154000e+03 1.104000e+01 3.505000e+00 1.270000e+01 6.328000e+02 8.503000e+01
45 5.397000e+03 1.112000e+01 3.567000e+00 1.282000e+01 6.192000e+02 8.503000e+01
46 5.884000e+03 1.123000e+01 3.647000e+00 1.297000e+01 6.050000e+02 8.503000e+01
47 6.370999e+03 1.127000e+01 3.672000e+00 1.302000e+01 6.006000e+02 8.503000e+01
#---------------------------------end of all inputs-----------------------------------------
===== end qssp input =====
===== begin qssp output =====
===== end qssp output =====
===== begin qssp error =====
===== end qssp error =====
qssp had a non-zero exit state: -11
qssp has been invoked as "fomosto_qssp2010"

I don't know what went wrong. Can you give me some suggestions?

Error: get_unit(): Bad internal unit KIND

I am working through the BEAT tutorial Example 3: Rectangular source. This time, I am using a Linux server (CentOS Linux 7).

I installed BEAT using the Anaconda installation method. Since I don't have super user privileges, I installed fomosto-psgrn-pscmp and fomosto-qssp with ./configure --prefix=/u/sarh0/rzinke/local and added /u/sarh0/rzinke/local/bin to my PATH.

I made it part way through Example 3. I successfully built the Green's functions for the geodetic data, but then for the seismic data I am getting an error.

$ beat build_gfs Laquila --datatypes='seismic' --execute

yields:

pyrocko.fomosto.qssp - INFO     Starting step 1 / 2, block 1 / 8
pyrocko.fomosto.qssp - INFO     Starting step 1 / 2, block 2 / 8
pyrocko.fomosto.qssp - INFO     Starting step 1 / 2, block 3 / 8
pyrocko.fomosto.qssp - INFO     Starting step 1 / 2, block 4 / 8
pyrocko.fomosto.qssp - WARNING  qssp emitted something via stderr: 

At line 72 of file qpgetinp.f
Internal Error: get_unit(): Bad internal unit KIND

pyrocko.fomosto.qssp - WARNING  qssp emitted something via stderr: 

At line 72 of file qpgetinp.f
Internal Error: get_unit(): Bad internal unit KIND

pyrocko.fomosto.qssp - WARNING  qssp emitted something via stderr: 

At line 72 of file qpgetinp.f
Internal Error: get_unit(): Bad internal unit KIND

Traceback (most recent call last):
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/bin/beat", line 33, in <module>
    sys.exit(load_entry_point('beat==1.0', 'console_scripts', 'beat')())
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
    globals()['command_' + command](args)
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1038, in command_build_gfs
    force=options.force)
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/heart.py", line 1709, in seis_construct_gf
    store_dir, nworkers=sf.nworkers, force=force)
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/pyrocko/fomosto/qssp.py", line 1141, in build
    step=step, iblock=iblock)
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/pyrocko/gf/builder.py", line 165, in build
    Interrupted, store.StoreError)):
  File "/u/sar-r0/rzinke/python/miniconda3/envs/beat/lib/python3.6/site-packages/pyrocko/parimap.py", line 139, in parimap
    raise exc
pyrocko.fomosto.qssp.QSSPError: ===== begin qssp input =====

The full output is longer, but it ends with

At line 72 of file qpgetinp.f
Internal Error: get_unit(): Bad internal unit KIND

Is this a Fortran compiler issue?

beat sample Laquila error

Hello,

I am trying to complete Example 3[https://hvasbath.github.io/beat/examples/Rectangular.html] from the tutorial.
I have completed part of the tutorial, however, when I attempt to run beat sample Laquila I get an error message that I have trouble interpreting. The message is quite long, and in addition to pasting it below, I have attached it in a text file if that is easier to read.

I have also included my config_geometry.yaml file if that helps.

I am using a Mac running Mojave 10.14.6.

Please let me know where to start fixing this, and how I can help troubleshoot!

beat sample Laquila Error.txt

config_geometry.txt

error when using multi fault--could not broadcast input array from shape (2,1) into shape (1,2) for "magnitude" attribute in lpoint

Hi @hvasbath
The developed version of beat works well for single fault, but I meet an error when I tried to use multi-faults. The error occurred when the command "beat summarize --stage_number=-1" and the error message is:

[lxrtest@dqcentos7insar BalochistanEQ2]$ beat summarize --stage_number=-1
beat         - INFO     Loading problem ...
config       - INFO     All hierarchicals ok!
config       - INFO     All hyperparameters ok!
config       - INFO     All priors ok!
models       - INFO     ... Initialising Geometry Optimizer ... 

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

geodetic     - INFO     Number of geodetic datasets: 2 
geodetic     - INFO     Number of geodetic data points: 919 
geodetic     - WARNING  Covariance estimation not implemented (yet)! Using imported covariances!
geodetic     - INFO     Initialising corrections ...
heart        - INFO     Not correcting Sen201410_201712asc for Euler Pole
heart        - INFO     Not correcting Sen201410_201712asc for Ramps
heart        - INFO     Not correcting Sen201410_201712des for Euler Pole
heart        - INFO     Not correcting Sen201410_201712des for Ramps
models       - INFO     ... Building model ...

models       - INFO     Initialised hyperparameter h_SAR with size 1 
models       - INFO     Optimization for 1 hyperparameters in total!
geodetic     - INFO     Geodetic optimization on: 
 depth, dip, east_shift, length, north_shift, rake, slip, strike, width
geodetic     - INFO     Evaluating config for Euler Pole corrections for datasets...
geodetic     - INFO     No Euler Pole correction!
geodetic     - INFO     Evaluating config for Ramps corrections for datasets...
geodetic     - INFO     Initialized 0 hierarchical parameters.
geodetic     - INFO     Applying corrections! ...
models       - INFO     Model building was successful! 

beat         - INFO     Summarizing stage(s): -1
beat         - INFO     Summarizing stage under: /data/lxrtest/insarlab/BEAT/BalochistanEQ2/geometry/stage_-1
backend      - INFO     Loading multitrace from /data/lxrtest/insarlab/BEAT/BalochistanEQ2/geometry/stage_-1
backend      - INFO     Adding derived variables magnitude to trace.
100%|████████████████████████████████████████████████████████████████████████████████████████████████████████████████| 1000/1000 [00:02<00:00, 405.79it/s]
Traceback (most recent call last):
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/bin/beat", line 11, in <module>
    load_entry_point('beat==1.1', 'console_scripts', 'beat')()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 2041, in main
    globals()['command_' + command](args)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/apps/beat.py", line 1105, in command_summarize
    rtrace.record_buffer()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/backend.py", line 333, in record_buffer
    self._write_data_to_file()
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/backend.py", line 728, in _write_data_to_file
    lpoint2file(fh, self.varnames, data, lpoint)
  File "/home/lxrtest/wclsolid/pyenv/versions/miniconda3-4.7.12/lib/python3.7/site-packages/beat-1.1-py3.7-linux-x86_64.egg/beat/backend.py", line 711, in lpoint2file
    data[names] = array
ValueError: could not broadcast input array from shape (2,1) into shape (1,2)

I tried to debug to find the reason for this error and it seems that when it tried to summarize "magnitude" attribute in lpoint into "data", the size of "magnitude" is (2, 1) but the size of other attributes in lpoint is (2,). The details are following:

(Pdb) names
'magnitude'
(Pdb) lpoint[-1]
array([[7.17854651],
       [7.26610615]])

(Pdb) lpoint[0]
array([0.49622183, 0.44478329])

Could you please give me some suggestions about how to change lpoint?
Thank you!

Warning: The source's sub-sources are further than lambda_min/4 apart

Hi, when I run "beat sample Laquila" in Scenario I, there are many WARNING said that "The source's sub-sources are further than lambda_min/4 apart". Does it matter?

beat.sampler.distributed.MPIError:

Hi all,

1) Running the second example with n_chains=8, I got:

(beat) lab@lab:~$ beat sample Laquila_dc
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors 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: 35
models - INFO ... Building model ...

models - INFO Optimization for 35 hyperparameters in total!
seismic - INFO Seismic optimization on:
depth, dip, duration, east_shift, magnitude, north_shift, rake, strike, time
seismic - INFO Retrieving seismic data-covariances with structure "variance" for any_P_0 ...
seismic - INFO Preparing data of "any_P_0" for optimization
seismic - INFO Initializing synthesizer for "any_P_0"
models - INFO Model building was successful!

models - INFO Using "bin" backend to store samples!
models - INFO ... Initiate Metropolis for Parallel Tempering...
proposal_distribution: Normal, tune_interval=10, n_chains=8

models.base - INFO ... Starting Parallel Tempering ...

distributed - INFO MPI is installed!
distributed - INFO Detected python3
distributed - INFO Done initialising mpi runner
distributed - INFO sampler directory: /home/lab/bayesiano/beat/beat/sampler/pt.py
Traceback (most recent call last):
File "/home/lab/anaconda3/envs/beat/bin/beat", line 11, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
globals()'command_' + command
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 775, in command_sample
sample(step, problem)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/base.py", line 196, in sample
record_worker_chains=pa.record_worker_chains)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/pt.py", line 836, in pt_sample
project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 226, in run_mpi_sampler
samplerdir, n_jobs=n_jobs, loglevel=loglevel, project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 158, in run
in the directory %s'''.lstrip() % (program, self.tempdir))
beat.sampler.distributed.MPIError: mpiexec has been invoked as "mpiexec"
in the directory /tmp/mpiexec-56qm7lyu
distributed - WARNING not removing temporary directory: /tmp/mpiexec-56qm7lyu

------------------------------------------------------------*-

2. Running again in --loglevel=debug:

(beat) lab@lab:~$ beat sample Laquila_dc
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors 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: 35
models - INFO ... Building model ...

models - INFO Optimization for 35 hyperparameters in total!
seismic - INFO Seismic optimization on:
depth, dip, duration, east_shift, magnitude, north_shift, rake, strike, time
seismic - INFO Retrieving seismic data-covariances with structure "variance" for any_P_0 ...
seismic - INFO Preparing data of "any_P_0" for optimization
seismic - INFO Initializing synthesizer for "any_P_0"
models - INFO Model building was successful!

models - INFO Using "bin" backend to store samples!
models - INFO ... Initiate Metropolis for Parallel Tempering...
proposal_distribution: Normal, tune_interval=10, n_chains=8

models.base - INFO ... Starting Parallel Tempering ...

distributed - INFO MPI is installed!
distributed - INFO Detected python3
distributed - INFO Done initialising mpi runner
distributed - INFO sampler directory: /home/lab/bayesiano/beat/beat/sampler/pt.py
Traceback (most recent call last):
File "/home/lab/anaconda3/envs/beat/bin/beat", line 11, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
globals()'command_' + command
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 775, in command_sample
sample(step, problem)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/base.py", line 196, in sample
record_worker_chains=pa.record_worker_chains)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/pt.py", line 836, in pt_sample
project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 226, in run_mpi_sampler
samplerdir, n_jobs=n_jobs, loglevel=loglevel, project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 158, in run
in the directory %s'''.lstrip() % (program, self.tempdir))
beat.sampler.distributed.MPIError: mpiexec has been invoked as "mpiexec"
in the directory /tmp/mpiexec-56qm7lyu
distributed - WARNING not removing temporary directory: /tmp/mpiexec-56qm7lyu
(beat) lab@lab:~$ beat sample Laquila_dc --loglevel=debug
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors ok!
models - INFO ... Initialising Geometry Optimizer ...

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

seismic - DEBUG Setting up seismic structure ...

utility - DEBUG Valid distance range: [30.000000, 90.000000]!
utility - DEBUG Distance of station TSUM: 61.349855 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station RCBR: 65.567308 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station RES: 54.143590 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KMI: 73.109289 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SSPA: 64.946397 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station BTDF: 89.437223 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ANWB: 67.439618 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station MBO: 38.167185 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station TEIG: 84.816714 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KBL: 43.633450 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station MSEY: 60.417579 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station DBIC: 38.986606 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station HIA: 68.244425 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station YAK: 64.054778 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station QIZ: 82.062195 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KBS: 36.698309 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station LSZ: 58.931762 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station TLY: 58.456908 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station INCN: 79.971426 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station DWPF: 74.925103 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station WMQ: 52.564826 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station AAK: 44.192670 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station XAN: 71.681441 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KMBO: 48.335702 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KURK: 44.334011 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ADYE: 36.168619 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SDV: 79.664513 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ARU: 31.995457 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SAML: 86.199701 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ALE: 46.405043 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station OTAV: 91.209244 [deg]
utility - DEBUG Outside defined distance range!
utility - DEBUG Distance of station ABKT: 34.203018 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station AKTO: 31.451507 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station BJT: 72.303754 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SSE: 81.349799 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station CHTO: 74.397519 [deg]
utility - DEBUG Inside defined distance range!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
models - INFO ... Building model ...

models - DEBUG models - INFO seismic - INFO depth, dip, seismic - INFO heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG seismic - INFO heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG heart - DEBUG a: 1238982169.0
b: 1238982174.0
c: 1238982224.0
d: 1238982239.0 Optimization for 1 sources
Optimization for 35 hyperparameters in total!
Seismic optimization on:
duration, east_shift, magnitude, north_shift, rake, strike, time
Retrieving seismic data-covariances with structure "variance" for any_P_0 ...
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 61.349855 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 65.567308 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 54.143590 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 73.109289 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 64.946397 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 89.437223 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 67.439618 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 38.167185 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 84.816714 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 43.633450 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 60.417579 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 38.986606 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 68.244425 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 64.054778 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 82.062195 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 36.698309 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 58.931762 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 58.456908 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 79.971426 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 74.925103 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 52.564826 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 44.192670 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 71.681441 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 48.335702 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 44.334011 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 36.168619 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 79.664513 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 31.995457 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 86.199701 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 46.405043 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 34.203018 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 31.451507 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 72.303754 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 81.349799 [deg]
Wavename not specified using first tabulated phase! any_P
Arrival time for wavename "any_P" distance 74.397519 [deg]
Preparing data of "any_P_0" for optimization
Arrival time for wavename "any_P" distance 61.349855 [deg]
Arrival time for wavename "any_P" distance 65.567308 [deg]
Arrival time for wavename "any_P" distance 54.143590 [deg]
Arrival time for wavename "any_P" distance 73.109289 [deg]
Arrival time for wavename "any_P" distance 64.946397 [deg]
Arrival time for wavename "any_P" distance 89.437223 [deg]
Arrival time for wavename "any_P" distance 67.439618 [deg]
Arrival time for wavename "any_P" distance 38.167185 [deg]
Arrival time for wavename "any_P" distance 84.816714 [deg]
Arrival time for wavename "any_P" distance 43.633450 [deg]
Arrival time for wavename "any_P" distance 60.417579 [deg]
Arrival time for wavename "any_P" distance 38.986606 [deg]
Arrival time for wavename "any_P" distance 68.244425 [deg]
Arrival time for wavename "any_P" distance 64.054778 [deg]
Arrival time for wavename "any_P" distance 82.062195 [deg]
Arrival time for wavename "any_P" distance 36.698309 [deg]
Arrival time for wavename "any_P" distance 58.931762 [deg]
Arrival time for wavename "any_P" distance 58.456908 [deg]
Arrival time for wavename "any_P" distance 79.971426 [deg]
Arrival time for wavename "any_P" distance 74.925103 [deg]
Arrival time for wavename "any_P" distance 52.564826 [deg]
Arrival time for wavename "any_P" distance 44.192670 [deg]
Arrival time for wavename "any_P" distance 71.681441 [deg]
Arrival time for wavename "any_P" distance 48.335702 [deg]
Arrival time for wavename "any_P" distance 44.334011 [deg]
Arrival time for wavename "any_P" distance 36.168619 [deg]
Arrival time for wavename "any_P" distance 79.664513 [deg]
Arrival time for wavename "any_P" distance 31.995457 [deg]
Arrival time for wavename "any_P" distance 86.199701 [deg]
Arrival time for wavename "any_P" distance 46.405043 [deg]
Arrival time for wavename "any_P" distance 34.203018 [deg]
Arrival time for wavename "any_P" distance 31.451507 [deg]
Arrival time for wavename "any_P" distance 72.303754 [deg]
Arrival time for wavename "any_P" distance 81.349799 [deg]
Arrival time for wavename "any_P" distance 74.397519 [deg]
Pre-processing data ...
Filtering, tapering, chopping ... trace_samples: 1641
taper times: --- !pf.CosTaper

heart - DEBUG trace: Trace (IU, TSUM, filt, Z)
timerange: 2009-04-06 01:32:44 - 2009-04-06 02:27:24
delta t: 2

heart - DEBUG chopped trace: Trace (IU, TSUM, filt, Z)
timerange: 2009-04-06 01:42:54 - 2009-04-06 01:43:42
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1716
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982197.0
b: 1238982202.0
c: 1238982252.0
d: 1238982267.0

heart - DEBUG trace: Trace (IU, RCBR, filt, Z)
timerange: 2009-04-06 01:33:12 - 2009-04-06 02:30:22
delta t: 2

heart - DEBUG chopped trace: Trace (IU, RCBR, filt, Z)
timerange: 2009-04-06 01:43:22 - 2009-04-06 01:44:10
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1508
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982119.0
b: 1238982124.0
c: 1238982174.0
d: 1238982189.0

heart - DEBUG trace: Trace (CN, RES, filt, Z)
timerange: 2009-04-06 01:31:52 - 2009-04-06 02:22:06
delta t: 2

heart - DEBUG chopped trace: Trace (CN, RES, filt, Z)
timerange: 2009-04-06 01:42:04 - 2009-04-06 01:42:52
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1853
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982245.0
b: 1238982250.0
c: 1238982300.0
d: 1238982315.0

heart - DEBUG trace: Trace (IC, KMI, filt, Z)
timerange: 2009-04-06 01:34:00 - 2009-04-06 02:35:44
delta t: 2

heart - DEBUG chopped trace: Trace (IC, KMI, filt, Z)
timerange: 2009-04-06 01:44:10 - 2009-04-06 01:44:58
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1702
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982193.0
b: 1238982198.0
c: 1238982248.0
d: 1238982263.0

heart - DEBUG trace: Trace (IU, SSPA, filt, Z)
timerange: 2009-04-06 01:33:08 - 2009-04-06 02:29:50
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SSPA, filt, Z)
timerange: 2009-04-06 01:43:18 - 2009-04-06 01:44:06
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 2169
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982331.0
b: 1238982336.0
c: 1238982386.0
d: 1238982401.0

heart - DEBUG trace: Trace (MS, BTDF, filt, Z)
timerange: 2009-04-06 01:35:26 - 2009-04-06 02:47:42
delta t: 2

heart - DEBUG chopped trace: Trace (MS, BTDF, filt, Z)
timerange: 2009-04-06 01:45:36 - 2009-04-06 01:46:24
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1748
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982209.0
b: 1238982214.0
c: 1238982264.0
d: 1238982279.0

heart - DEBUG trace: Trace (CU, ANWB, filt, Z)
timerange: 2009-04-06 01:33:24 - 2009-04-06 02:31:38
delta t: 2

heart - DEBUG chopped trace: Trace (CU, ANWB, filt, Z)
timerange: 2009-04-06 01:43:34 - 2009-04-06 01:44:22
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1229
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981993.0
b: 1238981998.0
c: 1238982048.0
d: 1238982063.0

heart - DEBUG trace: Trace (G, MBO, filt, Z)
timerange: 2009-04-06 01:29:48 - 2009-04-06 02:10:44
delta t: 2

heart - DEBUG chopped trace: Trace (G, MBO, filt, Z)
timerange: 2009-04-06 01:39:58 - 2009-04-06 01:40:46
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 2078
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982309.0
b: 1238982314.0
c: 1238982364.0
d: 1238982379.0

heart - DEBUG trace: Trace (IU, TEIG, filt, Z)
timerange: 2009-04-06 01:35:02 - 2009-04-06 02:44:16
delta t: 2

heart - DEBUG chopped trace: Trace (IU, TEIG, filt, Z)
timerange: 2009-04-06 01:45:14 - 2009-04-06 01:46:02
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1324
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982039.0
b: 1238982044.0
c: 1238982094.0
d: 1238982109.0

heart - DEBUG trace: Trace (IU, KBL, filt, Z)
timerange: 2009-04-06 01:30:32 - 2009-04-06 02:14:38
delta t: 2

heart - DEBUG chopped trace: Trace (IU, KBL, filt, Z)
timerange: 2009-04-06 01:40:44 - 2009-04-06 01:41:32
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1624
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982163.0
b: 1238982168.0
c: 1238982218.0
d: 1238982233.0

heart - DEBUG trace: Trace (II, MSEY, filt, Z)
timerange: 2009-04-06 01:32:38 - 2009-04-06 02:26:44
delta t: 2

heart - DEBUG chopped trace: Trace (II, MSEY, filt, Z)
timerange: 2009-04-06 01:42:48 - 2009-04-06 01:43:36
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1244
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982001.0
b: 1238982006.0
c: 1238982056.0
d: 1238982071.0

heart - DEBUG trace: Trace (GT, DBIC, filt, Z)
timerange: 2009-04-06 01:29:54 - 2009-04-06 02:11:20
delta t: 2

heart - DEBUG chopped trace: Trace (GT, DBIC, filt, Z)
timerange: 2009-04-06 01:40:06 - 2009-04-06 01:40:54
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1763
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982215.0
b: 1238982220.0
c: 1238982270.0
d: 1238982285.0

heart - DEBUG trace: Trace (IC, HIA, filt, Z)
timerange: 2009-04-06 01:33:30 - 2009-04-06 02:32:14
delta t: 2

heart - DEBUG chopped trace: Trace (IC, HIA, filt, Z)
timerange: 2009-04-06 01:43:40 - 2009-04-06 01:44:28
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1686
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982187.0
b: 1238982192.0
c: 1238982242.0
d: 1238982257.0

heart - DEBUG trace: Trace (IU, YAK, filt, Z)
timerange: 2009-04-06 01:33:02 - 2009-04-06 02:29:12
delta t: 2

heart - DEBUG chopped trace: Trace (IU, YAK, filt, Z)
timerange: 2009-04-06 01:43:12 - 2009-04-06 01:44:00
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 2026
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982295.0
b: 1238982300.0
c: 1238982350.0
d: 1238982365.0

heart - DEBUG trace: Trace (IC, QIZ, filt, Z)
timerange: 2009-04-06 01:34:48 - 2009-04-06 02:42:18
delta t: 2

heart - DEBUG chopped trace: Trace (IC, QIZ, filt, Z)
timerange: 2009-04-06 01:45:00 - 2009-04-06 01:45:48
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1201
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981981.0
b: 1238981986.0
c: 1238982036.0
d: 1238982051.0

heart - DEBUG trace: Trace (GE, KBS, filt, Z)
timerange: 2009-04-06 01:29:34 - 2009-04-06 02:09:34
delta t: 2

heart - DEBUG chopped trace: Trace (GE, KBS, filt, Z)
timerange: 2009-04-06 01:39:46 - 2009-04-06 01:40:34
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1598
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982153.0
b: 1238982158.0
c: 1238982208.0
d: 1238982223.0

heart - DEBUG trace: Trace (IU, LSZ, filt, Z)
timerange: 2009-04-06 01:32:28 - 2009-04-06 02:25:42
delta t: 2

heart - DEBUG chopped trace: Trace (IU, LSZ, filt, Z)
timerange: 2009-04-06 01:42:38 - 2009-04-06 01:43:26
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1586
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982149.0
b: 1238982154.0
c: 1238982204.0
d: 1238982219.0

heart - DEBUG trace: Trace (II, TLY, filt, Z)
timerange: 2009-04-06 01:32:24 - 2009-04-06 02:25:14
delta t: 2

heart - DEBUG chopped trace: Trace (II, TLY, filt, Z)
timerange: 2009-04-06 01:42:34 - 2009-04-06 01:43:22
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1984
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982283.0
b: 1238982288.0
c: 1238982338.0
d: 1238982353.0

heart - DEBUG trace: Trace (IU, INCN, filt, Z)
timerange: 2009-04-06 01:34:38 - 2009-04-06 02:40:44
delta t: 2

heart - DEBUG chopped trace: Trace (IU, INCN, filt, Z)
timerange: 2009-04-06 01:44:48 - 2009-04-06 01:45:36
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1888
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982255.0
b: 1238982260.0
c: 1238982310.0
d: 1238982325.0

heart - DEBUG trace: Trace (IU, DWPF, filt, Z)
timerange: 2009-04-06 01:34:08 - 2009-04-06 02:37:02
delta t: 2

heart - DEBUG chopped trace: Trace (IU, DWPF, filt, Z)
timerange: 2009-04-06 01:44:20 - 2009-04-06 01:45:08
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1481
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982107.0
b: 1238982112.0
c: 1238982162.0
d: 1238982177.0

heart - DEBUG trace: Trace (IC, WMQ, filt, Z)
timerange: 2009-04-06 01:31:42 - 2009-04-06 02:21:02
delta t: 2

heart - DEBUG chopped trace: Trace (IC, WMQ, filt, Z)
timerange: 2009-04-06 01:41:52 - 2009-04-06 01:42:40
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1334
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982043.0
b: 1238982048.0
c: 1238982098.0
d: 1238982113.0

heart - DEBUG trace: Trace (II, AAK, filt, Z)
timerange: 2009-04-06 01:30:36 - 2009-04-06 02:15:02
delta t: 2

heart - DEBUG chopped trace: Trace (II, AAK, filt, Z)
timerange: 2009-04-06 01:40:48 - 2009-04-06 01:41:36
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1827
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982235.0
b: 1238982240.0
c: 1238982290.0
d: 1238982305.0

heart - DEBUG trace: Trace (IC, XAN, filt, Z)
timerange: 2009-04-06 01:33:50 - 2009-04-06 02:34:42
delta t: 2

heart - DEBUG chopped trace: Trace (IC, XAN, filt, Z)
timerange: 2009-04-06 01:44:00 - 2009-04-06 01:44:48
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1408
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982075.0
b: 1238982080.0
c: 1238982130.0
d: 1238982145.0

heart - DEBUG trace: Trace (IU, KMBO, filt, Z)
timerange: 2009-04-06 01:31:10 - 2009-04-06 02:18:04
delta t: 2

heart - DEBUG chopped trace: Trace (IU, KMBO, filt, Z)
timerange: 2009-04-06 01:41:20 - 2009-04-06 01:42:08
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1335
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982045.0
b: 1238982050.0
c: 1238982100.0
d: 1238982115.0

heart - DEBUG trace: Trace (II, KURK, filt, Z)
timerange: 2009-04-06 01:30:38 - 2009-04-06 02:15:06
delta t: 2

heart - DEBUG chopped trace: Trace (II, KURK, filt, Z)
timerange: 2009-04-06 01:40:50 - 2009-04-06 01:41:38
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1195
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981977.0
b: 1238981982.0
c: 1238982032.0
d: 1238982047.0

heart - DEBUG trace: Trace (ZF, ADYE, filt, Z)
timerange: 2009-04-06 01:29:32 - 2009-04-06 02:09:20
delta t: 2

heart - DEBUG chopped trace: Trace (ZF, ADYE, filt, Z)
timerange: 2009-04-06 01:39:42 - 2009-04-06 01:40:30
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1978
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982281.0
b: 1238982286.0
c: 1238982336.0
d: 1238982351.0

heart - DEBUG trace: Trace (IU, SDV, filt, Z)
timerange: 2009-04-06 01:34:36 - 2009-04-06 02:40:30
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SDV, filt, Z)
timerange: 2009-04-06 01:44:46 - 2009-04-06 01:45:34
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1120
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981941.0
b: 1238981946.0
c: 1238981996.0
d: 1238982011.0

heart - DEBUG trace: Trace (II, ARU, filt, Z)
timerange: 2009-04-06 01:28:54 - 2009-04-06 02:06:12
delta t: 2

heart - DEBUG chopped trace: Trace (II, ARU, filt, Z)
timerange: 2009-04-06 01:39:06 - 2009-04-06 01:39:54
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 2105
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982315.0
b: 1238982320.0
c: 1238982370.0
d: 1238982385.0

heart - DEBUG trace: Trace (IU, SAML, filt, Z)
timerange: 2009-04-06 01:35:10 - 2009-04-06 02:45:18
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SAML, filt, Z)
timerange: 2009-04-06 01:45:20 - 2009-04-06 01:46:08
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1370
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982061.0
b: 1238982066.0
c: 1238982116.0
d: 1238982131.0

heart - DEBUG trace: Trace (II, ALE, filt, Z)
timerange: 2009-04-06 01:30:54 - 2009-04-06 02:16:32
delta t: 2

heart - DEBUG chopped trace: Trace (II, ALE, filt, Z)
timerange: 2009-04-06 01:41:06 - 2009-04-06 01:41:54
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1159
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981959.0
b: 1238981964.0
c: 1238982014.0
d: 1238982029.0

heart - DEBUG trace: Trace (II, ABKT, filt, Z)
timerange: 2009-04-06 01:29:14 - 2009-04-06 02:07:50
delta t: 2

heart - DEBUG chopped trace: Trace (II, ABKT, filt, Z)
timerange: 2009-04-06 01:39:24 - 2009-04-06 01:40:12
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1110
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981935.0
b: 1238981940.0
c: 1238981990.0
d: 1238982005.0

heart - DEBUG trace: Trace (KZ, AKTO, filt, Z)
timerange: 2009-04-06 01:28:50 - 2009-04-06 02:05:48
delta t: 2

heart - DEBUG chopped trace: Trace (KZ, AKTO, filt, Z)
timerange: 2009-04-06 01:39:00 - 2009-04-06 01:39:48
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1838
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982239.0
b: 1238982244.0
c: 1238982294.0
d: 1238982309.0

heart - DEBUG trace: Trace (IC, BJT, filt, Z)
timerange: 2009-04-06 01:33:54 - 2009-04-06 02:35:08
delta t: 2

heart - DEBUG chopped trace: Trace (IC, BJT, filt, Z)
timerange: 2009-04-06 01:44:04 - 2009-04-06 01:44:52
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 2012
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982291.0
b: 1238982296.0
c: 1238982346.0
d: 1238982361.0

heart - DEBUG trace: Trace (IC, SSE, filt, Z)
timerange: 2009-04-06 01:34:44 - 2009-04-06 02:41:46
delta t: 2

heart - DEBUG chopped trace: Trace (IC, SSE, filt, Z)
timerange: 2009-04-06 01:44:56 - 2009-04-06 01:45:44
delta t: 2

heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1879
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982251.0
b: 1238982256.0
c: 1238982306.0
d: 1238982321.0

heart - DEBUG trace: Trace (IU, CHTO, filt, Z)
timerange: 2009-04-06 01:34:06 - 2009-04-06 02:36:42
delta t: 2

heart - DEBUG chopped trace: Trace (IU, CHTO, filt, Z)
timerange: 2009-04-06 01:44:16 - 2009-04-06 01:45:04
delta t: 2

heart - DEBUG Returning chopped traces ...
seismic - INFO Initializing synthesizer for "any_P_0"
heart - DEBUG Synthetics generation time: 0.017754
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982169.0
b: 1238982174.0
c: 1238982224.0
d: 1238982239.0

heart - DEBUG trace: Trace (IU, TSUM, 0, Z)
timerange: 2009-04-06 01:42:38 - 2009-04-06 01:44:18
delta t: 2

heart - DEBUG chopped trace: Trace (IU, TSUM, 0, Z)
timerange: 2009-04-06 01:42:54 - 2009-04-06 01:43:42
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982197.0
b: 1238982202.0
c: 1238982252.0
d: 1238982267.0

heart - DEBUG trace: Trace (IU, RCBR, 0, Z)
timerange: 2009-04-06 01:43:06 - 2009-04-06 01:44:46
delta t: 2

heart - DEBUG chopped trace: Trace (IU, RCBR, 0, Z)
timerange: 2009-04-06 01:43:22 - 2009-04-06 01:44:10
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982119.0
b: 1238982124.0
c: 1238982174.0
d: 1238982189.0

heart - DEBUG trace: Trace (CN, RES, 0, Z)
timerange: 2009-04-06 01:41:48 - 2009-04-06 01:43:28
delta t: 2

heart - DEBUG chopped trace: Trace (CN, RES, 0, Z)
timerange: 2009-04-06 01:42:04 - 2009-04-06 01:42:52
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982245.0
b: 1238982250.0
c: 1238982300.0
d: 1238982315.0

heart - DEBUG trace: Trace (IC, KMI, 0, Z)
timerange: 2009-04-06 01:43:54 - 2009-04-06 01:45:34
delta t: 2

heart - DEBUG chopped trace: Trace (IC, KMI, 0, Z)
timerange: 2009-04-06 01:44:10 - 2009-04-06 01:44:58
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982193.0
b: 1238982198.0
c: 1238982248.0
d: 1238982263.0

heart - DEBUG trace: Trace (IU, SSPA, 0, Z)
timerange: 2009-04-06 01:43:02 - 2009-04-06 01:44:42
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SSPA, 0, Z)
timerange: 2009-04-06 01:43:18 - 2009-04-06 01:44:06
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982331.0
b: 1238982336.0
c: 1238982386.0
d: 1238982401.0

heart - DEBUG trace: Trace (MS, BTDF, 0, Z)
timerange: 2009-04-06 01:45:20 - 2009-04-06 01:47:00
delta t: 2

heart - DEBUG chopped trace: Trace (MS, BTDF, 0, Z)
timerange: 2009-04-06 01:45:36 - 2009-04-06 01:46:24
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982209.0
b: 1238982214.0
c: 1238982264.0
d: 1238982279.0

heart - DEBUG trace: Trace (CU, ANWB, 0, Z)
timerange: 2009-04-06 01:43:18 - 2009-04-06 01:44:58
delta t: 2

heart - DEBUG chopped trace: Trace (CU, ANWB, 0, Z)
timerange: 2009-04-06 01:43:34 - 2009-04-06 01:44:22
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981993.0
b: 1238981998.0
c: 1238982048.0
d: 1238982063.0

heart - DEBUG trace: Trace (G, MBO, 0, Z)
timerange: 2009-04-06 01:39:42 - 2009-04-06 01:41:22
delta t: 2

heart - DEBUG chopped trace: Trace (G, MBO, 0, Z)
timerange: 2009-04-06 01:39:58 - 2009-04-06 01:40:46
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982309.0
b: 1238982314.0
c: 1238982364.0
d: 1238982379.0

heart - DEBUG trace: Trace (IU, TEIG, 0, Z)
timerange: 2009-04-06 01:44:58 - 2009-04-06 01:46:38
delta t: 2

heart - DEBUG chopped trace: Trace (IU, TEIG, 0, Z)
timerange: 2009-04-06 01:45:14 - 2009-04-06 01:46:02
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982039.0
b: 1238982044.0
c: 1238982094.0
d: 1238982109.0

heart - DEBUG trace: Trace (IU, KBL, 0, Z)
timerange: 2009-04-06 01:40:28 - 2009-04-06 01:42:08
delta t: 2

heart - DEBUG chopped trace: Trace (IU, KBL, 0, Z)
timerange: 2009-04-06 01:40:44 - 2009-04-06 01:41:32
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982163.0
b: 1238982168.0
c: 1238982218.0
d: 1238982233.0

heart - DEBUG trace: Trace (II, MSEY, 0, Z)
timerange: 2009-04-06 01:42:32 - 2009-04-06 01:44:12
delta t: 2

heart - DEBUG chopped trace: Trace (II, MSEY, 0, Z)
timerange: 2009-04-06 01:42:48 - 2009-04-06 01:43:36
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982001.0
b: 1238982006.0
c: 1238982056.0
d: 1238982071.0

heart - DEBUG trace: Trace (GT, DBIC, 0, Z)
timerange: 2009-04-06 01:39:50 - 2009-04-06 01:41:30
delta t: 2

heart - DEBUG chopped trace: Trace (GT, DBIC, 0, Z)
timerange: 2009-04-06 01:40:06 - 2009-04-06 01:40:54
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982215.0
b: 1238982220.0
c: 1238982270.0
d: 1238982285.0

heart - DEBUG trace: Trace (IC, HIA, 0, Z)
timerange: 2009-04-06 01:43:24 - 2009-04-06 01:45:04
delta t: 2

heart - DEBUG chopped trace: Trace (IC, HIA, 0, Z)
timerange: 2009-04-06 01:43:40 - 2009-04-06 01:44:28
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982187.0
b: 1238982192.0
c: 1238982242.0
d: 1238982257.0

heart - DEBUG trace: Trace (IU, YAK, 0, Z)
timerange: 2009-04-06 01:42:56 - 2009-04-06 01:44:36
delta t: 2

heart - DEBUG chopped trace: Trace (IU, YAK, 0, Z)
timerange: 2009-04-06 01:43:12 - 2009-04-06 01:44:00
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982295.0
b: 1238982300.0
c: 1238982350.0
d: 1238982365.0

heart - DEBUG trace: Trace (IC, QIZ, 0, Z)
timerange: 2009-04-06 01:44:44 - 2009-04-06 01:46:24
delta t: 2

heart - DEBUG chopped trace: Trace (IC, QIZ, 0, Z)
timerange: 2009-04-06 01:45:00 - 2009-04-06 01:45:48
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981981.0
b: 1238981986.0
c: 1238982036.0
d: 1238982051.0

heart - DEBUG trace: Trace (GE, KBS, 0, Z)
timerange: 2009-04-06 01:39:30 - 2009-04-06 01:41:10
delta t: 2

heart - DEBUG chopped trace: Trace (GE, KBS, 0, Z)
timerange: 2009-04-06 01:39:46 - 2009-04-06 01:40:34
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982153.0
b: 1238982158.0
c: 1238982208.0
d: 1238982223.0

heart - DEBUG trace: Trace (IU, LSZ, 0, Z)
timerange: 2009-04-06 01:42:22 - 2009-04-06 01:44:02
delta t: 2

heart - DEBUG chopped trace: Trace (IU, LSZ, 0, Z)
timerange: 2009-04-06 01:42:38 - 2009-04-06 01:43:26
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982149.0
b: 1238982154.0
c: 1238982204.0
d: 1238982219.0

heart - DEBUG trace: Trace (II, TLY, 0, Z)
timerange: 2009-04-06 01:42:18 - 2009-04-06 01:43:58
delta t: 2

heart - DEBUG chopped trace: Trace (II, TLY, 0, Z)
timerange: 2009-04-06 01:42:34 - 2009-04-06 01:43:22
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982283.0
b: 1238982288.0
c: 1238982338.0
d: 1238982353.0

heart - DEBUG trace: Trace (IU, INCN, 0, Z)
timerange: 2009-04-06 01:44:32 - 2009-04-06 01:46:12
delta t: 2

heart - DEBUG chopped trace: Trace (IU, INCN, 0, Z)
timerange: 2009-04-06 01:44:48 - 2009-04-06 01:45:36
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982255.0
b: 1238982260.0
c: 1238982310.0
d: 1238982325.0

heart - DEBUG trace: Trace (IU, DWPF, 0, Z)
timerange: 2009-04-06 01:44:04 - 2009-04-06 01:45:44
delta t: 2

heart - DEBUG chopped trace: Trace (IU, DWPF, 0, Z)
timerange: 2009-04-06 01:44:20 - 2009-04-06 01:45:08
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982107.0
b: 1238982112.0
c: 1238982162.0
d: 1238982177.0

heart - DEBUG trace: Trace (IC, WMQ, 0, Z)
timerange: 2009-04-06 01:41:36 - 2009-04-06 01:43:16
delta t: 2

heart - DEBUG chopped trace: Trace (IC, WMQ, 0, Z)
timerange: 2009-04-06 01:41:52 - 2009-04-06 01:42:40
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982043.0
b: 1238982048.0
c: 1238982098.0
d: 1238982113.0

heart - DEBUG trace: Trace (II, AAK, 0, Z)
timerange: 2009-04-06 01:40:32 - 2009-04-06 01:42:12
delta t: 2

heart - DEBUG chopped trace: Trace (II, AAK, 0, Z)
timerange: 2009-04-06 01:40:48 - 2009-04-06 01:41:36
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982235.0
b: 1238982240.0
c: 1238982290.0
d: 1238982305.0

heart - DEBUG trace: Trace (IC, XAN, 0, Z)
timerange: 2009-04-06 01:43:44 - 2009-04-06 01:45:24
delta t: 2

heart - DEBUG chopped trace: Trace (IC, XAN, 0, Z)
timerange: 2009-04-06 01:44:00 - 2009-04-06 01:44:48
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982075.0
b: 1238982080.0
c: 1238982130.0
d: 1238982145.0

heart - DEBUG trace: Trace (IU, KMBO, 0, Z)
timerange: 2009-04-06 01:41:04 - 2009-04-06 01:42:44
delta t: 2

heart - DEBUG chopped trace: Trace (IU, KMBO, 0, Z)
timerange: 2009-04-06 01:41:20 - 2009-04-06 01:42:08
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982045.0
b: 1238982050.0
c: 1238982100.0
d: 1238982115.0

heart - DEBUG trace: Trace (II, KURK, 0, Z)
timerange: 2009-04-06 01:40:34 - 2009-04-06 01:42:14
delta t: 2

heart - DEBUG chopped trace: Trace (II, KURK, 0, Z)
timerange: 2009-04-06 01:40:50 - 2009-04-06 01:41:38
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981977.0
b: 1238981982.0
c: 1238982032.0
d: 1238982047.0

heart - DEBUG trace: Trace (ZF, ADYE, 0, Z)
timerange: 2009-04-06 01:39:26 - 2009-04-06 01:41:06
delta t: 2

heart - DEBUG chopped trace: Trace (ZF, ADYE, 0, Z)
timerange: 2009-04-06 01:39:42 - 2009-04-06 01:40:30
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982281.0
b: 1238982286.0
c: 1238982336.0
d: 1238982351.0

heart - DEBUG trace: Trace (IU, SDV, 0, Z)
timerange: 2009-04-06 01:44:30 - 2009-04-06 01:46:10
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SDV, 0, Z)
timerange: 2009-04-06 01:44:46 - 2009-04-06 01:45:34
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981941.0
b: 1238981946.0
c: 1238981996.0
d: 1238982011.0

heart - DEBUG trace: Trace (II, ARU, 0, Z)
timerange: 2009-04-06 01:38:50 - 2009-04-06 01:40:30
delta t: 2

heart - DEBUG chopped trace: Trace (II, ARU, 0, Z)
timerange: 2009-04-06 01:39:06 - 2009-04-06 01:39:54
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982315.0
b: 1238982320.0
c: 1238982370.0
d: 1238982385.0

heart - DEBUG trace: Trace (IU, SAML, 0, Z)
timerange: 2009-04-06 01:45:04 - 2009-04-06 01:46:44
delta t: 2

heart - DEBUG chopped trace: Trace (IU, SAML, 0, Z)
timerange: 2009-04-06 01:45:20 - 2009-04-06 01:46:08
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982061.0
b: 1238982066.0
c: 1238982116.0
d: 1238982131.0

heart - DEBUG trace: Trace (II, ALE, 0, Z)
timerange: 2009-04-06 01:40:50 - 2009-04-06 01:42:30
delta t: 2

heart - DEBUG chopped trace: Trace (II, ALE, 0, Z)
timerange: 2009-04-06 01:41:06 - 2009-04-06 01:41:54
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981959.0
b: 1238981964.0
c: 1238982014.0
d: 1238982029.0

heart - DEBUG trace: Trace (II, ABKT, 0, Z)
timerange: 2009-04-06 01:39:08 - 2009-04-06 01:40:48
delta t: 2

heart - DEBUG chopped trace: Trace (II, ABKT, 0, Z)
timerange: 2009-04-06 01:39:24 - 2009-04-06 01:40:12
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238981935.0
b: 1238981940.0
c: 1238981990.0
d: 1238982005.0

heart - DEBUG trace: Trace (KZ, AKTO, 0, Z)
timerange: 2009-04-06 01:38:44 - 2009-04-06 01:40:24
delta t: 2

heart - DEBUG chopped trace: Trace (KZ, AKTO, 0, Z)
timerange: 2009-04-06 01:39:00 - 2009-04-06 01:39:48
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982239.0
b: 1238982244.0
c: 1238982294.0
d: 1238982309.0

heart - DEBUG trace: Trace (IC, BJT, 0, Z)
timerange: 2009-04-06 01:43:48 - 2009-04-06 01:45:28
delta t: 2

heart - DEBUG chopped trace: Trace (IC, BJT, 0, Z)
timerange: 2009-04-06 01:44:04 - 2009-04-06 01:44:52
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982291.0
b: 1238982296.0
c: 1238982346.0
d: 1238982361.0

heart - DEBUG trace: Trace (IC, SSE, 0, Z)
timerange: 2009-04-06 01:44:40 - 2009-04-06 01:46:20
delta t: 2

heart - DEBUG chopped trace: Trace (IC, SSE, 0, Z)
timerange: 2009-04-06 01:44:56 - 2009-04-06 01:45:44
delta t: 2

heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982251.0
b: 1238982256.0
c: 1238982306.0
d: 1238982321.0

heart - DEBUG trace: Trace (IU, CHTO, 0, Z)
timerange: 2009-04-06 01:44:00 - 2009-04-06 01:45:40
delta t: 2

heart - DEBUG chopped trace: Trace (IU, CHTO, 0, Z)
timerange: 2009-04-06 01:44:16 - 2009-04-06 01:45:04
delta t: 2

heart - DEBUG Post-process time 0.023978
heart - DEBUG Stack traces time 0.000000
heart - DEBUG Returning...
seismic - DEBUG Teleseismic forward model on test model takes: 0.827361
models - INFO Model building was successful!

models - INFO Using "bin" backend to store samples!
models - INFO ... Initiate Metropolis for Parallel Tempering...
proposal_distribution: Normal, tune_interval=10, n_chains=8

models.base - INFO ... Starting Parallel Tempering ...

distributed - INFO MPI is installed!
distributed - INFO Detected python3
distributed - INFO Done initialising mpi runner
distributed - INFO sampler directory: /home/lab/bayesiano/beat/beat/sampler/pt.py
distributed - DEBUG Running mpi with arguments: -n 9 python3 /home/lab/bayesiano/beat/beat/sampler/pt.py debug /home/lab/Laquila_dc
Traceback (most recent call last):
File "/home/lab/anaconda3/envs/beat/bin/beat", line 11, in
load_entry_point('beat==1.0', 'console_scripts', 'beat')()
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 1686, in main
globals()'command_' + command
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/apps/beat.py", line 775, in command_sample
sample(step, problem)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/models/base.py", line 196, in sample
record_worker_chains=pa.record_worker_chains)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/pt.py", line 836, in pt_sample
project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 226, in run_mpi_sampler
samplerdir, n_jobs=n_jobs, loglevel=loglevel, project_dir=project_dir)
File "/home/lab/anaconda3/envs/beat/lib/python3.6/site-packages/beat-1.0-py3.6-linux-x86_64.egg/beat/sampler/distributed.py", line 158, in run
in the directory %s'''.lstrip() % (program, self.tempdir))
beat.sampler.distributed.MPIError: mpiexec has been invoked as "mpiexec"
in the directory /tmp/mpiexec-7bw3hmzo
distributed - WARNING not removing temporary directory: /tmp/mpiexec-7bw3hmzo

------------------------------------------------------------*-

3) Running:

mpiexec -n 9 python3 /home/lab/bayesiano/beat/beat/sampler/pt.py debug /home/lab/Laquila_dc

in the temporary directory: /tmp/mpiexec-7bw3hmzo

I got:

(beat) lab@lab:/tmp/mpiexec-7bw3hmzo$ mpiexec -n 9 python3 /home/lab/bayesiano/beat/beat/sampler/pt.py debug /home/lab/Laquila_dc

There are not enough slots available in the system to satisfy the 9
slots that were requested by the application:

python3

Either request fewer slots for your application, or make more slots
available for use.

A "slot" is the Open MPI term for an allocatable unit where we can
launch a process. The number of slots available are defined by the
environment in which Open MPI processes are run:

Hostfile, via "slots=N" clauses (N defaults to number of
processor cores if not provided)
The --host command line parameter, via a ":N" suffix on the
hostname (N defaults to 1 if not provided)
Resource manager (e.g., SLURM, PBS/Torque, LSF, etc.)
If none of a hostfile, the --host command line parameter, or an
RM is present, Open MPI defaults to the number of processor cores

In all the above cases, if you want Open MPI to default to the number
of hardware threads instead of the number of processor cores, use the
--use-hwthread-cpus option.

Alternatively, you can use the --oversubscribe option to ignore the
number of available slots when deciding the number of processes to
launch.

------------------------------------------------------------*-

4) running the line:

mpiexec -n 9 python3 /home/lab/bayesiano/beat/beat/sampler/pt.py debug /home/lab/Laquila_dc --use-hwthread-cpus

Got:

(beat) lab@lab:/tmp/mpiexec-7bw3hmzo$ mpiexec -n 9 python3 /home/lab/bayesiano/beat/beat/sampler/pt.py --use-hwthread-cpus

There are not enough slots available in the system to satisfy the 9
slots that were requested by the application:

python3

Either request fewer slots for your application, or make more slots
available for use.

A "slot" is the Open MPI term for an allocatable unit where we can
launch a process. The number of slots available are defined by the
environment in which Open MPI processes are run:

Hostfile, via "slots=N" clauses (N defaults to number of
processor cores if not provided)
The --host command line parameter, via a ":N" suffix on the
hostname (N defaults to 1 if not provided)
Resource manager (e.g., SLURM, PBS/Torque, LSF, etc.)
If none of a hostfile, the --host command line parameter, or an
RM is present, Open MPI defaults to the number of processor cores

In all the above cases, if you want Open MPI to default to the number
of hardware threads instead of the number of processor cores, use the
--use-hwthread-cpus option.

Alternatively, you can use the --oversubscribe option to ignore the
number of available slots when deciding the number of processes to
launch.

------------------------------------------------------------*-

5) Now, Running the second example with n_chains=5, I got:

(beat) lab@lab:~$ beat sample Laquila_dc --loglevel=debug
config - INFO All hierarchicals ok!
config - INFO All hyperparameters ok!
config - INFO All priors ok!
models - INFO ... Initialising Geometry Optimizer ...

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

seismic - DEBUG Setting up seismic structure ...

utility - DEBUG Valid distance range: [30.000000, 90.000000]!
utility - DEBUG Distance of station TSUM: 61.349855 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station RCBR: 65.567308 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station RES: 54.143590 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KMI: 73.109289 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SSPA: 64.946397 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station BTDF: 89.437223 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ANWB: 67.439618 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station MBO: 38.167185 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station TEIG: 84.816714 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KBL: 43.633450 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station MSEY: 60.417579 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station DBIC: 38.986606 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station HIA: 68.244425 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station YAK: 64.054778 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station QIZ: 82.062195 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KBS: 36.698309 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station LSZ: 58.931762 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station TLY: 58.456908 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station INCN: 79.971426 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station DWPF: 74.925103 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station WMQ: 52.564826 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station AAK: 44.192670 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station XAN: 71.681441 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KMBO: 48.335702 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station KURK: 44.334011 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ADYE: 36.168619 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SDV: 79.664513 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ARU: 31.995457 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SAML: 86.199701 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station ALE: 46.405043 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station OTAV: 91.209244 [deg]
utility - DEBUG Outside defined distance range!
utility - DEBUG Distance of station ABKT: 34.203018 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station AKTO: 31.451507 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station BJT: 72.303754 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station SSE: 81.349799 [deg]
utility - DEBUG Inside defined distance range!
utility - DEBUG Distance of station CHTO: 74.397519 [deg]
utility - DEBUG Inside defined distance range!
heart - INFO Consistent number of datasets and targets in any_P_0 wavemap!
heart - INFO Number of seismic datasets for wavemap: any_P_0: 35
models - INFO ... Building model ...

models - DEBUG Optimization for 1 sources
models - INFO Optimization for 35 hyperparameters in total!
seismic - INFO Seismic optimization on:
depth, dip, duration, east_shift, magnitude, north_shift, rake, strike, time
seismic - INFO Retrieving seismic data-covariances with structure "variance" for any_P_0 ...
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 61.349855 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 65.567308 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 54.143590 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 73.109289 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 64.946397 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 89.437223 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 67.439618 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 38.167185 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 84.816714 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 43.633450 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 60.417579 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 38.986606 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 68.244425 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 64.054778 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 82.062195 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 36.698309 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 58.931762 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 58.456908 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 79.971426 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 74.925103 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 52.564826 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 44.192670 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 71.681441 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 48.335702 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 44.334011 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 36.168619 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 79.664513 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 31.995457 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 86.199701 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 46.405043 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 34.203018 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 31.451507 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 72.303754 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 81.349799 [deg]
heart - DEBUG Wavename not specified using first tabulated phase! any_P
heart - DEBUG Arrival time for wavename "any_P" distance 74.397519 [deg]
seismic - INFO Preparing data of "any_P_0" for optimization
heart - DEBUG Arrival time for wavename "any_P" distance 61.349855 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 65.567308 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 54.143590 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 73.109289 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 64.946397 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 89.437223 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 67.439618 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 38.167185 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 84.816714 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 43.633450 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 60.417579 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 38.986606 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 68.244425 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 64.054778 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 82.062195 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 36.698309 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 58.931762 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 58.456908 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 79.971426 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 74.925103 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 52.564826 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 44.192670 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 71.681441 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 48.335702 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 44.334011 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 36.168619 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 79.664513 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 31.995457 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 86.199701 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 46.405043 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 34.203018 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 31.451507 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 72.303754 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 81.349799 [deg]
heart - DEBUG Arrival time for wavename "any_P" distance 74.397519 [deg]
heart - DEBUG Pre-processing data ...
heart - DEBUG Filtering, tapering, chopping ... trace_samples: 1641
heart - DEBUG taper times: --- !pf.CosTaper
a: 1238982169.0
b: 1238982174.0
c: 1238982224.0
d: 1238982239.0