b'\n \n \n
 
\n  
sfisaac2 (4.0)
index
user/zone/Misaac2.c
\n 2D Bending ray tracing in multi-layered media\n

\n \n \n \n \n \n
 
\n Synopsis
       sfisaac2 < refl.rsf > xrefl.rsf aniso=vti.rsf layer= velocity= xgradient= zgradient= xref= zref= velocity= xinitial= number= xs= xr= vstatus= min=(xx[0]xx[nr2+1])? xx[0]:xx[nr2+1] niter=100 debug=n tol=0.000001/v_inp[0] ns=2 ns2=nr2+2 ds=1 s0=0 order=3
\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
file aniso=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
bool debug=n [y/n]
\tDebug flag
\n
\n \n\n \n \n
\n  
float ds=1
\tStep increment
\n
\n \n\n \n \n
\n  
floats layer=
\tLayer sequence [nr2+1]
\n
\n \n\n \n \n
\n  
float max=(xx[0]>xx[nr2+1])? xx[0]:xx[nr2+1]
\tThe maximum boundary if not entered, set to max(xs,xr)
\n
\n \n\n \n \n
\n  
float min=(xx[0]
\tThe minimum boundary if not entered, set to min(xs,xr)
\n
\n \n\n \n \n
\n  
int niter=100
\tThe number of iterations
\n
\n \n\n \n \n
\n  
int ns=2
\tDimension of output reflection points (x,z)
\n
\n \n\n \n \n
\n  
int ns2=nr2+2
\tDimension of output reflection points (the number of points)
\n
\n \n\n \n \n
\n  
int number=
\tNumber of intersecting points [nr2]
\n
\n \n\n \n \n
\n  
int order=3
\tInterpolation order
\n
\n \n\n \n \n
\n  
float s0=0
\tStaring position
\n
\n \n\n \n \n
\n  
double tol=0.000001/v_inp[0]
\tAssign a default value for tolerance
\n
\n \n\n \n \n
\n  
floats velocity=
\tAssign velocity km/s [N-1]
\n
\n \n\n \n \n
\n  
int vstatus=
\tVelocity status (0 for constant v, 1 for gradient v, and 2 for VTI)
\n
\n \n\n \n \n
\n  
floats xgradient=
\tAssign x-gradient [N-1]
\n
\n \n\n \n \n
\n  
floats xinitial=
\t [nr2]
\n
\n \n\n \n \n
\n  
float xr=
\tReceiver
\n
\n \n\n \n \n
\n  
floats xref=
\tAssign x-reference point [N-1]
\n
\n \n\n \n \n
\n  
float xs=
\tSource
\n
\n \n\n \n \n
\n  
floats zgradient=
\tAssign z-gradient [N-1]
\n
\n \n\n \n \n
\n  
floats zref=
\tAssign z-reference point [N-1]
\n
\n \n
'