# History file 'cv-inv1.h' to perform a single iteration of
# simultaneous inversion for a velocity model and hypocentral
# coordinates.
# The history file should be run from the history file 'cv-inv.h'.
#
# Input files required
  chk.pl: "" "cv-mod.out"
  chk.pl: "" "cv-src.out"
  chk.pl: "" "rp01.out" # all files listed in 'cv-crtp.dat' and 'cv-crts.dat'
                        # are required, but only 'rp01.out' is checked for

# Initial, template and updated models
  MODEL='cv-mod.out' MODIN='cv-mod.out'  MODOUT='cv-mod.out'  NEGPAR=0

# Calculating matrices for inversion
# 557 P-arrivals + 940 S-arrivals = 1497 arrivals  (M2=1497 data points)
  PTS='cv-rec.dat'  PTS1='cv-src.out'
  ERRMUL=1.  # ERRMUL is not used by 'invtt.for' ver.5.70
  M1='m1.out'  M2='m2.out'  INVSRC='invsrc.out'
  GM1='gm1.out'  GM2='gm2.out'  GM3='dat0.out'  DM1='dm1.out'
  FTT='cv-arrp.dat'  CRTOUT='cv-crtp.dat'  DIST=10.0
  invtt:                                                                      #
  M1IN='m1.out'  M2IN='m2.out'
  FTT='cv-arrs.dat'  CRTOUT='cv-crts.dat'  DIST=10.0
  invtt:                                                                      #
  M1IN=' '  M2IN=' '

# Inverting the matrix of data variances
  MATIN1='dm1.out'                       MATOUT='dm2.out'
  MATFUN='inv'
  matfun:

# Initial standard deviation 'cv-dif.out' of arrival times
  MATIN1='dm2.out'   MATIN2='gm2.out'    MATOUT='dm2gm2.out'
  SYMMETRY=      MATT1=  MATT2=
  matmul:
  MATIN1='gm2.out'   MATIN2='dm2gm2.out' MATOUT='gm0.out'
  SYMMETRY='diag' MATT1=1 MATT2=
  matmul:
  MATIN1='gm0.out'                       MATOUT='cv-dif.ou'
  SPARSE=-1 MATFUN='sqrt'
  matfun:
  copy.pl: "cv-dif.ou@" "cv-dif.out"

# Controlling smoothness of the final model
  M1='m1mod.out'  M1LOC='m1.out'
  SOBOLEV='sob22.dat'
  SOBOLEV='cv-sob.dat'
  MODPAR='mod0.out'  MODL2=' '  MODSOB='modsob22.out'
  SOBW01=1.0
  SOBW02=1.7
  invsoft:                                                                    #

# Controlling the step of iterations
  SOBOLEV='sob11.dat'
  MODPAR=' '  MODL2=' '  MODSOB='modsob11.out'
  SOBW01=1.0
  SOBW02=1.7
  invsoft:                                                                    #
  SOBOLEV=' '
  MODPAR=' '  MODL2=' '  MODSOB='xstep.out'
  XDLOC=100.0
  invsoft:                                                                    #

# Matrix operations
  MATIN1='dm2.out'   MATIN2='gm1.out'    MATOUT='dm2gm1.out'
  SYMMETRY=      MATT1=  MATT2=1
  matmul:
  MATIN1='gm1.out'   MATIN2='dm2gm1.out' MATOUT='sm1d.out'
  SYMMETRY='sym' MATT1=  MATT2=
  matmul:
  MATIN1='dm2.out'   MATIN2='gm2.out'    MATOUT='gm3.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
  MATIN1='gm1.out'   MATIN2='gm3.out'    MATOUT='gm4d.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
# Controlling smoothness of the final model
  MATIN1='modsob22.out'   MATIN2='mod0.out'    MATOUT='gm4s.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:
# SOBMUL=0.1           16.52 19%
# SOBMUL=0.5  >50%     16.87 17%
# SOBMUL=1.    23%     17.05 16%
# SOBMUL=1.5   12%     17.15 16%
# SOBMUL=2.    11%     17.26 15%
# SOBMUL=5.     9%     17.91 12%
# SOBMUL=10.    7%     18.56  9%
# SOBMUL=20.    3%     19.26  5%
# SOBMUL=100.   0.1%   20.21  0.4%
# SOBMUL=1000.  0.002% 20.30  0%

# SOBMUL=1.    14%     16.35 19%
# SOBMUL=1.5    9%     16.61 18%
# SOBMUL=2.     6%     16.76 17%
# SOBMUL=3.     4%     16.96 16%
# SOBMUL=5.     1.4%   17.31 15%
  MATIN1='gm4d.out'  MATIN2='gm4s.out' MATOUT='gm4.out'
  COEF1=  COEF2=-25.
  matlin:
  MATIN1='sm1d.out'  MATIN2='modsob22.out' MATOUT='sm1s.out'
  COEF1=  COEF2=25.
  matlin:
# Limiting the step of iterations
  MATIN1='sm1s.out'  MATIN2='modsob11.out' MATOUT='sm1.out'
  COEF1=  COEF2=0.
  matlin:
  MATIN1='sm1.out'  MATIN2='xstep.out' MATOUT='sm2.out'
  COEF1=  COEF2=4.
  matlin:
# Continuing matrix operations
  MATIN1='sm2.out'                       MATOUT='sm3.out'
  matinv:
  MATIN1='sm3.out'   MATIN2='gm4.out'    MATOUT='gm5.out'
  SYMMETRY=' '   MATT1=  MATT2=
  matmul:

# Updating the model
  M1='m1mod.out'  MODNEW='gm5.out'
  modmod:                                                                     #

# Updating source coordinates
  M1MOD='m1mod.out'  M1='m1.out'  MODNEW='gm5.ou@'  INVSRC='invsrc.out'
  SRC='cv-src.out'  SRCNEW='cv-srcn.out'
  srcsrc:                                                                     #

# Distances of new sources from the reference vertical plane
  MODEL='cv-srf.dat'
  INTF='cv-srcn.out'  INTFOUT='cv-dist.out'  KSRFC=1  KOLUMN=4
  intf:                                                                       #

# Estimated standard deviation 'cv-dev.out' of arrival times
  MATIN1='gm1.out'   MATIN2='gm5.out'    MATOUT='gm8.out'
  SYMMETRY=      MATT1=1 MATT2=
  matmul:
  MATIN1='gm2.out'   MATIN2='gm8.out' MATOUT='gm9.out'
  COEF1=  COEF2=-1
  matlin:
  MATIN1='dm2.out'   MATIN2='gm9.out'    MATOUT='dm2gm9.out'
  SYMMETRY=       MATT1=  MATT2=
  matmul:
  MATIN1='gm9.out'   MATIN2='dm2gm9.out' MATOUT='gm0.out'
  SYMMETRY='diag' MATT1=1 MATT2=
  matmul:
  MATIN1='gm0.out'                       MATOUT='cv-dev.ou'
  SPARSE=-1 MATFUN='sqrt'
  matfun:
  copy.pl: "cv-dev.ou@" "cv-dev.out"