WIP update tutorials example1

hvasbath · Sep 17, 2019 · c652a56 · c652a56
1 parent 1ce20b5
commit c652a56
Show file tree

Hide file tree

Showing 12 changed files with 169 additions and 219 deletions.
diff --git a/data/examples/FullMT/config_geometry.yaml b/data/examples/FullMT/config_geometry.yaml
@@ -1,22 +1,22 @@
 --- !beat.BEATconfig
-name: MTtest_nont
+name: FullMT
 date: dummy
 event: !pf.Event
   lat: 29.07
   lon: 34.73
   time: 1995-11-22 04:15:26.200000
-  name: 112295A
   depth: 8000.0
+  name: 112295A
   magnitude: 7.20583885303153
   region: ARAB REPUBLIC OF EGYPT
   catalog: gCMT
   moment_tensor: !pf.MomentTensor
-    mnn: -3.255e+19
-    mee: 4.616e+19
-    mdd: -1.361e+19
-    mne: 4.493e+19
-    mnd: 3.887e+19
-    med: -9.1e+17
+    mnn: -0.43283071
+    mee: 0.65741974
+    mdd: -0.22458903
+    mne: 0.63839719
+    mnd: 0.50698292
+    med: 0.02063122
     strike1: 294.0906468639183
     dip1: 77.25911284933075
     rake1: -148.53160420242887
@@ -26,26 +26,27 @@ event: !pf.Event
     moment: 7.223676349339025e+19
     magnitude: 7.20583885303153
   duration: 22.0
-project_dir: /home/vasyurhm/BEATS/MTtest_nont
+  tags: []
+project_dir: /home/vasyurhm/BEATS/FullMT
 problem_config: !beat.ProblemConfig
   mode: geometry
   source_type: MTSource
-  stf_type: HalfSinusoid
+  stf_type: Triangular
   n_sources: 1
   datatypes: [seismic]
   hyperparameters:
     h_any_P_0_Z: !beat.heart.Parameter
       name: h_any_P_0_Z
       form: Uniform
       lower: [0.0]
-      upper: [5.0]
-      testvalue: [2.0]
+      upper: [4.0]
+      testvalue: [1.0]
     h_any_P_1_T: !beat.heart.Parameter
       name: h_any_P_1_T
       form: Uniform
       lower: [0.0]
-      upper: [6.0]
-      testvalue: [3.0]
+      upper: [4.0]
+      testvalue: [1.0]
   priors:
     depth: !beat.heart.Parameter
       name: depth
@@ -56,74 +57,74 @@ problem_config: !beat.ProblemConfig
     duration: !beat.heart.Parameter
       name: duration
       form: Uniform
-      lower: [0.0]
-      upper: [10.0]
-      testvalue: [5.0]
+      lower: [5.0]
+      upper: [30.0]
+      testvalue: [22.0]
     east_shift: !beat.heart.Parameter
       name: east_shift
       form: Uniform
-      lower: [-10.0]
-      upper: [30.0]
-      testvalue: [10.0]
+      lower: [-20.0]
+      upper: [20.0]
+      testvalue: [0.0]
     magnitude: !beat.heart.Parameter
       name: magnitude
       form: Uniform
-      lower: [4.5]
-      upper: [5.1]
-      testvalue: [4.8]
+      lower: [5.0]
+      upper: [6.0]
+      testvalue: [5.5]
     mdd: !beat.heart.Parameter
       name: mdd
       form: Uniform
       lower: [-1.4142135623730951]
       upper: [1.4142135623730951]
-      testvalue: [-0.08682409]
+      testvalue: [-0.22458903]
     med: !beat.heart.Parameter
       name: med
       form: Uniform
       lower: [-1.0]
       upper: [1.0]
-      testvalue: [-0.25767963]
+      testvalue: [0.02063122]
     mee: !beat.heart.Parameter
       name: mee
       form: Uniform
       lower: [-1.4142135623730951]
       upper: [1.4142135623730951]
-      testvalue: [-0.75868967]
+      testvalue: [0.65741974]
     mnd: !beat.heart.Parameter
       name: mnd
       form: Uniform
       lower: [-1.0]
       upper: [1.0]
-      testvalue: [0.14554675]
+      testvalue: [0.50698292]
     mne: !beat.heart.Parameter
       name: mne
       form: Uniform
       lower: [-1.0]
       upper: [1.0]
-      testvalue: [0.51322155]
+      testvalue: [0.63839719]
     mnn: !beat.heart.Parameter
       name: mnn
       form: Uniform
       lower: [-1.4142135623730951]
       upper: [1.4142135623730951]
