# History file 'sal-inv1.h' to fit the SEG/EAGE 3-D Salt Model by a smooth
# 3-D model.
# =====================================================================

# Input files required
  chk.pl: "data/sal/"  "sal-m1.dat"
  chk.pl: "model/"     "sob22.dat"
  chk.pl: "forms/"     "inv.cal"
  chk.pl: "forms/"     "sqrt.cal"
  chk.pl: "forms/"     "addsob.cal"
# Output of 'sal-vel.h'
  chk.pl: ""           "s-sa825.out"

# Initial model
# ~~~~~~~~~~~~~
  MODEL='sal-m1.dat'  MODIN='sal-m1.dat'

# Form of the files with matrices
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# FORMM='formatted'
  FORMM='unformatted'
# MATFORM='formatted'
  MATFORM='unformatted'

# Inversion of velocities
# ~~~~~~~~~~~~~~~~~~~~~~~
  ICLASS=2
  NEGPAR=1

# Data grid
  N1=28   N2=28   N3=9
  D1=500  D2=500  D3=500
  O1=0    O2=0    O3=0

  N1=23   N2=23   N3=11
  D1=600  D2=600  D3=400
  O1=300  O2=300

  N1=14   N2=14   N3=21
  D1=1000 D2=1000 D3=200
  O1=250  O2=250

  N1=18   N2=18   N3=17
  D1=800  D2=800  D3=250
  O1=-50  O2=-50  O3=0


# Calculating matrices for inversion
  M1='m1.out'  M2='m2.out'
  SOBOLEV='sob22.dat'  MODSOB='modsob.out'  SOBW01=1
  invsoft:
                       MODSOB=              SOBW01=
  GM1='gm1.out'  GM2='gm2.out'  GM3=' '  DM1='dm1.out'
  GRD='s-sa500.out'  INDFUN=1 MPAR=1 POWERM=-1 ERRMUL=84.   # 7056 points
  GRD='s-sa640.out'  INDFUN=1 MPAR=1 POWERM=-1 ERRMUL=76.28 # 5819 points
  GRD='s-sa000.out'  INDFUN=1 MPAR=1 POWERM=-1 ERRMUL=64.16 # 4116 points
  GRD='s-sa825.out'  INDFUN=1 MPAR=1 POWERM=-1 ERRMUL=74.22 # 5508 points
  invpts:
  GRD=                INDFUN=  MPAR=  POWERM=   ERRMUL=

# Matrix operations
  MATIN1='dm1.out'                       MATOUT='dm2.out'
  MATFUN='inv'
  matfun:
  MATIN1='dm2.out'   MATIN2='gm1.out'    MATOUT='dm2gm1.out'
  SYMMETRY=      MATT1=  MATT2=1
  matmul:
  MATIN1='gm1.out'   MATIN2='dm2gm1.out' MATOUT='sm1.out'
  SYMMETRY='sym' MATT1=  MATT2=
  matmul:
  MATIN1='sm1.out'   MATIN2='modsob.out' MATOUT='sm2.out'
  COEF1=  COEF2=2500000000
  matlin:
  MATIN1='sm2.out'                       MATOUT='sm3.out'
  matinv:
  MATIN1='dm2.out'   MATIN2='gm2.out'    MATOUT='gm3.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
  MATIN1='gm1.out'   MATIN2='gm3.out'    MATOUT='gm4.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
  MATIN1='sm3.out'   MATIN2='gm4.out'    MATOUT='gm5.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:

# Updating the model
  M1='m1.out'  MODNEW='gm5.out'
  MODOUT='sal-m1.out'
  modmod:

#  Sobolev norm of the slowness in the updated model
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MATIN1='modsob.out'    MATIN2='gm5.out'    MATOUT='gm6.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
  MATIN1='gm5.out'       MATIN2='gm6.out'    MATOUT='gm7.out'
  SYMMETRY='diag' MATT1=1 MATT2=0
  matmul:
  MATIN1='gm7.out'                       MATOUT='sal-sns.out'
  SPARSE=-1 MATFUN='sqrt'                MATFORM='formatted'
  matfun:

# List of output files
# ~~~~~~~~~~~~~~~~~~~~
# 'sal-m1.out' ... Inverted smooth SEG/EAGE 3-D Salt Model.
# 'sal-sns.ou@'... Sobolev norm of the slowness in the model.