# Calculating the amplitudes of Gaussian packets
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  N1=750      N2=96       N3=1

  GRD='timess.out'  NDIG=0
  binasc:

  CAL='mul.cal' GRD1='timess.out' GRD2='gpa-sm.dat' GRD3='timess.out'
  grdcal:

  echo.pl:   "7500*0"    ">time.out"
  copy.pl:   "time.out"  "times.out"
  append.pl: "times.out" "timess.out"
  append.pl: "times.out" "time.out"

#  grdte3:

  N1=750      N2=116       N3=1
  CAL='mul.cal' GRD1='times.out' GRD2='cos-win.out' GRD3='times.out'
#  grdcal:

# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Calculating the "initial parameters of GPs"
  NXTP=14688

  N1=36
  N2=34  D2=0.404040E-04  O2=-0.666667E-3
  N3=12  D3=0.272364

  AR0=0.25E-6
  GRDTER='fgby22.out'

  grdte:

  AR0=-0.25E-6
#  AR0=0.0
  GRDTER='fgbr22.out'

  grdte:

  AR0=1.0
#  AR3=0.138071
  AR3=0.276142
#  TI2=1200.0
  TI2=1575.0
  GRDTER='cos2.out'

#  grdte:

  CAL='mul.cal' GRD1='cos2.out' GRD2='cos2.out' GRD3='cos2.out'

#  grdcal:

  CAL='div.cal' GRD1='fgby22.out' GRD2='cos2.out' GRD3='fgby22.out'

#  grdcal:
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  N1NEW=1  N2NEW=1  N3NEW=1
  GNORM=-999
  GRD='fgbr22.out' GRDNEW='fgbrmin.out'

  grdnorm:

  GRD='fgby22.out' GRDNEW='fgbymin.out'

  grdnorm:

  GNORM=999
  GRD='fgby22.out' GRDNEW='fgbymax.out'

  grdnorm:

  GRD='fgbr22.out' GRDNEW='fgbrmax.out'

  grdnorm:
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  N1=750      N2=116       N3=1
  CREF=0.166667  VREF=0  VCIRC=4003.7
  GRD='times.out'  PS='times.ps'

#  grdps:
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Calculating the amplitudes of GPs
  NX=116      NT=750
  DX=25.0     DT=0.004
#  OX=425.0    OT=0.0

#  AN=0.5E-6  AK=15.0
#  AN=0.25E-6  AK=11.781

#  VSURF=2000.0
  VSURF=1500.0
  GWM=3.0

  FMIN=5.0
  FMAX=55.0
  FREF=35.0
#  PARK=0.5
#  PARN=0.5

#  NXTP=42042

  TIMESEC='times.out'
#  GPOUTR='ampr000.out'
#  GPOUTI='ampi000.out'
  SYNTSEC='synts.out'

  gpanal:
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  GRD='synts.out'  PS='synts.ps'

#  grdps:

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CAL='sub.cal' GRD1='synts.out' GRD2='times.out' GRD3='huh.out'

#  grdcal:

#  VCIRC=100.0
#  VCIRC=
  GRD='huh.out' PS='huh.ps'

#  grdps: