EXAMPLE: 3-C Earthquake Data 
============================

(1) See example script "run_polfre.cmd". It basically handles a list
of parameters which then are used to call the program.

(2) "plot_baz_dop.gmt" is the GMT script to plot the BAZ as function of
frequency and time. "pal.cpt" is the color palette. Output is "plot.ps"

(3) Execution sequence: run_polfre.cmd; plot_baz_dop.gmt ; gv plot.ps

==========================================================================

Files in directory:
--------------------

SAC data:
2007.031.03.15.55.3250.G.PEL..LHZ.D.SACd, 2007.031.03.15.55.3250.G.PEL..LHN.D.SACd
2007.031.03.15.55.3250.G.PEL..LHE.D.SACd

LHE.asc, LHN.asc, LHZ.asc: Amplitude normalized traces used for plotting
                           in the GMT script.
pal.cpt                  : Periodic color palette to plot BAZ with GMT. 
plot_baz_dop.gmt         : GMT script to plot results.
plot_baz_dop_trace.gmt   : GMT script to plot input data + results.
run_polfre.cmd           : Script to run polfre with parameters set to reproduce
                           the figures of this directory.
azi_dopm.asc             : Program output. Generated after execution of run_polfre.cmd
plot.ps plot.pdf         : Figures generated with plot_baz_dop.gmt
plot1.ps plot1.pdf       : Figures generated with plot_baz_dop_trace.gmt
Figures                  : Fig_PEL_baz_dop.pdf, Fig_PEL_baz_dop.ps
                           Fig_PEL_baz_dop_trace.pdf, Fig_PEL_baz_dop_trace.ps      

==========================================================================

In this example one can observe the R1 and R2 fundamental mode Rayleigh
waves and the microseismic noise. The surface wave dispersion and
BAZ can be verified using the information from the SAC-header. 

SAC> r 2007.031.03.15.55.3250.G.PEL..LHZ.D.SACd
SAC> lh stla stlo evla evlo baz gcarc dist
  
  
  FILE: 2007.031.03.15.55.3250.G.PEL..LHZ.D.SACd - 1
 ----------------------------------------------

      stla = -3.314600e+01
      stlo = -7.067500e+01
      evla = -2.959000e+01
      evlo = -1.779300e+02
       baz = 2.363943e+02
     gcarc = 8.709932e+01
      dist = 9.682092e+03

===========================================================================

Event information from USGS:
----------------------------

Magnitude 6.5 - KERMADEC ISLANDS, NEW ZEALAND

2007 January 31 03:15:55 UTC

Magnitude	6.5
Date-Time	Wednesday, January 31, 2007 at 03:15:55 (UTC)
= Coordinated Universal Time
Wednesday, January 31, 2007 at 3:15:55 PM 
= local time at epicenter
Time of Earthquake in other Time Zones
Location	29.593°S, 177.935°W
Depth	53.7 km (33.4 miles)
Region	KERMADEC ISLANDS, NEW ZEALAND
Distances	55 km (35 miles) SSW of Raoul Island, Kermadec Islands
Location Uncertainty	horizontal +/- 9.2 km (5.7 miles); depth +/- 4.5 km (2.8 miles)
Parameters	Nst=175, Nph=175, Dmin=37.9 km, Rmss=0.96 sec, Gp= 29°,
M-type=moment magnitude (Mw), Version=7
Source	USGS NEIC (WDCS-D)
Event ID	us2007yfak