# 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.

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

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

# Ray tracing
  MODEL='mar-mod.dat'  # Model
  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" "-"

# 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
  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: