b'\n \n \n
 
\n  
sfrays3 (4.0)
index
system/seismic/Mrays3.c
\n Ray tracing by a Runge-Kutta integrator in 3-D.\n

\n \n \n \n \n \n
 
\n Synopsis
       sfrays3 < vel.rsf > rays.rsf shotfile=shots.rsf anglefile=angles.rsf vel=y order=4 nt= dt= sym=y escvar=n zshot=o[0] yshot=o[1] + 0.5*(n[1]-1)*d[1] xshot=o[2] + 0.5*(n[2]-1)*d[2] na= nb= a0=0. amax=360. b0=0. bmax=180. > rays.rsf
Rays can be plotted with sfplotrays.
\n\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
float a0=0.
\tFirst azimuth angle in degrees (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
float amax=360.
\tMaximum azimuth angle in degrees (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
string anglefile=
\tfile with initial angles (auxiliary input file name)
\n
\n \n\n \n \n
\n  
float b0=0.
\tFirst inclination angle in degrees (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
float bmax=180.
\tMaximum inclination angle in degrees (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
float dt=
\tSampling in time
\n
\n \n\n \n \n
\n  
bool escvar=n [y/n]
\tIf y - output escape values, n - trajectories
\n
\n \n\n \n \n
\n  
int na=
\tNumber of azimuths (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
int nb=
\tNumber of inclinations (if anglefile is not specified)
\n
\n \n\n \n \n
\n  
int nt=
\tNumber of time steps
\n
\n \n\n \n \n
\n  
int order=4
\tInterpolation order
\n
\n \n\n \n \n
\n  
string shotfile=
\tfile with shot locations (auxiliary input file name)
\n
\n \n\n \n \n
\n  
bool sym=y [y/n]
\tif y, use symplectic integrator
\n
\n \n\n \n \n
\n  
bool vel=y [y/n]
\tIf y, input is velocity; if n, slowness
\n
\n \n\n \n \n
\n  
float xshot=o[2] + 0.5*(n[2]-1)*d[2]
\tshot location crossline (if shotfile is not specified)
\n
\n \n\n \n \n
\n  
float yshot=o[1] + 0.5*(n[1]-1)*d[1]
\tshot location inline (if shotfile is not specified)
\n
\n \n\n \n \n
\n  
float zshot=o[0]
\tshot location in depth (if shotfile is not specified)
\n
\n \n
'