!Switch
[rho, mu, lambda]: !Include Aqaba_rhomulambda.yaml
[bulkFriction]: !ConstantMap
map:
bulkFriction: 0.6
[plastCo]: !IdentityMap
components:
- !SphericalDomainFilter
radius: 1800
centre:
x: -1540
y: -22115
z: -7200
components: !ConstantMap
map:
plastCo: 1.0e20
- !EvalModel
parameters: [mu]
model: !Switch
[mu]: !Include Aqaba_rhomulambda.yaml
components: !FunctionMap
map:
plastCo: return 0.0001*mu;
[s_xx, s_yy, s_zz, s_xy, s_yz, s_xz]: !Include Aqaba_initial_stress.yaml
!LayeredModel
map: !AffineMap
matrix:
z1: [0.0, 0.0, 1.0]
translation:
z1: 0
interpolation: lower
parameters: [rho, mu, lambda]
nodes:
10000.000000: [2348.000000, 1.865899e+10, 2.138201e+10]
-2500.000000: [2348.000000, 1.865899e+10, 2.138201e+10]
-5000.000000: [2441.000000, 2.174986e+10, 2.311560e+10]
-7500.000000: [2441.000000, 2.174986e+10, 2.311560e+10]
-10000.000000: [2546.000000, 2.560064e+10, 2.722188e+10]
-12500.000000: [2705.000000, 3.228965e+10, 3.433229e+10]
-15000.000000: [2784.000000, 3.602053e+10, 3.828058e+10]
-17500.000000: [2784.000000, 3.602053e+10, 3.828058e+10]
-20000.000000: [2833.000000, 3.844907e+10, 4.085213e+10]
-22500.000000: [2888.000000, 4.133043e+10, 4.394226e+10]
-25000.000000: [2939.000000, 4.410815e+10, 4.688483e+10]
-27500.000000: [3029.000000, 4.929139e+10, 5.235072e+10]
-30000.000000: [3107.000000, 5.413101e+10, 5.753690e+10]
-32500.000000: [3119.000000, 5.488824e+10, 5.839891e+10]
-35000.000000: [3119.000000, 5.488824e+10, 5.839891e+10]
-37500.000000: [3119.000000, 5.488824e+10, 5.839891e+10]
-400000.000000: [3144.000000, 5.652157e+10, 6.000752e+10]
On the other hand, when I compute explicitly PlasticCo beforehand, then SeisSol does not crash...
!Switch
[rho, mu, lambda]: !Include Aqaba_rhomulambda.yaml
[bulkFriction]: !ConstantMap
map:
bulkFriction: 0.6
[plastCo]: !IdentityMap
components:
- !SphericalDomainFilter
radius: 1800
centre:
x: -1540
y: -22115
z: -7200
components: !ConstantMap
map:
plastCo: 1.0e20
- !LayeredModel
map: !AffineMap
matrix:
z1: [0.0, 0.0, 1.0]
translation:
z1: 0
interpolation: lower
parameters: [plastCo]
nodes:
10000.000000: [1.865899e+06]
-2500.000000: [1.865899e+06]
-5000.000000: [2.174986e+06]
-7500.000000: [2.174986e+06]
-10000.000000: [2.560064e+06]
-12500.000000: [3.228965e+06]
-15000.000000: [3.602053e+06]
-17500.000000: [3.602053e+06]
-20000.000000: [3.844907e+06]
-22500.000000: [4.133043e+06]
-25000.000000: [4.410815e+06]
-27500.000000: [4.929139e+06]
-30000.000000: [5.413101e+06]
-32500.000000: [5.488824e+06]
-35000.000000: [5.488824e+06]
-37500.000000: [5.488824e+06]
-400000.000000: [5.652157e+06]
[s_xx, s_yy, s_zz, s_xy, s_yz, s_xz]: !Include Aqaba_initial_stress.yaml
This is also a problem we experienced with the Ridgecrest setup.