You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardexpand all lines: README.md
+1-1
Original file line number
Diff line number
Diff line change
@@ -1,7 +1,7 @@
1
1
# BforStrain
2
2
These codes compute strain rates on a trianglular 2D mesh using a geodetic velocity field. The method is based on the paper by K. Johnson (2024) available at https://doi.org/10.1029/2023JB027472 . Body-forces in a thin elastic plate are used to compute a velocity and strain rate field computed at the centroids of the triangles in the mesh. The inversion solves for the distribution of body forces (at the triangle nodes) that best-fits the geodetic velocity field. You can specify a range of smoothing values (called beta) that minimize the magnitude of the body forces. There is also an optional second step that further minimizes strain rates below a specified threshold.
3
3
4
-
Paramters can be specified in params.py. The data used in Johnson (2024) from Zeng (2022) https://doi.org/10.1785/0220220180 is currently loaded, along with creeping fault traces from the 2023 US National Seismic Hazard Model.
4
+
Parameters can be specified in params.py. The data used in Johnson (2024) from Zeng (2022) https://doi.org/10.1785/0220220180 is currently loaded, along with creeping fault traces from the 2023 US National Seismic Hazard Model.
5
5
6
6
To run the demo code, use the code below to download the repository and initialize a new environment
0 commit comments