# History file '98-ss.h' to compute synthetic seismograms in model 98
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required
  chk.pl: "data/98/" "98-crt2.dat"
  chk.pl: "data/98/" "98-mod.dat"
  chk.pl: "data/98/" "98-src.dat"
  chk.pl: "data/98/" "98-rec.dat"
  chk.pl: "data/98/" "98-rpa2.dat"
  chk.pl: "crt/"     "writ.dat"
  chk.pl: "crt/"     "source.dat"
  chk.pl: "forms/"   "hsv.dat"
  chk.pl: "forms/"   "echo.pl"
# Output file of '98-mod.h' used here as an input file
  chk.pl:  ""         "98-srf.wrl"  # VRML
# chk.pl: ""         "98-srf.mx"    # GOCAD

# Model
  MODEL='98-mod.dat'

# Source and receivers
  SRC='98-src.dat'  REC='98-rec.dat'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='98-crt2.dat'
  CODE='98-crt2.dat'
  RPAR='98-rpa2.dat'
  WRIT='writ.dat'
  CRTLOG='98-crt2.out'

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

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

# Data describing the filtration of the source time function
  FMIN=4.  FLOW=20.  FHIGH=184.  FMAX=200. # (cosine filter)

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

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

# Data to control seismogram plotting (program SP)
  KODESP=2    SPCHRH=0.40
  SPTMIN=0.5  SPTMAX=3.0  SPTLEN=25.0  SPTDIV=5    SPTSUB=5
  SPXMIN=0.8  SPXMAX=4.2  SPXLEN=15.3  SPXDIV=-17  SPXSUB=1
  SPYMIN=0.8  SPYMAX=4.2
  NORMSP=1    SPAMP=-2500

# Data to control ray parameters plotting (program RPPLOT)
  ITHOM=0 ISUC=1  HRBAS=0.2 HRTWO=0.3 HRAUX=0.1 IWAVES=4
  echo.pl: "' '   ' '   's01i.out' 't04.out'" ">crtout.tmp"
  CRTOUT='crtout.tmp' RPPLOT='98-plot1.ps'

# Programs to calculate Green function
  crt:        # two-point ray tracing
  rpplot:     # ray parameters plotting
  CRTOUT=
  green:      # Green function

# Frequency-domain seismograms
  SS='98-ssf.gse'  SP1='98-ssf1.ps'  SP2='98-ssf2.ps'  SP3='98-ssf3.ps'
  greenss:    # response function
  ss:         # synthetic seismograms
  sp:         # seismogram plotting

# Time-domain seismograms
  SIGGSE='98-sig.gse' SS=' '
  ss:         # source time function
  SS='98-sst.gse'  SP1='98-sst1.ps'  SP2='98-sst2.ps'  SP3='98-sst3.ps'
  greenss:    # synthetic seismograms
  sp:         # seismogram plotting

# Comparison of the seismograms
  SS='98-sst.gse'   KOLOR=2   # time-domain seismograms (red)
  SS1='98-ssf.gse'  KOLOR1=1  # frequency-domain seismograms (black)
  SP1='98-ssb1.ps'  SP2='98-ssb2.ps'  SP3='98-ssb3.ps'
  sp:         # seismogram plotting

# Exponentially scaled seismograms
  NORMSP=0  SPAMP=-1  SPEXP=5
  SP1='98-sse1.ps'  SP2='98-sse2.ps'  SP3='98-sse3.ps'
  sp:         # seismogram plotting

# Rewriting the two-point rays into a formatted file
  CRTOUT='crtout.tmp'
  echo.pl:"'r01.out'   ' '   'r01i.out' ' '" ">crtout.tmp"
  RAYS='98-ray1.out'                   crtray: '98-ss.h' /
  echo.pl:"'r02.out'   ' '   'r02i.out' ' '" ">crtout.tmp"
  RAYS='98-ray2.out'                   crtray: '98-ss.h' /
  echo.pl:"'r03.out'   ' '   'r03i.out' ' '" ">crtout.tmp"
  RAYS='98-ray3.out'                   crtray: '98-ss.h' /
  echo.pl:"'r04.out'   ' '   'r04i.out' ' '" ">crtout.tmp"
  RAYS='98-ray4.out'                   crtray: '98-ss.h' /
  CRTOUT=

# Adding two-point rays into the VRML or GOCAD file with surfaces
  VRML='VRML2'  WRL='98-srf.wrl'  WRLOUT='98-ray.wrl'  # VRML
# VRML='GOCAD'  WRL='98-srf.mx'   WRLOUT='98-ray.mx'   # GOCAD
# Setting rays blue
  KOLLIN=0  R=0.00  G=0.00  B=1.00
# Example of parameters for colouring rays according to travel time
# KOLLIN=4  COLORS='hsv.dat'  VPER=1.0  VREF=0.0  CREF=0.166667
  LIN='98-ray1.out'
  NAME='Rays1'  PROPERTIES=' '  # just for GOCAD
# NAME='Rays1'  PROPERTIES='T'  # just for GOCAD
  linwrl:
  WRL='98-ray.wrl'  # VRML
# WRL='98-ray.mx'   # GOCAD
  LIN='98-ray2.out'
  NAME='Rays2'  PROPERTIES=' '  # just for GOCAD
# NAME='Rays2'  PROPERTIES='T'  # just for GOCAD
  linwrl:
  LIN='98-ray3.out'
  NAME='Rays3'  PROPERTIES=' '  # just for GOCAD
# NAME='Rays3'  PROPERTIES='T'  # just for GOCAD
  linwrl:
  LIN='98-ray4.out'
  NAME='Rays4'  PROPERTIES=' '  # just for GOCAD
# NAME='Rays4'  PROPERTIES='T'  # just for GOCAD
  linwrl:

# Adding source and receivers into the VRML or GOCAD file
  PTS='98-src.dat'
  KOLPTS=0  R=1.00  G=0.00  B=1.00
  NAME='Source'  PROPERTIES=' '  # just for GOCAD
  ptswrl:
  PTS='98-rec.dat'
  KOLPTS=0  R=1.00  G=1.00  B=1.00  SHIFT3=0.001
  NAME='Receivers'  PROPERTIES=' '  # just for GOCAD
  ptswrl:

#=======================================================================
# List of output PostScript figures
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# '98-plot1.ps' Ray parameters of the fourth elementary wave.
# '98-sig.ps'   Sketch of source time function (output of SS).

#  Synthetic seismograms
# '98-ssf1.ps' Frequency-domain, 1st component.
# '98-ssf2.ps' Frequency-domain, 2nd component.
# '98-ssf3.ps' Frequency-domain, 3rd component.
# '98-sst1.ps' Time-domain, 1st component.
# '98-sst2.ps' Time-domain, 2nd component.
# '98-sst3.ps' Time-domain, 3rd component.
# '98-ssb1.ps' Both (comparison), 1st component.
# '98-ssb2.ps' Both (comparison), 2nd component.
# '98-ssb3.ps' Both (comparison), 3rd component.
# '98-sse1.ps' Exponential scale, 1st component.
# '98-sse2.ps' Exponential scale, 2nd component.
# '98-sse3.ps' Exponential scale, 3rd component.

# Output VRML or GOCAD file containing 3-D sculptures of rays
# and interfaces
# 98-ray.wrl or 98-ray.mx