# Calculating velocity differences in models P1 and P1I
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
# ~~~~~~~~~~~~~~~~~~~~
  chk.pl: "forms/"   "sub.cal"
  chk.pl: "data/p1/" "p1-vp.dat"
# output of 'p1-inv.h'
  chk.pl: ""         "p1-vs.out"
  chk.pl: ""         "p1-mod.dat"
  chk.pl: ""         "p1-mods.dat"
# output of 'p1i-inv.h'
  chk.pl: ""         "p1i-mod.dat"

# Plot of the gridded velocity data
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  VREF=5.45  CREF=0.1666667  VCIRC=-1.5
  HSIZE=24 YSIGN=-1 ROTATE=-90 VOFFSET=13
  O1=0.        D1=0.1     N1=474
  O3=0.        D3=0.1     N3=61
  GRD='p1-vp.dat'  PS='p1-vp.ps'
  grdps:
  VREF=3.14  CREF=0.096225   VCIRC=-0.866
  GRD='p1-vs.out'  PS='p1-vs.ps'
  grdps:

# Plot of the gridded velocities in the smooth model
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='p1-mod.dat'
  VREF=5.45  CREF=0.1666667  VCIRC=-1.5
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  MPAR=1   VEL='p1-vep.out'
  grid:
  GRD='p1-vep.out'  PS='p1-vep.ps'
  grdps:
  MODEL='p1-mods.dat'
  VREF=3.14  CREF=0.096225   VCIRC=-0.866
  MPAR=2   VEL='p1-ves.out'
  grid:
  GRD='p1-ves.out'  PS='p1-ves.ps'
  grdps:

# Plot of the gridded velocities in the model with interfaces
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='p1i-mod.dat'
  VREF=5.45  CREF=0.1666667  VCIRC=-1.5
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  MPAR=1   VEL='p1i-vep.out'
  grid:
  GRD='p1i-vep.out'  PS='p1i-vep.ps'
  grdps:
  VREF=3.14  CREF=0.096225   VCIRC=-0.866
  MPAR=2   VEL='p1i-ves.out'
  grid:
  GRD='p1i-ves.out'  PS='p1i-ves.ps'
  grdps:
  O3=          D3=        N3=


# Subtracting the velocities
# ~~~~~~~~~~~~~~~~~~~~~~~~~~

# Smooth model - model with intefaces
  CAL='sub.cal'
  N1=474  N2=61  N3=
  GRD1='p1-vep.out'  GRD2='p1i-vep.out' GRD3='p1-difvp.out'
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.4
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1-difvp.out'  PS='p1-difvp.ps'
  grdps:
  N1=474  N2=61  N3=
  GRD1='p1-ves.out'  GRD2='p1i-ves.out' GRD3='p1-difvs.out'
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.23 #=0.4/sqrt(3)
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1-difvs.out'  PS='p1-difvs.ps'
  grdps:

# Smooth model - data
  GRD1='p1-vep.out'  GRD2='p1-vp.dat' GRD3='p1-dvp.out'
  N1=474  N2=61  N3=
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.4
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1-dvp.out'  PS='p1-dvp.ps'
  grdps:

  GRD1='p1-ves.out'  GRD2='p1-vs.out' GRD3='p1-dvs.out'
  N1=474  N2=61  N3=
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.23
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1-dvs.out'  PS='p1-dvs.ps'
  grdps:

# Model with interfaces - data
  GRD1='p1i-vep.out'  GRD2='p1-vp.dat' GRD3='p1i-dvp.out'
  N1=474  N2=61  N3=
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.4
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1i-dvp.out'  PS='p1i-dvp.ps'
  grdps:
  GRD1='p1i-ves.out'  GRD2='p1-vs.out' GRD3='p1i-dvs.out'
  N1=474  N2=61  N3=
  grdcal:
  GRD1=              GRD2=              GRD3=
  VREF=0      CREF=0.1666667  VCIRC=0.23
  O1=0.        D1=0.1     N1=474
  O2=          D2=        N2=
  O3=0.        D3=0.1     N3=61
  GRD='p1i-dvs.out'  PS='p1i-dvs.ps'
  grdps: