# History file 'm2d-ss.h' generates the Born seismogram from the RF
# format file computed by history file 'm2d-born.h' and visualizes it
# together with the ray-theory seismogram computed by history file
# 'm2d-ss2.h'.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/born/" "m2d-rec.dat"
# output of 'm2d-born.h'
  chk.pl: ""           "rf-born.out"
# output of 'm2d-ss2.h'
  chk.pl: ""           "m2d-rss.gse"

# Data for ss.for
# 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

# The input RF file
  RF='rf-born.out'

# Data for sp.for

# Data to control seismogram plotting (program SP) in model M1
  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...Born approximation'
  SPTEXT4='RED...Ray theory'

## Data to control seismogram plotting (program SP) in model M2
#  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...Born approximation'
#  SPTEXT4='RED...Ray theory'
#
## Data to control seismogram plotting (program SP) in model M3
#  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...Born approximation'
#  SPTEXT4='RED...Ray theory'

# Other input data
  REC='m2d-rec.dat'

# Frequency-domain seismograms
  SS='m2d-ss.gse'    # born seismogram
  SS1='m2d-rss.gse'  # ray-theory seismogram
  SP1='m2d-ss1.ps'  SP2='m2d-ss2.ps'  SP3='m2d-ss3.ps'
  KOLOR=1  KOLOR1=2
  SPTDEC=1
  ss:         # synthetic seismograms
  sp:         # seismogram plotting