# History file 'p1-grid.h' computes the elastic parameters at the
# gridpoints in both unperturbed (background) and perturbed model and
# then it computes the perturbations of these elastic parameters.
# The history file also generates several pictures.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/p1i/" "pert2.cal"
  chk.pl: "forms/"    "sub.cal"
  chk.pl: "data/p1i/" "p1-mod1c.dat"
  chk.pl: "data/p1i/" "p1-mod2.dat"

# Computation of the perturbations of the elastic parameters at
# the gridpoints in the smaller grid for the Born approximation
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Grid dimensions
  N1=2500     N2=1   N3=650
  D1=0.005    D2=0.  D3=0.005
  O1=21.0025  O2=0.  O3=2.0025

# Perturbed model
  MODEL='p1-mod1c.dat'
  LFREE=1
  #modchk:

# Discretization of the elastic parameters at the gridpoints
# in the perturbed model.
# P-wave velocity
  MPAR=1
  VEL='p1-vp1.out'
  grid:
# S-wave velocity
  ICB=
  MPAR=2
  VEL='p1-vs1.out'
  grid:
# Density
  MPAR=3
  VEL='p1-rho1.out'
  grid:

# Unperturbed model
  MODEL='p1-mod2.dat'
  LFREE=1
  #modchk:

# Discretization of the elastic parameters at the gridpoints
# in the unperturbed model
# P-wave velocity
  MPAR=1
  VEL='p1-vp2.out'
  grid:
# S-wave velocity
  MPAR=2
  VEL='p1-vs2.out'
  grid:
# Density
  MPAR=3
  VEL='p1-rho2.out'
  grid:

# Computation of the perturbations of the elastic parameters
# at the gridpoints
  CAL='pert2.cal'
  GRD1='p1-vp1.out'
  GRD2='p1-vp2.out'
  GRD3='p1-vpper.out'
  GRD4='p1-vs1.out'
  GRD5='p1-vs2.out'
  GRD6='p1-vsper.out'
  GRD7='p1-rho1.out'
  GRD8='p1-rho2.out'
  GRD9='p1-rhper.out'
  grdcal:

# Generating pictures
# ~~~~~~~~~~~~~~~~~~~

# Grid dimensions
  N1=4730   N2=1   N3=600
  D1=0.01   D2=0.  D3=0.01
  O1=0.005  O2=0.  O3=0.005

# Perturbed model
  MODEL='p1-mod1c.dat'
  LFREE=1
  #modchk:

# Discretization of the elastic parameters at the gridpoints
# in the perturbed model.
# P-wave velocity
  ICB='p1-icb.out'
  MPAR=1
  VEL='p1-vp1.out'
  grid:

# Unperturbed model
  MODEL='p1-mod2.dat'
  LFREE=1
  #modchk:

# Discretization of the elastic parameters at the gridpoints
# in the unperturbed model.
# P-wave velocity
  ICB=' ' #reseting to default
  MPAR=1
  VEL='p1-vp2.out'
  grid:

# Plotting P-wave velocity perturbation.
  GRD1='p1-vp1.out'
  GRD2='p1-vp2.out'
  GRD3='p1-vpdif.out'
  CAL='sub.cal'
  grdcal:

  SHOWPAGE=1
  GRD='p1-vpdif.out'
  PS='p1-vpdif.ps'
  VREF=0
  CREF=0.1666
  VCIRC=-0.02
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  grdps:

# Plotting indeces of the blocks.
  SHOWPAGE=0
  GRD='p1-icb.out'
  PS='p1-icb.ps'
  VREF=8
  CREF=0.6666
  VCIRC=-16
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  VMIN=7.9
  VMAX=8.1
  R=1
  G=1
  B=0
  grdps:

# Plotting P-wave velocity in the unperturbed model.
  SHOWPAGE=0
  GRD='p1-vp2.out'
  PS='p1-vep.ps'
  VREF=5.45
  CREF=0.1666667
  VCIRC=-1.5
  HSIZE=23.65 YSIGN=-1 ROTATE=-90 VOFFSET=13.35  HOFFSET=
  VMIN=-999999  #reseting VMIN, VMAX, R, G, B to default
  VMAX=999999000000.
  R=0.80
  G=0.80
  B=0.80
  grdps: