# Calculation of the discretization steps for the decomposition 
# of the wavefield into Gaussian packets
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/mar"  "mgb-grd.h"
  chk.pl: "forms/"    "append.pl"
  chk.pl: "forms/"    "del.pl"
 
# Deleting old file with maximum and minimum values of R0, Y0
  del.pl: "ryminmax.out"    
  del.pl: "mgpstep.out"    
 
# Determination of maximum and minimum values of the optimized
# and smoothed initial parameters of Gaussian beams R0 (GBR22), 
# Y0 (GBY22)
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Name of the file with minimum and maximum values of R0, Y0
  RYMM="ryminmax.out"
# Grid parameters of initial Hamiltonian hypersurface (XR,TR,P)
  N1=26  N2=48  N3=36
  D1=1.  D2=1.  D3=1. 
  O1=1.  O2=1.  O3=1.
# R0*R0/Y0+Y0
  CAL='gb-ry0m.cal'
  GRD1='r0m.out' GRD2='y0m.out' GRD3='ry0m.out'
  grdcal:
# Grid of output files with minimum a maximum values of R0, Y0
  N1NEW=1  N2NEW=1  N3NEW=1
# Minimum of R0
  GNORM=-999 GRD='r0m.out' GRDNEW='ryminmax.tmp'
  grdnorm:
  append.pl: "ryminmax.out" "ryminmax.tmp"
# Maximum of R0
  GNORM=999  GRD='r0m.out' GRDNEW='ryminmax.tmp'
  grdnorm:
  append.pl: "ryminmax.out" "ryminmax.tmp"
# Minimum of Y0
  GNORM=-999 GRD='ry0m.out' GRDNEW='ryminmax.tmp'
  grdnorm:
  append.pl: "ryminmax.out" "ryminmax.tmp"
# Maximum of Y0
  GNORM=999  GRD='ry0m.out' GRDNEW='ryminmax.tmp'
  grdnorm:
  append.pl: "ryminmax.out" "ryminmax.tmp"
 
# Discretization steps of the decomposed wavefield
  NX=96 DX=25. OX=0.  NT=725 DT=0.004 OT=0.
# Minimum velocity in the depth of the receivers
  VSURF=1850.
# Frequencies of the Gaussian packets
  FMIN=5.0 FMAX=55.0 FREF=35.0
# Parameters
  PKSI=0.8  PNY=0.5  KAPPA=1.253314137  PKSIT=0.7  IVER=1
  GPSTEP='mgpstep.out'  GBGRD='gb-grd.tmp' GPSTEPI='mgbstep.ini'
  FGBR22='r0m.out'  FGBY22='y0m.out'
  DSRC=25. OXRINI=425. 
  gpstep:

# Calculating new gridded values of optimized GB parameter R0 (GBR22), 
# Y0 (GBY22)
  append.pl: "gb-grd.tmp" "mgb-grd.h"
  go.pl: "gb-grd.tmp" "gp-grd.out"