Synopsis
		sfefd3dmt < Fvp.rsf vs=Fvs.rsf rho=Frho.rsf > Frvz.rsf rvx=Frvx.rsf rvy=Frvy.rsf wavx=Fwavx.rsf wavy=Fwavy.rsf wavz=Fwavz.rsf M= verb=n ifwfd=n nb=30 nt= dt= fm=20.0 ft=0 jsnap=1 ct=0.01 sx=nxpad/2 sy=nypad/2 sz=nzpad/2 free=n

By Yangkang Chen, 2020
Revised in July, 2021
Currently still in a draft version

DEMO:
https://github.com/chenyk1990/tutorials/blob/main/demo/efd3dmt/SConstruct

Parameters
		floats M= [9]   float ct=0.01 for absorbing boundary   float dt= time sampling interval   float fm=20.0 dominant freq of Ricker wavelet   bool free=n [y/n] if free surface   int ft=0 first recorded time   bool ifwfd=n [y/n] if output wavefield   int jsnap=1 interval for snapshots   int nb=30 thickness of ABC layer   int nt= number of time steps   file rho= auxiliary input file name   file rvx= auxiliary output file name   file rvy= auxiliary output file name   int sx=nxpad/2   int sy=nypad/2   int sz=nzpad/2   bool verb=n [y/n] verbosity   file vs= auxiliary input file name   file wavx= auxiliary output file name   file wavy= auxiliary output file name   file wavz= auxiliary output file name