# 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
  N1=0  N2=1  N3=1  M1='m2.out'
  CAL='inv.cal' GRD1='dm1.out' GRD2='dm2.out'
  grdcal:                                                                     #

# Initial standard deviation 'cv-dif.out' of arrival times
  M1=' '  M2='m2.out'  GM1='gm2.out'  DM1='dm2.out'  SM1='gm0.out'
  gmdmgmt:                                                                    #
  N1=0  N2=1  N3=1  M1=' '
  CAL='sqrt.cal' GRD1='gm0.out' GRD2='cv-dif.out' GRD3=
  grdcal:                                                                     #

# 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
  M1='m1.out'  M2='m2.out'  GM1='gm1.out'  DM1='dm2.out'  SM1='sm1d.out'
  gmdmgmt:                                                                    #
  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='gm4d.out'
  gmgm:                                                                       #
# Controlling smoothness of the final model
  M1='m1.out'  M2=' '   SM1='modsob22.out'  GM1='mod0.out'  GM2='gm4s.out'
  smgm:                                                                       #
  N1=0  N2=0  N3=1  M1='m1.out'  M2=' '
  CAL='subsob.cal'  GRD1='gm4d.out'  GRD2='gm4s.out'  GRD3='gm4.out'
# 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%
  grdcal:                                                                     #
  N1=0  N2=0  N3=1  M1='m1.out'  M2='m1.out'
  CAL='addsob.cal'  GRD1='sm1d.out'  GRD2='modsob22.out'  GRD3='sm1s.out'
  grdcal:                                                                     #
# Limiting the step of iterations
  SOBMUL=0.0
  CAL='addsob.cal'  GRD1='sm1s.out'  GRD2='modsob11.out'  GRD3='sm1.out'
  grdcal:                                                                     #
  SOBMUL=2.0
  CAL='addsob.cal'  GRD1='sm1.out'  GRD2='xstep.out'  GRD3='sm2.out'
  grdcal:                                                                     #
# Continuing matrix operations
  M1='m1.out'               SM1='sm2.out'  SM2='sm3.out'
  sminv:                                                                      #
  M1='m1.out'  M2=' '       SM1='sm3.out'  GM1='gm4.out'  GM2='gm5.out'
  smgm:                                                                       #

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

# Updating source coordinates
  M1MOD='m1mod.out'  M1='m1.out'  MODNEW='gm5.out'  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
  M1=' '  M2='m1.out'  M3='m2.out'
  GM1='gm5.out'  GM2='gm1.out'  GM3='gm8.out'
  gmgm:                                                                       #
  M3=' '
  N1=0  N2=1  N3=1  M1='m2.out'
  CAL='sub.cal' GRD1='gm2.out' GRD2='gm8.out' GRD3='gm9.out'
  grdcal:                                                                     #
  M1=' '  M2='m2.out'  GM1='gm9.out'  DM1='dm2.out'  SM1='gm0.out'
  gmdmgmt:                                                                    #
  N1=0  N2=1  N3=1  M1=' '
  CAL='sqrt.cal' GRD1='gm0.out' GRD2='cv-dev.out' GRD3=
  grdcal:                                                                     #