# History file 'sc-ss.h' for computation of coupling-ray
# theory S-wave synthetic seismograms in models SC1_I, SC1_II and ORT.
# The seismograms are compared with results calculated by FORT
# and by Fourier Method.

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

# Coupling ray theory synthetic seismograms - anisotropic common ray
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

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

# Data for ray tracing
  SRC='sc1-src.dat'     REC='sc1-rec.dat'     SOURCE='force.dat'
  DCRT='sc1-crt.dat'    CODE='sc1-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.000250  NFFT=2048
  OF=       NF=        DF=


# Ray tracing and synthetic seismograms - model SC1_I
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='sc1-mod.dat'
# Common anisotropic ray
  CRTANI=1  RPAR='sc1-rpaa.dat' CRTLOG='sc1-crta.out'
  SS='sc1-ssc.gse'
  crt:
  green:
  greenss:
  ss:
# Common anisotropic ray, second-order perturbation of travel time
  QIRAY=1
  SS='sc1-sscs.gse'
  crt:
  green:
  greenss:
  ss:
# Prevailing-frequency approximation
  SS='sc1-sscf.gse'
  SINGLF=50.
  green:
  SINGLF=
  greenss:
  ss:
  QIRAY=
# SH and SV rays, prevailing-frequency approximation
  CRTLOG='sc1-crt1.out'
  KSWAVE=1  QIRAY=0  CRTANI=1
  TIA1=0.8660254  TIA2=0.5  TIA3=0.0  DEGREE=2
  crt:
  SINGLF=50.
  GREEN='sc1-ss01.out'
  green:
  KSWAVE=2  CRTLOG='sc1-crt2.out'
  crt:
  GREEN='sc1-ss02.out'
  green:
  GREEN1='sc1-ss01.out'
  GREEN2='sc1-ss02.out'
  GREEN='sc1-ss12.out'
  greenti:
  SS='sc1-ss12.gse'
  greenss:
  ss:
  QIRAY=
  KSWAVE=
  TIA1=  TIA2=  TIA3=  DEGREE=

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

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

# FM plus coupling ray theory plus prevailing-frequency approximation
# seismograms
#
  KOLOR=2 red   KOLOR1=3 green   KOLOR2=1 black
  SS='sc1-sscf.gse'  # prevailing-frequency approximation
  SS1='sc1-sscs.gse' # coupling ray theory
  SS2='sc1-x-fm.gse' SPAMP2=100000000000000.     # FM
  SP1='sc1-sf1.ps'   SP2=' '            SP3=' '
  sp:
  SS2='sc1-y-fm.gse'
  SP1=' '            SP2='sc1-sf2.ps'   SP3=' '
  sp:
  SS2='sc1-z-fm.gse'
  SP1=' '            SP2='           '  SP3='sc1-sf3.ps'
  sp:
  KOLOR= KOLOR1= KOLOR2= KOLOR3=
  SS=    SS1=    SS2=    SS3=
  SPAMP= SPAMP1= SPAMP2= SPAMP3=
# Only FM plus prevailing-frequency approximation seismograms
 #  KOLOR=2 (red) KOLOR1=1 (black) KOLOR2=
 #  SS='sc1-sscf.gse'
 #  SS2=               SPAMP2=
# Only FM plus SH and SV prevailing-frequency approximation seismograms
 #  KOLOR=2 (red) KOLOR1=1 (black) KOLOR2=
 #  SS='sc1-ss12.gse'
 #  SS2=               SPAMP2=
 #  SS1='sc1-x-fm.gse' SPAMP1=100000000000000.
 #  SP1='sc1-sf1.ps'   SP2=' '            SP3=' '
 #  sp:
 #  SS1='sc1-y-fm.gse'
 #  SP1=' '            SP2='sc1-sf2.ps'   SP3=' '
 #  sp:
 #  SS1='sc1-z-fm.gse'
 #  SP1=' '            SP2='           '  SP3='sc1-sf3.ps'
 #  sp:
 #  KOLOR= KOLOR1= KOLOR2= KOLOR3=
 #  SS=    SS1=    SS2=    SS3=
 #  SPAMP= SPAMP1= SPAMP2= SPAMP3=

