close depths requires larger kmax

[Dengda98]

这一点在大震中距情况下更明显，会出现很多伪影，kmax不够大会导致DCM或PTAM这些收敛算法结果偏差较大。
此时大震中距也需要精度更高的SAFIM

[Dengda98]

一个简单的脚本可以复现这个现象，当源点和场点接近时，会需要更大的kmax来消除伪影。而当震中距较大时，时窗会比较大，则虚频率变小，核函数在极点附近变得更尖锐，这要求更小的积分步长（尽管SAFIM能很好地解决这个问题，但这也要求更严格的精度），进一步形成伪影。这是造成伪影形成的多个因素。

#!/bin/bash 

set -euo pipefail

seq -f "%.2f" 0.5 0.01 1.5 > dists
distarr=($(cat dists | tr '\n' ','))
distarr=${distarr%,} 

depsrc=0.1
deprcv=0

nt=400
dt=0.02 

modname="mod"
out="GRN"

cat > $modname <<EOF
0.01    1.5   1.0   1.78
EOF

# compute Green's Functions
grt greenfn -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/${deprcv} -Rdists -Gv -K+s5 # 只生成垂直力源的格林函数

gmt sac ${out}/*/VFZ.sac -R0/8/0.4/1.6 -JX10c/10c -BWSen -Baf -Ek -M0.5c -W0p,red -pdf test