From bb80e0d2ba0b025e76f809104230b7844d25b8fa Mon Sep 17 00:00:00 2001 From: Thomas-Ulrich Date: Tue, 27 Jan 2026 14:28:42 +0100 Subject: [PATCH 1/3] fix typo --- sulawesi/Sulawesi_faultLSW.yaml | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/sulawesi/Sulawesi_faultLSW.yaml b/sulawesi/Sulawesi_faultLSW.yaml index 6aa8f6e..5bf7993 100644 --- a/sulawesi/Sulawesi_faultLSW.yaml +++ b/sulawesi/Sulawesi_faultLSW.yaml @@ -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 From a7159396b453717c4b203118841ee6a038c26e2b Mon Sep 17 00:00:00 2001 From: Thomas-Ulrich Date: Tue, 27 Jan 2026 14:49:27 +0100 Subject: [PATCH 2/3] recenter receivers --- sulawesi/GPSfiltered_proj_200m_8e6.dat | 61 ++++++++++++-------------- 1 file changed, 28 insertions(+), 33 deletions(-) diff --git a/sulawesi/GPSfiltered_proj_200m_8e6.dat b/sulawesi/GPSfiltered_proj_200m_8e6.dat index 6a52b85..a2c7674 100644 --- a/sulawesi/GPSfiltered_proj_200m_8e6.dat +++ b/sulawesi/GPSfiltered_proj_200m_8e6.dat @@ -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 From a098b26fbb9e36a512939e08916fffe640531d9c Mon Sep 17 00:00:00 2001 From: Thomas-Ulrich Date: Tue, 27 Jan 2026 14:51:06 +0100 Subject: [PATCH 3/3] update parameters --- sulawesi/faultreceiver_hypocenter.dat | 2 +- sulawesi/parametersLSW.par | 6 +++--- sulawesi/parametersRS.par | 5 +++-- 3 files changed, 7 insertions(+), 6 deletions(-) diff --git a/sulawesi/faultreceiver_hypocenter.dat b/sulawesi/faultreceiver_hypocenter.dat index f3c4c95..5771897 100644 --- a/sulawesi/faultreceiver_hypocenter.dat +++ b/sulawesi/faultreceiver_hypocenter.dat @@ -1 +1 @@ -17690.222084 1475585.606433 -9991.459379 +17690.222084 55585.6 -9991.459379 diff --git a/sulawesi/parametersLSW.par b/sulawesi/parametersLSW.par index f2a5292..315f624 100644 --- a/sulawesi/parametersLSW.par +++ b/sulawesi/parametersLSW.par @@ -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' / @@ -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 @@ -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. diff --git a/sulawesi/parametersRS.par b/sulawesi/parametersRS.par index 753e173..e9c202f 100644 --- a/sulawesi/parametersRS.par +++ b/sulawesi/parametersRS.par @@ -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' @@ -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.