In this exercise you will apply the tdmt_invc_iso code to a California earthquake.

SAC data files 1998* polezero response files *zp are provided. MODEL_gil7 defines the velocity model used
to compute Green's functions.

Station	Distance(km)	azimuth(deg from north)
BKS	142.		331.5
CMB	171.		33.6
HOPS	286.		330.8
JRSC	99.6		316.5
KCC	201.		71.1
MHC	67.		346.4
ORV	311.		359.4
PKD	122.		137.0
SAO	2.0		56.8
YBH	563.		349.3

The station and event latitude and longitudes are already included in the sac files so that they may be
rotated into radial and transverse components.

The processing steps are as follows:
1. using sac read each component separately to remove instrument response.
	a. read 19980812141000.BKS.BHE
	b. rmean
	c. transfer from polezero s BKS_EVBB.zp
	d. mul 100.    			<converts to cm/s>
	e. write tmpe
	f. repeat for N and Z components
	g. read tmpe tmpn
	h. rot to gcp
	i. read more tmpz
	j. int
	k. bp co 0.02 0.05 p 2
	l. interpolate delta 1.0
	m. write tmp2 tmp1 tmp3		<write radial, transverse and vertical components>
2. sac2helm date=BKS.data		<combine the tmp[1,2,3] files in to a single ascii data file>
3. repeat 1 and 2 for other stations
4. compute Green's functions
	a. usage: FKRPROG < MODEL_gil7
	b. usage: wvint9
	c. usage: run_fkrsortiso
5. filter Green's functions
	a. usage: run_filtsyniso gil7_142d8.disp gil7_142d8
6. run tdmt_invc_iso
	a. usage: tdmt_invc_iso
	b. obtain a solution using only BKS
	c. next add data for CMB and obtain a moment tensor solution
	d. next add data for KCC and obtain a moment tensor solution
	e. finally improve the fit by adjusting the zcor values
	f. add more stations if you desire
	g. examine the sensitivity of solution to various combinations of stations
	h. compute Green's functions for different depths and obtain the best fitting depth