-      testvalue: [0.84551376]
+      testvalue: [-0.43283071]
     north_shift: !beat.heart.Parameter
       name: north_shift
       form: Uniform
-      lower: [0.0]
-      upper: [40.0]
-      testvalue: [20.0]
+      lower: [-20.0]
+      upper: [20.0]
+      testvalue: [0.0]
     time: !beat.heart.Parameter
       name: time
       form: Uniform
-      lower: [-6.0]
-      upper: [4.0]
-      testvalue: [-2.7]
+      lower: [-20.0]
+      upper: [0.0]
+      testvalue: [-11.0]
 seismic_config: !beat.SeismicConfig
   datadir: ./
   noise_estimator: !beat.SeismicNoiseAnalyserConfig
     structure: variance
-    pre_arrival_time: 3.0
+    pre_arrival_time: 4.0
   pre_stack_cut: true
   station_corrections: false
   waveforms:
@@ -136,32 +137,16 @@ seismic_config: !beat.SeismicConfig
       lower_corner: 0.01
       upper_corner: 0.1
       order: 3
-    distances: [0.0, 90.0]
+    distances: [0.0, 10.0]
     interpolation: multilinear
     arrival_taper: !beat.heart.ArrivalTaper
       a: -30.0
       b: -20.0
       c: 230.0
       d: 250.0
-  - !beat.WaveformFitConfig
-    include: true
-    preprocess_data: true
-    name: any_P
-    channels: [T]
-    filterer: !beat.heart.Filter
-      lower_corner: 0.01
-      upper_corner: 0.1
-      order: 3
-    distances: [0.0, 90.0]
-    interpolation: multilinear
-    arrival_taper: !beat.heart.ArrivalTaper
-      a: -30.0
-      b: -20.0
-      c: 230.0
-      d: 250.0
-  dataset_specific_residual_noise_estimation: false
+  dataset_specific_residual_noise_estimation: true
   gf_config: !beat.SeismicGFConfig
-    store_superdir: /home/vasyurhm/BEATS/MTtest
+    store_superdir: /home/vasyurhm/GF
     reference_model_idx: 0
     n_variations: [0, 1]
     earth_model_name: ak135-f-average.m
@@ -186,39 +171,40 @@ seismic_config: !beat.SeismicConfig
     reference_location: !beat.heart.ReferenceLocation
       lat: 29.07
       lon: 34.73
-      elevation: 0.0
       depth: 0.0
       station: AqabaMT
     code: qseis
-    sample_rate: 1.0
+    sample_rate: 2.0
     rm_gfs: true
 sampler_config: !beat.SamplerConfig
-  name: PT
+  name: SMC
+  backend: bin
   progressbar: true
   buffer_size: 1000
-  parameters: !beat.ParallelTemperingConfig
-    tune_interval: 10
-    proposal_dist: MultivariateNormal
+  buffer_thinning: 50
+  parameters: !beat.SMCConfig
+    tune_interval: 50
     check_bnd: true
-    rm_flag: false
-    n_samples: 50000
-    n_chains: 8
-    swap_interval: [10, 30]
-    beta_tune_interval: 1000
-    n_chains_posterior: 3
-    resample: false
-    thin: 1
-    burn: 0.6
+    rm_flag: true
+    n_jobs: 25
+    n_steps: 300
+    n_chains: 2000
+    coef_variation: 1.0
+    stage: 0
+    proposal_dist: MultivariateNormal
+    update_covariances: false
 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: 1
+    n_jobs: 4
     n_steps: 25000
     n_chains: 20
     thin: 5

diff --git a/data/examples/FullMT/seismic_data.pkl b/data/examples/FullMT/seismic_data.pkl
diff --git a/data/examples/Mogi/config_geometry.yaml b/data/examples/Mogi/config_geometry.yaml
diff --git a/docs/_static/example1/FullMT_corr_hist_max.png b/docs/_static/example1/FullMT_corr_hist_max.png
diff --git a/docs/_static/example1/FullMT_data.png b/docs/_static/example1/FullMT_data.png
diff --git a/docs/_static/example1/FullMT_data_old.png b/docs/_static/example1/FullMT_data_old.png
diff --git a/docs/_static/example1/FullMT_stage_-1_max.png b/docs/_static/example1/FullMT_stage_-1_max.png
diff --git a/docs/_static/example1/FullMT_waveforms_max.png b/docs/_static/example1/FullMT_waveforms_max.png
diff --git a/docs/_static/example1/FullMT_waveforms_ref.png b/docs/_static/example1/FullMT_waveforms_ref.png
diff --git a/docs/_static/example1/FullMT_windowed.png b/docs/_static/example1/FullMT_windowed.png