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

# Convert decomposed wavefield
  N1=750 N2=96 N3=1 N1NEW=725 GRD='ss.out' NDIG=0
  binasc:

# Discretization steps of the decomposed wavefield
  NX=96 DX=25. OX=0.  NT=725  DT=0.004 OT=0.
# Parameters
 #PKSI=0.8  PNY=0.5  KAPPA=1.253314137
  GWM=3.0     RELAMP=0.001
#
  N1=725 N2=96  N1NEW=1 N2NEW=1 
  GRD='ss.out'  GRDNEW='grdnew.out'  GNORM=2
  grdnorm:

# Calculating the amplitudes of Gaussian packets
  SS='ss.out'
  gpanal:
  RELERR=0.05
  gpsynt: