# History file 'p1-ss2.h' computes the ray-theory seismograms
# in model p1-8-10%.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "forms/"    "echo.pl"
  chk.pl: "forms/"    "copy.pl"
  chk.pl: "forms/"    "append.pl"
  chk.pl: "crt/"      "source.dat"
  chk.pl: "data/p1i/" "p1-mod1c.dat"
  chk.pl: "data/p1i/" "p1-src.dat"
  chk.pl: "data/p1i/" "p1-recb.dat"
  chk.pl: "data/p1i/" "p1i-crt.dat"
  chk.pl: "data/p1i/" "p1i-codb.dat"
  chk.pl: "data/p1i/" "p1i-rpab.dat"
  chk.pl: "data/p1i/" "p1-writ.dat"
  chk.pl: "data/p1i/" "p1-bnd.dat"

# Data for ray tracing
  MODEL='p1-mod1c.dat'
  SRC='p1-src.dat'  REC='p1-recb.dat'
  DCRT='p1i-crt.dat'
  CODE='p1i-codb.dat'
  RPAR='p1i-rpab.dat'
  WRIT='p1-writ.dat'
  CRTLOG='p1i-log.out'

# Kind of ray parameters (for details see 'init.for')
  INIPAR=3

# Reference frequency for the attenuation
  TIMUL=1  FREF=10

# Data describing the source time function
  KSIG=1 # (Gabor signal)
  SIGT=0.  SIGF=10.  SIGW=4.  SIGPH=0.
  SIGA=258.174 # (258.174=12*pi*rho*v**3)

# Data describing the filtration of the source time function
  FMIN=1.  FLOW=2.  FHIGH=19.  FMAX=20. # (cosine filter)

# Time step and time interval for the Fast Fourier Transform
  DT=0.001  NFFT=16384

# Data to control optional plotting of program SS
# SIGPLOT='p1-sig.ps'   MPTS=100

# GSE file with synthetic seismograms
  SS='p1-rss.gse'

# Data to control seismogram plotting (program SP)
  SPTMIN=1.5    SPTMAX=5.4    SPXLEN=17  SPTDIV=3.9  SPTSUB=10
  SPTLEN=25
  NORMSP=1   SPAMP=5.5e6
  SPTEXT3='BLACK...Ray theory'
  SP1='p1-rss1.ps'  SP2='p1-rss2.ps'  SP3='p1-rss3.ps'

# Ray tracing
# ~~~~~~~~~~~
  CRTSCRO='swrg'
  CRTSCRO='swr'
  CRTANI=0
  crt:

# Plotting the ray
# ~~~~~~~~~~~~~~~~~
  CRTOUT='crtout.tmp'
  echo.pl: "'r06.out' * 'r06i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r06.tmp'
  crtray:

  PICDAT='picdat.tmp' PICTURE='p1i-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: "'p1i-r06.tmp'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'POINTS'                                  " ">>picdat.tmp"
  echo.pl: "'p1-recb.dat'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  pictures:
  copy.pl:      "p1-icb.ps" "p1i-v06.ps"
  append.pl:                "p1i-v06.ps" "p1i-ray.ps"


  CRTOUT='crtout.tmp'
  echo.pl: " *      's01.out' 's01i.out' /      " ">  crtout.tmp"
  KALL=-2
  crtpts:  # To determine the numbers of the reflected waves.
           # Only wave number 6 incidents into the receiver.
  MODEL=
  green:
  greenss:
  ss:
  sp: