# History file 'sal-ray1.h' computes two point rays in the smoothed
# SEG/EAGE Salt Model, plots computed rays and synthetic seismograms.
# Source is located in the centre of 21x21 square of receivers
# above the salt crest.
# The time of computation is approximately 1 hour
# on a PC equipped with Athlon XP 1700+ and with Linux.
# ======================================================================
#
# Input files required
  chk.pl: "data/sal/" "sal-mod2.dat"
  chk.pl: "data/sal/" "sal-src1.dat"
  chk.pl: "data/sal/" "sal-rec1.dat"
  chk.pl: "data/sal/" "sal-rec1.plt"
  chk.pl: "data/sal/" "sal-dcr1.dat"
  chk.pl: "data/sal/" "sal-cod1.dat"
  chk.pl: "data/sal/" "sal-rpa1.dat"
  chk.pl: "data/sal/" "sal-wri1.dat"
  chk.pl: "crt/"      "sourcex.dat"
  chk.pl: "forms/"    "append.pl"
  chk.pl: "forms/"    "del.pl"
  chk.pl: "forms/"    "echo.pl"

# Input data for complete ray tracing (see crtin.for)
  MODEL='sal-mod2.dat' SRC='sal-src1.dat'    REC='sal-rec1.dat'
  DCRT='sal-dcr1.dat'  CODE='sal-cod1.dat'   RPAR='sal-rpa1.dat'
  WRIT='sal-wri1.dat'  CRTLOG='sal-log.out'  INIPAR=3

# File containing the complex-valued seismic moment
  SOURCE=sourcex.dat # Explosive isotropic source

# Data describing the source time function
  KSIG=1 # Gabor signal
  SIGT=0.  SIGF=30.  SIGW=4.  SIGPH=0.  SIGA=1.

# Data describing the filtration of the source time function
  FMIN=1.  FLOW=5.   FHIGH=56.   FMAX=60.  # (cosine filter)

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

# Data to control seismogram plotting (for details see 'sp.for')
  KODESP=2     SPCHRH=0.40
  SPTMIN=1.0   SPTMAX=4.0     SPTLEN=25.0  SPTDIV=6    SPTSUB=5
  SPXMIN=4160. SPXMAX=11200.  SPXLEN=15.3  SPXDIV=-22  SPXSUB=1
  SPYMIN=7640. SPYMAX=7640.
  NORMSP=0     SPAMP=1.

# Ray tracing of refracted and reflected waves
  crt:

# Ray parameters plotting
  ISANG=1   ITHOM=0    ISUC=1    IRAUX=0 IRBAS=0
  HRBAS=0.2 HRTWO=0.3  HRAUX=0.1 IWAVES=1
  CRTOUT='crtout.tmp'
  echo.pl: "' '   ' '   's011i.out' 't01.out'" ">crtout.tmp"
  RPPLOT='rpcode1.ps'  IWAVES=1
  rpplot:
  echo.pl: "' '   ' '   's012i.out' 't02.out'" ">crtout.tmp"
  RPPLOT='rpcode2.ps'  IWAVES=2
  rpplot:
  echo.pl: "' '   ' '   's013i.out' 't03.out'" ">crtout.tmp"
  RPPLOT='rpcode3.ps'  IWAVES=3
  rpplot:
  echo.pl: "' '   ' '   's014i.out' 't04.out'" ">crtout.tmp"
  RPPLOT='rpcode4.ps'  IWAVES=4
  rpplot:
  echo.pl: "' '   ' '   's015i.out' 't05.out'" ">crtout.tmp"
  RPPLOT='rpcode5.ps'  IWAVES=5
  rpplot:
  echo.pl: "' '   ' '   's016i.out' 't06.out'" ">crtout.tmp"
  RPPLOT='rpcode6.ps'  IWAVES=6
  rpplot:
  echo.pl: "' '   ' '   's017i.out' 't07.out'" ">crtout.tmp"
  RPPLOT='rpcode7.ps'  IWAVES=7
  rpplot:
  echo.pl: "' '   ' '   's018i.out' 't08.out'" ">crtout.tmp"
  RPPLOT='rpcode8.ps'  IWAVES=8
  rpplot:

# Deleting old files
  del.pl:'s01.out'
  del.pl:'s01i.out'
  del.pl:'sal-ray1.mx'

# Computed quantities at the given surface and corresponding
# quantities at the initial points of rays for all calculated
# elementary waves are concatenated.
  append.pl: 's01.out'  's011.out'
  append.pl: 's01.out'  's012.out'
  append.pl: 's01.out'  's013.out'
  append.pl: 's01.out'  's014.out'
  append.pl: 's01.out'  's015.out'
  append.pl: 's01.out'  's016.out'
  append.pl: 's01.out'  's017.out'
  append.pl: 's01.out'  's018.out'
  append.pl: 's01i.out' 's011i.out'
  append.pl: 's01i.out' 's012i.out'
  append.pl: 's01i.out' 's013i.out'
  append.pl: 's01i.out' 's014i.out'
  append.pl: 's01i.out' 's015i.out'
  append.pl: 's01i.out' 's016i.out'
  append.pl: 's01i.out' 's017i.out'
  append.pl: 's01i.out' 's018i.out'

# Frequency-domain seismograms in selected receivers
  CRTOUT=' '
  green:      # Green function
  greenss:    # response function
  ss:         # synthetic seismograms
  REC='sal-rec1.plt'
  SP1='sal-sf1.ps' SP2='sal-sf2.ps' SP3='sal-sf3.ps'
  sp:         # seismogram plotting

