Synopsis
		sfisaac3 < refl.rsf > xrefl.rsf aniso=vti.rsf layer= velocity= xgradient= ygradient= zgradient= xref= yref= zref= xinitial= yinitial= xinitial= number= xs= ys= xr= yr= vstatus= order=3 xmin= ymin= xmax= ymax= niter= debug=n tol=0.000001/v_inp[0] ns=3 ns2=nr3+2 ds=1 s0=0 order=3

Parameters
		file aniso= auxiliary input file name   bool debug=n [y/n] Debug flag   float ds=1 Step increment   floats layer= Layer sequence [nr3+1]   int niter= The number of iterations   int ns=3 Dimension of output reflection points (x,y,z)   int ns2=nr3+2 Dimension of output reflection points (the number of points)   int number= Number of reflectors   int order=3 Interpolation order   float s0=0 Staring position   double tol=0.000001/v_inp[0] Assign a default value for tolerance   floats velocity= Assign velocity km/s [N-1]   int vstatus= Velocity status (0 for constant v, 1 for gradient v, and 2 for VTI)   floats xgradient= Assign x-gradient [N-1]   floats xinitial= [nr3]   float xmax=   float xmin=   float xr= x-Receiver   floats xref= Assign x-reference point [N-1]   float xs= x-Source   floats ygradient= Assign y-gradient [N-1]   floats yinitial= [nr3]   float ymax=   float ymin=   float yr= y-Receiver   floats yref= Assign y-reference point [N-1]   float ys= y-Source   floats zgradient= Assign z-gradient [N-1]   floats zref= Assign z-reference point [N-1]