# Calculating the amplitudes of Gaussian packets
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~   

# Convert decomposed wavefield (725 samples per trace are significant)
  N1=750 N2=96  N1NEW=725
  GRD='ss.out'  NDIG=0
  binasc:

# 10 zero traces are added to begin and end of decomposed wavefield
  echo.pl:   "7250*0"   ">time.out"
  copy.pl:   "time.out" "ssd.out"
  append.pl: "ssd.out"  "ss.out"
  append.pl: "ssd.out"  "time.out"
 
# Input files
  SSD='ssd.out'        GPSTEP='mgpstep.out'
  FGBR22='fgbr22.out'  FGBY22='fgby22.out'
# Discretization steps of the decomposed wavefield
  NX=116 DX=25. OX=0.  NT=725 DT=0.004 OT=0.
# Parameters
  PARKSI=0.8  PARNY=0.5  KAPPA=1.253314137
  GWM=3.0     RELAMP=0.001
# Output files
  SSC='ssc.out'
# Calculating the amplitudes of Gaussian packets
  gpanal: