Merge pull request #740 from baagaard-usgs/fix-examples-fault-output-fields

baagaard-usgs · web-flow · commit b3bd749f924f · 2024-06-05T21:50:54.000-06:00
FIX: Add fault traction change to output in examples.
diff --git a/docs/user/components/faults/FaultCohesiveKin.md b/docs/user/components/faults/FaultCohesiveKin.md
@@ -63,7 +63,7 @@ Example of setting `FaultCohesiveKin` Pyre properties and facilities in a parame
 label = fault
 edge = fault_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Two earthquakes with different slip time functions.
 eq_ruptures = [quake10, quake50]
diff --git a/docs/user/examples/reverse-2d/step05-onefault.md b/docs/user/examples/reverse-2d/step05-onefault.md
@@ -53,7 +53,7 @@ label = fault
 label_value = 20
 edge = fault_end
 edge_value = 21
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.UniformDB
diff --git a/docs/user/examples/reverse-2d/step06-twofaults-elastic.md b/docs/user/examples/reverse-2d/step06-twofaults-elastic.md
@@ -66,7 +66,7 @@ label_value = 20
 edge = fault_end
 edge_value = 21
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.UniformDB
@@ -80,7 +80,7 @@ label_value = 22
 edge = splay_end
 edge_value = 23
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.splay.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.UniformDB
diff --git a/docs/user/examples/strikeslip-2d/common-information.md b/docs/user/examples/strikeslip-2d/common-information.md
@@ -98,8 +98,8 @@ interfaces = [fault]
 label = fault
 label_value = 20
 
