# History file 'm2d-ss2.h' computes the ray-theory seismograms in
# models M1, M2 and M3 for the configuration of the source and receiver
# suitable for 2D computations.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/born/" "m1-mod.dat"
#  chk.pl: "data/born/" "m2-mod.dat"
#  chk.pl: "data/born/" "m3-mod.dat"
  chk.pl: "data/born/" "m2d-src.dat"
  chk.pl: "data/born/" "m2d-rec.dat"
  chk.pl: "data/born/" "m2d-crt2.dat"
  chk.pl: "data/born/" "m2d-rpa2.dat"
  chk.pl: "crt/"       "writ.dat"
  chk.pl: "crt/"       "source.dat"
  chk.pl: "forms/"     "hsv.dat"
  chk.pl: "forms/"     "echo.pl"

# Model

# Model M1-2D
  MODEL='m1-mod.dat'
  LFREE=1
  # modchk:
## Model M2-2D
#  MODEL='m2-mod.dat'
#  LFREE=1
#  # modchk:
## Model M3-2D
#  MODEL='m3-mod.dat'
#  LFREE=1
#  # modchk:

# Source and receivers
  SRC='m2d-src.dat'  REC='m2d-rec.dat'
  SOURCE='source.dat'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='m2d-crt2.dat'
  CODE='m2d-crt2.dat'
  RPAR='m2d-rpa2.dat'
  WRIT='writ.dat'
  CRTLOG='m2d-crt2.out'

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

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

# Data describing the filtration of the source time function
  FMIN=1.  FLOW=2.  FHIGH=19.  FMAX=20. # (cosine filter)

# Time step and time interval for the Fast Fourier Transform
  DT=0.0025  NFFT=8192

# Data to control optional plotting of program SS
  SIGPLOT='m2d-sig.ps'  MPTS=100

# Data to control seismogram plotting (program SP) in M1-2D
  KODESP=0    SPCHRH=0.40
  SPTMIN=3.2  SPTMAX=3.6  SPTLEN=25.0  SPTDIV=4    SPTSUB=10
  SPXMIN=0.0  SPXMAX=2.0  SPXLEN=17.0  SPXDIV=1  SPXSUB=1
  NORMSP=1    SPAMP=-9.5e8
  SPTEXT3='BLACK...Ray theory'
## Data to control seismogram plotting (program SP) in M2-2D
#  KODESP=0    SPCHRH=0.40
#  SPTMIN=2.9  SPTMAX=3.3  SPTLEN=25.0  SPTDIV=4    SPTSUB=10
#  SPXMIN=0.0  SPXMAX=2.0  SPXLEN=17.0  SPXDIV=1  SPXSUB=1
#  NORMSP=1    SPAMP=-4.e8
#  SPTEXT3='BLACK...Ray theory'
## Data to control seismogram plotting (program SP) in M3-2D
#  KODESP=0    SPCHRH=0.40
#  SPTMIN=2.8  SPTMAX=3.6  SPTLEN=25.0  SPTDIV=8    SPTSUB=10
#  SPXMIN=0.0  SPXMAX=2.0  SPXLEN=17.0  SPXDIV=1  SPXSUB=1
#  NORMSP=1    SPAMP=-6.5e8
#  SPTEXT3='BLACK...Ray theory'

# Data to control ray parameters plotting (program RPPLOT)
  ITHOM=0 ISUC=1  HRBAS=0.2 HRTWO=0.3 HRAUX=0.1 IWAVES=1
  echo.pl: "' '   ' '   's01i.out' 't01.out'" ">crtout.tmp"
  CRTOUT='crtout.tmp' RPPLOT='m2d-plot1.ps'

# Programs to calculate the Green function
  crt:        # two-point ray tracing
  #rpplot:     # ray parameters plotting
  CRTOUT=
  green:      # Green function

# Frequency-domain seismograms
  SS='m2d-rss.gse'  SP1='m2d-rss1.ps'  SP2='m2d-rss2.ps'  SP3='m2d-rss3.ps'
  greenss:    # response function
  ss:         # synthetic seismograms
  sp:         # seismogram plotting