view release on metacpan or  search on metacpan

lib/App/SeismicUnixGui/sunix/model/susynvxz.pm  view on Meta::CPAN

 susynvxz >outfile [optional parameters]				



 Required Parameters:							

 <vfile		file containing velocities v[nx][nz]		

 nx=			number of x samples (2nd dimension)		

 nz=			number of z samples (1st dimension)		

 Optional Parameters:							

 nxb=nx		band centered at midpoint			

 nxd=1			skipped number of midponits			

 dx=100		x sampling interval (m)				

 fx=0.0		first x sample					

 dz=100		z sampling interval (m)				

 nt=101		number of time samples				

 dt=0.04		time sampling interval (sec)			

 ft=0.0		first time (sec)				

 nxo=1		 	number of offsets				

 dxo=50		offset sampling interval (m)			

 fxo=0.0		first offset (m)				

 nxm=101		number of midpoints				

 dxm=50		midpoint sampling interval (m)			

 fxm=0.0		first midpoint (m)				

 fpeak=0.2/dt		peak frequency of symmetric Ricker wavelet (Hz)	

 ref="1:1,2;4,2"	reflector(s):  "amplitude:x1,z1;x2,z2;x3,z3;...

 smooth=0		=1 for smooth (piecewise cubic spline) reflectors

 ls=0			=1 for line source; default is point source	

 tmin=10.0*dt		minimum time of interest (sec)			

 ndpfz=5		number of diffractors per Fresnel zone		

 verbose=0		=1 to print some useful information		



 Notes:								

 This algorithm is based on formula (58) in Geo. Pros. 34, 686-703,	

 by N. Bleistein.							



 Offsets are signed - may be positive or negative.			", 

 Traveltime and amplitude are calculated by finite differences which	

 is done only in part of midpoints; in the skiped midpoint, interpolation

 is used to calculate traveltime and amplitude.			", 



 More than one ref (reflector) may be specified.			

 Note that reflectors are encoded as quoted strings, with an optional	

 reflector amplitude: preceding the x,z coordinates of each reflector.	

 Default amplitude is 1.0 if amplitude: part of the string is omitted.	







   CWP:  Zhenyue Liu, 07/20/92

	Many subroutines borrowed from Dave Hale's program: SUSYNLV



 Trace header fields set: trid, counit, ns, dt, delrt,

				tracl. tracr,

				cdp, cdpt, d2, f2, offset, sx, gx



=head2 User's notes (Juan Lorenzo)
untested

=cut


=head2 CHANGES and their DATES

=cut

use Moose;
our $VERSION = '0.0.1';


=head2 Import packages

=cut