# History file 'p1-mtt.h' computes the Green functions from the source
# and from the receiver, slowness vectors, travel times and other
# quantities at the gridpoints.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
 #chk.pl: "forms/"    "echo.pl"
 #chk.pl: "forms/"    "copy.pl"
 #chk.pl: "forms/"    "append.pl"
 #chk.pl: "data/p1i/" "p1-mod2.dat"
 #chk.pl: "data/p1i/" "p1-recb.dat"
 #chk.pl: "data/p1i/" "p1-src.dat"
 #chk.pl: "data/p1i/" "p1-codep.dat"
 #chk.pl: "data/p1i/" "p1-crt-g.dat"
 #chk.pl: "data/p1i/" "p1-crt-u.dat"
 #chk.pl: "data/p1i/" "p1-bnd.dat"
 #chk.pl: "crt/"      "writall.dat"
 #chk.pl: " "         "p1-vep.ps"     # output of p1-grid.h


# Computation of the Green function from the receiver
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Unperturbed model
  MODEL='p1-mod2.dat'

# Source of the Green function
  SRC='p1-recb.dat'

# Data filenames for CRT (used also in the next computation)
  WRIT='writall.dat'  # all rays are required by 'mtt.for'
  CRTLOG='p1-crt.out'

# Data filenames for CRT
  CODE='p1-codep.dat'
  DCRT='p1-crt-g.dat'
  RPAR='p1-crt-g.dat'

# Kind of ray parameters
  INIPAR=3

# Ray tracing
  CRTSCRO='swr'
  crt:

# Grid dimensions
  N1=2500     N2=1   N3=650
  D1=0.005    D2=0.  D3=0.005
  O1=21.0025  O2=0.  O3=2.0025

# Interpolation within ray tubes
  NUM='p1-num-g.out'
  MTT='p1-mtt-g.out'
  MP1='p1-mp1-g.out'  MP2='p1-mp2-g.out'  MP3='p1-mp3-g.out'
  MTT22='p1-mt22g.out'
  AMPR11='p1-ar11g.out'  AMPI11='p1-ai11g.out'
  AMPR21='p1-ar21g.out'  AMPI21='p1-ai21g.out'
  AMPR31='p1-ar31g.out'  AMPI31='p1-ai31g.out'
  AMPR12='p1-ar12g.out'  AMPI12='p1-ai12g.out'
  AMPR22='p1-ar22g.out'  AMPI22='p1-ai22g.out'
  AMPR32='p1-ar32g.out'  AMPI32='p1-ai32g.out'
  AMPR13='p1-ar13g.out'  AMPI13='p1-ai13g.out'
  AMPR23='p1-ar23g.out'  AMPI23='p1-ai23g.out'
  AMPR33='p1-ar33g.out'  AMPI33='p1-ai33g.out'
  mtt:

# Plotting several pictures with the rays
#
# Plotting the rays
  CRTOUT='crtout.tmp'
  echo.pl: "'r01.out' * 'r01i.out' /  " ">  crtout.tmp"
  RAYS='p1-r01.tmp'
  KALL=1
  crtray:

  PICDAT='picdat.tmp' PICTURE='p1-ray.ps'
  echo.pl: "'PROJECTION'                              " "> picdat.tmp"
  echo.pl: "16.35  27.2    0. -0.5  0. 0. -0.5 0. /   " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                              " ">>picdat.tmp"
  echo.pl: "'MSZSF' 1.01 / Marker size scale factor.  " ">>picdat.tmp"
  echo.pl: "'ITC'   0    / No texts drawn.            " ">>picdat.tmp"
  echo.pl: "'LWIDTH'   0.01    /                      " ">>picdat.tmp"
  echo.pl: "/                                         " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1-bnd.dat'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1-r01.tmp'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'POINTS'                                  " ">>picdat.tmp"
  echo.pl: "'p1-recb.dat'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  pictures:

