sfkirmod_newton (2.0-git)
index
user/zone/Mkirmod_newton.c
Kirchhoff 2-D/2.5-D modeling in layered media with bending ray tracing.

 
Synopsis
        sfkirmod_newton < modl.rsf > data.rsf curv=curv.rsf refl=refl.rsf picks=picks.rsf slopes=slopes.rsf aniso=vti.rsf velocity= xgradient= zgradient= xref= zref= lin=n adj=n absoff=n nt= dt=0.004 t0=0. ns=nx s0=x0 ds=dx nh=nx h0=0. dh=dx verb=n r0=1. r0=1. debug=n fwdxini=n vstatus= niter=500 tol=0.00001 order=3 cmp=n freq=0.2/dt rgrad= dip=

 
Parameters
       
 
bool absoff=n [y/n]
y - h0 is not in shot coordinate system
 
bool adj=n [y/n]
adjoint flag
 
file aniso=
auxiliary input file name
 
bool cmp=n [y/n]
compute CMP instead of shot gathers
 
file curv=
auxiliary input file name
 
bool debug=n [y/n]
debug flag
 
float dh=dx
offset increment
 
string dip=
reflector dip file
 
float ds=dx
shot/midpoint increment
 
float dt=0.004
time sampling
 
float freq=0.2/dt
peak frequency for Ricker wavelet
 
bool fwdxini=n [y/n]
use the result of previous iteration to be the xinitial of the next one
 
float h0=0.
first offset
 
bool lin=n [y/n]
if linear operator
 
int nh=nx
number of offsets
 
int niter=500
The number of iterations
 
int ns=nx
number of shots (midpoints if cmp=y)
 
int nt=
time samples
 
int order=3
Interpolation order
 
string picks=
auxiliary output file name
 
float r0=1.
normal reflectivity (if constant)
 
string refl=
auxiliary input file name
 
string rgrad=
AVO gradient file (B/A)
 
float s0=x0
first shot (midpoint if cmp=y)
 
string slopes=
auxiliary output file name
 
float t0=0.
time origin
 
double tol=0.00001
Assign a default value for tolerance
 
floats velocity=
Assign velocity km/s [nc]
 
bool verb=n [y/n]
verbosity flag
 
int vstatus=
Velocity status (0 for constant v,1 for gradient v, and 2 for vti)
 
floats xgradient=
[nc]
 
floats xref=
Assign x-reference point [nc]
 
floats zgradient=
[nc]
 
floats zref=
Assign z-reference point [nc]

 
Used In
       

 
RSF
        su/rsfmodel