# Frequency-domain individual seismograms of elementary waves
# in selected receivers
  REC='sal-rec1.dat'
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r01.out' 's011.out ' 's011i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s1.gse'
  ss:         # synthetic seismograms

# 2nd elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r02.out' 's012.out ' 's012i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s2.gse'
  ss:         # synthetic seismograms

# 3rd elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r03.out' 's013.out ' 's013i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s3.gse'
  ss:         # synthetic seismograms

# 4th elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r04.out' 's014.out ' 's014i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s4.gse'
  ss:         # synthetic seismograms

# 5th elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r05.out' 's015.out ' 's015i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s5.gse'
  ss:         # synthetic seismograms

# 6th elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r06.out' 's016.out ' 's016i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s6.gse'
  ss:         # synthetic seismograms

# 7th elementary wave
  CRTOUT='crtout.tmp' SS=' '
  echo.pl:"'r07.out' 's017.out ' 's017i.out'" ">crtout.tmp"
  green:      # Green function
  greenss:    # response function
  SS='sal-s7.gse'
  ss:         # synthetic seismograms

# There are no rays at the selected receivers for
# 8th elementary wave

# All elementary waves seismogram plotting
  REC='sal-rec1.plt'
  SP1='sal-sf11.ps'  SP2='sal-sf12.ps'  SP3='sal-sf13.ps'
  SS='sal-s1.gse'  KOLOR=2
  SS1='sal-s2.gse' KOLOR1=3
  SS2='sal-s3.gse' KOLOR2=3
  SS3='sal-s4.gse' KOLOR3=4
  SS4='sal-s5.gse' KOLOR4=4
  SS5='sal-s6.gse' KOLOR5=1
  SS6='sal-s7.gse' KOLOR6=1
  sp:         # seismogram plotting

# Rewriting the two-point rays into a formatted file
  REC='sal-rec1.dat' CRTOUT='crtout.tmp'
  echo.pl:"'r01.out'   ' '   'r01i.out' ' '" ">crtout.tmp"
  RAYS='sal-r1.out'  crtray:
  echo.pl:"'r02.out'   ' '   'r02i.out' ' '" ">crtout.tmp"
  RAYS='sal-r2.out'  crtray:
  echo.pl:"'r03.out'   ' '   'r03i.out' ' '" ">crtout.tmp"
  RAYS='sal-r3.out'  crtray:
  echo.pl:"'r04.out'   ' '   'r04i.out' ' '" ">crtout.tmp"
  RAYS='sal-r4.out'  crtray:
  echo.pl:"'r05.out'   ' '   'r05i.out' ' '" ">crtout.tmp"
  RAYS='sal-r5.out'  crtray:
  echo.pl:"'r06.out'   ' '   'r06i.out' ' '" ">crtout.tmp"
  RAYS='sal-r6.out'  crtray:
  echo.pl:"'r07.out'   ' '   'r07i.out' ' '" ">crtout.tmp"
  RAYS='sal-r7.out'  crtray:
  echo.pl:"'r08.out'   ' '   'r08i.out' ' '" ">crtout.tmp"
  RAYS='sal-r8.out'  crtray:

# Adding two-point rays into GOCAD file
  VRML='GOCAD' WRL=' ' WRLOUT='sal-tmp.mx'  # GOCAD
# Setting blue colour for rays
  KOLLIN=0  R=0.00  G=0.00  B=1.00
  NAME='Rays'  PROPERTIES=' '   # just for GOCAD
  LIN='sal-r1.out'
  NAME='Rays1-refracted'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r2.out'
  NAME='Rays2-top-salt-1'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r3.out'
  NAME='Rays3-top-salt-2'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r4.out'
  NAME='Rays4-bottom-salt-1'  PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r5.out'
  NAME='Rays5-bottom-salt-2'  PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r6.out'
  NAME='Rays6-bottom-1'       PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r7.out'
  NAME='Rays7-bottom-2'       PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  LIN='sal-r8.out'
  NAME='Rays8-bottom-3'       PROPERTIES=' '  # just for GOCAD
  linwrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD

# Adding source and receivers into GOCAD file
  PTS='sal-src1.dat'
  KOLPTS=0  R=1.00  G=0.00  B=1.00
  NAME='Source-1'        PROPERTIES=' ' ASIZE=6 # just for GOCAD
  ptswrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  PTS='sal-rec1.dat'
  KOLPTS=0  R=1.00  G=1.00  B=1.00
  NAME='Receivers-1'     PROPERTIES=' ' ASIZE=3 # just for GOCAD
  ptswrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD
  PTS='sal-rec1.plt'
  KOLPTS=0  R=1.00  G=0.00  B=0.00 SHIFT3=-0.1
  NAME='Receivers-1-plt' PROPERTIES=' ' ASIZE=5 # just for GOCAD
  ptswrl:
  append.pl: 'sal-ray1.mx' 'sal-tmp.mx'   # GOCAD

# To save disk space delete temporary output files
  del.pl: '*.out'

# Important output files
# ~~~~~~~~~~~~~~~~~~~~~~
# 'rpcode*.ps'  ... PostScript files with ray parameters
# 'sal-sf*.ps'  ... PostScript files with seismograms
# 'sal-ray1.mx' ... GOCAD file with two-point rays,
#                   shot and receivers