sftkirinv (2.0-git)
2-D least-squares Kirchhoff pre-stack time migration with different regul..

        sftkirinv < in.rsf vel=vel.rsf > out.rsf offset=offset.rsf fdip=fdip.rsf ompchunk=1 ompnth=0 verb=n half=y amp=y ncdp=ncmp dcdp=dcmp cdp0=cmp0 reg=0 antialias=1.0 apt=ncmp rho=1.-1./nt niter=5 nw=3 err=
regularization (or preconditioning) operator:
reg=0: no regularization;
reg=1: regularization => first derivative along offset axis;
reg=2: precondition => causual integration along offset axis;
reg=3: precondition => triangle smoother along offset axis,
reg=4: precondition => local slope constraints along t-x plane and smoothing along offset axis

bool amp=y [y/n]
if y, use amplitue factor
float antialias=1.0
float apt=ncmp
migration aperture
float cdp0=cmp0

float dcdp=dcmp

string err=
output file for error
string fdip=
auxiliary input file name
bool half=y [y/n]
if y, the third axis is half-offset instead of full offset
int ncdp=ncmp

int niter=5
number of iterations
int nw=3

string offset=
auxiliary input file name
int ompchunk=1
OpenMP data chunk size
int ompnth=0
OpenMP available threads
int reg=0
regularization type
float rho=1.-1./nt
Leaky integration constant
file vel=
auxiliary input file name
bool verb=n [y/n]
verbosity flag