# History file 'sal-inv3.h' to fit the SEG/EAGE 3-D Salt Model by a model
# with interfaces. Input data for the salt interface are taken from the
# GOCAD surfaces, data from file 'fault_a.ts' are also used.
# =====================================================================

# Input files required
 #chk.pl: "data/sal/" "sal-m3.dat"
 #chk.pl: "model/"     "sob22.dat"
 #chk.pl: "forms/"     "inv.cal"
 #chk.pl: "forms/"     "addsob.cal"
# Output of 'sal-ts.h'
  chk.pl: ""           "top.pts"
  chk.pl: ""           "hrz_1.pts"
  chk.pl: ""           "gp2.pts"
  chk.pl: ""           "bottom.pts"
  chk.pl: ""           "fault_a.pts"
# Output of 'sal-salt.h'
  chk.pl: ""           "salt-to1.pts"
  chk.pl: ""           "salt-to2.pts"
  chk.pl: ""           "salt-bo1.pts"

# Initial model and its consistency check
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='sal-m3.dat'  MODIN='sal-m3.dat'
  N1=20  N2=20  N3=10  LFREE=1  MODLOG='modchk.out'
 #modchk:

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

# Inversion of interfaces
# ~~~~~~~~~~~~~~~~~~~~~~~
  ICLASS=1            # 646 model parameters for surfaces

# Calculating matrices for inversion
  M1='m1.out'  M2='m2.out'
  SOBOLEV='sob22.dat'  MODSOB='modsob.out'  SOBW00=1
  invsoft:
                       MODSOB=              SOBW00=
  GM1='gm1.out'  GM2='gm2.out'  GM3=' '  DM1='dm1.out'
  MPAR=0
  PTS='top.pts'      INDFUN=1  ERRMUL=2.000  invpts:  #     4 points
  M2IN='m2.out'                                       #
  PTS='hrz_1.pts'    INDFUN=2  ERRMUL=24.52  invpts:  #   601 points
  PTS='gp2.pts'      INDFUN=3  ERRMUL=42.19  invpts:  #  1780 points
  PTS='bottom.pts'   INDFUN=4  ERRMUL=2.000  invpts:  #     4 points
  PTS='salt-t1.pts'  INDFUN=5  ERRMUL=28.16  invpts:  #   793 points
  PTS='salt-b1.pts'  INDFUN=6  ERRMUL=24.02  invpts:  #   577 points
  PTS='fault_a.pts'  INDFUN=7  ERRMUL=31.91  invpts:  #  1018 points
  PTS='salt-t2.pts'  INDFUN=8  ERRMUL=14.04  invpts:  #   197 points


  M2IN=' '  PTS=' '  INDFUN=   ERRMUL=                #       total

# Matrix operations
  N1=0  N2=1  N3=1  M1='m2.out'
  CAL='inv.cal' GRD1='dm1.out' GRD2='dm2.out'
  grdcal:
  M1='m1.out'  M2='m2.out'  GM1='gm1.out'  DM1='dm2.out'  SM1='sm1.out'
  gmdmgmt:
  N1=0  N2=0  N3=1  M1='m1.out' M2='m1.out'
  CAL='addsob.cal'  GRD1='sm1.out'  GRD2='modsob.out'  GRD3='sm2.out'
  SOBMUL=20000 8000
  grdcal:
  SOBMUL=
  M1='m1.out'               SM1='sm2.out'  SM2='sm3.out'
  sminv:
  M1='m2.out'  M2=' '       DM1='dm2.out'  GM1='gm2.out'  GM2='gm3.out'
  dmgm:
  M1='m1.out'  M2='m2.out'  GM1='gm1.out'  GM2='gm3.out'  GM3='gm4.out'
  gmgm:
  M1='m1.out'  M2=' '       SM1='sm3.out'  GM1='gm4.out'  GM2='gm5.out'
  smgm:

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


