-
This is a trimmed version of Quaycle.jl with only elastic part. Full building, testing and docs building are not included here.
-
Citation: Wei, M., & Pengcheng, S. (2021). Synchronization of earthquake cycles of adjacent segments on oceanic transform faults revealed by numerical simulation in the framework of rate-and-state friction. Journal of Geophysical Research: Solid Earth, 126, e2020JB020231. https://doi.org/10.1029/2020JB020231
This repo constains the sample scripts for reproducing simulations and the outputs of catalogues.
-
scripts:
s01-domain.jl
: create a fault domains02-greensfunc.jl
: compute dislocation-stress Green's functions03-parameters.jl
: set fault parameterss04-solve.jl
: solve modelsscanfunc.jl
: functions for computing event catalogues from raw output datasolve.job
: an example sbatch script for solving each model using one node (multi-threading within node, array job)
-
outputs (for catalogue only):
All model outputs are in HDF5.
ic{#}/
denotes different initial conditions.extend/
contains extended simulations (up to 1200 or 1800 years, end-year is appened in those output names).Catalog output names are like
resf-{parameter group}{value}.h5
. LF denote left half fault, RF right half fault. Fields in the catalogue outputs are:t
: time step in secondti
: time step in second (downsampled if_stride != 1
)maxva
: max velocity in m/s (LF)maxvb
: max velocity in m/s (RF)maxvf
: max velocity in m/s (whole fault)mwa
: moment magnitude (LF)mwb
: moment magnitude (RF)ixba
: event start time (LF), correspond toti
ixbb
: event start time (RF), correspond toti
ixea
: event end index (LF), correspond toti
ixeb
: event end index (RF), correspond toti
riax
: event initial rutpure along strike position (m) (LF)ribx
: event initial rutpure along strike position (m) (RF)riaz
: event initial rutpure down dip position (m) (LF)ribz
: event initial rutpure down dip strike position (m) (RF)
Slip ratio output names are like
srf-{parameter group}{value}.h5
. Fields in the catalogue outputs are:sr
: seismic ratioar
: afterslip ratio
{output type}-s-{parameter group}{value}.h5
denote single VW model.ODEs solution output names are like
otf-{parameter group}{value}.h5
(not included in this repo). They contains three fields:t
: time step (s)v
: velocity (m/s)θ
: state variable (s)δ
: displacement (m)