b'\n \n \n
 
\n  
sfmpilrrtm_ts (4.0)
index
user/zhiguang/Mmpilrrtm_ts.c
\n One-step lowrank RTM with time-shift imaging condition \n

\n \n \n \n \n \n
 
\n Synopsis
       sfmpilrrtm_ts Fimg2=Fimg2.rsf Fsrc=Fsrc.rsf Fpadvel=Fvel.rsf Ffwf=Ffwf.rsf Fbwf=Fbwf.rsf taper=0 thresh=0.92 wantwf=n verb=n pad1=1 nb= srctrunc=0.4 rectx=2 rectz=2 repeat=2 scalet=1 snap=100 snapshot=0 nds= gpz= spx= spz= rnx= ndr=1 nr0=0 ntau= dtau= tau0= ncut=0 path1= path2=
\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
file Fbwf=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Ffwf=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fimg2=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fpadvel=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
file Fsrc=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
float dtau=
\tinterval of time-shift
\n
\n \n\n \n \n
\n  
int gpz=
\tdepth of geophone
\n
\n \n\n \n \n
\n  
int nb=
\tpadded boundary width
\n
\n \n\n \n \n
\n  
int ncut=0
\tnumber of cutting samples for generating non-negative source wavelet
\n
\n \n\n \n \n
\n  
int ndr=1
\treceiver interval
\n
\n \n\n \n \n
\n  
int nds=
\tsource interval in number of dx
\n
\n \n\n \n \n
\n  
int nr0=0
\treceiver starting point in rnx
\n
\n \n\n \n \n
\n  
int ntau=
\tnumber of time-shift
\n
\n \n\n \n \n
\n  
int pad1=1
\tpadding factor on the first axis
\n
\n \n\n \n \n
\n  
string path1=
\tpath of left matrices \'./mat/left\'
\n
\n \n\n \n \n
\n  
string path2=
\tpath of right matrices \'./mat/left\'
\n
\n \n\n \n \n
\n  
int rectx=2
\tsource smoothing in x-direction
\n
\n \n\n \n \n
\n  
int rectz=2
\tsource smoothing in z-direction
\n
\n \n\n \n \n
\n  
int repeat=2
\trepeat numbers of source smoothing
\n
\n \n\n \n \n
\n  
int rnx=
\tcoverage area of one shot
\n
\n \n\n \n \n
\n  
int scalet=1
\twavefield storage interval
\n
\n \n\n \n \n
\n  
int snap=100
\twavefield output interval when wantwf=y
\n
\n \n\n \n \n
\n  
int snapshot=0
\tprint out the wavefield snapshots of this shot
\n
\n \n\n \n \n
\n  
int spx=
\thorizontal location of source
\n
\n \n\n \n \n
\n  
int spz=
\tdepth of source
\n
\n \n\n \n \n
\n  
float srctrunc=0.4
\tsource truncation
\n
\n \n\n \n \n
\n  
int taper=0
\tif not 0, tapering in the frequency domain
\n
\n \n\n \n \n
\n  
float tau0=
\torigin of time-shift
\n
\n \n\n \n \n
\n  
float thresh=0.92
\ttapering threshold
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]
\tverbosity flag
\n
\n \n\n \n \n
\n  
bool wantwf=n [y/n]
\tif true, output wavefield of a certain (snapshot=) shot
\n
\n \n
'