# Least squares migration of Cascadia Dataset 
# This makefile only contains the rule that test one node (OBS A)
# the python codes (under pycode/) actually runs the whole survey

#------------------------------------------------------------------
#----- Velocity Model ---------------------------------------------

Input/velCas2.H:
	cp ../VelModel/Data/velCas2D.H Input/velCas2.H 

#------------------------------------------------------------------
#----- Wavelet ----------------------------------------------------


$I/spike.H:
	Wavelet n1=2048 d1=0.001792 fhigh=100 wavelet=spike tdelay=0.4  | Scale scale_to=1. > $@
	< $@ Spectra | Graph title='Spike wavelet' n1=400 | Tube &
	
	
$I/spkBand40.H: $I/spike.H
	< $I/spike.H Window n1=1000  |  Bandpass flo=5. fhi=30. nphi=8 phase=0 > $@
	< $@ Graph title="zero-phase Bandpass 0-5 30-40"  out=t1.v >/dev/null 
	< $@ Spectra | Graph  title='spectra of bandpass 0-5 30-40' n1=100 out=t2.v >/dev/null
	vp_SideBySideAniso t1.v t2.v > t3.v
	tube t3.v &	

$I/spkBand80.H: $I/spike.H
	< $I/spike.H Window n1=1000  |  Bandpass flo=5. fhi=65. nphi=8 phase=0 > $@
	< $@ Graph title="zero-phase Bandpass 0-5 65-80"  out=t1.v >/dev/null 
	< $@ Spectra | Graph  title='spectra of bandpass 0-5 65-80' n1=200 out=t2.v >/dev/null
	vp_SideBySideAniso t1.v t2.v > t3.v
	tube t3.v &

#------------------------------------------------------------------
#----- image file -------------------------------

Data/mInterpI%.H:
	Interp < 6shots/L2I/m/mCasITotaliter$*.H n1out=160 o1out=1200 d1out=5 > t1.H
	Window3d <t1.H n2=800 f2=200 > $@
	Rm t1.H

Data/mInterpJ%.H:
	Interp < 6shots/L2J/m/mCasJTotaliter$*.H n1out=160 o1out=1200 d1out=5 > t1.H
	Window3d <t1.H n2=800 f2=200 > $@
	Rm t1.H

#------------------------------------------------------------------
#----- Forward Testing --------------------------------------------

Data/ForCas3I.H:
	bin/ForwardOneDir.x velin=Input/velCas2D.H par=Par/ForCasA.par Gx=8231.  Gz=1310. bmask=2500. smooth=0 lpfil=0 VDer=0 gauss=0 MWfile='dummy.H' DW=1 DWfile=syndat/DatWgt1/tmWgtdn3.H wfldsou=souwfld/G/swCasG3.H reflfile=6shots/L2I/m/mCasITotaliter10.H outtag=$@ >/dev/null
	


