# History file 'kis-ss.h' for computation of coupling-ray
# theory S-wave synthetic seismograms in model KISS.
# The seismograms are compared with results calculated by FORT
# and by Fourier Method.

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: "data/qi/" "kis-mod.dat"
  chk.pl: "data/qi/" "qi-src.dat"
  chk.pl: "data/qi/" "qi-rec0.dat"
  chk.pl: "data/qi/" "qi-crt.dat"
  chk.pl: "data/qi/" "qi-rpaa.dat"
  chk.pl: "data/qi/" "force.dat"
  chk.pl: "crt/"     "writ.dat"
# GSE files for FORT - results of "data/qi/fort/fort2gse.h"
  chk.pl: ""         "kis-x-fo.gse"
  chk.pl: ""         "kis-y-fo.gse"
  chk.pl: ""         "kis-z-fo.gse"
# GSE files for FM - results of "data/qi/fm/fm2gse.h"
  chk.pl: ""         "kis-x-fm.gse"
  chk.pl: ""         "kis-y-fm.gse"
  chk.pl: ""         "kis-z-fm.gse"

# Data to control seismogram plotting (program SP)
  KODESP=3     SPCHRH=0.25
  SPTMIN=0.40  SPTMAX=0.55  SPTLEN=6.00   SPTDIV=3   SPTSUB=5
  SPXMIN=0.00  SPXMAX=0.60  SPXLEN=15.0   SPXDIV=6   SPXSUB=5
  SPTTEXT='Time (s)'  SPXTEXT='Depth (km)'           SPTDEC=2
  NORMSP=1     SPAMP=1.00
  CALCOPS='0.50 setlinewidth'
  GSEWIDTH=8 (writing GSE files in greenss.for)

# Coupling ray theory synthetic seismograms - anisotropic common ray
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  SRC='qi-src.dat'      REC='qi-rec0.dat'     SOURCE='force.dat'
  DCRT='qi-crt.dat'     CODE='qi-crt.dat'     WRIT='writ.dat'
  INIPAR=3

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

# Data describing the filtration of the source time function
  FMIN=0.0  FLOW=5.0  FHIGH=60.  FMAX=100.  # (cosine filter)

# Time and frequency step and interval for the Fast Fourier Transform
# DT=0.000488  NFFT=2048
# OF=0.0    NF=101
# DF=1.000576332
  DT=0.002000  NFFT=2048
  OF=       NF=        DF=

  MODEL='kis-mod.dat' CRTANI=1  RPAR='qi-rpaa.dat'  CRTLOG='qi-crta.out'
  SS='kis-ssc.gse'
  crt:
  green:
  greenss:
  ss:
  QIRAY=1
  SS='kis-sscs.gse'
  crt:
  green:
  greenss:
  ss:
  SS='kis-sscf.gse'
  crt:
  SINGLF=50.
  green:
  SINGLF=
  greenss:
  ss:
  QIRAY=

# FM plus FORT plus both coupling ray theory seismograms

  KOLOR=3 (green) KOLOR1=2 (red)    KOLOR2=4 (blue)    KOLOR3=1 (black)
  SS='kis-ssc.gse'
  SS1='kis-x-fo.gse'
  SS2='kis-sscs.gse'
  SS3='kis-x-fm.gse' SPAMP3=100000000000000.
  SP1='kis-ss1.ps'   SP2=' '            SP3=' '
  sp:
  SS1='kis-y-fo.gse'
  SS3='kis-y-fm.gse'
  SP1=' '            SP2='kis-ss2.ps'   SP3=' '
  sp:
  SS1='kis-z-fo.gse'
  SS3='kis-z-fm.gse'
  SP1=' '            SP2='           '  SP3='kis-ss3.ps'
  sp:
  KOLOR= KOLOR1= KOLOR2= KOLOR3=
  SS=    SS1=    SS2=    SS3=
  SPAMP= SPAMP1= SPAMP2= SPAMP3=

