PARDIR  = ./Par
P	= ${PARDIR}

# SEAM_EARTH_MODEL
VEL_DATA = /data/seam/earthmodel/Vp_xyz_10m.bin
RHO_DATA = /data/seam/earthmodel/Rho_xyz_10m.bin

wavelet_25f_1ms.H:
	Wavelet n1=8192 d1=0.0005 wavelet=ricker2 phase=0 fund=20 tdelay=0.075 > tmp.H
	Math file1=tmp.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm tmp.H

%_spec.H: %.H
	Rtoc < $*.H > cmplx$*.H
	Pad n2=512 < cmplx$*.H | Ft3d sign1=0 sign2=1 center1=0 center2=1 | Cabs > $@
	rm cmplx$*.H

### Run Modeling #########################################################

vel.H: ${VEL_DATA}
	echo 	o1=0 o2=0 o3=0  >>$@
	echo	n1=3501 n2=4001 n3=1501  >>$@
	echo	d1=10 d2=10 d3=10  >>$@
	echo	data_format=""xdr_float""  >>$@
	echo	label1="x(m)" label2="y(m)" label3="z(m)" >>$@
	echo	esize= 4  >>$@
	echo in="$<" >> $@

%_wind.H: %.H
	Window3d f1=1100 n1=500 f2=2220 n2=500 n3=600 < $<  > $*2.H
	#f1 was 1000 for gp257 images
	cp $*2.H $@
	swapbytes < /scratch/chrisl/$*2.H@ > /scratch/chrisl/$@@
	echo in="/scratch/chrisl/$@@" >> $@
	rm -rf $*2.H

vel3d.H: vel_wind.H
	Transp plane=12 reshape=2,3,4 < $< > $@
	echo "o1=0 d1=10. o2=0 d2=10. o3=0 d3=10." >> $@

survey3d.H: vel3d.H vel3d_sea.H wavelet_25f_1ms.H model3d.x ${P}/surv3d.p
	model3d.x < vel3d.H par=${P}/surv3d.p wavelet=wavelet_25f_1ms.H > su_tran.H
	model3d.x < vel3d_sea.H par=${P}/surv3d.p wavelet=wavelet_25f_1ms.H > su_tran_sea.H
	Math file1='su_tran.H' file2='su_tran_sea.H' exp='file1-file2' > surv_tran.H
	Transp plane=12 reshape=2,3,4,5 < surv_tran.H > $@
	rm su_tran.H

shot3d_basic.H: vel3d_basic.H wavelet_25f_1ms.H model3d.x ${P}/par3d_basic.p
	Window3d n1=5000 < wavelet_25f_1ms.H > wave2.H
	model3d.x < vel3d_basic.H par=${P}/par3d_basic.p wavelet=wave2.H > s_basic.H
	Transp plane=12 reshape=2,3,4,5 < s_basic.H > $@

shot3d_basic_sea.H: vel3d_basic_sea.H wavelet_25f_1ms.H model3d.x ${P}/par3d_basic.p
	Window3d n1=5000 < wavelet_25f_1ms.H > wave2.H
	model3d.x < vel3d_basic_sea.H par=${P}/par3d_basic.p wavelet=wave2.H > s_basic_sea.H
	Transp plane=12 reshape=2,3,4,5 < s_basic_sea.H > $@

shot3d_basic2.H: vel3d_basic.H wavelet_25f_1ms.H model3d.x ${P}/par3d_basic2.p
	Window3d n1=5000 < wavelet_25f_1ms.H > wave2.H
	model3d.x < vel3d_basic.H par=${P}/par3d_basic2.p wavelet=wave2.H > s_basic2.H
	Transp plane=12 reshape=2,3,4,5 < s_basic2.H > $@

shot3d_basic_sea2.H: vel3d_basic_sea.H wavelet_25f_1ms.H model3d.x ${P}/par3d_basic2.p
	Window3d n1=5000 < wavelet_25f_1ms.H > wave2.H
	model3d.x < vel3d_basic_sea.H par=${P}/par3d_basic2.p wavelet=wave2.H > s_basic_sea2.H
	Transp plane=12 reshape=2,3,4,5 < s_basic_sea2.H > $@

vel3d_basic.H:
	Vel vc=1500 n1=200 d1=10 o1=0 n2=200 d2=10 o2=0 n3=200 d3=10 o3=0 	\
		xa1=0  xb1=4000   z1=400   alfa1=20    vr1=2100    const1=1 > $@

vel3d_basic_sea.H:
	Vel vc=1500 n1=200 d1=10 o1=0 n2=200 d2=10 o2=0 n3=200 d3=10 o3=0 > $@

