# History file 'm2d-iwf.h' computes the Incident WaveField
# (the wavefield incidents at the grid) in models M1, M2, M3
# for the 2D ray-based Born approximation.
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

# Input files required
  chk.pl: "data/born/" "m2d-iwf.cal"
# output of 'm2d-mtt.h'
  chk.pl: ""           "m2d-p1s.out"
  chk.pl: ""           "m2d-p2s.out"
  chk.pl: ""           "m2d-p3s.out"
  chk.pl: ""           "m2d-r11u.out"
  chk.pl: ""           "m2d-r12u.out"
  chk.pl: ""           "m2d-r13u.out"
  chk.pl: ""           "m2d-r21u.out"
  chk.pl: ""           "m2d-r22u.out"
  chk.pl: ""           "m2d-r23u.out"
  chk.pl: ""           "m2d-r31u.out"
  chk.pl: ""           "m2d-r32u.out"
  chk.pl: ""           "m2d-r33u.out"
  chk.pl: ""           "m2d-i11u.out"
  chk.pl: ""           "m2d-i12u.out"
  chk.pl: ""           "m2d-i13u.out"
  chk.pl: ""           "m2d-i21u.out"
  chk.pl: ""           "m2d-i22u.out"
  chk.pl: ""           "m2d-i23u.out"
  chk.pl: ""           "m2d-i31u.out"
  chk.pl: ""           "m2d-i32u.out"
  chk.pl: ""           "m2d-i33u.out"
  chk.pl: ""           "m2d-numu.out"

# 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

# Moment tensor
  SM11=1.    SM12=0.    SM13=1.
  SM21=0.    SM22=1.    SM23=0.
  SM31=-1.   SM32=0.    SM33=1.

# Incident wavefield

# 1-st component, real part
  CAL='m2d-iwf.cal'
  GRD1='m2d-p1s.out'
  GRD2='m2d-p2s.out'
  GRD3='m2d-p3s.out'
  GRD4='m2d-r11u.out'
  GRD5='m2d-r12u.out'
  GRD6='m2d-r13u.out'
  GRD7='m2d-aur1.out'
  GRDCALNUM='m2d-numu.out'
  grdcal:

# 2-nd component, real part
  GRD4='m2d-r21u.out'
  GRD5='m2d-r22u.out'
  GRD6='m2d-r23u.out'
  GRD7='m2d-aur2.out'
  grdcal:

# 3-rd component, real part
  GRD4='m2d-r31u.out'
  GRD5='m2d-r32u.out'
  GRD6='m2d-r33u.out'
  GRD7='m2d-aur3.out'
  grdcal:

# 1-st component, imaginary part
  GRD4='m2d-i11u.out'
  GRD5='m2d-i12u.out'
  GRD6='m2d-i13u.out'
  GRD7='m2d-aui1.out'
  grdcal:

# 2-nd component, imaginary part
  GRD4='m2d-i21u.out'
  GRD5='m2d-i22u.out'
  GRD6='m2d-i23u.out'
  GRD7='m2d-aui2.out'
  grdcal:

# 3-rd component, imaginary part
  GRD4='m2d-i31u.out'
  GRD5='m2d-i32u.out'
  GRD6='m2d-i33u.out'
  GRD7='m2d-aui3.out'
  grdcal: