%procedure for generating sensitivity kernels for each period
% need phase velocity
% all routines written by Yingie Yang, 2004

1. run createsac.f to create a single frequency wave
 (compile with $SACLIB/sacio.a)
	output file name 50.sac (50 means period)

2.  look at file in sac
     decide on cut window near center of sine wave
    (depends on way data was windowed,
     usually ~200s for short periods, 
     may change for long periods)	

3. run sac macro "bpcut2.m"
      cut1, cut2 = narrow cut
      input filename 50.sac
      output filename 50.cut800t1000.sac	


4.  run the follwoing in SAc

fft 50.cut800t1000.sac
wsp
r   50.cut800t1000.sac.am
   (look at spectrum, see where peak is, 
   and where first lobe starts and ends)

5. run getdatafromsac.f
     (compile with $SACLIB/sacio.a)
   input filename 50.cut800t1000.am.sac
   output filename 50.cut800t1000.am.sac.dat
  	

  look at 50.cut800t1000.am.sac.dat file.
  find highest value at spectrum peak
  find smallest value before and after this
  remove lines above and below this range	
  In the other words, this part is to retain the main lobe 

6. run sensitivity.f  like "sensitivity < sensinp"
     (normal compiling)
  input phase velocity
  input spectral file: 50.cut800t1000.am.sac.dat
 input smoothed file name  sens50s65km.dat 
 input scalelen ( about 65 km )

done!