# Plotting P-wave velocity in the unperturbed model with the rays
  copy.pl:    "p1-vep.ps" "p1-v-g.ps"
  append.pl:              "p1-v-g.ps" "p1-ray.ps"

# Plotting indices of the blocks with the rays
  copy.pl:    "p1-icb.ps" "p1-cov-g.ps"
  append.pl:              "p1-cov-g.ps" "p1-ray.ps"

# Computing the numbers of arrivals and travel times for plotting
# N1=4730   N2=1   N3=600
# D1=0.01   D2=0.  D3=0.01
  O1=0.005  O2=0.  O3=0.005
  N1=2365   N2=1   N3=300
  D1=0.02   D2=0.  D3=0.02

  NUM='p1-numgg.out'
  MTT='p1-mttg.out'
  MHIST='p1-histg.out'
  MP1=' '  MP2=' '  MP3=' '
  MTT22=' '
  AMPR11=' '  AMPI11=' '
  AMPR21=' '  AMPI21=' '
  AMPR31=' '  AMPI31=' '
  AMPR12=' '  AMPI12=' '
  AMPR22=' '  AMPI22=' '
  AMPR32=' '  AMPI32=' '
  AMPR13=' '  AMPI13=' '
  AMPR23=' '  AMPI23=' '
  AMPR33=' '  AMPI33=' '
  mtt:
  MHIST=

# Plotting the numbers of arrivals with the rays
  SHOWPAGE=0
  GRD='p1-numgg.out'
  PS='p1-num-g.ps'
  VPLUS=1
  VREF=0
  CREF=0.1666
  VCIRC=
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  grdps:
  append.pl:               "p1-num-g.ps" "p1-ray.ps"

# Plotting travel times at the gridpoints,
# sorted according to their value
  VPLUS= VCIRC=1  VREF=0  CREF=0.166667
  MGRD='p1-mttg.out'  GRD='p1-ttg00.out' PS='p1-ttg01.ps'
  mgrd:
  grdps:
  N4=

# Plotting grid travel times sorted according to the histories
  HIST='p1-histg.out' GRD='p1-thg00.out' PS='p1-thg01.ps'
  mttgrd>mttgrdg.out:
  grdps:
  N4=

# Computation of the Green function from the source
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Source of the Green function
  SRC='p1-src.dat'

# Data filenames for CRT
  CODE='p1-codep.dat'
  DCRT='p1-crt-u.dat'
  RPAR='p1-crt-u.dat'

# Ray tracing
  crt:

# Grid dimensions
  N1=2500     N2=1   N3=650
  D1=0.005    D2=0.  D3=0.005
  O1=21.0025  O2=0.  O3=2.0025

# Interpolation within ray tubes
  NUM='p1-num-u.out'
  MTT='p1-mtt-u.out'
  MP1='p1-mp1-u.out'  MP2='p1-mp2-u.out'  MP3='p1-mp3-u.out'
  MTT22='p1-mt22u.out'
  AMPR11='p1-ar11u.out'  AMPI11='p1-ai11u.out'
  AMPR21='p1-ar21u.out'  AMPI21='p1-ai21u.out'
  AMPR31='p1-ar31u.out'  AMPI31='p1-ai31u.out'
  AMPR12='p1-ar12u.out'  AMPI12='p1-ai12u.out'
  AMPR22='p1-ar22u.out'  AMPI22='p1-ai22u.out'
  AMPR32='p1-ar32u.out'  AMPI32='p1-ai32u.out'
  AMPR13='p1-ar13u.out'  AMPI13='p1-ai13u.out'
  AMPR23='p1-ar23u.out'  AMPI23='p1-ai23u.out'
  AMPR33='p1-ar33u.out'  AMPI33='p1-ai33u.out'
  MP4='p1-p1s.out'  MP5='p1-p2s.out'  MP6='p1-p3s.out'
 #mtt:

