# Synthetic seismograms in 2-D model P1I
# ======================================

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: "data/p1/" "p1i-mod.dat"
  chk.pl: "data/p1/" "p1-src.dat"
  chk.pl: "data/p1/" "p1-rec.dat"
  chk.pl: "data/p1/" "p1i-crt.dat"
  chk.pl: "data/p1/" "p1i-cod.dat"
  chk.pl: "data/p1/" "p1i-rpa.dat"
  chk.pl: "data/p1/" "p1-writ.dat"
  chk.pl: "data/p1/" "p1-bnd.dat"
  chk.pl: "crt/"     "source.dat"
  chk.pl: "forms/"   "echo.pl"
  chk.pl: "forms/"   "copy.pl"
  chk.pl: "forms/"   "append.pl"

# Data for ray tracing
  MODEL='p1i-mod.dat'
  SRC='p1-src.dat'
  REC='p1-rec.dat'
  DCRT='p1i-crt.dat'
  CODE='p1i-cod.dat'
  RPAR='p1i-rpa.dat'
  WRIT='p1-writ.dat'
  CRTLOG='p1i-log.out'
# Kind of ray parameters (for details see 'init.for')
  INIPAR=3

# Reference frequency for the attenuation
  TIMUL=1  FREF=100

# Data describing the source time function
  KSIG=1 # (Gabor signal)
  SIGT=0.  SIGF=100.  SIGW=4.  SIGPH=0.
  SIGA=258.174 # (258.174=12*pi*rho*v**3)

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

# GSE file with synthetic seismograms
  SS='p1i-ss.gse'

# Data to control seismogram plotting (program SP)
  SP1='p1i-ss1.ps'  SP2=' '           SP3='p1i-ss3.ps'
  KODESP=   SPCHRH=0.15
  SPTMIN=       SPTMAX=6.     SPTDIV=6.  SPTLEN=12.5
  SPXLEN=19
  NORMSP=1  SPAMP=1000
  CALCOPS='0.7 setlinewidth'

# Ray tracing
# ~~~~~~~~~~~
  CRTSCRO='swrg'
  CRTSCRO='swr'
  crt:

# Plotting the rays
# ~~~~~~~~~~~~~~~~~
  CRTOUT='crtout.tmp'
  echo.pl: "'r01.out' * 'r01i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r01.out'
  crtray:
  echo.pl: "'r02.out' * 'r02i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r02.out'
  crtray:
  echo.pl: "'r03.out' * 'r03i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r03.out'
  crtray:
  echo.pl: "'r04.out' * 'r04i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r04.out'
  crtray:
  echo.pl: "'r05.out' * 'r05i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r05.out'
  crtray:
  echo.pl: "'r06.out' * 'r06i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r06.out'
  crtray:
  echo.pl: "'r07.out' * 'r07i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r07.out'
  crtray:
  echo.pl: "'r08.out' * 'r08i.out' /  " ">  crtout.tmp"
   RAYS='p1i-r08.out'
  crtray:

  O1=0.        D1=0.1     N1=474
  O3=0.        D3=0.1     N3=61
  O2=          D2=        N2=
  O1=0.        D1=0.01    N1=4731
  O3=0.        D3=0.01    N3=601
  O1=0.        D1=0.05    N1=947
  O3=0.        D3=0.05    N3=121
  MPAR=1   VEL='p1i-vep.out'  ICB=  'p1i-icb.out'
  grid:
  VREF=5.45  CREF=0.1666667  VCIRC=-1.5
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  GRD='p1i-icb.out'  PS='p1i-icb.ps' SHOWPAGE=0
  GRD='p1i-vep.out'  PS='p1i-vep.ps' SHOWPAGE=0
  grdps:

  PICDAT='picdat.tmp' PICTURE='p1i-ray.ps'
  echo.pl: "'PROJECTION'                              " "> picdat.tmp"
  echo.pl: "16.35  27.2    0. -0.5  0. 0. -0.5 0. /   " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                              " ">>picdat.tmp"
  echo.pl: "'MTYPE' 2    / Marker type                " ">>picdat.tmp"
  echo.pl: "'MSZSF' 0.50 / Marker size scale factor.  " ">>picdat.tmp"
  echo.pl: "'ITC'   0    / No texts drawn.            " ">>picdat.tmp"
  echo.pl: "'LWIDTH'   0.01    /                      " ">>picdat.tmp"
  echo.pl: "/                                         " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1-bnd.dat'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'LINES'                                   " ">>picdat.tmp"
  echo.pl: "'p1i-ray.out'                             " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"
  echo.pl: "'ATTRIBUTES'                              " ">>picdat.tmp"
  echo.pl: "'MTYPE' 2    / Marker type                " ">>picdat.tmp"
  echo.pl: "'MSZSF' 0.20 / Marker size scale factor.  " ">>picdat.tmp"
  echo.pl: "'LWIDTH'   1.50    /                      " ">>picdat.tmp"
  echo.pl: "/                                         " ">>picdat.tmp"
  echo.pl: "'POINTS'                                  " ">>picdat.tmp"
  echo.pl: "'p1-rec.dat'                              " ">>picdat.tmp"
  echo.pl: "                                          " ">>picdat.tmp"

  copy.pl:                 "p1i-r01.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v01.ps"
  append.pl:               "p1i-v01.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r02.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v02.ps"
  append.pl:               "p1i-v02.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r03.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v03.ps"
  append.pl:               "p1i-v03.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r04.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v04.ps"
  append.pl:               "p1i-v04.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r05.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v05.ps"
  append.pl:               "p1i-v05.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r06.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v06.ps"
  append.pl:               "p1i-v06.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r07.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v07.ps"
  append.pl:               "p1i-v07.ps" "p1i-ray.ps"
#
  copy.pl:                 "p1i-r08.out" "p1i-ray.out"
  pictures:
  copy.pl:    "p1i-vep.ps" "p1i-v08.ps"
  append.pl:               "p1i-v08.ps" "p1i-ray.ps"

# Calculating seismograms
# ~~~~~~~~~~~~~~~~~~~~~~~
  echo.pl: " *  's01.out' 's01i.out' /  " ">  crtout.tmp"
  green:
  greenss:
  ss:
  sp: