Add ability to refine meshes used in output to better capture higher order discretization

.vscode/c_cpp_properties.json

@@ -23,7 +23,7 @@
             ],
             "cStandard": "c11",
             "cppStandard": "c++14",
-            "intelliSenseMode": "macos-clang-x64"
+            "intelliSenseMode": "macos-clang-arm64"
         },
         {
             "name": "linux-gcc",

CHANGES.md

@@ -7,16 +7,19 @@ The version numbers are in the form `MAJOR.MINOR.PATCH`, where major releases in
 
 ## Version 4.2.0
 
+* **Changed**
+  * Updated `examples/strikeslip-2d` Steps 4-7 to use a more realistic slip distribution and mesh refinement in output.
+  * Updated to PETSc 3.22.0
+  * Switch CI from Azure Pipelines to GitHub Actions.
 * **Added**
   * Default filenames for progress monitor and parameters file are set from the simulation name like
   the other output files.
   * Consistency check for material properties and scales used in nondimensionalization (currently just the shear modulus).
   * Add section to User Guide on troubleshooting solver issues.
-* **Changed**
-  * Switch CI from Azure Pipelines to GitHub Actions.
+  * Allow output on a finer resolution mesh than used in the simulation to facilitate accurate representation of fields with a basis order of 2 or greater.
 * **Fixed**
   * Fix inconsistency in normal direction on fault surfaces. Orientation was correct but direction was flipped at some locations. This affected local slip direction and the resulting deformation close to the fault. This bug fix was not in version 4.1.3.
-  * Update autoconf macros for numpy for compatibility with location of include files in numpy version 2.x.
+  * Update autoconf numpy macros for compatibility with location of include files in numpy version 2.x.
 
 ## Version 4.1.3 (2024/07/31)
 

docs/Makefile.am

@@ -473,13 +473,15 @@ dist_noinst_DATA = \
 	user/examples/strikeslip-2d/figs/step04-solution.jpg \
 	user/examples/strikeslip-2d/figs/step05-diagram.pdf \
 	user/examples/strikeslip-2d/figs/step05-diagram.svg \
+	user/examples/strikeslip-2d/figs/step05_greensfns-impulses.pdf \
+	user/examples/strikeslip-2d/figs/step05_greensfns-impulses.svg \
 	user/examples/strikeslip-2d/figs/step05-solution.jpg \
-	user/examples/strikeslip-2d/figs/step06-solution.pdf \
-	user/examples/strikeslip-2d/figs/step06-solution.svg \
-	user/examples/strikeslip-2d/figs/step07a-solution.pdf \
-	user/examples/strikeslip-2d/figs/step07a-solution.svg \
-	user/examples/strikeslip-2d/figs/step07b-solution.pdf \
-	user/examples/strikeslip-2d/figs/step07b-solution.svg \
+	user/examples/strikeslip-2d/figs/step06_inversion-results.pdf \
+	user/examples/strikeslip-2d/figs/step06_inversion-results.svg \
+	user/examples/strikeslip-2d/figs/step07a_catmip-results.pdf \
+	user/examples/strikeslip-2d/figs/step07a_catmip-results.svg \
+	user/examples/strikeslip-2d/figs/step07b_cfcatmip-results.pdf \
+	user/examples/strikeslip-2d/figs/step07b_cfcatmip-results.svg \
 	user/examples/subduction-2d/index.md \
 	user/examples/subduction-2d/common-information.md \
 	user/examples/subduction-2d/exercises.md \

docs/user/components/meshio/OutputObserver.md

@@ -21,4 +21,8 @@ Abstract base class for managing output of solution information.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 

docs/user/components/meshio/OutputPhysics.md

@@ -33,6 +33,10 @@ Implements `OutputObserver`.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 
 ## Example
 

docs/user/components/meshio/OutputSoln.md

@@ -26,4 +26,8 @@ Implements `OutputObserver`.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 

docs/user/components/meshio/OutputSolnBoundary.md

@@ -29,14 +29,18 @@ Implements `OutputSoln`.
 * `label`=\<str\>: Name of label identifier for external boundary.
   - **default value**: ''
   - **current value**: '', from {default}
-  - **validator**: <function validateLabel at 0x11f366af0>
+  - **validator**: <function validateLabel at 0x10816d800>
 * `label_value`=\<int\>: Value of label identifier for external boundary (tag of physical group in Gmsh files).
   - **default value**: 1
   - **current value**: 1, from {default}
 * `output_basis_order`=\<int\>: Basis order for output.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 
 ## Example
 

docs/user/components/meshio/OutputSolnDomain.md

@@ -30,6 +30,10 @@ Implements `OutputSoln`.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 
 ## Example
 

docs/user/components/meshio/OutputSolnPoints.md

@@ -36,6 +36,10 @@ Implements `OutputSoln`.
   - **default value**: 1
   - **current value**: 1, from {default}
   - **validator**: (in [0, 1])
+* `refine_levels`=\<int\>: Number of mesh refinement levels for output.
+  - **default value**: 0
+  - **current value**: 0, from {default}
+  - **validator**: (greater than or equal to 0)
 
 ## Example