# Two point ray tracing in smoothed SEG/EAGE Salt Model
# followed by plotting of the computed rays. Source and
# receivers positions are similar to the Phase A 
# SEG Salt Model project.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
 #chk.pl: "data/sal/" "sal-dcr1.dat"
 #chk.pl: "data/sal/" "sal-rp1.dat"
 #chk.pl: "data/sal/" "sal-cod1.dat"
 #chk.pl: "data/sal/" "sal-mod2.dat"
 #chk.pl: "data/sal/" "sal-src1.dat"
 #chk.pl: "data/sal/" "sal-rec1.dat"
 #chk.pl: "data/sal/" "sal-rec1.plt"
 #chk.pl: "crt/"      "writ.dat"
 #chk.pl: "crt/"      "sourcex.dat"
 #chk.pl: "forms/"    "echo.pl"
 #chk.pl: "wrl/"      "sal-srf.wrl"  # VRML
 #chk.pl: "wrl/"      "sal-srf.mx"   # GOCAD
# For optional computation of individual seismograms
# of elementary waves
 #chk.pl: "data/sal/" "sal-wri1.dat"

# Model, source and receivers
  MODEL='sal-mod2.dat' SRC='sal-src1.dat' REC='sal-rec1.dat'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='sal-dcr1.dat'
  RPAR='sal-rp1.dat'
  CODE='sal-cod1.dat' 
  WRIT='writ.dat' 
  CRTLOG='sal-log1.out'

# Kind of ray parameters (for details see 'crtin.for')
  INIPAR=3

# File containing the complex-valued seismic moment
  SOURCE=sourcex.dat # Explosive isotropic source

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

# Data describing the filtration of the source time function
  FMIN=1.  FLOW=5.   FHIGH=56.   FMAX=60.  # (cosine filter)

# Time step and time interval for the Fast Fourier Transform
  DT=0.001  NFFT=4096

# Data to control seismogram plotting (for details see 'sp.for')
  KODESP=2     SPCHRH=0.40
  SPTMIN=1.0   SPTMAX=4.0     SPTLEN=25.0  SPTDIV=6    SPTSUB=5
  SPXMIN=4160. SPXMAX=11200.  SPXLEN=15.3  SPXDIV=-22  SPXSUB=1
  SPYMIN=7640. SPYMAX=7640.
  NORMSP=0     SPAMP=1.

# Ray tracing of refracted and reflected waves
  crt:

# Frequency-domain seismograms in selected receivers
  green:      # Green function
  greenss:    # response function
  ss:         # synthetic seismograms 
  REC='sal-rec1.plt' 
  SP1='sal-sf1.ps' SP2='sal-sf2.ps' SP3='sal-sf3.ps'
  sp:         # seismogram plotting

# Optional computation of individual seismograms
# of elementary waves
# Ray tracing of refracted and reflected waves
 #REC='sal-rec1.dat' WRIT='sal-wri1.dat' 
 #crt:

# Frequency-domain seismograms in selected receivers
# 1st elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r01.out' 's011.out ' 's011i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s1.gse'
 #ss:         # synthetic seismograms 

# 2nd elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r02.out' 's012.out ' 's012i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s2.gse'
 #ss:         # synthetic seismograms 

# 3rd elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r03.out' 's013.out ' 's013i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s3.gse'
 #ss:         # synthetic seismograms 

# 4th elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r04.out' 's014.out ' 's014i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s4.gse'
 #ss:         # synthetic seismograms 

# 5th elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r05.out' 's015.out ' 's015i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s5.gse'
 #ss:         # synthetic seismograms 

# 6th elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r06.out' 's016.out ' 's016i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s6.gse'
 #ss:         # synthetic seismograms 

# 7th elementary wave
 #CRTOUT='crtout.tmp' SS=' '
 #echo.pl:"'r07.out' 's017.out ' 's017i.out'" ">crtout.tmp"
 #green:      # Green function
 #greenss:    # response function
 #SS='sal-s7.gse'
 #ss:         # synthetic seismograms 

# There are no rays at the selected receivers for 
# 8th elementary wave

