Program sfefd3dmt | 3D 8-th order elastic wave propagation with sponge ABC and moment tensor source
[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
[COMMENTS]
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
[DIRECTORY]
user/chenyk