# Subtracting the seismograms

  DT=0.
  SS1='kis-x-fm.gse' COEF1=-100000000000000.
  SS2='kis-x-fo.gse' COEF2=1.
  SS='kis-fo1.gse'
  gsecal:
  SS1='kis-y-fm.gse'  SS2='kis-y-fo.gse'  SS='kis-fo2.gse'   gsecal:
  SS1='kis-z-fm.gse'  SS2='kis-z-fo.gse'  SS='kis-fo3.gse'   gsecal:
  SS1='kis-x-fm.gse'  SS2='kis-ssc.gse'   SS='kis-cr1.gse'   gsecal:
  SS1='kis-y-fm.gse'  SS2='kis-ssc.gse'   SS='kis-cr2.gse'   gsecal:
  SS1='kis-z-fm.gse'  SS2='kis-ssc.gse'   SS='kis-cr3.gse'   gsecal:
  SS1='kis-x-fm.gse'  SS2='kis-sscs.gse'  SS='kis-crs1.gse'  gsecal:
  SS1='kis-y-fm.gse'  SS2='kis-sscs.gse'  SS='kis-crs2.gse'  gsecal:
  SS1='kis-z-fm.gse'  SS2='kis-sscs.gse'  SS='kis-crs3.gse'  gsecal:
  SPAMP=2.
  KOLOR=3 (green) KOLOR1=2 (red)    KOLOR2=4 (blue)    KOLOR3=
  SS='kis-cr1.gse'
  SS1='kis-fo1.gse'
  SS2='kis-crs1.gse'
  SP1='kis-sd1.ps'   SP2=' '            SP3=' '
  sp:
  SS='kis-cr2.gse'
  SS1='kis-fo2.gse'
  SS2='kis-crs2.gse'
  SP1=' '            SP2='kis-sd2.ps'   SP3=' '
  sp:
  SS='kis-cr3.gse'
  SS1='kis-fo3.gse'
  SS2='kis-crs3.gse'
  SP1=' '            SP2=' '            SP3='kis-sd3.ps'
  sp:
  KOLOR= KOLOR1= KOLOR2= KOLOR3=
  SS=    SS1=    SS2=    SS3=
  SPAMP= SPAMP1= SPAMP2= SPAMP3=

# FM plus coupling ray theory plus single-frequency approximation seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  KOLOR=2 red   KOLOR1=3 green   KOLOR2=1 black
  SS='kis-sscf.gse'  # single-frequency approximation
  SS1='kis-sscs.gse' # coupling ray theory
  SS2='kis-x-fm.gse' SPAMP2=100000000000000.     # FM
  SP1='kis-sf1.ps'   SP2=' '            SP3=' '
  sp:
  SS2='kis-y-fm.gse'
  SP1=' '            SP2='kis-sf2.ps'   SP3=' '
  sp:
  SS2='kis-z-fm.gse'
  SP1=' '            SP2='           '  SP3='kis-sf3.ps'
  sp:

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# 'kis-ss1.ps','kis-ss2.ps','kis-ss3.ps' ... Synthetic seismograms
#       compared with FORT and with FM.
# 'kis-sd1.ps','kis-sd2.ps','kis-sd3.ps' ... Differences of the
#       seismograms from the FM seismograms.
#       Colors of seismograms are
#         Black - Fourier Method
#         Red - First Order Ray Tracing
#         Green - Coupling Ray Theory along anisotropic common rays
#         Blue - CRT with second-order travel time perturbation
# 'kis-sf1.ps','kis-sf2.ps','kis-sf3.ps' ... Single-frequency approximation
#       seismograms compared with standard coupling ray theory and with FM.
#       Colors of seismograms are
#         Black - Fourier Method
#         Red - Single-frequency approximation of the Coupling Ray Theory
#               along anisotropic common rays
#         Green - CRT with second-order travel time perturbation