# All elementary waves seismogram plotting
 #REC='sal-rec1.plt' 
 #SP1='sal-sf11.ps'  SP2='sal-sf12.ps'  SP3='sal-sf13.ps'
 #SS='sal-s1.gse'  KOLOR=2
 #SS1='sal-s2.gse' KOLOR1=3
 #SS2='sal-s3.gse' KOLOR2=3
 #SS3='sal-s4.gse' KOLOR3=4
 #SS4='sal-s5.gse' KOLOR4=4
 #SS5='sal-s6.gse' KOLOR5=1
 #SS6='sal-s7.gse' KOLOR6=1
 #sp:         # seismogram plotting

# Rewriting the two-point rays into a formatted file
  REC='sal-rec1.dat' CRTOUT='crtout.tmp'
  echo.pl:"'r01.out'   ' '   'r01i.out' ' '" ">crtout.tmp"
  RAYS='sal-r1.out'  crtray:
  echo.pl:"'r02.out'   ' '   'r02i.out' ' '" ">crtout.tmp"
  RAYS='sal-r2.out'  crtray:
  echo.pl:"'r03.out'   ' '   'r03i.out' ' '" ">crtout.tmp"
  RAYS='sal-r3.out'  crtray:
  echo.pl:"'r04.out'   ' '   'r04i.out' ' '" ">crtout.tmp"
  RAYS='sal-r4.out'  crtray:
  echo.pl:"'r05.out'   ' '   'r05i.out' ' '" ">crtout.tmp"
  RAYS='sal-r5.out'  crtray:
  echo.pl:"'r06.out'   ' '   'r06i.out' ' '" ">crtout.tmp"
  RAYS='sal-r6.out'  crtray:
  echo.pl:"'r07.out'   ' '   'r07i.out' ' '" ">crtout.tmp"
  RAYS='sal-r7.out'  crtray:
  echo.pl:"'r08.out'   ' '   'r08i.out' ' '" ">crtout.tmp"
  RAYS='sal-r8.out'  crtray:

# Adding two-point rays into the VRML or GOCAD file with surfaces
  VRML='VRML2'  WRL='sal-srf.wrl'  WRLOUT='sal-ray1.wrl'  # VRML
 #VRML='GOCAD'  WRL='sal-srf.mx'   WRLOUT='sal-ray1.mx'   # GOCAD
# Setting blue colour for rays
  KOLLIN=0  R=0.00  G=0.00  B=1.00
  NAME='Rays'  PROPERTIES=' '   # just for GOCAD
  LIN='sal-r1.out'
  NAME='Rays1-refracted'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  WRL='sal-ray1.wrl'  # VRML
 #WRL='sal-ray1.mx'   # GOCAD
  LIN='sal-r2.out'
  NAME='Rays2-top-salt-1'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r3.out'
  NAME='Rays3-top-salt-2'     PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r4.out'
  NAME='Rays4-bottom-salt-1'  PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r5.out'
  NAME='Rays5-bottom-salt-2'  PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r6.out'
  NAME='Rays6-bottom-1'       PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r7.out'
  NAME='Rays7-bottom-2'       PROPERTIES=' '  # just for GOCAD
  linwrl:
  LIN='sal-r8.out'
  NAME='Rays8-bottom-3'       PROPERTIES=' '  # just for GOCAD
  linwrl:

# Adding source and receivers into the VRML or GOCAD file
  PTS='sal-src1.dat'
  KOLPTS=0  R=1.00  G=0.00  B=1.00
  NAME='Source'      PROPERTIES=' '  # just for GOCAD
  ptswrl:
  PTS='sal-rec1.dat'
  KOLPTS=0  R=1.00  G=1.00  B=1.00
  NAME='Receivers'   PROPERTIES=' '  # just for GOCAD
  ptswrl:
  PTS='sal-rec1.plt'
  KOLPTS=0  R=1.00  G=0.00  B=0.00 SHIFT3=-0.1
  NAME='Receivers-1' PROPERTIES=' '  # just for GOCAD
  ptswrl: