sfvweks3d (2.0-git)
index
user/jsun/Mvweks3d.c
3D visco-elastic time-domain pseudo-spectral (k-space) modeling using shared-memory parallel FFT

 
Synopsis
        sfvweks3d < Fwav.rsf ccc=Fccc.rsf den=Fden.rsf sou=Fsou.rsf rec=Frec.rsf wfl=Fwfl.rsf > Fdat.rsf qqq=Fqqq.rsf verb=n snap=n free=n ssou=0 dabc=n opot=n back=n kspace=n qmod=0 nbell=5 jdata=1 jsnap=nt rfreq=1000 fcut=500 pcut=0.2 nqz=sf_n(az) nqx=sf_n(ax) nqy=sf_n(ay) oqz=sf_o(az) oqx=sf_o(ax) oqy=sf_o(ay)
sou wavelet (nx,ny,nc,nt)
rec data (nx,ny,nc,nt)
sou geometry (nc,nx,ny)
rec geometry (nc,nx,ny)

 
Parameters
       
 
bool back=n [y/n]
backward extrapolation flag (for rtm)
 
file ccc=
auxiliary input file name
 
bool dabc=n [y/n]
absorbing BC
 
file den=
auxiliary input file name
 
float fcut=500
cutoff frequency for Q-compensation
 
bool free=n [y/n]
free surface flag
 
int jdata=1

 
int jsnap=nt

 
bool kspace=n [y/n]
k-space method (ps) flag
 
int nbell=5
bell size
 
int nqx=sf_n(ax)

 
int nqy=sf_n(ay)

 
int nqz=sf_n(az)

 
bool opot=n [y/n]
output potentials -> 1*scalar, 3*vector potentials
 
float oqx=sf_o(ax)

 
float oqy=sf_o(ay)

 
float oqz=sf_o(az)

 
float pcut=0.2
pcut/2 is tapered portion w.r.t. 1
 
int qmod=0
q modeling switch
 
file qqq=
auxiliary input file name
 
file rec=
auxiliary input file name
 
float rfreq=1000
reference frequency for constant-Q
 
bool snap=n [y/n]
wavefield snapshots flag
 
file sou=
auxiliary input file name
 
int ssou=0
0 -> acceleration source; 1 -> stress source; 2 -> displacement source
 
bool verb=n [y/n]
verbosity flag
 
file wfl=
auxiliary output file name