Multi-dimensional earthquake cycle simulation code based on Boundary Element Method with H-matrices.
What's new in Version 2025.10.0
- Merged fluid versions (main_fv.f90 and main_fv3dp.d90) to main_LH.f90 with more flexible parameters and boundary conditions
- Simulations using analytical solutions of fluid pressure for point injection source
- Improved output (now the arrays are NOT reordered in output)
Future extensions
- Fault opening
- GPU optimization
Features
-
2D and 3D planar/nonplanar faults in a full/half space
-
Rectangular and triangular meshes for 3D problems (require .stl mesh file for triangular mesh)
-
Rate and state friction law (aging law, slip law, cut-off velocity model, flash heating, and modified CNS model)
-
Fault-zone linear and nonlinear viscous flow
-
Quasi-dynamic approximation using radiation-damping term
-
Along-fault fluid pressure diffusion with permeability evolution
-
Many examples for input and parameter files
-
Hybrid MPI/open-MP parallelization
See Documentation for more information.
This software is freely available under the MIT license. If you write a paper using this code, please cite the following paper:
So Ozawa, Akihiro Ida, Tetsuya Hoshino, Ryosuke Ando (2023), "Large-scale earthquake sequence simulations of 3D geometrically complex faults using the boundary element method accelerated by lattice H-matrices", Geophysical Journal International,232 (3), 1471-1481 https://doi.org/10.1093/gji/ggac386
If you write a paper using this code with fluid pressure evolution, please cite the following paper:
So Ozawa, Yuyun Yang, Eric M. Dunham (2024), "Fault Valve Instability: A mechanism for slow slip events", Journal of Geophysical Research: Solid Earth, 129, https://doi.org/10.1029/2024JB029165,