# History file 'u2d-inv.h' to run 'invsoft.for' and 'invtt.for'
# with demo data in 2-D model UNCONFORMITY
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Input files required:
 #chk.pl: "data/u2d/" "u2d-mod.dat"
 #chk.pl: "data/u2d/" "u2d-src.dat"
 #chk.pl: "data/u2d/" "u2d-rec.dat"
 #chk.pl: "data/u2d/" "u2d-rpa.dat"
 #chk.pl: "data/u2d/" "u2d-wri.dat"
 #chk.pl: "data/u2d/" "u2d-pts.dat"
 #chk.pl: "data/u2d/" "u2d-ftt.dat"
 #chk.pl: "data/u2d/" "u2d-crt1.dat"
 #chk.pl: "model/"    "sob22.dat"

# Initial model (model UNCONFORMITY by Cormier and Mellen 1984)
  MODEL='u2d-mod.dat'

# Calculating L2 and Sobolev scalar products of basis functions
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# Data describing Sobolev scalar products of basis functions
  SOBOLEV='sob22.dat'  # Minimizing second derivatives
  M1='u2d-m1.out'      # Output number of basis functions

# L2 scalar products of surface, velocity and density basis functions
  MODL2='u2d-l2.out'  MODSOB=' '  SOBW00=1  SOBW01=1  SOBW02=1  SOBW03=1
  invsoft: 'u2d-inv.h' /

# Sobolev scalar products of surface basis functions
  MODL2=' '  MODSOB='u2d-sobs.out'  SOBW00=1  SOBW01=0  SOBW02=0  SOBW03=0
  invsoft: 'u2d-inv.h' /

# Sobolev scalar products of velocity basis functions
  MODL2=' '  MODSOB='u2d-sobv.out'  SOBW00=0  SOBW01=1  SOBW02=1  SOBW03=0
  invsoft: 'u2d-inv.h' /

# Ray tracing
# ~~~~~~~~~~~
# Source and receivers
  SRC='u2d-src.dat'  REC='u2d-rec.dat'

# Data filenames for CRT (for details see 'crt.for')
  DCRT='u2d-crt1.dat'
  CODE='u2d-crt1.dat'
  RPAR='u2d-rpa.dat'
  WRIT='u2d-wri.dat'
  CRTLOG='u2d-log.out'

# Kind of ray parameters (for details see 'init.for')
  INIPAR=2  # geographic-like spherical coordinates

# 2-D two-point ray tracing
  crt: 'u2d-inv.h' /

# Travel-time derivatives with respect to basis functions
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  CRTOUT=' '           # default names of output files of program CRT
  PTS='u2d-pts.dat'    # coordinates of shot and receiver points
  FTT='u2d-ftt.dat'    # travel times to be inverted
  DIST=0.0001          # maximum distance between the point and rays
  VPOWER=1             # velocity power for the travel-time check sum
  M1=' '               # output file with number of model coefficients
  M2='u2d-m2.out'      # output file with number of values to invert
  GM1='u2d-gm1.out'    # matrix of the derivatives of travel times
  GM2='u2d-gm2.out'    # vector of travel-time residuals
  GM3='u2d-gm3.out'    # vector of travel times in the model
  DM1='u2d-dm1.out'    # diagonal matrix of travel-time variances
  INVLOG='u2d-itt.out' # output log file

# Generating system of equations for travel-time inversion
  invtt: 'u2d-inv.h' /

# ======================================================================

# Output files of 'invsoft.for':
# 'u2d-m1.out'   Number of basis functions
# 'u2d-l2.out'   L2 scalar products of all basis functions
# 'u2d-sobs.out' Sobolev scalar products of surface basis functions
# 'u2d-sobv.out' Sobolev scalar products of velocity basis functions

# Temporary files (may be deleted):
# 'r0*.out'
# 's0*.out'

# Output files of 'invtt.for':
# 'u2d-m2.out'
# 'u2d-gm1.out'
# 'u2d-gm2.out'
# 'u2d-gm3.out'
# 'u2d-dm1.out'
# 'u2d-itt.out'