sfeikods (4.0 Meikonal.c 7107 2011-04-10 02:04:14Z ivlad)
index
user/llisiw/Meikods.c
Fast marching with source perturbation.

 
Synopsis
        sfeikods < vel.rsf > time.rsf shotfile=shots.rsf tdl1=tdl1.rsf tds1=tds1.rsf tdl2=tdl2.rsf tds2=tds2.rsf vel=y order=2 sweep=n br1=d1 br2=d2 br3=d3 plane1=n plane2=n plane3=n b1=plane[2]? n1: (int) (br1/d1+0.5) b2=plane[1]? n2: (int) (br2/d2+0.5) b3=plane[0]? n3: (int) (br3/d3+0.5) zshot=0. yshot=o2 + 0.5*(n2-1)*d2 xshot=o3 + 0.5*(n3-1)*d3 l=1

 
Parameters
       
 
int b1=plane[2]? n1: (int) (br1/d1+0.5)

 
int b2=plane[1]? n2: (int) (br2/d2+0.5)

 
int b3=plane[0]? n3: (int) (br3/d3+0.5)
Constant-velocity box around the source (in samples)
 
float br1=d1

 
float br2=d2

 
float br3=d3
Constant-velocity box around the source (in physical dimensions)
 
int l=1
source perturbation direction
 
int order=2 [1,2]
Accuracy order
 
bool plane1=n [y/n]

 
bool plane2=n [y/n]

 
bool plane3=n [y/n]
plane-wave source
 
string shotfile=
File with shot locations (n2=number of shots, n1=3) (auxiliary input file name)
 
bool sweep=n [y/n]
if y, use fast sweeping instead of fast marching
 
string tdl1=
auxiliary output file name
 
string tdl2=
auxiliary output file name
 
string tds1=
auxiliary output file name
 
string tds2=
auxiliary output file name
 
bool vel=y [y/n]
if y, the input is velocity; n, slowness squared
 
float xshot=o3 + 0.5*(n3-1)*d3

 
float yshot=o2 + 0.5*(n2-1)*d2

 
float zshot=0.
Shot location (used if no shotfile)

 
Used In
       

 
GALLERY
        gradient/fakirmig
marmousi/recurkir
overthrust/kirch

 
GEO384S
        hw6/kirchhoff

 
SEP
        fat/migration

 
TCCS
        eikods/check
eikods/marm
eikods/migration
t2dweak/alaska
t2dweak/field
time2depth2/beinew