# History file 'm2d-mtt.h' computes the perturbations of the elastic
# parameters, the Green functions from the source and from the receiver,
# slowness vectors, travel times and other quantities at the gridpoints.
# Designed for 2D computations in models M1, M2 and M3.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/born/" "pert.cal"
  chk.pl: "data/born/" "m1-mod.dat"
#  chk.pl: "data/born/" "m2-mod.dat"
#  chk.pl: "data/born/" "m3-mod.dat"
  chk.pl: "data/born/" "m-bg-mod.dat"
  chk.pl: "data/born/" "m2d-src.dat"
  chk.pl: "data/born/" "m2d-rec.dat"
  chk.pl: "data/born/" "m2d-cod.dat"
  chk.pl: "data/born/" "m12d-crt.dat"
#  chk.pl: "data/born/" "m22d-crt.dat"
#  chk.pl: "data/born/" "m32d-crt.dat"
  chk.pl: "crt/"       "writall.dat"

# Grid dimensions for 2D computations in model M1
  N1=100   N2=1   N3=400
  D1=0.1   D2=0.  D3=0.025
  O1=0.05  O2=5.  O3=10.0125

## Grid dimensions for 2D computations in model M2
#  N1=343        N2=1   N3=800
#  D1=0.043723   D2=0.  D3=0.017489
#  O1=-4.978139  O2=5.  O3=6.017489
#
## Grid dimensions for 2D computations in model M3
#  N1=600        N2=1   N3=3200
#  D1=0.033333   D2=0.  D3=0.00625
#  O1=-4.983333  O2=5.  O3=0.003125

# Discretization of the elastic parameters and computation of
# the perturbations of elastic parameters
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Perturbed model (model M1)
  MODEL='m1-mod.dat'
  LFREE=1
  #modchk:

## Perturbed model (model M2)
#  MODEL='m2-mod.dat'
#  LFREE=1
#  #modchk:
#
# Perturbed model (model M3)
#  MODEL='m3-mod.dat'
#  LFREE=1
#  #modchk:


# Discretization of the elastic parameters at the gridpoints
# in the perturbed model.
# P-wave velocity
  MPAR=1
  VEL='m2d-vp1.out'
  grid:
# S-wave velocity
  MPAR=2
  VEL='m2d-vs1.out'
  grid:
# Density
  MPAR=3
  VEL='m2d-rh1.out'
  grid:

# Background (unperturbed) model
  MODEL='m-bg-mod.dat'
  LFREE=1
  #modchk:

# Discretization of the elastic parameters at the gridpoints
# in the unperturbed model.
# P-wave velocity
  MPAR=1
  VEL='m2d-vp2.out'
  grid:
# S-wave velocity
  MPAR=2
  VEL='m2d-vs2.out'
  grid:
# Density
  MPAR=3
  VEL='m2d-rh2.out'
  grid:

# Computation of the perturbations of the elastic parameters at
# the gridpoints.
  CAL='pert.cal'

  GRD1='m2d-vp1.out'
  GRD2='m2d-vp2.out'
  GRD3='m2d-vpp.out'

  GRD4='m2d-vs1.out'
  GRD5='m2d-vs2.out'
  GRD6='m2d-vsp.out'

  GRD7='m2d-rh1.out'
  GRD8='m2d-rh2.out'
  GRD9='m2d-rhp.out'
  grdcal:

# Computation of the Green function from the receiver.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Source of the Green function
  SRC='m2d-rec.dat'

# Data filenames for CRT (used also in the next computation)
  CODE='m2d-cod.dat'
  WRIT='writall.dat'  # all rays are required by 'mtt.for'

# Data filenames for CRT

# Data filenames for CRT in model M1, 2D computation
  DCRT='m12d-crt.dat'
  RPAR='m12d-crt.dat'
  CRTLOG='m2d-crtg.out'

## Data filenames for CRT in model M2, 2D computation
#  DCRT='m22d-crt.dat'
#  RPAR='m22d-crt.dat'
#  CRTLOG='m2d-crtg.out'
#
## Data filenames for CRT in model M3, 2D computation
#  DCRT='m32d-crt.dat'
#  RPAR='m32d-crt.dat'
#  CRTLOG='m2d-crtg.out'

# Kind of ray parameters (used also in the next computation)
  INIPAR=3

# Ray tracing
  crt:

# MTT output files
  NUM='m2d-numg.out'
  MTT='m2d-ttg.out'
  MP1='m2d-p1g.out'  MP2='m2d-p2g.out'  MP3='m2d-p3g.out'
  MTT22='m2d-tt2g.out' #for the correction term in 2D Born approximation
  AMPR11='m2d-r11g.out'  AMPI11='m2d-i11g.out'
  AMPR21='m2d-r21g.out'  AMPI21='m2d-i21g.out'
  AMPR31='m2d-r31g.out'  AMPI31='m2d-i31g.out'
  AMPR12='m2d-r12g.out'  AMPI12='m2d-i12g.out'
  AMPR22='m2d-r22g.out'  AMPI22='m2d-i22g.out'
  AMPR32='m2d-r32g.out'  AMPI32='m2d-i32g.out'
  AMPR13='m2d-r13g.out'  AMPI13='m2d-i13g.out'
  AMPR23='m2d-r23g.out'  AMPI23='m2d-i23g.out'
  AMPR33='m2d-r33g.out'  AMPI33='m2d-i33g.out'
  mtt:

# Computation of the Green function from the source.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Source of the Green function
  SRC='m2d-src.dat'

# Data filenames for CRT
  CRTLOG='m2d-crtu.out'

# Ray tracing
  crt:

# MTT output files
  NUM='m2d-numu.out'
  MTT='m2d-ttu.out'
  MP1='m2d-p1u.out'  MP2='m2d-p2u.out'  MP3='m2d-p3u.out'
  MTT22='m2d-tt2u.out' #for the correction term in 2D Born approximation
  AMPR11='m2d-r11u.out'  AMPI11='m2d-i11u.out'
  AMPR21='m2d-r21u.out'  AMPI21='m2d-i21u.out'
  AMPR31='m2d-r31u.out'  AMPI31='m2d-i31u.out'
  AMPR12='m2d-r12u.out'  AMPI12='m2d-i12u.out'
  AMPR22='m2d-r22u.out'  AMPI22='m2d-i22u.out'
  AMPR32='m2d-r32u.out'  AMPI32='m2d-i32u.out'
  AMPR13='m2d-r13u.out'  AMPI13='m2d-i13u.out'
  AMPR23='m2d-r23u.out'  AMPI23='m2d-i23u.out'
  AMPR33='m2d-r33u.out'  AMPI33='m2d-i33u.out'
  MP4='m2d-p1s.out'  MP5='m2d-p2s.out'  MP6='m2d-p3s.out'
  mtt: