This folder has the group velocity data used in our paper. (Ma and Masters 2014, A new global Rayleigh and Love wave group velocity dataset for constraining lithosphere properties)

(1) file names
There is a separate file for each frequency and wave type (gv=Rayleigh and gvl=Love).  For example, ph10.abs means Rayleigh wave at 10mHz.

BH data are also added to our LH dataset. Measurements on BH data have file names starting with gvbh or gvlbh.

(2) content of each file
An example of reading these files is

        read(7,887) kk,ky,kd,kh,km,ssk,th,ph,d0
  887   format(i5,i6,i4,2i3,f7.3,3f9.3)
        
        do k=1,kk
          read(7,888) nsta,tbest,terr,
     +      scale,scale_se,pol,ravg,time,azisv,distsv,snr,ttmx,ktyp
  888       format(a4,11f10.3,2x,a2)
        enddo


The meaning of each variable is
kk:  number of records in this cluster
ky,kd,kh,km,ssk:  event origin time (PDE catalog) (year, day, hour, minute and second)
th,ph,d0:  PDE location (colatitude,longitude and depth)

nsta:  station name
tbest,terr:  travel time anomaly (dt) and its error  (in seconds)
scale,scale_se:  amplitude variation and its error (can be ignored)
pol:  centroid frequency (Shapiro and Singh 1999)
ravg:  average cross-correlation coefficient
time:  1D reference arrival time
azisv,distsv:  azimuth and distance (in degrees)
snr:  signal to noise ratio
ttmx:  arrival time of the peak of the envelope
ktyp:  network code

One of the most useful information is the average velocity (km/sec) between each source-station pair, and it can be estimated by using
distsv*111.1949/(tbest+time)


If you have any questions or comments or find any bugs, feel free to email me (Zhitu Ma z1ma@ucsd.edu).