first run figure5aStaggeredFDTra.m, figure5bStaggeredFD.m to get seismic records, then run figure6acompareSeimogramsVx,figure6bCompareSeis_recordVx; figure6cCompareSeis_recordTxx, figure6dcompareSeimogramsTxx to get the figures in figure 5 and figure 6.
first run figure8aStaggeredFDTra.m, figure8bStaggeredFD.m to get seismic records, then run figure8cCompareSeis_recordVx, figure8dcompareSeimogramsVx; figure9abCompareSeis_recordTxx, figure9cdCompareSeimogramsTxx to get the figures in figure 8 and figure 9.
first run Figure11aStaggeredFDTra.m,figure11bStaggeredFD.m,figure11cKspace to get seismic records, then run figure11compareSeimogramsVz.m, Figure11CompareSeis_recordVz.m; figure12compareSeimogramsTxx.m, figure12CompareSeis_recordTxx to get the figures in figure 11 and figure 12.
This is only for the convenient of the reviewers.
the codes require the CREWES MATLAB toolbox for boundary condition, the link is https://www.crewes.org/ResearchLinks/FreeSoftware/