SUTIVEL -  SU Transversely Isotropic velocity table builder			computes vnmo or vphase as a function of Thomsen's parameters and	theta and optionally interpolate to constant increments in slowness Optional Parameters:							 a=2500.		alpha (vertical p velocity)			 b=1250.		beta (vertical sv velocity)			 e=.20			epsilon (horiz p-wave anisotropy)		 d=.10			delta (strange parameter)			 maxangle=90.0		max angle in degrees				 nangle=9001		number of angles to compute			 verbose=0		set to 1 to see full listing			 np=8001		number of slowness values to output		 option=1		1=output vnmo(p) (result used for TI DMO)				2=output vnmo(theta) in degrees						3=output vnmo(theta) in radians						4=output vphase(p)							5=output vphase(theta) in degrees					6=output vphase(theta) in radians					7=output first derivative vphase(p)					8=output first derivative vphase(theta) in degrees			9=output first derivative vphase(theta) in radians			10=output second derivative vphase(p)					11=output second derivative vphase(theta) in degrees			12=output second derivative vphase(theta) in radians			13=( 1/vnmo(0)^2 -1/vnmo(theta)^2 )/p^2 test vs theta			   (result should be zero for all theta for d=e)			14=return vnmo(p) for weak anisotropy		 normalize=0		=1 means scale vnmo by cosine and scale vphase by 			    1/sqrt(1+2*e*sin(theta)*sin(theta)			 		   (only useful for vphase when d=e for constant				result)								=0 means output vnmo or vphase unnormalized	 Output on standard output is ascii text with:				 line   1: number of values						 line   2: abscissa increment (p or theta increment, always starts at zero) line 3-n: one value per line						 Author: (visitor to CSM form Mobil) John E. Anderson, Spring 1994