defaults: extend physical bounds for lat, lon to WGS84

heart: GeodeticDataset.setup_corrections() check for los_vector
hvasbath · Apr 16, 2024 · 3b3f50f · 3b3f50f
1 parent b03d2b4
commit 3b3f50f
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 11 deletions.
diff --git a/beat/defaults.py b/beat/defaults.py
@@ -184,11 +184,11 @@ def __getitem__(self, k):
         default_bounds=(-200.0, 200.0),
         unit=u_nanostrain,
     ),
-    "pole_lat": Bounds(
-        physical_bounds=(-90.0, 90.0), default_bounds=(0, 1), unit=u_deg
+    "lat": Bounds(
+        physical_bounds=(-90.0, 90.0), default_bounds=(30.0, 30.5), unit=u_deg
     ),
-    "pole_lon": Bounds(
-        physical_bounds=(-180.0, 180.0), default_bounds=(0, 1), unit=u_deg
+    "lon": Bounds(
+        physical_bounds=(-180.0, 180.0), default_bounds=(30.0, 30.5), unit=u_deg
     ),
     "omega": Bounds(
         physical_bounds=(-10.0, 10.0), default_bounds=(0.5, 0.6), unit=u_deg_myr
@@ -239,20 +239,14 @@ def __getitem__(self, k):
         physical_bounds=(0.0, 1.0), default_bounds=(0.0, 1.0), unit=u_hyp
     ),
     "hypers": Bounds(
-        physical_bounds=(-4.0, 10.0), default_bounds=(-2.0, 6.0), unit=u_hyp
+        physical_bounds=(-10.0, 10.0), default_bounds=(-2.0, 6.0), unit=u_hyp
     ),
     "ramp": Bounds(
         physical_bounds=(-0.005, 0.005), default_bounds=(-0.005, 0.005), unit=u_rad
     ),
     "offset": Bounds(
         physical_bounds=(-0.05, 0.05), default_bounds=(-0.05, 0.05), unit=u_m
     ),
-    "lat": Bounds(
-        physical_bounds=(30.0, 30.5), default_bounds=(30.0, 30.5), unit=u_deg
-    ),
-    "lon": Bounds(
-        physical_bounds=(30.0, 30.5), default_bounds=(30.0, 30.5), unit=u_deg
-    ),
     "traction": Bounds(physical_bounds=(0, 1000), default_bounds=(0, 50), unit=u_pa),
     "strike_traction": Bounds(
         physical_bounds=(-15000, 15000), default_bounds=(-50, 50), unit=u_pa

diff --git a/beat/heart.py b/beat/heart.py
@@ -1099,6 +1099,9 @@ def setup_corrections(self, event, correction_configs):
 
                 data_mask = self.get_data_mask(corr_conf)
 
+                if self.los_vector is None:
+                    self.update_los_vector()
+
                 corr.setup_correction(
                     locy=locy,
                     locx=locx,