# 
# History file 'kis-pfa.h' for computation of S-wave travel times and
# synthetic seismograms in model KIS using the prevailing-frequency
# approximation of the coupling-ray theory.
# Comparing synthetic seismograms calculated by program CRTPFA
# with synthetic seismograms calculated by program GREEN.
# Calculating coupling-ray-theory travel times of S1 and S2 waves
# and their differences.

# 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/" "kis-crt.dat"
  chk.pl: "data/qi/" "kis-rpaa.dat"
  chk.pl: "data/qi/" "force.dat"
  chk.pl: "crt/"     "writall.dat"

# 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)'
  NORMSP=1     SPAMP=1.00
  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.002000  NFFT=2048
  OF=       NF=        DF=

# Coupling ray theory synthetic seismograms along two-point rays,
# anisotropic common ray,
# second-order perturbations of travel time,
# prevailing-frequency approximation
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  SRC='qi-src.dat'      REC='qi-rec0.dat'     SOURCE='force.dat'
  DCRT='kis-crt.dat'    CODE='kis-crt.dat'    WRIT='writall.dat'
  INIPAR=3
  MODEL='kis-mod.dat' CRTANI=1  RPAR='kis-rpaa.dat' CRTLOG='kis-crta.out'
  QIRAY=1
  CRTSCRO='SWR'
  crt:
#  RPPLOT='g.ps'
#  rpplot:
  SS='kis-sscf.gse'
  SINGLF=50.
  green:
  greenss:
  ss:
  SP1='kis-ss1f.ps'  SP2='kis-ss2f.ps'  SP3='kis-ss3f.ps'
  sp:
  QIRAY=
# Coupling ray theory along all successful rays
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTOUT='crtout.tmp'  CRTNEW='crtnew.tmp'
  echo.pl: "'r01.out'     's01.out'     'r01i.out'     't01.out'    / " "> crtout.tmp"
  echo.pl: "'r01-pfa.out' 's01-pfa.out' 'r01i-pfa.out' 't01.out'    / " "> crtnew.tmp"
  crtpfa:
  FTRGL1='t01-ok.out'  FTRGL2='t01-ko.out'
  CRTOUT='crtnew.tmp'
  pfatubes:
# Green function from the quantities already computed along two-point rays
  CRTOUT='crtnew.tmp'  MODEL=' '
  echo.pl: "'r01-pfa.out' 's01-pfa.out' 'r01i-pfa.out' 't01-ok.out' / " "> crtnew.tmp"
  SS='kis-sscp.gse'
  green:
  greenss:
  ss:
  SP1='kis-ss1p.ps'  SP2='kis-ss2p.ps'  SP3='kis-ss3p.ps'
  sp:
# Comparison of the seismograms
  SS='kis-sscf.gse' SRC=   REC=
  SS1='kis-sscp.gse'
  SP1='kis-s1fp.ps'  SP2='kis-s2fp.ps'  SP3='kis-s3fp.ps'
  sp:
# Green function by interpolation within ray tubes
  PTS='qi-rec0.dat'  MTTPTS='mtt.out'
  COLUMN01=X1  COLUMN02=X2  COLUMN03=X3  COLUMN04=MX4  COLUMN05=MX5
  COLUMN06=MX6  COLUMN07=MTT  COLUMN08=MTI  COLUMN09=MP1  COLUMN10=MP2
  COLUMN11=MP3  COLUMN12=MP4  COLUMN13=MP5  COLUMN14=MP6  COLUMN15=AMPR11
  COLUMN16=AMPI11  COLUMN17=AMPR21  COLUMN18=AMPI21  COLUMN19=AMPR31
  COLUMN20=AMPI31  COLUMN21=AMPR12  COLUMN22=AMPI12  COLUMN23=AMPR22
  COLUMN24=AMPI22  COLUMN25=AMPR32  COLUMN26=AMPI32  COLUMN27=AMPR13
  COLUMN28=AMPI13  COLUMN29=AMPR23  COLUMN30=AMPI23  COLUMN31=AMPR33
  COLUMN32=AMPI33
  mtt:
#  porovnej:
  PTS=               MTTPTS=

