-- Notes On Reproducibility --

All of the figures are CR.
The following python files will contains all the rules and references to codes for running.
However, the Input (original raw data) is absence. While SEP has permission for the field data
to be published, there is no permission for its release.

# #######################################################
To generate the pre-processing (Figure 5):
*Absence Input: original segy files used in pyExtract.py

Run the following in sequence


1. python pycode/pyExtract.py  # Convert from segy to sep format
2. python pycode/pyTMall.py    # Zero out noisy traces
3. python pycode/pyPEFtm.py    # Apply gapped decon
4. python pycode/pySDE.py      # Interpolation and bandpass, convert to zerophase
5. python pycode/pySDFG.py     
6. make dataCom.pdf            # Generate figure

# #######################################################
To generate the migration and LSM result (Figure 6) :
* Absence Input: Velocity model for Cascadia dataset (Input/velCas2Dnew.H)

Run the following in sequence
# Preprocessing must be done before inversion


# Calculate the source side wavefields
1. make Input/spkBand80.H     # Create input wavelet
2. python pycode/pySWB.py     # calculate incident wflds
3. python pycode/pySWCD.py    
4. python pycode/pySWE.py
5. python pycode/pySWFG.py


# generate data space weighting file for inversion
6. python pycode/pyMuteWgt.py

# Run inversion (*** Long Compute Time ***)
7. python pyinv/pysolverRegI.py   # For LSRTM-HM inversion (will need 3 nodes)
8. python pyinv/pysolverRegJ.py   # For LSRTM-M inversion (will need 3 nodes)

# Generate Figures
9. make RTMinv.pdf



# #######################################################
To generate the Zoom image (Figure 7)

1. make Zoom.pdf



# #######################################################
To generate the Data residual plot (Figure 9)


1. make Data/ForCas3I.H
2. make DatspaceResid.pdf