## Auxiliary calculations for debugging purposes
## Values at the receivers
#  crtpts:
## Pictures of the ray parameter domain and of the reference surface
#  CRTOUT=  CALCOPS='0.01 setlinewidth'
#  IRBAS=  IRTWO=  IRAUX=  ITHOM=  ISANG=1 ISHP=  ISUC=1
#  PLIM1=0.  PLIM2=1.  PLIM3=0.  PLIM4=1.
#  RPPLOT='g.ps'
#  rpplot:
#  IRBAS=  IRTWO=  IRAUX=  ITHOM=  ISANG=0 ISHP=  ISUC=
#  PLIM1=0.  PLIM2=1.5  PLIM3=-1. PLIM4=1.
#  RPPLOT='x.ps'
#  rpplot:


# Ray tracing and synthetic seismograms - model SC1_II
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='sc2-mod.dat'
# Common anisotropic ray
  CRTANI=1  RPAR='sc1-rpaa.dat' CRTLOG='sc2-crta.out'
  SRC='sc1-src.dat'    REC='sc1-rec.dat'
  DCRT='sc2-crt.dat'   CODE='sc2-crt.dat'
  SS='sc2-ssc.gse'
  crt:
  green:
  greenss:
  ss:
# Common anisotropic ray, second-order perturbation of travel time
  QIRAY=1
  SS='sc2-sscs.gse'
  crt:
  green:
  greenss:
  ss:
# Prevailing-frequency approximation
  SS='sc2-sscf.gse'
  SINGLF=50.
  green:
  SINGLF=
  greenss:
  ss:
  QIRAY=
# SH and SV rays, prevailing-frequency approximation
  TIA1=0.866   TIA2=0.5  TIA3=1.  DEGREE=2  QIRAY=0
  KSWAVE=1  CRTLOG='sc2-crt1.out'
  crt:
  SINGLF=50.
  GREEN='sc2-ss01.out'
  green:
  KSWAVE=2  CRTLOG='sc2-crt2.out'
  crt:
  GREEN='sc2-ss02.out'
  green:
  GREEN1='sc2-ss01.out'
  GREEN2='sc2-ss02.out'
  GREEN='sc2-ss12.out'
  greenti:
  SS='sc2-ss12.gse'
  greenss:
  ss:
  QIRAY=  KSWAVE=
  TIA1=   TIA2=     TIA3=    DEGREE=

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

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

# FM plus coupling ray theory plus prevailing-frequency approximation
# seismograms
#
  KOLOR=2 red   KOLOR1=3 green   KOLOR2=1 black
  SS='sc2-sscf.gse'  # prevailing-frequency approximation
  SS1='sc2-sscs.gse' # coupling ray theory
  SS2='sc2-x-fm.gse' SPAMP2=100000000000000.     # FM
  SP1='sc2-sf1.ps'   SP2=' '            SP3=' '
  sp:
  SS2='sc2-y-fm.gse'
  SP1=' '            SP2='sc2-sf2.ps'   SP3=' '
  sp:
  SS2='sc2-z-fm.gse'
  SP1=' '            SP2='           '  SP3='sc2-sf3.ps'
  sp:
  KOLOR= KOLOR1= KOLOR2= KOLOR3=
  SS=    SS1=    SS2=    SS3=
  SPAMP= SPAMP1= SPAMP2= SPAMP3=
# Only FM plus prevailing-frequency approximation seismograms
 #  KOLOR=2 (red) KOLOR1=1 (black) KOLOR2=
 #  SS='sc2-sscf.gse'
 #  SS2=               SPAMP2=
