sfpremig (2.0-git)
index
user/jsun/Mpremig.c
Pseudo-spectral pre-stack source-receiver source independent diffraction imaging

 
Synopsis
        sfpremig < Fi.rsf > Fo.rsf vel=Fv.rsf dat_v=Fd_v.rsf dat_2=Fi1.rsf snaps=snaps.rsf spx= spz= f0= t0= A= snap=0 cmplx=y pad1=1 abc=n roll=n born=n diff=n justrec=n nbt= nbb=nbt nbl=nbt nbr=nbt ct= cb=ct cl=ct cr=ct verb=n ps=n adj=n vref=1500 nt= dt= gpl=-1 gpl_v=-1 src=0 n_srcs=1 shtbgn= shtend= shtint= which=0 gpx=-1 gpz=-1 gpx_v=-1 gpz_v=-1 offset=0 split=1

 
Parameters
       
 
floats A=
wavelet amplitude [n_srcs]
 
bool abc=n [y/n]
absorbing flag
 
bool adj=n [y/n]
use pseudo-spectral
 
bool born=n [y/n]
born modeling flag
 
float cb=ct

 
float cl=ct

 
bool cmplx=y [y/n]
use complex fft
 
float cr=ct

 
float ct=

 
file dat_2=
auxiliary input file name
 
string dat_v=
auxiliary input file name
 
bool diff=n [y/n]
diffraction imaging flag
 
float dt=

 
floats f0=
wavelet peak freq [n_srcs]
 
int gpl=-1
geophone length
 
int gpl_v=-1
geophone height
 
int gpx=-1
geophone position x
 
int gpx_v=-1
geophone position x
 
int gpz=-1
geophone position z
 
int gpz_v=-1
geophone position z
 
bool justrec=n [y/n]
just need full waveform record (no born or rtdm)
 
int n_srcs=1
source type
 
int nbb=nbt

 
int nbl=nbt

 
int nbr=nbt

 
int nbt=

 
int nt=

 
int offset=0
nearest offset
 
int pad1=1
padding factor on the first axis
 
bool ps=n [y/n]
use pseudo-spectral
 
bool roll=n [y/n]
if n, receiver is independent of source location and gpl=nx
 
int shtbgn=

 
int shtend=

 
int shtint=

 
int snap=0
interval for snapshots
 
file snaps=
auxiliary output file name
 
int split=1
receiver split
 
ints spx=
shot position x [n_srcs]
 
ints spz=
shot position z [n_srcs]
 
int src=0
source type
 
floats t0=
wavelet time lag [n_srcs]
 
file vel=
auxiliary input file name
 
bool verb=n [y/n]
verbosity
 
float vref=1500
reference velocity (default using water)
 
int which=0