-# Output `slip` on the fault.
-observers.observer.data_fields = [slip]
+# Output `slip` and the change in fault traction on the fault.
+observers.observer.data_fields = [slip, traction_change]
 ```
 
 All simulations in this directory use Dirichlet boundary conditions on the -x and +x boundaries with zero displacements,
diff --git a/docs/user/examples/subduction-2d/step01-coseismic.md b/docs/user/examples/subduction-2d/step01-coseismic.md
@@ -34,7 +34,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.SimpleDB
diff --git a/docs/user/examples/subduction-2d/step02-interseismic.md b/docs/user/examples/subduction-2d/step02-interseismic.md
@@ -46,7 +46,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault_slabtop.eq_ruptures]
 rupture = pylith.faults.KinSrcConstRate
@@ -64,7 +64,7 @@ label_value = 22
 edge = fault_slabbot_edge
 edge_value = 32
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault_slabbot.eq_ruptures]
 rupture = pylith.faults.KinSrcConstRate
diff --git a/docs/user/examples/subduction-2d/step03-eqcycle.md b/docs/user/examples/subduction-2d/step03-eqcycle.md
@@ -34,7 +34,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 eq_ruptures = [creep, earthquake]
 
diff --git a/examples/meshing-cubit/surface-nurbs/merge-surfs/pylithapp.cfg b/examples/meshing-cubit/surface-nurbs/merge-surfs/pylithapp.cfg
@@ -104,7 +104,7 @@ label = fault
 edge = fault_edge
 
 observers.observer.writer.filename = output/faultslip-fault.h5
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.UniformDB
diff --git a/examples/reverse-2d/step05a_onefault.cfg b/examples/reverse-2d/step05a_onefault.cfg
@@ -75,7 +75,7 @@ label = fault
 label_value = 20
 edge = fault_end
 edge_value = 21
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 # We have uniform slip, so we use a `UniformDB` spatial database.
diff --git a/examples/reverse-2d/step05b_onefault.cfg b/examples/reverse-2d/step05b_onefault.cfg
@@ -82,7 +82,7 @@ label = fault
 label_value = 20
 edge = fault_end
 edge_value = 21
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 # We have uniform slip, so we use a `UniformDB` spatial database.
diff --git a/examples/reverse-2d/step05c_onefault.cfg b/examples/reverse-2d/step05c_onefault.cfg
@@ -96,7 +96,7 @@ label = fault
 label_value = 20
 edge = fault_end
 edge_value = 21
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 # We have uniform slip, so we use a `UniformDB` spatial database.
diff --git a/examples/reverse-2d/step06_twofaults_elastic.cfg b/examples/reverse-2d/step06_twofaults_elastic.cfg
@@ -94,7 +94,7 @@ label_value = 20
 edge = fault_end
 edge_value = 21
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 # Prescribe 2 meters of reverse (right-lateral) slip.
@@ -110,7 +110,7 @@ label_value = 22
 edge = splay_end
 edge_value = 23
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.splay.eq_ruptures.rupture]
 # Prescribe 1 meter of reverse (right-lateral) slip.
diff --git a/examples/reverse-2d/step07_twofaults_maxwell.cfg b/examples/reverse-2d/step07_twofaults_maxwell.cfg
@@ -102,7 +102,7 @@ label_value = 20
 edge = fault_end
 edge_value = 21
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 origin_time = 0.0*year
@@ -117,7 +117,7 @@ label_value = 22
 edge = splay_end
 edge_value = 23
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.splay.eq_ruptures.rupture]
 origin_time = 39.99*year
diff --git a/examples/reverse-2d/step08_twofaults_powerlaw.cfg b/examples/reverse-2d/step08_twofaults_powerlaw.cfg
@@ -126,7 +126,7 @@ label_value = 20
 edge = fault_end
 edge_value = 21
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 origin_time = 0.0*year
@@ -141,7 +141,7 @@ label_value = 22
 edge = splay_end
 edge_value = 23
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.splay.eq_ruptures.rupture]
 origin_time = 39.99*year
diff --git a/examples/strikeslip-2d/pylithapp.cfg b/examples/strikeslip-2d/pylithapp.cfg
@@ -138,8 +138,8 @@ interfaces = [fault]
 label = fault
 label_value = 20
 
-# Output `slip` on the fault.
-observers.observer.data_fields = [slip]
+# Output `slip` and the change in tractions on the fault.
+observers.observer.data_fields = [slip, traction_change]
 
 
 # ----------------------------------------------------------------------
diff --git a/examples/strikeslip-2d/step05_greensfns.cfg b/examples/strikeslip-2d/step05_greensfns.cfg
@@ -93,6 +93,9 @@ fault = pylith.faults.FaultCohesiveImpulses
 # Limit Green's functions to left-lateral slip (dof=1)
 impulse_dof = [1]
 
+# Output `slip` on the fault.
+observers.observer.data_fields = [slip]
+
 # We use a spatial database to specify the region of the fault where we want
 # slip impulses for Green's functions.
 db_auxiliary_field = spatialdata.spatialdb.SimpleDB
diff --git a/examples/subduction-2d/step01_coseismic.cfg b/examples/subduction-2d/step01_coseismic.cfg
@@ -45,7 +45,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 # We specify slip as a function of depth, so we use a `SimpleDB` with
diff --git a/examples/subduction-2d/step02_interseismic.cfg b/examples/subduction-2d/step02_interseismic.cfg
@@ -56,7 +56,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault_slabtop.eq_ruptures]
 # Switch to constant slip rate time function.
@@ -80,7 +80,7 @@ label_value = 22
 edge = fault_slabbot_edge
 edge_value = 32
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Switch to constant slip rate time function.
 [pylithapp.problem.interfaces.fault_slabbot.eq_ruptures]
diff --git a/examples/subduction-2d/step03_eqcycle.cfg b/examples/subduction-2d/step03_eqcycle.cfg
@@ -57,7 +57,7 @@ label_value = 21
 edge = fault_slabtop_edge
 edge_value = 31
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Set earthquake sources to an array consisting of creep and 1 earthquake.
 eq_ruptures = [creep, earthquake]
@@ -96,7 +96,7 @@ label_value = 22
 edge = fault_slabbot_edge
 edge_value = 32
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Switch to constant slip rate time function.
 [pylithapp.problem.interfaces.fault_slabbot.eq_ruptures]
diff --git a/examples/subduction-3d/step02_coseismic.cfg b/examples/subduction-3d/step02_coseismic.cfg
@@ -78,7 +78,7 @@ interfaces = [fault_slabtop]
 label = fault_slabtop_patch
 edge = fault_slabtop_patch_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # We prescribe uniform slip, so we use ia `UniformDB`.
 # The earthquake rupture occurs at t=10 years.
diff --git a/examples/subduction-3d/step03_interseismic.cfg b/examples/subduction-3d/step03_interseismic.cfg
@@ -87,7 +87,7 @@ edge = fault_slabbot_edge
 # along-strike and up-dip directions.
 ref_dir_1 = [-0.1, 0, 0.9]
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Use the constant slip rate time function.
 eq_ruptures.rupture = pylith.faults.KinSrcConstRate
@@ -104,7 +104,7 @@ db_auxiliary_field.data = [+2.0*cm/year, -4.0*cm/year, 0.0*cm/year, 0.0*year]
 label = fault_slabtop
 edge = fault_slabtop_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Use the constant slip rate time function.
 eq_ruptures.rupture = pylith.faults.KinSrcConstRate
diff --git a/examples/subduction-3d/step04_eqcycle.cfg b/examples/subduction-3d/step04_eqcycle.cfg
@@ -89,7 +89,7 @@ edge = fault_slabbot_edge
 # along-strike and up-dip directions.
 ref_dir_1 = [-0.1,0,0.9]
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Use the constant slip rate time function.
 eq_ruptures.rupture = pylith.faults.KinSrcConstRate
@@ -106,7 +106,7 @@ db_auxiliary_field.data = [+2.0*cm/year, -4.0*cm/year, 0.0*cm/year, 0.0*year]
 label = fault_slabtop
 edge = fault_slabtop_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Define 3 earthquake ruptures with different origin times.
 eq_ruptures = [creep, eq1, eq2]
@@ -141,7 +141,7 @@ db_auxiliary_field.query_type = linear
 label = fault_splay
 edge = fault_splay_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # Origin time for splay fault earthquake.
 eq_ruptures.rupture.origin_time = 249.999*year
diff --git a/examples/subduction-3d/step06_slowslip.cfg b/examples/subduction-3d/step06_slowslip.cfg
@@ -108,7 +108,7 @@ interfaces = [fault_slabtop]
 label = fault_slabtop_patch
 edge = fault_slabtop_patch_edge
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 # We use a time history slip function.
 [pylithapp.problem.interfaces.fault_slabtop.eq_ruptures]
diff --git a/pylith/faults/FaultCohesiveKin.py b/pylith/faults/FaultCohesiveKin.py
@@ -37,7 +37,7 @@ class FaultCohesiveKin(FaultCohesive, ModuleFaultCohesiveKin):
             label = fault
             edge = fault_edge
 
-            observers.observer.data_fields = [slip]
+            observers.observer.data_fields = [slip, traction_change]
 
             # Two earthquakes with different slip time functions.
             eq_ruptures = [quake10, quake50]
diff --git a/tests/fullscale/cornercases/faults-2d/pylithapp.cfg b/tests/fullscale/cornercases/faults-2d/pylithapp.cfg
@@ -84,7 +84,7 @@ interfaces = [fault]
 
 [pylithapp.problem.interfaces.fault]
 label = fault
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 
 # ----------------------------------------------------------------------
diff --git a/tests/fullscale/linearelasticity/faults-2d/twoblocks_ic.cfg b/tests/fullscale/linearelasticity/faults-2d/twoblocks_ic.cfg
@@ -29,7 +29,7 @@ interfaces = [fault]
 label = fault_xmid
 label_value = 20
 
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
 db_auxiliary_field = spatialdata.spatialdb.UniformDB
diff --git a/tests/manual/2d/faults_2d_buried/faulttest.cfg b/tests/manual/2d/faults_2d_buried/faulttest.cfg
@@ -30,7 +30,7 @@ id = 10
 label = fault
 edge = fault_edge
 observers.observer.writer.filename = faultedge_test-fault.h5
-observers.observer.data_fields = [slip]
+observers.observer.data_fields = [slip, traction_change]
 observers.observer.info_fields = []
 
 [pylithapp.problem.interfaces.fault.eq_ruptures.rupture]
