Programs to compute phase velocities (between a source and a station).

Step 1 : Compute spectra using the script run.ida, which calls the 
program idaah

example to compute spectrum for station BKS
type:
run.ida bks BKS
 
--------------
idaah <<!
CALIFORNIA 07/13/91 			/* title of event (needs file DISTCALI)
ida$1r					/* output file will be idabksr
lpz91194.ah				/* input data file in old ah format	
400. 10.				/* period range for computation
4.0 2.6 4.1 3.3 4.0 3.4 3.9 3.30    /* group velocity Umax Umin, for R1,R2,R3,R4
1					/* iwrit=1 to write time series
1					/* component  (Z =1)
1					/* number of channels read
VLPZ					/* names of channels
1               /*   NUMBER OF STATIONS
$2    1 	/* name of station (here BKS) and which trains (here only R1)
! 
----------------

Step 2: compute phase velocities using script run.vit

type:

run.vit bks

--------------
vitphas97 <<!
1					/*calculate phase velocities = 1 */
0					/* antipodal data = 1
$1phasr2a				/* output file (will contain phase velo)
/home/u1/barbara/mineig/b1066ms		/* model file name */
M1066B					/* model title */
ida$1r					/* spectra filename
3					/* spheroidal = 3 toroidal =2 */
1					/* iopt =1 for surface waves */
200 42.182 -125.641			/* noe (not used) lat and long of epicenter */
10. 300. 60. 500. 60. 500. 60. 500.	/* period ranges for R1, R2, R3, R4 */
0.01 0.01 0.01 0.01			/* left over from old version??? */
vitmend		/* reference phase velocity file */
315 78 -178 15.6 10.4 0.21 /* cmt: strike,dip,slip,depth,tau,moment in 10**21N_m
0				/* mode number (0= fundamental mode)
/home/u1/barbara/mineig/out66s		/* eigenfunction file */
256.0				/* reference period
!