该测试已考虑到 #215 的问题,手动提高了每个频率的积分上限。
但当震中距过小时,Bessel函数振荡缓慢,低频段的积分上限不够大,且积分步长不够小,形成信号尾端的低频假象。应该手动/自动再适当微调 k0。
#!/bin/bash
rm halfspace* -rf
cat > mod0 <<EOF
1.0 0.7 0.3 1.3 80 20
EOF
grt greenfn -Mmod0 -N2000/0.001 -D0/0 -Ohalfspace1 -R0.2 -Gh -K+s5 -Cd # -K+s5+k520 -L30
# grt greenfn -Mmod0 -N2000/0.1+a -D0/0 -Ohalfspace1 -R20 -Gh -K+s5 -Cd
sac <<EOF
r halfspace1/*/*;
p1
saveimg test1.pdf
q
EOF
pdf2svg test1.pdf test1.svg