# Only FM plus SH and SV prevailing-frequency approximation seismograms
 #  KOLOR=2 (red) KOLOR1=1 (black) KOLOR2=
 #  SS='sc2-ss12.gse'
 #  SS2=               SPAMP2=
 #  SS1='sc2-x-fm.gse' SPAMP1=100000000000000.
 #  SP1='sc2-sf1.ps'   SP2=' '            SP3=' '
 #  sp:
 #  SS1='sc2-y-fm.gse'
 #  SP1=' '            SP2='sc2-sf2.ps'   SP3=' '
 #  sp:
 #  SS1='sc2-z-fm.gse'
 #  SP1=' '            SP2='           '  SP3='sc2-sf3.ps'
 #  sp:
 #  KOLOR= KOLOR1= KOLOR2= KOLOR3=
 #  SS=    SS1=    SS2=    SS3=
 #  SPAMP= SPAMP1= SPAMP2= SPAMP3=

# Ray tracing and synthetic seismograms - model ORT
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='ort-mod.dat'  CRTANI=1  RPAR='ort-rpaa.dat' CRTLOG='ort-crta.out'
  SRC='sc1-src.dat'    REC='ort-rec.dat'
  DCRT='ort-crt.dat'   CODE='ort-crt.dat'
  SPTMIN=0.65  SPTMAX=1.15  SPTLEN=6.0   SPTDIV=5   SPTSUB=10
  SPXMIN=0.00  SPXMAX=1.60  SPXLEN=15.0  SPXDIV=8   SPXSUB=5
  DT=0.002000  NFFT=2048
  OF=       NF=        DF=
  SS='ort-ssc.gse'
  crt:
  green:
  greenss:
  ss:
  QIRAY=1
  SS='ort-sscs.gse'
  crt:
  green:
  greenss:
  ss:
  SS='ort-sscf.gse'
  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)
  SPAMP=0.15
  SS='ort-ssc.gse'
  SS1='ort-x-fo.gse'
  SS2='ort-sscs.gse'
  SS3='ort-x-fm.gse' SPAMP3=15000000000000.
  SP1='ort-ss1.ps'   SP2=' '            SP3=' '
  sp:
  SS1='ort-y-fo.gse'
  SS3='ort-y-fm.gse'
  SP1=' '            SP2='ort-ss2.ps'   SP3=' '
  sp:
  SS1='ort-z-fo.gse'
  SS3='ort-z-fm.gse'
  SP1=' '            SP2='           '  SP3='ort-ss3.ps'
  sp:
  KOLOR= KOLOR1= KOLOR2= KOLOR3=
  SS=    SS1=    SS2=    SS3=
  SPAMP= SPAMP1= SPAMP2= SPAMP3=

# Subtracting the seismograms

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

# FM plus coupling ray theory plus prevailing-frequency approximation
# seismograms

  KOLOR=2 red   KOLOR1=3 green   KOLOR2=1 black
  SPAMP=0.15
  SS='ort-sscf.gse'  # prevailing-frequency approximation
  SS1='ort-sscs.gse' # coupling ray theory
  SS2='ort-x-fm.gse' SPAMP2=15000000000000.      # FM
  SP1='ort-sf1.ps'   SP2=' '            SP3=' '
  sp:
  SS2='ort-y-fm.gse'
  SP1=' '            SP2='ort-sf2.ps'   SP3=' '
  sp:
  SS2='ort-z-fm.gse'
  SP1=' '            SP2='           '  SP3='ort-sf3.ps'
  sp:

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# '???-ss1.ps','???-ss2.ps','???-ss3.ps' ... Synthetic seismograms
#       compared with FORT and with FM.
# '???-sd1.ps','???-sd2.ps','???-sd3.ps' ... Differences of the
#       seismograms from the FM seismograms.
#       "???" identifies the model
#         sc1 ... model SC1_I
#         sc2 ... model SC1_II
#         ort ... model ORT
# '???-sf1.ps','???-sf2.ps','???-sf3.ps' ... Prevailing-frequency
#       approximation seismograms compared with standard coupling ray
#       theory and with FM.
# Colors of seismograms '*ss*.ps' 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
# Colors of seismograms '*sf*.ps' are
#         Black - Fourier Method
#         Red - Single Frequency approximation of coupling ray theory
#         Green - CRT with second-order travel time perturbation