# Indices of complex blocks relevant to the velocities
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  N1=68   N2=68   N3=20  D1=200  D2=200  D3=200  O1=0    O2=0    O3=40
  N1=46   N2=46   N3=14  D1=300  D2=300  D3=300  O1=0    O2=0    O3=100
  N1=34   N2=34   N3=20  D1=400  D2=400  D3=200  O1=200  O2=200  O3=40
  N1=23   N2=23   N3=21  D1=600  D2=600  D3=200  O1=200  O2=200  O3=100
  N1=17   N2=17   N3=20  D1=800  D2=800  D3=200  O1=400  O2=400  O3=40
  N1=11   N2=11   N3=20  D1=1340 D2=1340 D3=200  O1=40   O2=40   O3=40
  N1=17   N2=17   N3=39  D1=800  D2=800  D3=100  O1=400  O2=400  O3=40

  ISRF=3  VEL='sal-s3.out'  grid:
  ISRF=4  VEL='sal-s4.out'  grid:
  ISRF=
  CAL='sal-icb.cal'
  GRD2='sal-s3.out' GRD3='sal-s4.out' GRD4='sal-icb.out'
  SALTV1=1000  SALTV2=1501  SALTV5=4449  SALTD3=152.4  SALTD4=1
  GRD1='s-va200.out'
  GRD1='s-vc200.out'
  GRD1='s-va200r.out'
  GRD1='s-vc200r.out'
  GRD1='s-va300.out'
  GRD1='s-va420.out'
  GRD1='s-va620.out'
  GRD1='s-va820.out'
  GRD1='s-va134.out'
  GRD1='s-vc810r.out'
  grdcal:


# Inversion of velocities
# ~~~~~~~~~~~~~~~~~~~~~~~
  MODEL='sal-m3.tmp'    MODIN='sal-m3.tmp'
  ICLASS=2
  NEGPAR=1

# 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'
  GRDICB='sal-icb.out' MPAR=1 POWERM=-1
  GRD='s-sa200.out'   ERRMUL=304.11 #92480 points
  GRD='s-sc200.out'   ERRMUL=304.11 #92480 points
  GRD='s-sa200r.out'  ERRMUL=304.11 #92480 points
  GRD='s-sc200r.out'  ERRMUL=304.11 #92480 points
  GRD='s-sa300.out'   ERRMUL=172.12 #29624 points
  GRD='s-sa420.out'   ERRMUL=152.05 #23120 points
  GRD='s-sa620.out'   ERRMUL=105.40 #11109 points
  GRD='s-sa820.out'   ERRMUL=76.026 #5780  points
  GRD='s-sa134.out'   ERRMUL=50.408 #2541  points
  GRD='s-sc810r.out'  ERRMUL=106.16 #11271 points
  invpts:
  GRD= GRDICB= MPAR= POWERM= ERRMUL=

# Matrix operations
  N1=0  N2=1  N3=1  M1='m2.out'
  CAL='inv.cal' GRD1='dm1.out' GRD2='dm2.out'
  grdcal:
  M1='m1.out'  M2='m2.out'  GM1='gm1.out'  DM1='dm2.out'  SM1='sm1.out'
  gmdmgmt:

  N1=0  N2=0  N3=1  M1='m1.out' M2='m1.out'
  CAL='addsob.cal'  GRD1='sm1.out'  GRD2='modsob.out'  GRD3='sm2.out'
  SOBMUL=5000 2500 10000 25000 0 50000 200000 100000 1000 75000 150000
  grdcal:
  SOBMUL=
  M1='m1.out'               SM1='sm2.out'  SM2='sm3.out'
  sminv:
  M1='m2.out'  M2=' '       DM1='dm2.out'  GM1='gm2.out'  GM2='gm3.out'
  dmgm:
  M1='m1.out'  M2='m2.out'  GM1='gm1.out'  GM2='gm3.out'  GM3='gm4.out'
  gmgm:
  M1='m1.out'  M2=' '       SM1='sm3.out'  GM1='gm4.out'  GM2='gm5.out'
  smgm:

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

## Sobolev norm of the slowness in the updated model
##~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  M1='m1.out'  M2=' '  SM1='modsob.out' GM1='gm5.out'  GM2='gm6.out'
  smgm:
  M1=' '  M2='m1.out'  GM1='gm5.out'  GM2='gm6.out'  GM3='gm7.out'
  gmgm:
  N1=0  N2=1  N3=1  M1=' '
  CAL='sqrt.cal' GRD1='gm7.out' GRD2='s-sns.out' FORMMW='formatted'
  grdcal:

# List of output files
# ~~~~~~~~~~~~~~~~~~~~
# 'sal-m3.tmp' ... Inverted model after inversion of interfaces.
# 'sal-m3.out' ... Inverted model (after inversions of interfaces and
#             of velocities).