update commits in hdf5_loop branch from v2.0

NikkhooWalter2015/Makefile.debug

@@ -17,3 +17,4 @@ clean:
 	rm -f *.o debug_nikkhoo
 
 .PHONY: clean
+

NikkhooWalter2015/debug_nikkhoo.f90

@@ -29,9 +29,10 @@ program debug_nikkhoo
   logical :: casep_log, casen_log, casez_log
 
   ! Initialize test data (single point for debugging)
-  x = [-1.0_DP/3.0_DP]
-  y = [-1.0_DP/3.0_DP]
-  z = [-14.0_DP/3.0_DP]
+
+  x = [7.0_DP]
+  y = [-1.0_DP]
+  z = [-5.0_DP]
 
   p1 = [-1.0_DP, -1.0_DP, -5.0_DP]
   p2 = [1.0_DP, -1.0_DP, -5.0_DP]
@@ -143,7 +144,7 @@ program debug_nikkhoo
   write(*,*) ''
 
   ! Determine configuration
-  call trimode_finder(y_td, z_td, x_td, p1_td, p2_td, p3_td, trimode)
+  call trimode_finder(x_td, y_td, z_td, p1_td, p2_td, p3_td, trimode)
 
   casep_log = (trimode == 1)
   casen_log = (trimode == -1)

NikkhooWalter2015/test_nikkhoo.f90

@@ -9,18 +9,27 @@ program test_nikkhoo
   implicit none
 
   ! Test parameters
-  integer, parameter :: n_points = 7
+  integer, parameter :: n_points = 1
   integer :: i
   real(DP), dimension(n_points) :: x, y, z
   real(DP), dimension(3) :: p1, p2, p3
   real(DP) :: ss, ds, ts, mu, lambda
   real(DP), dimension(n_points, 6) :: stress, strain
 
   ! Initialize test data
-  x = [-1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP/3.0_DP, 7.0_DP, -7.0_DP, -1.0_DP, -1.0_DP]
-  y = [-1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP, -1.0_DP, -3.0_DP, 3.0_DP]
-  z = [-3.0_DP, -14.0_DP/3.0_DP, -6.0_DP, -5.0_DP, -5.0_DP, -6.0_DP, -3.0_DP]
+  ! Original test points
+  !x = [-1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP/3.0_DP, 7.0_DP, -7.0_DP, -1.0_DP, -1.0_DP, &
+  !     3.0_DP, -3.0_DP, -1.0_DP, -1.0_DP, 1.0_DP, -1.0_DP, -1.0_DP, 1.0_DP]
+  !y = [-1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP/3.0_DP, -1.0_DP, -1.0_DP, -3.0_DP, 3.0_DP, &
+  !     -3.0_DP, 3.0_DP, -1.0_DP, 1.0_DP, -1.0_DP, -1.0_DP, 1.0_DP, -1.0_DP]
+  !z = [-3.0_DP, -14.0_DP/3.0_DP, -6.0_DP, -5.0_DP, -5.0_DP, -6.0_DP, -3.0_DP, &
+  !     -6.0_DP, -3.0_DP, -1.0_DP, -1.0_DP, -1.0_DP, -8.0_DP, -8.0_DP, -8.0_DP]
 
+
+  x = [-1.0_DP/3.0_DP]
+  y = [-1.0_DP/3.0_DP]
+  z = [-14.0_DP/3.0_DP]
+
   p1 = [-1.0_DP, -1.0_DP, -5.0_DP]
   p2 = [1.0_DP, -1.0_DP, -5.0_DP]
   p3 = [-1.0_DP, 1.0_DP, -4.0_DP]

README.md

@@ -1,3 +1,20 @@
+### SEAS Benchmark BP6
+
+The Seismological Grand Challenge in Subduction Megathrust Earthquake Modeling (SEAS) Benchmark Problem 6 (BP6) is a rate-and-state, quasi-dynamic earthquake cycle problem on a 3D dipping fault. This repository includes an example setup and data for BP6 under `example2/`.
+
+- **Objective**: Validate and compare numerical implementations for quasi-dynamic cycles with rate-and-state friction on a geometrically complex fault.
+- **What’s included**: Input files, derived profiles, and a reference sketch of the BP6 configuration.
+- **Where to start**: See `example2/` for inputs and the mesh, and `src/` for the solver modules.
+
+#### BP6 Configuration Sketch
+
+![BP6 configuration sketch](example2/bp6-sketch.png)
+
+The sketch illustrates the dipping fault geometry, nucleation region, and boundary conditions typically used for BP6. Use it alongside the files under `example2/` to reproduce the configuration and compare with benchmark outputs.
+
+
+
+
 # TriBIE
 
 **TriBIE** (Triangular Boundary Integral Equation) is a high-performance Fortran90 parallel computing framework for simulating earthquake cycles, slow slip events, and aseismic transients on complex 3D fault geometries. The code employs hybrid MPI+OpenMP parallelization with advanced computational optimizations including SIMD vectorization and dynamic load balancing to efficiently model rate-and-state friction physics on triangular fault meshes embedded in elastic half-space media.