Skip to content

Commit c18909d

Browse files
Thomas-Ulrichfabian-kutscheraAangniu
authored
Preparation Quakeworx workshop, including symbolic links (#55)
* Preparation Quakeworx workshop * fix file format * add missing file * add another missing file * add new example, cdb_tpv23 * add the input files for cdb_tpv23 * update cdb_tpv23 notebook * add jupyter notebook for tpv13, with simple fault output * add link to the existing images in the notebook for tpv13 * adapt output dir for quakeworx * update notebook for cdb_tpv23 * further polishing of the cdb_tpv23 notebook * polish the modiefied tpv13_qwx notebook * another round of edits for tpv13 quakeworx notebook * clean up Kaikoura and remove RS * Revert "clean up Kaikoura and remove RS" This reverts commit f4a1de5. * update notebooks * add gitignore * smaller visualization updates * change file paths to those on the App * minor edit to description Kaikoura * change the output dir to those that we use for training * small changes to text cells * Add a link to download the input files to ReadMe. * replace by symbolic links --------- Co-authored-by: fabian-kutschera <fkutschera@ucsd.edu> Co-authored-by: aangniu <aangniu96@gmail.com> Co-authored-by: Aangniu <38255610+Aangniu@users.noreply.github.com>
1 parent 1b2bad0 commit c18909d

31 files changed

+2251
-0
lines changed

‎.DS_Store‎

8 KB
Binary file not shown.

‎.gitignore‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
*.ipynb_checkpoints/

‎quakeworx/README.md‎

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,5 @@
1+
# [Quakeworx](https://quakeworx.org)
2+
3+
The files in this directory have been adjusted compared to the "standard" SeisSol training repository to be readily used for a science gateway for seismic simulations: [Quakeworx Kick-off: Advancing Earthquake Science and Cybertraining in Seismology](https://quakeworx.org/events/quakeworx-kick-off-workshop)
4+
5+
The input files for the Gateway can be downloaded [here](https://syncandshare.lrz.de/getlink/fiVsNGoWemY1u7To7Gh1pg/quakeworx.zip) as well if you do not have access to Git.

‎quakeworx/cdb_tpv23/cdb_tpv23.ipynb‎

Lines changed: 680 additions & 0 deletions
Large diffs are not rendered by default.

‎quakeworx/cdb_tpv23/off_fault_damage.png‎

33.8 KB
Loading

‎quakeworx/cdb_tpv23/parameters_qwx.par‎

Lines changed: 111 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,111 @@
1+
&equations
2+
!yaml file defining spatial dependance of material properties
3+
MaterialFileName = 'tpv23_damMaterial.yaml'
4+
/
5+
6+
&IniCondition
7+
/
8+
9+
&Boundaries
10+
BC_fs = 1 ! enable free surface boundaries
11+
BC_dr = 1 ! enable fault boundaries
12+
BC_of = 1 ! enable absorbing boundaries
13+
/
14+
15+
&DynamicRupture
16+
FL = 16 ! Friction law
17+
!0: none, 16:LSW, 103: RS with strong velocity weakening
18+
!yaml file defining spatial dependance of fault properties
19+
ModelFileName = 'tpv23_fault.yaml'
20+
21+
!characteristic time of the forced nucleation
22+
t_0 = 0.5
23+
24+
!reference vector for defining strike and dip direction
25+
XRef = 0.0 ! Reference point
26+
YRef = -1.0
27+
ZRef = 0.0
28+
refPointMethod = 1
29+
30+
RF_output_on = 1 ! Rupture front ascii output
31+
OutputPointType = 5 ! Type (0: no output, 3: ascii file, 4: paraview file, 5: 3+4)
32+
SlipRateOutputType=1 ! 0: (smoother) slip rate output evaluated from the difference between the velocity on both side of the fault
33+
! 1: slip rate output evaluated from the fault tractions and the failure criterion (less smooth but usually more accurate where the rupture front is well developped)
34+
/
35+
36+
!see: https://seissol.readthedocs.io/en/latest/fault-output.html
37+
! parameterize paraview file output
38+
&Elementwise
39+
printIntervalCriterion = 2 ! 1=iteration, 2=time
40+
printtimeinterval_sec = 1.0 ! Time interval at which output will be written
41+
OutputMask = 1 1 1 0 1 1 1 1 1 1 1 ! turn on and off fault outputs
42+
refinement_strategy = 2
43+
refinement = 1
44+
/
45+
46+
! parameterize ascii fault file outputs
47+
&Pickpoint
48+
printtimeinterval = 1 ! Index of printed info at timesteps
49+
OutputMask = 1 1 1 1 ! turn on and off fault outputs
50+
nOutpoints = 14
51+
PPFileName = 'tpv23_3km_faultreceivers.dat'
52+
/
53+
54+
&SourceType
55+
/
56+
57+
&SpongeLayer
58+
/
59+
60+
&MeshNml
61+
MeshFile = 'stepover23_3km_qwx.puml.h5'
62+
meshgenerator = 'PUML' ! Name of meshgenerator (Gambit3D-fast, Netcdf or PUML)
63+
/
64+
65+
&Discretization
66+
CFL = 0.99 ! CFL number (<=1.0)
67+
FixTimeStep = 5 ! Manualy chosen minimum time
68+
ClusteredLTS = 1 ! 1 for Global time stepping, 2,3,5,... Local time stepping (advised value 2)
69+
!ClusteredLTS defines the multi-rate for the time steps of the clusters 2 for Local time stepping
70+
/
71+
72+
&Output
73+
FaultOutputFlag = 1 ! DR output (add this line only if DR is active)
74+
OutputFile = 'outputs/cdb_tpv23'
75+
Format = 6 ! Format (10= no output, 6=hdf5 output)
76+
! |stress |vel
77+
iOutputMask = 1 1 1 1 1 1 1 1 1
78+
! |strain |eta
79+
iPlasticityMask = 0 0 0 0 0 0 1
80+
TimeInterval = 5. ! Index of printed info at time
81+
refinement = 1
82+
83+
! Free surface output
84+
SurfaceOutput = 1
85+
SurfaceOutputRefinement = 2
86+
SurfaceOutputInterval = 5.0
87+
88+
printIntervalCriterion = 2 ! Criterion for index of printed info: 1=timesteps,2=time,3=timesteps+time
89+
90+
pickdt = 0.005 ! Pickpoint Sampling
91+
pickDtType = 1 ! Pickpoint Type
92+
RFileName = 'tpv23_receivers.dat' ! Record Points in extra file
93+
94+
xdmfWriterBackend = 'posix' ! (optional) The backend used in fault, wavefield,
95+
! and free-surface output. The HDF5 backend is only supported when SeisSol is compiled with
96+
! HDF5 support.
97+
98+
EnergyOutput = 1 ! Computation of energy, written in csv file
99+
EnergyTerminalOutput = 1 ! Write energy to standard output
100+
EnergyOutputInterval = 0.5
101+
/
102+
103+
&AbortCriteria
104+
EndTime = 14
105+
/
106+
107+
&Analysis
108+
/
109+
110+
&Debugging
111+
/

‎quakeworx/cdb_tpv23/stepover23_3km_qwx.puml.h5‎

9.08 MB
Binary file not shown.

‎quakeworx/cdb_tpv23/tpv23.png‎

224 KB
Loading

‎quakeworx/cdb_tpv23/tpv23_3km_faultreceivers.dat‎

Lines changed: 29 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,29 @@
1+
0.0e3 -3.0e3 0.0e3
2+
5.0e3 -3.0e3 0.0e3
3+
20.0e3 -3.0e3 0.0e3
4+
4.0e3 -3.0e3 -5.0e3
5+
4.0e3 -3.0e3 -5.1e3
6+
4.0e3 -3.0e3 -5.5e3
7+
5.0e3 -3.0e3 -5.0e3
8+
5.0e3 -3.0e3 -5.1e3
9+
5.0e3 -3.0e3 -5.5e3
10+
0.0e3 -3.0e3 -10.0e3
11+
5.0e3 -3.0e3 -10.0e3
12+
6.0e3 -3.0e3 -10.0e3
13+
6.0e3 -3.0e3 -9.0e3
14+
6.0e3 -3.0e3 -8.0e3
15+
6.0e3 -3.0e3 -7.0e3
16+
6.0e3 -3.0e3 -6.0e3
17+
6.0e3 -3.0e3 -5.0e3
18+
10.0e3 -3.0e3 -10.0e3
19+
20.0e3 -3.0e3 -10.0e3
20+
7.0e3 -3.0e3 -11.0e3
21+
7.0e3 -3.0e3 -10.0e3
22+
7.0e3 -3.0e3 -9.0e3
23+
7.5e3 -3.0e3 -11.0e3
24+
7.5e3 -3.0e3 -10.0e3
25+
7.5e3 -3.0e3 -9.0e3
26+
12.5e3 -3.0e3 -11.0e3
27+
12.5e3 -3.0e3 -10.0e3
28+
12.5e3 -3.0e3 -9.0e3
29+
12.5e3 -3.0e3 -7.0e3

‎quakeworx/cdb_tpv23/tpv23_damMaterial.yaml‎

Lines changed: 46 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,46 @@
1+
!Switch
2+
[rho,mu0,lambda0,xi0,mu,lambda,gamma,epsxx_alpha,epsyy_alpha,epszz_alpha,epsxy_alpha,epsyz_alpha,epszx_alpha]: !ConstantMap
3+
map:
4+
rho: 2760
5+
mu0: 32038120320
6+
lambda0: 32043759360
7+
xi0: 0.75
8+
mu: 32038120320
9+
lambda: 32043759360
10+
gamma: 0.00e10
11+
epsxx_alpha: -4.6921e-5
12+
epsyy_alpha: -5.9238e-4
13+
epszz_alpha: -4.6921e-5
14+
epsxy_alpha: 4.5787e-4
15+
epsyz_alpha: 0.0
16+
epszx_alpha: 0.0
17+
[gammaR,Cd]: !FunctionMap
18+
map:
19+
gammaR: |
20+
depthVar = 1000.0;
21+
dGam = 0.0;
22+
uGam = 37200000000.0;
23+
if (z >= -depthVar) {
24+
return uGam - dGam*(depthVar + z)/depthVar;
25+
}
26+
return uGam;
27+
Cd: |
28+
rx = sqrt(x*x);
29+
rz = sqrt((z+7500)*(z+7500));
30+
rx1 = sqrt((x+7500)*(x+7500));
31+
rx2 = sqrt((x-7500)*(x-7500));
32+
r = max(rx, rz);
33+
depthCd = 1000.0;
34+
CdUni = 0.000005;
35+
CdDropTo = 0.0000005;
36+
if (z > -depthCd) {
37+
return CdUni - (CdUni - CdDropTo)*(depthCd+z)/depthCd;
38+
}
39+
if (max(rx1,rz) <=1800) {
40+
return CdUni;
41+
}
42+
if (max(rx2,rz) <=1800) {
43+
return CdUni;
44+
}
45+
return CdUni;
46+

0 commit comments

Comments
 (0)