# Calculation of geometrical spreading from gridded slowness vectors
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required:
  chk.pl: "data/vgr/" "vgr-velq.dat"
  chk.pl: "data/vgr/" "vgr-srca.dat"
  chk.pl: "data/vgr/" "vgr-ttta.dat"
  chk.pl: "data/vgr/" "vgr-srcb.dat"
  chk.pl: "data/vgr/" "vgr-tttb.dat"
  chk.pl: "data/vgr/" "q2.cal"
  chk.pl: "data/vgr/" "vgr-mod.dat"
  chk.pl: "data/vgr/" "vgr-srcq.dat"
  chk.pl: "data/vgr/" "vgr-crt.dat"
  chk.pl: "crt/"      "codep.dat"
  chk.pl: "crt/"      "writall.dat"

# Grid dimensions
  N1=202    N2=202
  D1=0.005  D2=0.005
  O1=-.005  O2=0.000
# Data to calculate Q2
  VSRC=1.   DSRC=0.014142136  QMIN=-999999.  QMAX=999999.
# Data to plot Q2
  VMIN=-5.  VMAX=5.  VCIRC=1.  VREF=0.  CREF=0.166667

# Travel-time and slowness-vector calculation
  NFSMAX=-1           # Second-order method
  NET='vgr-ttta.dat'  # List of input and output filenames, source A
  net:                # Program to run
  NET='vgr-tttb.dat'  # List of input and output filenames, source B
  net:                # Program to run

# Geometrical spreading
  CAL='q2.cal'
  GRD1='tta.out'  GRD2='p1a.out'  GRD3='p2a.out'
  GRD4='ttb.out'  GRD5='p1b.out'  GRD6='p2b.out'  GRD7='vgr-q.out'
  grdcal:
  GRD='vgr-q.out'  PS='vgr-q.ps'
  grdps:


# Ray tracing
  MODEL='vgr-mod.dat'  # Velocity model
  SRC='vgr-srcq.dat'   # Source coordinates
  DCRT='vgr-crt.dat'
  RPAR='vgr-crt.dat'
  CODE='codep.dat'     # Direct P wave
  WRIT='writall.dat'   # Writing all rays into output files
  CRTLOG='vgr-log.out'
  INIPAR=-3            # Ray parameters = azimuth. equidist. projection
  crt:

# New grid
  N2=  N3=202
  D2=  D3=0.005
  O2=  O3=0.000

# Interpolation of the geometrical spreading in ray cells
  NUM='num.tmp'  MTT=' '  MPQ13='qmtt.tmp'
  mtt:

# Plotting the geometrical spreading
  MGRD='qmtt.tmp'  GRD='vgr-qmtt.out'  PS='vgr-qmtt.ps'
  mgrd:
  grdps: