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

 
Synopsis
        sfeweks3d < Fwav.rsf ccc=Fccc.rsf den=Fden.rsf sou=Fsou.rsf rec=Frec.rsf wfl=Fwfl.rsf > Fdat.rsf verb=n snap=n free=n ssou=0 dabc=n opot=n back=n kspace=n nbell=5 jdata=1 jsnap=nt 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
 
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)

 
file rec=
auxiliary input file name
 
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