docs: update Laquila example3 RectangularSource

hvasbath · Apr 10, 2024 · 7393774 · 7393774
1 parent 9152cf5
commit 7393774
Show file tree

Hide file tree

Showing 3 changed files with 143 additions and 65 deletions.
diff --git a/beat/backend.py b/beat/backend.py
@@ -84,9 +84,8 @@ def _create_flat_names_summary(varname, shape):
     if not shape or sum(shape) == 1:
         return [varname]
 
-    labels = (num.ravel(xs).tolist() for xs in num.indices(shape))
-    labels = (map(str, [xs]) for xs in labels)
-    return ["{}{}".format(varname, "".join(idxs)) for idxs in zip(*labels)]
+    labels = num.indices(shape).ravel().tolist()
+    return ["{}[{}]".format(varname, idxs) for idxs in labels]
 
 
 def _create_shape(flat_names):

diff --git a/data/examples/Laquila/config_geometry.yaml b/data/examples/Laquila/config_geometry.yaml
@@ -4,12 +4,12 @@ date: '2009-04-06'
 event: !pf.Event
   lat: 42.29
   lon: 13.35
-  time: 2009-04-06 01:32:49.190000
-  name: 200904060132A
+  time: '2009-04-06 01:32:49.190000057'
   depth: 12000.0
+  name: '200904060132A'
   magnitude: 6.343080192483292
-  region: CENTRAL ITALY
-  catalog: gCMT
+  region: 'CENTRAL ITALY'
+  catalog: 'gCMT'
   moment_tensor: !pf.MomentTensor
     mnn: 1.43e+18
     mee: 1.87e+18
@@ -29,111 +29,175 @@ event: !pf.Event
 project_dir: /home/vasyurhm/BEATS/LaquilaJointPonlyUPDATE
 problem_config: !beat.ProblemConfig
   mode: geometry
-  source_types: [RectangularSource]
+  source_types:
+  - RectangularSource
   stf_type: HalfSinusoid
   decimation_factors:
     geodetic: 4
     seismic: 1
-  n_sources: [1]
-  datatypes: [geodetic, seismic]
+  n_sources:
+  - 1
+  datatypes:
+  - geodetic
+  - seismic
   hyperparameters:
     h_SAR: !beat.heart.Parameter
       name: h_SAR
       form: Uniform
-      lower: [-1.0]
-      upper: [5.0]
-      testvalue: [2.0]
+      lower:
+      - -1.0
+      upper:
+      - 5.0
+      testvalue:
+      - 2.0
     h_any_P_0_Z: !beat.heart.Parameter
       name: h_any_P_0_Z
       form: Uniform
-      lower: [0.0]
-      upper: [3.0]
-      testvalue: [0.5]
+      lower:
+      - 0.0
+      upper:
+      - 3.0
+      testvalue:
+      - 0.5
   priors:
     depth: !beat.heart.Parameter
       name: depth
       form: Uniform
-      lower: [0.0]
-      upper: [5.0]
-      testvalue: [2.95]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      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]
+      lower:
+      - 5.0
+      upper:
+      - 20.0
+      testvalue:
+      - 10.0
 geodetic_config: !beat.GeodeticConfig
-  datadir: /home/vasyurhm/BEATS/RectangularStatic/insar
-  names: [Laquila_dscxn, Laquila_ascxn]
-  types: [SAR]
-  calc_data_cov: true
+  types:
+    SAR: !beat.SARDatasetConfig
+      datadir: /home/vasyurhm/BEATS/RectangularStatic/insar
+      names:
+      - Laquila_dscxn
+      - Laquila_ascxn
+  noise_estimator: !beat.GeodeticNoiseAnalyserConfig
+    structure: import
+    max_dist_perc: 0.2
   interpolation: multilinear
-  fit_plane: true
+  corrections_config: !beat.GeodeticCorrectionsConfig
+    euler_poles:
+    - !beat.EulerPoleConfig
+      enabled: false
+    ramp: !beat.RampConfig
+      enabled: true
+      dataset_names:
+      - Laquila_dscxn
+      - Laquila_ascxn
+    strain_rates:
+    - !beat.StrainRateConfig
+      enabled: false
   dataset_specific_residual_noise_estimation: true
   gf_config: !beat.GeodeticGFConfig
-    store_superdir: /home/vasyurhm/GF/Laquila
     reference_model_idx: 0
-    n_variations: [0, 1]
+    n_variations:
+    - 0
+    - 1
+    store_superdir: /home/vasyurhm/GF/Laquila
     earth_model_name: ak135-f-continental.m
     nworkers: 4
     use_crust2: false
@@ -146,6 +210,7 @@ geodetic_config: !beat.GeodeticConfig
     error_depth: 0.1
     error_velocities: 0.1
     depth_limit_variation: 600.0
+    version: ''
     code: psgrn
     sample_rate: 1.1574074074074073e-05
     sampling_interval: 1.0
@@ -163,24 +228,33 @@ seismic_config: !beat.SeismicConfig
     include: true
     preprocess_data: true
     name: any_P
-    channels: [Z]
+    arrivals_marker_path: ./phase_markers.txt
+    quantity: displacement
+    channels:
+    - Z
     filterer:
     - !beat.heart.Filter
-      lower_corner: 0.001
+      lower_corner: 0.01
       upper_corner: 0.1
       order: 4
-    distances: [30.0, 90.0]
+      stepwise: true
+    distances:
+    - 30.0
+    - 90.0
     interpolation: multilinear
     arrival_taper: !beat.heart.ArrivalTaper
       a: -15.0
       b: -10.0
       c: 40.0
       d: 55.0
+    domain: time
   dataset_specific_residual_noise_estimation: false
   gf_config: !beat.SeismicGFConfig
-    store_superdir: /home/vasyurhm/BEATS/GF
     reference_model_idx: 0
-    n_variations: [0, 1]
+    n_variations:
+    - 0
+    - 1
+    store_superdir: /home/vasyurhm/BEATS/GF
     earth_model_name: ak135-f-continental.m
     nworkers: 4
     use_crust2: false
@@ -193,13 +267,16 @@ seismic_config: !beat.SeismicConfig
     error_depth: 0.1
     error_velocities: 0.1
     depth_limit_variation: 600.0
+    version: ''
     code: qssp
     sample_rate: 0.5
     rm_gfs: false
 sampler_config: !beat.SamplerConfig
   name: SMC
+  backend: bin
   progressbar: true
   buffer_size: 5000
+  buffer_thinning: 10
   parameters: !beat.SMCConfig
     tune_interval: 10
     check_bnd: true
@@ -210,11 +287,13 @@ sampler_config: !beat.SamplerConfig
     coef_variation: 1.0
     stage: 0
     proposal_dist: MultivariateCauchy
-    update_covariances: true
+    update_covariances: false
 hyper_sampler_config: !beat.SamplerConfig
   name: Metropolis
+  backend: bin
   progressbar: true
   buffer_size: 5000
+  buffer_thinning: 1000
   parameters: !beat.MetropolisConfig
     tune_interval: 50
     proposal_dist: Normal

diff --git a/data/examples/Laquila/geodetic_data.pkl b/data/examples/Laquila/geodetic_data.pkl