v1.3.2 - stable version - bp5 and fault roughness

This version contains the following major updates:

1. All adjustable parameters are in the user_defined_parameters.py for different cases.
2. Utilities to use the code are simplified to create.newcase, case.setup, and case.submit.
3. netCDF format is adopted for output.
4. Python-based post-processing utilities are developed to convert netCDF results to plots and movies.
5. EQquasi can stop at any time slice and restart from that time slice.
6. EQquasi is verified against BP5 of SEAS. The example code compset bp-5-2000 is tested.
7. Rough fault geometry is implemented and works for one earthquake cycle.
8. Normal stress evolution is implemented following Shi and Day (2013).

Code improvement.

Spins off case specific setup to input files. Code improvement. Tested on bp5-2000m with MUMPS. AZTEC performs poorly.

AZTEC incorporated

The parallel iterative solver AZTEC is implemented. The code is tested on SEAS BP5.

Version 1.2.1

Major Achievements:

• Combining the first cycle and later cycles calculated by EQquasi in one code;
• Simplification of code structure. No variable transference.
• Simulation of BP5-QD (quasi-dynamic with a radiated damping term).

Version: M-1.2.1
Author: Dunyu Liu ([email protected])
Branch: Master Parent version: M-1.2.0

Changes:

1. Remove slip rate caps in faulting.f90 when doing the Newton method.
2. Add the dipping dimension for off-fault station output.
3. Add dymax cap for element size along the y direction in the globalvar.f90.
4. Post processing matlab scripts are updated.
5. [Verification] The code is verified in SCEC SEAS BP5 with 2000 m, 1000 m, and 500 m resolutions.

1_0_0

Major achievments:

• EQsimu used in Liu et al. (2020) with MUMPS;
• EQdyna is modified from version 4.1.3.

Version: M-1.0.0
Author: DL
Branch: Master Parent version: Root

Changes:

• New: EQquasi is developed to simulate fully dynamic earthquake cycles on a bent fault.
• New: MUMPS, a parallel direct solver, is implemented to handle the heavy computing loads.
• New: The non-co-seismic deformation in the first and later cycles are handled in two codes due to initial conditions set in the Q0/ and loaded from input files in Qi/.
• New: libs and headfiles of MUMPS are put in ./library/
• New: License, REAME.md, Chagelog.md are created.

1_2_0

Major Achievements:

• Combining the first cycle and later cycles calculated by EQquasi in one code;
• Simplification of code structure. No variable transference.
• Simulation of BP5-QD (quasi-dynamic with a radiated damping term).

Version: M-1.2.0
Author: DL
Branch: Master Parent version: M-1.1.0

Problems:

• The first earthquake starts quickly and yields reasonable ‘rupture’ time contour. However, later cycles show abrupt change of maximum slip rate in the transition zone of frictional parameter a. A 1e-5 m/s will suddenly jump out of a sequence of 1e-9 m/s slip rates. Will model sizes affect it? The old models have 120 by 60 by 120 in x, y, z directions. Now it is changed to 160 by 70 by 50 instead.
  [Solved] Because the critical slip distance L is increased to 0.13/0.14m in BP5 compared to the 0.04 m in BP4, so the previous formula to calculate time step dtev = ksi * L/maxsliprate yields too rough time step and solution blows up. Ksi is reduced to 0.05 from 0.2 in the current model and it works. [Continued] Ksi is reduced to 0.025 based on equations in Lapusta et al. (2000). A matlab script is developed to process it.

Changes:

• New: A license is created.
• New: Use one code to simulate all earthquake cycles. For the first cycle, initial conditions are assigned. For later ones, initial conditions are imported from cplot_EQquasi.txt.
• New: Significant simplification of the code structure. No variable is transferred between major subroutines.
• New: files are renamed.
• New: library_bound.f90 is created. It contains two subroutines bound_load and bound_ft_ku to handle far_field loading and calculate KU on faults for use in subroutine faulting.
• New: library_output.f90 is created. It contains several subroutines output_onfault_st, output_offfault_st, output_onfault_transfer, output_timedy, output_globaldat, and output_prof. Future outputs could be implemented here.
• New: Read_Input_Files.f90 is copied from EQdyna 3D version 5.2.0 to import station locations and to read in icstart from currentcycle.txt, which is created by the batch script. Other systematic parameters should be implemented in a later time.
• New: crs.f90 is created to handle the CRS format of stiffness structure construction.
• New: elemcal.f90 is created to handle the construction of the stiffness matrix and mass array in the CRS format.
• New: ./script is created for MATLAB scripts for post processing.
  Improvement: on-screen information printing is improved.
• [Verification] SCEC SEAS benchmark BP-5-QD.