# History file 'mgp-an.h' to calculate
# the amplitudes of Gaussian packets
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# see  Zacek, K. (2006):  Decomposition of the wave field into optimized
#      Gaussian packets. Stud. Geophys. Geod., 50, 367-380.
#
# This example shows calculation for shot 77.
# Program gpsynt in mgb-an1.h needs MRAM=12000000 in file "ram.inc".
#
# Input files required
  chk.pl: "forms/"    "sub.cal"
  chk.pl: "forms/"    "echo.pl"
  chk.pl: "forms/"    "copy.pl"
  chk.pl: "forms/"    "append.pl"
  chk.pl: "data/mar/" "mgp-an1.h"
# Output of 'mgb-opt.h'
  chk.pl: " "         "r0m077.out"
  chk.pl: " "         "y0m077.out"
  chk.pl: " "         "step077.out"
# Input files 'seis-001.h@' to 'seis-240.h@' can be found at
# http://sw3d.mff.cuni.cz/software/marmousi/marmousi.htm
  chk.pl: ""          "seis-077.h@"

# Parameters for plotting seismograms
  N1=725  N2=96   N3=1  CREF=0.166667  VREF=0.
  ROTATE=-90  HSIZE=22.86  VSIZE=15.24

  echo.pl: "IN='seis-077.h@'" ">gpa.tmp"
  echo.pl: "GPAMPR='ampr077.out' GPAMPI='ampi077.out'" ">>gpa.tmp"
  echo.pl: "FGBR22='r0m077.out' FGBY22='y0m077.out' GPSTEP='step077.out'" ">>gpa.tmp"
  echo.pl: "GPA0R='a0r077.out' GPA0I='a0i077.out' SSC='ssc.out'" ">>gpa.tmp"
  append.pl: "gpa.tmp" "mgp-an1.h"
  go.pl: "gpa.tmp" "gpa.out"
  GRD='ss.out'  PS='ss077.ps'  grdps:
  GRD='ssc.out'  PS='ssc077.ps'  grdps:
  CAL='sub.cal' GRD1='ssc.out' GRD2='ss.out' GRD3='dif.out'  grdcal:
  GRD='dif.out' PS='dif077.ps'    grdps:

# Important output files
# ~~~~~~~~~~~~~~~~~~~~~~
# 'amprxxx.out'... Amplitudes of Gaussian packets - real part
# 'ampixxx.out'... Amplitudes of Gaussian packets - imaginary part
# 'a0rxxx.out'...  Amplitudes after replacing of the smallest amplitudes
#                  of Gaussian packets by zeroes - real part
# 'a0ixxx.out'...  Amplitudes after replacing of the smallest amplitudes
#                  of Gaussian packets by zeroes - imaginary part
# 'ssxxx.out'... Seismograms to be decomposed
# 'ssxxx.ps'...  Plots of seismograms to be decomposed
# 'ssc.out'... Seismograms after composition
# 'sscxxx.ps'...  Plots of seismograms after composition
# 'dif.out'...   Differences between seismograms
# 'difxxx.ps'...    Plots of differences between seismograms