FDMOD –声波方程有限差分正演模拟（二维）
格式：
fdmod <vfile >wfile nx= nz= tmax= xs= zs= [optional parameters] 必需的参数：
<vfile 速度输入文件（包含速度值v[nx][nz]）
>wfile 波场输出文件（包含每个时间步的波场值wave[nx][nz]）nx= x采样点个数(第二维)
nz= z采样点个数(第一维)
xs= 炮点x坐标
dxs= 炮点x坐标间隔
zs= 炮点z坐标
dzs= 炮点z坐标间隔
ns= 炮点个数
tmax= 最大记录时间
可选参数：
nt=1+tmax/dt 时间采样点数(dt决定结果的稳定度)
mt=1 波场输出时间切片的时间步长间隔
dx=1.0 x采样间隔
fx=0.0 x起始值
dz=1.0 z采样间隔
fz=0.0 z起始值
fmax = vmin/(10.0*h) 震源子波的最高频率
fpeak=0.5*fmax 雷克子波的峰值频率
dfile= 密度输入文件（包含密度值d[nx][nz]）
vsx= 垂直测线的x坐标
hsz= 水平测线的z坐标
rsx= 水平测线的起始检波器x坐标
rlen= 水平测线长度
rivl= 水平测线检波器采样间隔
vsfile= 垂直测线的输出文件data[nz][nt]
hsfile= 水平测线的输出文件data[nx][nt]
ssfile= 震源点检波器的输出文件data[nt]
verbose=0 =1 显示输出信息=2 更多输出信息
abs=1,1,1,1 模型的顶,底,左,右使用吸收边界条件
=0,1,1,1 顶部使用自由边界条件
PML 参数:
pml_max=1000.0 完全匹配边界条件参数
pml_thick=0 完全匹配半层厚度(0 = 不使用PML边界条件)
注意:
程序使用普通显式二维差分方法。