diff --git a/docs/source/Gallery/ex13/plot.py b/docs/source/Gallery/ex13/plot.py index d90b5a77..66b37d66 100644 --- a/docs/source/Gallery/ex13/plot.py +++ b/docs/source/Gallery/ex13/plot.py @@ -3,7 +3,7 @@ from scipy.io import netcdf_file def plot_disp(fLst): - fig, axs = plt.subplots(2, 3, figsize=(12, 6)) + fig, axs = plt.subplots(3, 3, figsize=(12, 8)) for i, f in enumerate(fLst): ax3 = axs[i] chs = ['E', 'N', 'Z'] @@ -18,11 +18,12 @@ def plot_disp(fLst): axs[0, 0].set_ylabel("Strike-slip", fontsize=12) axs[1, 0].set_ylabel("Dip-slip", fontsize=12) + axs[2, 0].set_ylabel("Tensile", fontsize=12) return fig, axs def plot_stress(fLst): - fig, axs = plt.subplots(2, 2, figsize=(9,7), gridspec_kw=dict(hspace=0.2)) - for i, f in enumerate([f1, f2]): + fig, axs = plt.subplots(2, 3, figsize=(12, 7), gridspec_kw=dict(hspace=0.2)) + for i, f in enumerate(fLst): ax = axs[0, i] Y = f.variables['north'][:] ax.plot(Y, f.variables['stress_ZZ'][:,0], 'bo-', ms=3, label='ZZ') @@ -38,33 +39,37 @@ def plot_stress(fLst): axs[0, 0].set_title("Strike-slip", fontsize=12) axs[0, 1].set_title("Dip-slip", fontsize=12) + axs[0, 2].set_title("Tensile", fontsize=12) axs[0, 0].set_ylabel("Normal", fontsize=12) axs[1, 0].set_ylabel("Shear", fontsize=12) - axs[0, 1].legend(bbox_to_anchor=(1.3, 1)) - axs[1, 1].legend(bbox_to_anchor=(1.3, 1)) + axs[0, 2].legend(bbox_to_anchor=(1.3, 1)) + axs[1, 2].legend(bbox_to_anchor=(1.3, 1)) return fig, axs # ====================================================================== with netcdf_file("stsyn_1_stkslip.nc", mmap=False) as f1, \ - netcdf_file("stsyn_1_dipslip.nc", mmap=False) as f2: - fig, axs = plot_disp([f1, f2]) + netcdf_file("stsyn_1_dipslip.nc", mmap=False) as f2, \ + netcdf_file("stsyn_1_tensile.nc", mmap=False) as f3: + fig, axs = plot_disp([f1, f2, f3]) fig.savefig("disp1.svg", bbox_inches='tight') - fig, axs = plot_stress([f1, f2]) + fig, axs = plot_stress([f1, f2, f3]) fig.savefig("stress1.svg", bbox_inches='tight') with netcdf_file("stsyn_2_stkslip.nc", mmap=False) as f1, \ - netcdf_file("stsyn_2_dipslip.nc", mmap=False) as f2: - fig, axs = plot_disp([f1, f2]) + netcdf_file("stsyn_2_dipslip.nc", mmap=False) as f2, \ + netcdf_file("stsyn_2_tensile.nc", mmap=False) as f3: + fig, axs = plot_disp([f1, f2, f3]) fig.savefig("disp2.svg", bbox_inches='tight') - fig, axs = plot_stress([f1, f2]) + fig, axs = plot_stress([f1, f2, f3]) fig.savefig("stress2.svg", bbox_inches='tight') with netcdf_file("stsyn_3_stkslip.nc", mmap=False) as f1, \ - netcdf_file("stsyn_3_dipslip.nc", mmap=False) as f2: - fig, axs = plot_disp([f1, f2]) + netcdf_file("stsyn_3_dipslip.nc", mmap=False) as f2, \ + netcdf_file("stsyn_3_tensile.nc", mmap=False) as f3: + fig, axs = plot_disp([f1, f2, f3]) fig.savefig("disp3.svg", bbox_inches='tight') - fig, axs = plot_stress([f1, f2]) + fig, axs = plot_stress([f1, f2, f3]) fig.savefig("stress3.svg", bbox_inches='tight') \ No newline at end of file diff --git a/docs/source/Gallery/ex13/run.sh b/docs/source/Gallery/ex13/run.sh index 2e5c3d2c..b6624d18 100755 --- a/docs/source/Gallery/ex13/run.sh +++ b/docs/source/Gallery/ex13/run.sh @@ -13,6 +13,9 @@ grt static stress stsyn_1_stkslip.nc grt static syn -M90/70/90 -Su1e6 -e -N -Gstgrn1.nc -Ostsyn_1_dipslip.nc grt static stress stsyn_1_dipslip.nc +grt static syn -M90/70 -Su1e6 -e -N -Gstgrn1.nc -Ostsyn_1_tensile.nc +grt static stress stsyn_1_tensile.nc + grt static greenfn -M$mod -D10/14 -Y3/3/1 -X-10/10/0.5 -e -Ostgrn2.nc @@ -22,6 +25,8 @@ grt static stress stsyn_2_stkslip.nc grt static syn -M90/70/90 -Su1e6 -e -N -Gstgrn2.nc -Ostsyn_2_dipslip.nc grt static stress stsyn_2_dipslip.nc +grt static syn -M90/70 -Su1e6 -e -N -Gstgrn2.nc -Ostsyn_2_tensile.nc +grt static stress stsyn_2_tensile.nc mod="mod" grt static greenfn -M$mod -D5/0 -Y3/3/1 -X-10/10/0.5 -e -Ostgrn3.nc @@ -32,6 +37,8 @@ grt static stress stsyn_3_stkslip.nc grt static syn -M90/70/90 -Su1e6 -e -N -Gstgrn3.nc -Ostsyn_3_dipslip.nc grt static stress stsyn_3_dipslip.nc +grt static syn -M90/70 -Su1e6 -e -N -Gstgrn3.nc -Ostsyn_3_tensile.nc +grt static stress stsyn_3_tensile.nc python plot.py diff --git a/docs/source/Gallery/ex14/run.sh b/docs/source/Gallery/ex14/run.sh old mode 100644 new mode 100755 diff --git a/docs/source/Tutorial/dynamic/run/run.py b/docs/source/Tutorial/dynamic/run/run.py index dd880c9f..6ad3d1ea 100755 --- a/docs/source/Tutorial/dynamic/run/run.py +++ b/docs/source/Tutorial/dynamic/run/run.py @@ -171,6 +171,23 @@ def plot_int_dif(stsyn:Stream, stsyn_int:Stream, stsyn_dif:Stream, comp:str, out # ----------------------------------------------------------------------------------- +# ----------------------------------------------------------------------------------- +# BEGIN SYN TS +# 接之前的代码 +idx = 2 +stgrn = stgrnLst[idx] # 选择格林函数 + +stsyn = pygrt.utils.gen_syn_from_gf_TS(stgrn, M0=1e24, strike=33, dip=50, az=30) +print(stsyn) +# 3 Trace(s) in Stream: +# .SYN..Z | 1970-01-01T00:00:00.000000Z - 1970-01-01T00:00:09.980000Z | 50.0 Hz, 500 samples +# .SYN..R | 1970-01-01T00:00:00.000000Z - 1970-01-01T00:00:09.980000Z | 50.0 Hz, 500 samples +# .SYN..T | 1970-01-01T00:00:00.000000Z - 1970-01-01T00:00:09.980000Z | 50.0 Hz, 500 samples +plot_syn(stsyn, "syn_ts.svg") +# END SYN TS +# ----------------------------------------------------------------------------------- + + # ----------------------------------------------------------------------------------- # BEGIN SYN MT idx = 2 diff --git a/docs/source/Tutorial/dynamic/run/run.sh b/docs/source/Tutorial/dynamic/run/run.sh index e4bffc46..d9418148 100755 --- a/docs/source/Tutorial/dynamic/run/run.sh +++ b/docs/source/Tutorial/dynamic/run/run.sh @@ -48,6 +48,13 @@ grt syn -GGRN/milrow_2_0_10 -S1e24 -A30 -M33/50/120 -Osyn_dc # END SYN DC # ----------------------------------------------------------------------------------- +# ----------------------------------------------------------------------------------- +# BEGIN SYN TS +# 合成结果在 syn_ts/ 目录下，以SAC格式保存。 +grt syn -GGRN/milrow_2_0_10 -S1e24 -A30 -M33/50 -Osyn_ts +# END SYN TS +# ----------------------------------------------------------------------------------- + # ----------------------------------------------------------------------------------- # BEGIN SYN MT # 合成结果在 syn_mt/ 目录下，以SAC格式保存。 diff --git a/docs/source/Tutorial/dynamic/syn.rst b/docs/source/Tutorial/dynamic/syn.rst index a8ad7fbe..0c1f833c 100644 --- a/docs/source/Tutorial/dynamic/syn.rst +++ b/docs/source/Tutorial/dynamic/syn.rst @@ -126,6 +126,31 @@ Python中合成动态位移的主函数为 :func:`gen_syn_from_gf_*() M.active = true; { diff --git a/pygrt/C_extension/src/static/grt_static_syn.c b/pygrt/C_extension/src/static/grt_static_syn.c index 5b39b466..66fd3510 100644 --- a/pygrt/C_extension/src/static/grt_static_syn.c +++ b/pygrt/C_extension/src/static/grt_static_syn.c @@ -209,7 +209,7 @@ static void getopt_from_command(GRT_MODULE_CTRL *Ctrl, int argc, char **argv){ }; break; - // 剪切震源， 张位错源 + // 剪切震源， 张裂源 case 'M': Ctrl->M.active = true; { diff --git a/pygrt/utils.py b/pygrt/utils.py index ad2afbe2..65cc6359 100755 --- a/pygrt/utils.py +++ b/pygrt/utils.py @@ -4,7 +4,7 @@ :date: 2024-07-24 该文件包含一些数据处理操作上的补充: - 1、剪切源、张位错源、单力源、爆炸源、矩张量源 通过格林函数合成理论地震图的函数\n + 1、剪切源、张裂源、单力源、爆炸源、矩张量源 通过格林函数合成理论地震图的函数\n 2、Stream类型的时域卷积、微分、积分 (基于numpy和scipy) \n 3、Stream类型写到本地sac文件，自定义名称 \n 4、读取波数积分和峰谷平均法过程文件 \n @@ -96,7 +96,7 @@ def _gen_syn_from_gf(st:Stream, calc_upar:bool, compute_type:GRT_SYN_TYPE, M0:fl srcName = ["EX", "VF", "HF", "DD", "DS", "SS"] allchs = [tr.stats.channel for tr in st] - # 为张位错计算 Vp/Vs + # 为张裂计算 Vp/Vs src_va = st[0].stats.sac['user6'] src_vb = st[0].stats.sac['user7'] VpVs_ratio = src_va / src_vb @@ -168,7 +168,7 @@ def _gen_syn_from_static_gf(grn:dict, calc_upar:bool, compute_type:GRT_SYN_TYPE, srcName = ["EX", "VF", "HF", "DD", "DS", "SS"] allchs = list(grn.keys()) - # 为张位错计算 Vp/Vs + # 为张裂计算 Vp/Vs src_va = grn['_src_va'] src_vb = grn['_src_vb'] VpVs_ratio = src_va / src_vb diff --git a/test/_compare_c_py/compare.py b/test/_compare_c_py/compare.py index 7efbc583..a2bd130a 100644 --- a/test/_compare_c_py/compare.py +++ b/test/_compare_c_py/compare.py @@ -198,6 +198,17 @@ def plot_static(resDct:dict, c_prefix:str): AVGRERR.append(compare3(strotation, f"syn_dc{suffix}/rotation_", ZNE=ZNE, dim2=True)) AVGRERR.append(compare3(ststress, f"syn_dc{suffix}/stress_", ZNE=ZNE, dim2=True)) + st = pygrt.utils.gen_syn_from_gf_TS(st_grn, S, stk, dip, az, ZNE=ZNE, calc_upar=True) + sigs = pygrt.sigs.gen_parabola_wave(0.6, dt) + pygrt.utils.stream_convolve(st, sigs) + AVGRERR.append(compare3(st, f"syn_ts{suffix}/", ZNE=ZNE)) + ststrain = pygrt.utils.compute_strain(st) + strotation = pygrt.utils.compute_rotation(st) + ststress = pygrt.utils.compute_stress(st) + AVGRERR.append(compare3(ststrain, f"syn_ts{suffix}/strain_", ZNE=ZNE, dim2=True)) + AVGRERR.append(compare3(strotation, f"syn_ts{suffix}/rotation_", ZNE=ZNE, dim2=True)) + AVGRERR.append(compare3(ststress, f"syn_ts{suffix}/stress_", ZNE=ZNE, dim2=True)) + st = pygrt.utils.gen_syn_from_gf_MT(st_grn, S, [M11,M12,M13,M22,M23,M33], az, ZNE=ZNE, calc_upar=True) sigs = pygrt.sigs.gen_ricker_wave(3, dt) pygrt.utils.stream_convolve(st, sigs) @@ -248,6 +259,15 @@ def plot_static(resDct:dict, c_prefix:str): update_dict(static_syn, strotation, "rotation_") AVGRERR2.append(static_compare3(static_syn, f"static/stsyn_dc{suffix}.nc")) + static_syn = pygrt.utils.gen_syn_from_gf_TS(static_grn, S, stk, dip, ZNE=ZNE, calc_upar=True) + ststrain = pygrt.utils.compute_strain(static_syn) + strotation = pygrt.utils.compute_rotation(static_syn) + ststress = pygrt.utils.compute_stress(static_syn) + update_dict(static_syn, ststrain, "strain_") + update_dict(static_syn, ststress, "stress_") + update_dict(static_syn, strotation, "rotation_") + AVGRERR2.append(static_compare3(static_syn, f"static/stsyn_ts{suffix}.nc")) + static_syn = pygrt.utils.gen_syn_from_gf_MT(static_grn, S, [M11,M12,M13,M22,M23,M33], ZNE=ZNE, calc_upar=True) ststrain = pygrt.utils.compute_strain(static_syn) strotation = pygrt.utils.compute_rotation(static_syn) diff --git a/test/_compare_c_py/test_compare_c_py.sh b/test/_compare_c_py/test_compare_c_py.sh index e74e72ed..c489b8d8 100755 --- a/test/_compare_c_py/test_compare_c_py.sh +++ b/test/_compare_c_py/test_compare_c_py.sh @@ -45,6 +45,11 @@ grt strain syn_dc$N grt rotation syn_dc$N grt stress syn_dc$N +grt syn -G$G -Osyn_ts$N -A$az -S$S -M$stk/$dip -Dp/0.6 -e $N +grt strain syn_ts$N +grt rotation syn_ts$N +grt stress syn_ts$N + M11=1 M12=-2 M13=-5 @@ -89,6 +94,11 @@ grt static strain stsyn_dc$N.nc grt static rotation stsyn_dc$N.nc grt static stress stsyn_dc$N.nc +grt static syn -S$S -M$stk/$dip -e $N -Gstgrn.nc -Ostsyn_ts$N.nc +grt static strain stsyn_ts$N.nc +grt static rotation stsyn_ts$N.nc +grt static stress stsyn_ts$N.nc + grt static syn -S$S -T$M11/$M12/$M13/$M22/$M23/$M33 -e $N -Gstgrn.nc -Ostsyn_mt$N.nc grt static strain stsyn_mt$N.nc grt static rotation stsyn_mt$N.nc diff --git a/test/_correct_full_wave/_Ref/syn_ts/R.sac b/test/_correct_full_wave/_Ref/syn_ts/R.sac new file mode 100644 index 00000000..023d75d0 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/R.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/T.sac b/test/_correct_full_wave/_Ref/syn_ts/T.sac new file mode 100644 index 00000000..b342e695 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/T.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/Z.sac b/test/_correct_full_wave/_Ref/syn_ts/Z.sac new file mode 100644 index 00000000..b9d7c0aa Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/Z.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/rR.sac b/test/_correct_full_wave/_Ref/syn_ts/rR.sac new file mode 100644 index 00000000..e1067618 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/rR.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/rT.sac b/test/_correct_full_wave/_Ref/syn_ts/rT.sac new file mode 100644 index 00000000..9a1a731b Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/rT.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/rZ.sac b/test/_correct_full_wave/_Ref/syn_ts/rZ.sac new file mode 100644 index 00000000..0f2c391f Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/rZ.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/sig.sac b/test/_correct_full_wave/_Ref/syn_ts/sig.sac new file mode 100644 index 00000000..5897aeb3 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/sig.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/tR.sac b/test/_correct_full_wave/_Ref/syn_ts/tR.sac new file mode 100644 index 00000000..f3121ef4 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/tR.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/tT.sac b/test/_correct_full_wave/_Ref/syn_ts/tT.sac new file mode 100644 index 00000000..b8c7478e Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/tT.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/tZ.sac b/test/_correct_full_wave/_Ref/syn_ts/tZ.sac new file mode 100644 index 00000000..bab873d1 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/tZ.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/zR.sac b/test/_correct_full_wave/_Ref/syn_ts/zR.sac new file mode 100644 index 00000000..e39fed87 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/zR.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/zT.sac b/test/_correct_full_wave/_Ref/syn_ts/zT.sac new file mode 100644 index 00000000..e6b782d1 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/zT.sac differ diff --git a/test/_correct_full_wave/_Ref/syn_ts/zZ.sac b/test/_correct_full_wave/_Ref/syn_ts/zZ.sac new file mode 100644 index 00000000..9c5bfde3 Binary files /dev/null and b/test/_correct_full_wave/_Ref/syn_ts/zZ.sac differ diff --git a/test/_correct_full_wave/correctness.sh b/test/_correct_full_wave/correctness.sh index 73ea4d80..5143740d 100755 --- a/test/_correct_full_wave/correctness.sh +++ b/test/_correct_full_wave/correctness.sh @@ -24,12 +24,14 @@ python ../compare_sac.py "GRN/seafloor_5_1.1_100/*" "_Ref/seafloor_5_1.1_100 grt syn -GGRN/milrow_2_3_10 -A22 -S1e20 -Dt/0.2/0.2/0.4 -e -Osyn_ex grt syn -GGRN/milrow_2_3_10 -A22 -S1e20 -F2/-1/4 -Dt/0.1/0.3/0.6 -e -Osyn_sf grt syn -GGRN/milrow_2_3_10 -A22 -S1e20 -M77/88/111 -Dp/0.6 -e -Osyn_dc +grt syn -GGRN/milrow_2_3_10 -A22 -S1e20 -M77/88 -Dp/0.6 -e -Osyn_ts grt syn -GGRN/milrow_2_3_10 -A22 -S1e20 -T1/-2/-5/0.5/3/1.2 -Dr/3 -e -Osyn_mt # also check the time functions python ../compare_sac.py "syn_ex/*" "_Ref/syn_ex/*" python ../compare_sac.py "syn_sf/*" "_Ref/syn_sf/*" python ../compare_sac.py "syn_dc/*" "_Ref/syn_dc/*" +python ../compare_sac.py "syn_ts/*" "_Ref/syn_ts/*" python ../compare_sac.py "syn_mt/*" "_Ref/syn_mt/*" # rotate to ZNE diff --git a/test/_correct_static_disp/_Ref/stsyn_ts.nc b/test/_correct_static_disp/_Ref/stsyn_ts.nc new file mode 100644 index 00000000..07165bfc Binary files /dev/null and b/test/_correct_static_disp/_Ref/stsyn_ts.nc differ diff --git a/test/_correct_static_disp/correctness.sh b/test/_correct_static_disp/correctness.sh index 1862674a..1ec3c7a4 100755 --- a/test/_correct_static_disp/correctness.sh +++ b/test/_correct_static_disp/correctness.sh @@ -12,6 +12,7 @@ grt static greenfn -M../milrow -D2/3 -X-3.1/3.1/0.2 -Y-2.1/2.1/0.2 -Cp -e -Ostgr grt static syn -S1e20 -e -Gstgrn.nc -Ostsyn_ex.nc grt static syn -S1e20 -e -F2/-1/4 -Gstgrn.nc -Ostsyn_sf.nc grt static syn -S1e20 -e -M77/88/111 -Gstgrn.nc -Ostsyn_dc.nc +grt static syn -S1e20 -e -M77/88 -Gstgrn.nc -Ostsyn_ts.nc grt static syn -S1e20 -e -T1/-2/-5/0.5/3/1.2 -Gstgrn.nc -Ostsyn_mt.nc # rotate to ZNE @@ -32,6 +33,7 @@ python ../compare_nc.py stgrn.nc _Ref/stgrn.nc python ../compare_nc.py stsyn_ex.nc _Ref/stsyn_ex.nc python ../compare_nc.py stsyn_sf.nc _Ref/stsyn_sf.nc python ../compare_nc.py stsyn_dc.nc _Ref/stsyn_dc.nc +python ../compare_nc.py stsyn_ts.nc _Ref/stsyn_ts.nc python ../compare_nc.py stsyn_mt.nc _Ref/stsyn_mt.nc python ../compare_nc.py stsyn_mt_ZNE.nc _Ref/stsyn_mt_ZNE.nc diff --git a/test/static_syn/test_static_syn.py b/test/static_syn/test_static_syn.py index d15af83e..e5d6aef6 100644 --- a/test/static_syn/test_static_syn.py +++ b/test/static_syn/test_static_syn.py @@ -16,6 +16,7 @@ stsyn = pygrt.utils.gen_syn_from_gf_EX(stgrn, 1e20, 22) stsyn = pygrt.utils.gen_syn_from_gf_SF(stgrn, 1e16, fN=-1, fE=2, fZ=-4, az=22) stsyn = pygrt.utils.gen_syn_from_gf_DC(stgrn, 1e20, strike=33, dip=44, rake=55, az=22) +stsyn = pygrt.utils.gen_syn_from_gf_TS(stgrn, 1e20, strike=33, dip=44, az=22) stsyn = pygrt.utils.gen_syn_from_gf_MT(stgrn, 1e20, [1, -2, -5, 0.5, 3, 1.2], az=22) diff --git a/test/static_syn/test_static_syn.sh b/test/static_syn/test_static_syn.sh index b5c9f438..5afed0ca 100755 --- a/test/static_syn/test_static_syn.sh +++ b/test/static_syn/test_static_syn.sh @@ -11,6 +11,7 @@ grt static syn -S1e20 -Gstgrn.nc -Ostsyn.nc grt static syn -S1e20 -F2/-1/4 -Gstgrn.nc -Ostsyn.nc grt static syn -S1e20 -M77/88/111 -Gstgrn.nc -Ostsyn.nc grt static syn -Su1e6 -M77/88/111 -Gstgrn.nc -Ostsyn.nc +grt static syn -Su1e6 -M77/88 -Gstgrn.nc -Ostsyn.nc grt static syn -S1e20 -T1/-2/-5/0.5/3/1.2 -Gstgrn.nc -Ostsyn.nc grt static syn -S1e20 -F2/-1/4 -e -Gstgrn.nc -Ostsyn.nc diff --git a/test/syn/test_syn.py b/test/syn/test_syn.py index 9bf9ea27..880bef45 100644 --- a/test/syn/test_syn.py +++ b/test/syn/test_syn.py @@ -18,6 +18,7 @@ stsyn = pygrt.utils.gen_syn_from_gf_EX(stgrn, 1e20, 22) stsyn = pygrt.utils.gen_syn_from_gf_SF(stgrn, 1e16, fN=-1, fE=2, fZ=-4, az=22) stsyn = pygrt.utils.gen_syn_from_gf_DC(stgrn, 1e20, strike=33, dip=44, rake=55, az=22) +stsyn = pygrt.utils.gen_syn_from_gf_TS(stgrn, 1e20, strike=33, dip=44, az=22) stsyn = pygrt.utils.gen_syn_from_gf_MT(stgrn, 1e20, [1, -2, -5, 0.5, 3, 1.2], az=22) diff --git a/test/syn/test_syn.sh b/test/syn/test_syn.sh index 60960841..1e8b7904 100755 --- a/test/syn/test_syn.sh +++ b/test/syn/test_syn.sh @@ -10,6 +10,7 @@ grt syn -GGRN/milrow_2_0_10 -A22 -S1e20 -Osyn grt syn -GGRN/milrow_2_0_10 -A22 -S1e16 -F-1/2/-4 -Osyn grt syn -GGRN/milrow_2_0_10 -A22 -S1e20 -M33/44/55 -Osyn grt syn -GGRN/milrow_2_0_10 -A22 -Su1e10 -M33/44/55 -Osyn +grt syn -GGRN/milrow_2_0_10 -A22 -Su1e10 -M33/44 -Osyn grt syn -GGRN/milrow_2_0_10 -A22 -S1e20 -T1/-2/-5/0.5/3/1.2 -Osyn grt syn -GGRN/milrow_2_0_10 -A22 -S1e20 -Dp/0.6 -Osyn