# Plotting several pictures with the rays
# Plotting the rays
  CRTOUT='crtout.tmp'
  echo.pl: "'r01.out' * 'r01i.out' /  " ">  crtout.tmp"
  RAYS='p1-r01.tmp'
  KALL=1
  crtray:

  PICDAT='picdat.tmp' PICTURE='p1-ray.ps'
  echo.pl: "'PROJECTION'                              " "> picdat.tmp"
  echo.pl: "16.35  27.2    0. -0.5  0. 0. -0.5 0. /   " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                              " ">>picdat.tmp"
  echo.pl: "'MSZSF' 1.01 / Marker size scale factor.  " ">>picdat.tmp"
  echo.pl: "'ITC'   0    / No texts drawn.            " ">>picdat.tmp"
  echo.pl: "'LWIDTH'   0.01    /                      " ">>picdat.tmp"
  echo.pl: "/                                         " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1-bnd.dat'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1-r01.tmp'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'POINTS'                                  " ">>picdat.tmp"
  echo.pl: "'p1-recb.dat'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  pictures:

# Plotting P-wave velocity in the unperturbed model with the rays
  copy.pl:    "p1-vep.ps" "p1-v-u.ps"
  append.pl:              "p1-v-u.ps" "p1-ray.ps"

# Plotting indices of the blocks with the rays
  copy.pl:    "p1-icb.ps" "p1-cov-u.ps"
  append.pl:              "p1-cov-u.ps" "p1-ray.ps"

# Computing the numbers of arrivals and travel times for plotting
# N1=4730   N2=1   N3=600
# D1=0.01   D2=0.  D3=0.01
  O1=0.005  O2=0.  O3=0.005
  N1=2365   N2=1   N3=300
  D1=0.02   D2=0.  D3=0.02

  NUM='p1-numuu.out'
  MTT='p1-mttu.out'
  MHIST='p1-histu.out'
  MP1=' '  MP2=' '  MP3=' '
  MTT22=' '
  AMPR11=' '  AMPI11=' '
  AMPR21=' '  AMPI21=' '
  AMPR31=' '  AMPI31=' '
  AMPR12=' '  AMPI12=' '
  AMPR22=' '  AMPI22=' '
  AMPR32=' '  AMPI32=' '
  AMPR13=' '  AMPI13=' '
  AMPR23=' '  AMPI23=' '
  AMPR33=' '  AMPI33=' '
  MP4=' '  MP5=' '  MP6=' '
  mtt:
  MHIST=

# Plotting the numbers of arrivals with the rays
  SHOWPAGE=0
  GRD='p1-numuu.out'
  PS='p1-num-u.ps'
  VPLUS=1
  VREF=0
  CREF=0.1666
  VCIRC=
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  grdps:
  append.pl:               "p1-num-u.ps" "p1-ray.ps"

# Plotting travel times at the gridpoints,
# sorted according to their value
  VPLUS= VCIRC=1  VREF=0  CREF=0.166667
  MGRD='p1-mttu.out'  GRD='p1-ttu00.out' PS='p1-ttu01.ps'
  mgrd:
  grdps:
  N4=

# Plotting grid travel times sorted according to the histories
  HIST='p1-histu.out' GRD='p1-thu00.out' PS='p1-thu01.ps'
  mttgrd>mttgrdh.out:
  grdps:
  N4=

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# p1*g.out... interpolated Green function from the receiver
# p1*u.out... interpolated Green function from the source
# PostScript pictures
# p1-v-g.ps.. velocity in the unperturbed model with rays from receiver
# p1-cov-g.ps indices of blocks with rays from receiver
# p1-num-g.ps numbers of arrivals with rays from receiver
# p1-ttg??.ps travel times with rays from receiver, sorted according the value
# p1-thg??.ps travel times with rays from receiver, sorted according the history
# p1-v-u.ps.. velocity in the unperturbed model with rays from source
# p1-cov-u.ps indices of blocks with rays from source
# p1-num-u.ps numbers of arrivals with rays from source
# p1-ttu??.ps travel times with rays from source, sorted according the value
# p1-thu??.ps travel times with rays from source, sorted according the history