# Coupling ray theory along all rays
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTOUT='crtout.tmp'  CRTNEW='crtnew.tmp'
  echo.pl: "'r01.out'     ' '           'r01i.out'     't01.out'    / " "> crtout.tmp"
  echo.pl: "'r01-pfa.out' ' '           'r01i-pfa.out' 't01.out'    / " "> crtnew.tmp"
  MODEL='kis-mod.dat'
  crtpfa:
  FTRGL1='t01-ok.out'  FTRGL2='t01-ko.out'
  CRTOUT='crtnew.tmp'
  pfatubes:
  CRTOUT='crtnew.tmp'  MODEL=' '
#  echo.pl: "'r01-pfa.out' ' '           'r01i-pfa.out' 't01-ko.out' / " "> crtnew.tmp"
#  RPPLOT='gnew-ko.ps'
#  rpplot:
#  echo.pl: "'r01-pfa.out' ' '           'r01i-pfa.out' 't01-ok.out' / " "> crtnew.tmp"
#  RPPLOT='gnew-ok.ps'
#  rpplot:
# Interpolation of coupling-ray-theory travel times within ray cells
  echo.pl: "'r01-pfa.out' ' '           'r01i-pfa.out' 't01-ok.out' / " "> crtnew.tmp"
  O1=50    O2=50  O3=0
  D1=0.001 D2=    D3=0.001
  N1=1001  N2=1   N3=601
  mtt:
# Plotting numbers of arrivals at the gridpoints
  VPLUS=2  CREF=0.166667  YSIGN=-1
  GRD='mtt-num.out' PS='kis-num.ps'
  grdps:
# Plotting travel times at the gridpoints, sorted according to their value
  VPLUS=   VCIRC=   VREF=0  CREF=0.166667 VSIGN=1.2
  NUM='mtt-num.out' MGRD='mtt-tt.out' GRD='kis-t00.out'
  mgrd:
  GRD='kis-t00.out' PS='kis-t01.ps'
  grdps:
  N4=
# Calculating and plotting travel-time differences between S1 and S2 waves
  IMGRD=1  GRD='kis-t01.out'
  mgrd:
  echo.pl: "DIF=$1-$2 " "> reldife.cal"
  echo.pl: "SUM=$1+$2 " ">>reldife.cal"
  echo.pl: "$3=DIF/SUM" ">>reldife.cal"
  CAL='reldife.cal' GRD1='kis-t02.out' GRD2='kis-t01.out' GRD3='kis-tdif.out'
  grdcal:
  GRD1= GRD2= GRD3=
  VPLUS=  VCIRC=0.021 VREF=0  CREF=0.166667  VSIGN=
  GRD='kis-tdif.out'  PS='kis-tdif.ps'
  grdps:

# List of the output files
# ~~~~~~~~~~~~~~~~~~~~~~~~
# 'kis-ss1?.ps','kis-ss2?.ps','kis-ss3?.ps' ... Synthetic seismograms
#       where ? identifies the calculation
#             f ... prevailing-frequency approximation of the coupling ray
#                   theory, the S1 and S2 Green tensors are calculated by
#                   program GREEN considering anisotropic model along the
#                   reference common anisotropic S-wave rays
#             p ... prevailing-frequency approximation of the coupling ray
#                   theory, the S1 and S2 Green tensors are calculated using
#                   program CRTPFA along the two-point rays and then program
#                   GREEN
# 'kis-s1fp.ps','kis-s2fp.ps','kis-s3fp.ps' ... Comparison of the seismograms.
#       Colors of seismograms are
#         Black - Seismograms calculated by GREEN
#         Red - Seismograms calculated by CRTPFA and GREEN
# 'kis-t01.ps','kis-t02.ps','kis-tdif.ps' ... Interpolated travel times
#       of the S1 and S2 waves and travel-time differences.
#       Figures used in Klimes & Bulant (2013):
#       Interpolation of the coupling-ray-theory S-wave Green tensor
#       within ray cells.
#       In: Seismic Waves in Complex 3-D Structures, Report 23,
#       Dep. Geophys., Charles Univ., Prague, pp. 203-218, fig. 5,
#       online at "http://sw3d.cz".
#