reposition plate to bottom of overburden is the plate elevation (depth) is less that the overburden thickness

Author: benk-mira

simpeg_drivers/plate_simulation/models/options.py

@@ -9,6 +9,8 @@
 # '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 
 
+from logging import getLogger
+
 import numpy as np
 from geoapps_utils.modelling.plates import PlateModel
 from geoapps_utils.utils.locations import topo_drape_elevation
@@ -18,13 +20,15 @@
     BaseModel,
     ConfigDict,
     Field,
-    field_validator,
     model_validator,
 )
 
 from simpeg_drivers.options import Deprecated
 
 
+logger = getLogger(__name__)
+
+
 class PlateOptions(BaseModel):
     """
     Parameters describing an anomalous plate.
@@ -112,3 +116,14 @@ class ModelOptions(BaseModel):
     background: float
     overburden_options: OverburdenOptions
     plate_options: PlateOptions
+
+    @model_validator(mode="after")
+    def plate_top_below_overburden(self):
+        if self.plate_options.geometry.elevation < self.overburden_options.thickness:
+            logger.warning(
+                "Overburden thickness exceeds the plate depth.  Adjusting"
+                "plate to bottom of overburden to preserve plate's geometry."
+            )
+            self.plate_options.geometry.elevation = self.overburden_options.thickness
+
+        return self

tests/plate_simulation/runtest/driver_test.py

@@ -8,6 +8,8 @@
 #                                                                                   '
 # '''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
 
+import logging
+
 from geoh5py.groups import SimPEGGroup
 from geoh5py.ui_json import InputFile
 
@@ -33,7 +35,7 @@
 
 
 # pylint: disable=too-many-statements
-def test_plate_simulation_params_from_input_file(tmp_path):
+def test_plate_simulation_params_from_input_file(tmp_path, caplog):
     opts = SyntheticsComponentsOptions(
         method="gravity",
         survey=SurveyOptions(n_stations=8, n_lines=8),
@@ -87,7 +89,9 @@ def test_plate_simulation_params_from_input_file(tmp_path):
         ifile.data["spacing"] = 10.0
         ifile.data["elevation"] = 20
 
-    params = PlateSimulationOptions.build(ifile)
+    with caplog.at_level(logging.WARNING):
+        params = PlateSimulationOptions.build(ifile)
+    assert "Overburden thickness exceeds the plate depth" in caplog.text
     assert isinstance(params.simulation, SimPEGGroup)
 
     simulation_parameters = params.simulation_parameters()
@@ -124,4 +128,5 @@ def test_plate_simulation_params_from_input_file(tmp_path):
 
     assert params.model.plate_options.number == 9
     assert params.model.plate_options.spacing == 10.0
-    assert params.model.plate_options.geometry.elevation == 20.0
+    # reset by validator
+    assert params.model.plate_options.geometry.elevation == 50.0