Skip to content

Thomas/bug fix typo #59

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
Thomas-Ulrich wants to merge 3 commits into
base: main
Choose a base branch
from
Draft

Thomas/bug fix typo #59

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
bb80e0d
fix typo
Thomas-Ulrich Jan 27, 2026
a715939
recenter receivers
Thomas-Ulrich Jan 27, 2026
a098b26
update parameters
Thomas-Ulrich Jan 27, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
61 changes: 28 additions & 33 deletions sulawesi/GPSfiltered_proj_200m_8e6.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,33 +1,28 @@
-3.5682315562e+04 1.2947634529e+06 2.3364754055e+01
1.4891031383e+04 1.3997871541e+06 5.3135325532e+01
1.9676785970e+04 1.4041030344e+06 -2.4679119913e+00
-1.2948099617e+04 1.4033134304e+06 -3.0757851155e+00
-7.6071372120e+03 1.4120552909e+06 -2.7815998650e+01
1.8540652680e+03 1.4216831598e+06 8.9499813056e+01
1.6942437298e+04 1.4213640950e+06 -1.8722838289e+01
-2.4879512780e+03 1.4195797035e+06 -1.9140854129e+01
4.8133998940e+03 1.4261977328e+06 -1.2959927733e+01
1.8231357573e+04 1.3995674563e+06 -9.1962562115e+00
1.4891031383e+04 1.3997871541e+06 5.3135325532e+01
2.2573966510e+04 1.3986845421e+06 5.3961913222e+01
1.2765025073e+05 1.5500883932e+06 1.9185553624e+02
9.4286419900e+03 1.4148280000e+06 1.4683454617e+01
2.1452790630e+04 1.4179296746e+06 1.0294644244e+02
8.1342568951e+04 1.5533007705e+06 4.6987668277e+01
1.8113136798e+04 1.4161587132e+06 -4.5656602629e+00
1.9671013439e+04 1.4187032806e+06 4.1189303023e+01
-2.7047341036e+04 1.3421252689e+06 7.5147612435e+01
3.2363403567e+04 1.4201455474e+06 7.4825818669e+02
3.4812869614e+04 1.4200357522e+06 8.5581435311e+02
2.7241561816e+04 1.4208074501e+06 4.0082389229e+02
8.0688740855e+04 1.3886431063e+06 4.3728419345e+01
6.8241902166e+04 1.3425238165e+06 1.1536835130e+03
2.9802460920e+04 1.4205871047e+06 6.8846728074e+02
2.2573966510e+04 1.3986845421e+06 5.3961913222e+01
3.4189114360e+04 1.3327711258e+06 7.0854827580e+02
3.3576539974e+04 1.4725709206e+06 -4.6693286695e+00
3.7930148273e+04 1.4204791590e+06 5.2945514929e+02
1.3369555128e+05 1.2890331143e+06 2.7526418650e+02
4.3607873742e+04 1.4215869069e+06 1.2837494666e+00
4.6602567175e+04 1.5140456184e+06 -5.4256807095e+01
4.3607873742e+04 1.4215869069e+06 1.2837494666e+00
-3.5682315562e+04 -125237 2.3364754055e+01
1.4891031383e+04 -20212.8 5.3135325532e+01
1.9676785970e+04 -15897 -2.4679119913e+00
-1.2948099617e+04 -16686.6 -3.0757851155e+00
-7.6071372120e+03 -7944.71 -2.7815998650e+01
1.8540652680e+03 1683.16 8.9499813056e+01
1.6942437298e+04 1364.09 -1.8722838289e+01
-2.4879512780e+03 -420.296 -1.9140854129e+01
4.8133998940e+03 6197.73 -1.2959927733e+01
1.8231357573e+04 -20432.5 -9.1962562115e+00
1.4891031383e+04 -20212.8 5.3135325532e+01
2.2573966510e+04 -21315.5 5.3961913222e+01
9.4286419900e+03 -5172 1.4683454617e+01
2.1452790630e+04 -2070.33 1.0294644244e+02
1.8113136798e+04 -3841.29 -4.5656602629e+00
1.9671013439e+04 -1296.72 4.1189303023e+01
-2.7047341036e+04 -77874.7 7.5147612435e+01
3.2363403567e+04 145.547 7.4825818669e+02
3.4812869614e+04 35.7522 8.5581435311e+02
2.7241561816e+04 807.45 4.0082389229e+02
8.0688740855e+04 -31356.9 4.3728419345e+01
2.2573966510e+04 -21315.5 5.3961913222e+01
3.4189114360e+04 -87228.9 7.0854827580e+02
3.3576539974e+04 52570.9 -4.6693286695e+00
1.3369555128e+05 -130967 2.7526418650e+02
4.3607873742e+04 1586.91 1.2837494666e+00
4.6602567175e+04 94045.6 -5.4256807095e+01
4.3607873742e+04 1586.91 1.2837494666e+00
2 changes: 1 addition & 1 deletion sulawesi/Sulawesi_faultLSW.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@
r_crit = 4000.0;
Vs = 3464.0;
if (r <= r_crit) {
return r/(0.7*Vs)+(0.081.0*r_crit/(0.7.0*Vs))*(1.0/(1.0-pow(r/r_crit, 2.0))-1.0);
return r/(0.7*Vs)+(0.081*r_crit/(0.7*Vs))*(1.0/(1.0-pow(r/r_crit, 2.0))-1.0);
}
return 1000000000.0;
[s_xx, s_yy, s_zz, s_xy, s_yz, s_xz]: !Include Sulawesi_initial_stressLSW.yaml
Expand Down
2 changes: 1 addition & 1 deletion sulawesi/faultreceiver_hypocenter.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1 +1 @@
17690.222084 1475585.606433 -9991.459379
17690.222084 55585.6 -9991.459379
6 changes: 3 additions & 3 deletions sulawesi/parametersLSW.par
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,9 +45,8 @@ refinement = 1

