# Ray tracing in a smoothed Marmousi model for a given point source
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# This history file calculates the gridded values corresponding to
# the point source specified in file 'mar-src.out'.
# This history file is designed to be run from history file
# 'mar-migr.h', which runs 'srp.for' to generate source coordinates
# 'mar-src.out' and then executes this history file for the source.
# The history file may also generate several PostScript figures.

# Form of the temporary files with intermediate results
  FORM='UNFORMATTED'

# Data required
  chk.pl: "data/mar/" "mar-mod.dat"
  chk.pl: "data/mar/" "mar-crt.dat"
#  chk.pl: "data/mar/" "mar-mpic.dat"   # Just for plotting
  chk.pl: "forms/"    "copy.pl"
  chk.pl: "crt/"      "codep.dat"
  chk.pl: "crt/"      "writall.dat"
  chk.pl: ""          "mar-src.out"

# Model
  MODEL='mar-mod.dat'

# Ray tracing
  SRC='mar-src.out'    # File with the coordinates of the initial point
  REC=' '              # No receivers in this example
  DCRT='mar-crt.dat'   # Numerical parameters for complete ray tracing
  CODE='codep.dat'     # Ray codes (direct P wave)
  INIPAR=2             # Geographic-like spherical ray parameters
  RPAR='mar-crt.dat'   # Ray-parameter domain and control
  WRIT='writall.dat'   # Names and extent of output files with rays
  CRTLOG='mar-log.out' # Output log file of ray tracing
  crt:                 # Running program crt.for with the above data

# Information about the current source
  copy.pl: "mar-src.out" "-"

# Grid for interpolation within ray cells
  N1=188  N2=575  D1=16  D2=16  O1=8  O2=8

# Interpolation within ray cells
  AMP2D3=1.  AMPKI1=25.
  NUM='mar-num.out'                      # Numbers of arrivals
  MTT='mar-mtt.out'  AMP='mar-mamp.out'  # Travel times, amplitudes
  AMPKI='mar-maki.out'                   # Amplitudes for Kirchhoff integral
  MTTSORT='AMP'  MTTORDER=-1.            # Sorting according to amplitudes
#  For plotting
#  MTTSORT=       MTTORDER=              # Sorting according to traveltimes
  mtt:

# Converting multi-valued grid into single-valued grids
  MGRD='mar-mtt.out'  GRD='mar-tt.out'
  mgrd:
  MGRD='mar-mamp.out' GRD='mar-amp.out'
  mgrd:
  MGRD='mar-maki.out' GRD='mar-aki.out'
  mgrd:
  N4=

## Generating several PostScript plots (disabled by default)
#
## Plotting rays
#  LIN='mar-bnd.out'  # Filename for the edges of the model box
#  bndlin:
#  NQ=3 KALL=1 RAYS='mar-rayt.out'
#  crtray:
#  PICDAT='mar-mpic.dat'  PICTURE='mar-ray.ps'  # Filenames to plot rays
#  CALCOPS='0.00 setlinewidth'
#  pictures:
#
## Plotting gridded travel times
#  HSIZE=7.50  VSIZE=23.00
#  VREF=0 VCIRC=1  SHOWPAGE=0
#  NUM='mar-num.out' MGRD='mar-mtt.out' GRD='mar-t00.out' PS='mar-t00.ps'
#  mgrd:
#  grdps:
#  N4=             SHOWPAGE=
#  copy.pl: "mar-t00.ps" "mar-t01r.ps"
#  append.pl:            "mar-t01r.ps" "mar-ray.ps"
#  copy.pl: "mar-t01.ps" "mar-t02r.ps"
#  append.pl:            "mar-t02r.ps" "mar-ray.ps"
#  copy.pl: "mar-t02.ps" "mar-t03r.ps"
#  append.pl:            "mar-t03r.ps" "mar-ray.ps"
#
## Plotting gridded numbers of arrivals
#  VCIRC= VPLUS=1  VREF=0 CREF=0.1666667
#  GRD='mar-num.out'  PS='mar-num.ps'
#  SHOWPAGE=0
#  grdps:
#  copy.pl: "mar-num.ps" "mar-numr.ps"
#  append.pl:            "mar-numr.ps" "mar-ray.ps"
#
## Calculating gridded slowness
#  MPAR=1  ICB=' '  VEL='mar-vel.out'
#  grid:
#  CAL='inv.cal'  GRD1='mar-vel.out'  GRD2='mar-u.out'
#  grdcal:
#
## Plotting gridded slowness
#  VCIRC=.000581818  VPLUS= VREF=.000666667  CREF=0.666667  # 1/1500 blue,
## 1/1755 cyan, 1/2115 green, 1/2661 yellow, 1/3587 red, 1/5500 magenta
#  GRD='mar-u.out'  PS='mar-u.ps'
#  SHOWPAGE=0
#  grdps:
#  copy.pl: "mar-u.ps" "mar-slor.ps"
#  append.pl:          "mar-slor.ps" "mar-ray.ps"

# Important output files
# ~~~~~~~~~~~~~~~~~~~~~~
# 'mar-num.out'... Grid with numbers of arrivals
# 'mar-mtt.out'... Multivalued grid with traveltimes
# 'mar-mamp.out'.. Multivalued grid with amplitudes
# 'mar-maki.out'.. Multivalued grid with amplitudes for Kirchhoff integral
# 'mar-mtt.out'... Singlevalued grids with traveltimes
# 'mar-mamp.out'.. Singlevalued grids with amplitudes
# 'mar-maki.out'.. Singlevalued grids with amplitudes for Kirchhoff integral

# Output PostScript plots (disabled by default)
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# 'mar-slor.ps'... Plot of the slowness together with rays
# 'mar-numr.ps'... Plot of the numbers of arrivals together with rays
# 'mar-t**r.ps'... Plot of the travel times together with rays