#
# History file 'sc-ss-an.h' for computation of anisotropic ray
# theory S-wave synthetic seismograms in models SC1_I, SC1_II and ORT
# using the program package ANRAY.
# The seismograms are compared with results calculated by Fourier
# Method.

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: "anray/data/qi" "anraysc1.qi"
  chk.pl: "anray/data/qi" "anraysc2.qi"
  chk.pl: "anray/data/qi" "anrayort.qi"
  chk.pl: "anray/data/qi" "fresan.scv"
  chk.pl: "anray/data/qi" "fresan.sct"
  chk.pl: "anray/data/qi" "fresan.scr"
  chk.pl: "anray/data/qi" "synfan.sc"
# 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"

# 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 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=

# =====================================================================
# Model SC1_I
# =====================================================================

# Anisotropic ray theory synthetic seismograms - package ANRAY
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  REC=
  SPXMAX=1.60  SPXDIV=8
  SPTMAX=0.70  SPTDIV=7   SPTSUB=5
  anray:   'anraysc1.qi' 1* 'lu1.out' 'lu2.out' /

# Vertical component
  fresan:  'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  LU8='lu8.out'
  SS='sc1-ss3a.gse'
  anraygse:
  SS='sc1-z-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc1-ss3a.gse' SPAMP1=1.              KOLOR1=4
  SP1=' '           SP2=' '           SP3='sc1-ss3a.ps'
  sp:

# Transverse component
  fresan:  'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='sc1-ss2a.gse'
  anraygse:
  SS='sc1-y-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc1-ss2a.gse' SPAMP1=1.              KOLOR1=4
  SP1=' '           SP2='sc1-ss2a.ps' SP3=' '
  sp:

# Radial component
  fresan:  'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='sc1-ss1a.gse'
  anraygse:
  SS='sc1-x-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc1-ss1a.gse' SPAMP1=1.              KOLOR1=4
  SP1='sc1-ss1a.ps' SP2=' '           SP3=' '
  sp:
  SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1=
  SPXMAX=1.40  SPXDIV=7
  SPTMAX=0.65  SPTDIV=3   SPTSUB=10

# =====================================================================
# Model SC1_II
# =====================================================================

# Anisotropic ray theory synthetic seismograms - package ANRAY
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  REC=
  SPXMAX=1.60  SPXDIV=8
  SPTMAX=0.70  SPTDIV=7   SPTSUB=5
  anray:   'anraysc2.qi' 1* 'lu1.out' 'lu2.out' /

# Vertical component
  fresan:  'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  LU8='lu8.out'
  SS='sc2-ss3a.gse'
  anraygse:
  SS='sc2-z-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc2-ss3a.gse' SPAMP1=1.              KOLOR1=4
  SP1=' '           SP2=' '           SP3='sc2-ss3a.ps'
  sp:

# Transverse component
  fresan:  'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='sc2-ss2a.gse'
  anraygse:
  SS='sc2-y-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc2-ss2a.gse' SPAMP1=1.              KOLOR1=4
  SP1=' '           SP2='sc2-ss2a.ps' SP3=' '
  sp:

# Radial component
  fresan:  'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='sc2-ss1a.gse'
  anraygse:
  SS='sc2-x-fm.gse'  SPAMP=100000000000000. KOLOR=1
  SS1='sc2-ss1a.gse' SPAMP1=1.              KOLOR1=4
  SP1='sc2-ss1a.ps' SP2=' '           SP3=' '
  sp:
  SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1=
  SPXMAX=1.40  SPXDIV=7
  SPTMAX=0.65  SPTDIV=3   SPTSUB=10

# =====================================================================
# Model ORT
# =====================================================================

# Anisotropic ray theory synthetic seismograms - package ANRAY
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  REC=
  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
  anray:   'anrayort.qi' 1* 'lu1.out' 'lu2.out' /

# Vertical component
  fresan:  'fresan.scv' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  LU8='lu8.out'
  SS='ort-ss3a.gse'
  anraygse:
  SS='ort-z-fm.gse'  SPAMP=15000000000000.  KOLOR=1
  SS1='ort-ss3a.gse' SPAMP1=0.15            KOLOR1=4
  SP1=' '           SP2=' '           SP3='ort-ss3a.ps'
  sp:

# Transverse component
  fresan:  'fresan.sct' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='ort-ss2a.gse'
  anraygse:
  SS='ort-y-fm.gse'  SPAMP=15000000000000.  KOLOR=1
  SS1='ort-ss2a.gse' SPAMP1=0.15            KOLOR1=4
  SP1=' '           SP2='ort-ss2a.ps' SP3=' '
  sp:

# Radial component
  fresan:  'fresan.scr' 1* 'lu2.out' 'lu7.out' ' ' /
  synfan:  'synfan.sc'  1* 'lu7.out' 'lu8.out' /
  SS='ort-ss1a.gse'
  anraygse:
  SS='ort-x-fm.gse'  SPAMP=15000000000000.  KOLOR=1
  SS1='ort-ss1a.gse' SPAMP1=0.15            KOLOR1=4
  SP1='ort-ss1a.ps' SP2=' '           SP3=' '
  sp:
  SPAMP=1. SPAMP1= KOLOR= KOLOR1= SS= SS1=
  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

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# '???-ss1a.ps','???-ss2a.ps','???-ss3a.ps' ... Anisotropic ray theory
#       synthetic seismograms compared with FM.
# "???" in the names of files identifies the model
#         sc1 ... model SC1_I
#         sc2 ... model SC1_II
#         ort ... model ORT
# Colors of seismograms '*ss?a.ps' are
#         Black - Fourier Method
#         Blue - Anisotropic ray theory
# Figures correspond to Psencik, Farra & Tessmer (2011):
# Comparison of the FORT approximation of the coupling ray theory with
# the Fourier pseudospectral method
# In: Seismic Waves in Complex 3-D Structures, Report 21,
# Dep. Geophys., Charles Univ., Prague, pp. 129-165.
#