! parameterize ascii fault file outputs
&Pickpoint
printtimeinterval = 1 ! Index of printed info at timesteps
printtimeinterval = 30 ! Index of printed info at timesteps
OutputMask = 1 1 1 1 1 1 1 1 1 1 1 ! turn on and off fault outputs
nOutpoints = 1
PPFileName = 'faultreceiver_hypocenter.dat'
/

Expand All @@ -70,7 +69,7 @@ ClusteredLTS = 2 ! 1 for Global time stepping, 2,3

&Output
OutputFile = 'output/SulawesiLSW'
Format = 10 ! Format (10= no output, 6=hdf5 output)
wavefieldoutput = 1
! |stress |vel
iOutputMask = 0 0 0 0 0 0 1 1 1
! |strain |eta
Expand All @@ -79,6 +78,7 @@ TimeInterval = 5. ! Index of printed info at time
refinement = 2

! off-fault ascii receivers
ReceiverOutput = 1
RFileName = 'GPSfiltered_proj_200m_8e6.dat' ! Record Points in extra file
pickdt = 0.01 ! Pickpoint Sampling
! (Optional) Synchronization point for receivers.
Expand Down
5 changes: 3 additions & 2 deletions sulawesi/parametersRS.par
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -50,7 +50,7 @@ refinement = 1

! parameterize ascii fault file outputs
&Pickpoint
printtimeinterval = 1 ! Index of printed info at timesteps
printtimeinterval = 30 ! Index of printed info at timesteps
OutputMask = 1 1 1 1 1 1 1 1 1 1 1 ! turn on and off fault outputs
nOutpoints = 1
PPFileName = 'faultreceiver_hypocenter.dat'
Expand All @@ -75,13 +75,14 @@ ClusteredLTS = 2 ! 1 for Global time stepping, 2,3

&Output
OutputFile = 'output/Sulawesi'
Format = 10 ! Format (10= no output, 6=hdf5 output)
wavefieldoutput = 1
! |stress |vel
iOutputMask = 0 0 0 0 0 0 1 1 1
TimeInterval = 5. ! Index of printed info at time
refinement = 2

! off-fault ascii receivers
ReceiverOutput = 1
RFileName = 'GPSfiltered_proj_200m_8e6.dat' ! Record Points in extra file
pickdt = 0.01 ! Pickpoint Sampling
! (Optional) Synchronization point for receivers.
Expand Down
Loading