###################3D
3D.vel :  
	Vel vc=1500 n1=51 d1=20 o1=0 n2=201 d2=20 o2=-2000 n3=201 d3=20 o3=-2000 	\
		xa1=-1500  xb1=1500   z1=500   alfa1=0    vr1=2000    const1=1 >$@.H datapath=$(PWD)/

3D.vel.exx:
	Window3d < /data/exxon-311/Vel/exxon-velocity.H min1=1000 j1=2 max1=4000 min2=13000 max2=18000 min3=2000 max3=7000 > $@; echo o1=0 >> $@

REF = reflectivity.x 
3D.ref.H: 3D.vel 
	cd ../../../Src; make $(BIN)/Triangle.x; cd $(PWD)
	$(REF)  vel=$< | Triangle.x wind3=2 wind2=2 wind1=2 >3D.ref$*
	Transp < 3D.ref$* | Transp plane=23 > 3D.ref$*.H 
	Transp < 3D.vel$* | Transp plane=23 | Smooth rect3=3 > 3D.vel$*.H

3D.green: greens_mpi 
	cd PSPI; make $G; cd ../
	#Wavelet n1=1200 d1=0.004 wavelet=ricker2 domain=time tdelay=0.04 fund=25 fhigh=50 | Transf f_min1=3 f_max=60 f_max1=45 | Transp plane=35> wav1-3D-wave.H
	Wavelet n1=1200 d1=0.004 wavelet=ricker2 domain=time tdelay=0.04 fund=25 fhigh=50 | Transf f_min1=3 f_max=45 f_max1=30 | Transp plane=35> wav1-3D-wave.H
	make 3D.ref.H 
	greens_mpi.x vel=3D.vel.H greens=$@.H wavelet=wav1-3D-wave.H  report=green.3D.log wfld_d=$@-wfld_d.H 		\
	nx_modl=    dx_modl=    ox_modl=    oy_modl=  dy_modl= ny_modl= nz_modl=    dz_modl=  oz_modl=0			\
	ox_surf=15000   dx_surf=    nx_surf=1   oy_surf=0 dy_surv= ny_surf=1  jw=1  nodes=15 nxpad_beg= nxpad_end=	\
	sign=+1 vnum=2 taper_pct=0.1 taper_type=2 tmin=-0.2 tmax=0.2							\
	only_green=n taper=y spike=n  wfld_d_out=y lloyd=n ff loop=w allfreq=n encode=n  random_constant=n 		\
	code_it=n ffd=y verb=1 jw=1 center_sou=y 
	mv $@.H data/ ;	touch $@

3DGreen = $(PBS_O_WORKDIR)/3DGreen/ 
3D.green.par: greens_mpi 
	cd PSPI; make $G; cd ../
	Wavelet n1=1200 d1=0.004 wavelet=ricker2 domain=time tdelay=0.04 fund=25 fhigh=50 | Transf f_min1=3 f_max=50 f_max1=30 | Transp plane=35 | Window3d squeeze=n  > wav1-3D-wave.H
	echo .datapath=$(3DGreen) > .datapath
	mkdir -p $(3DGreen)
	cp wav1-3D-wave.H $(3DGreen) 
	cp 3D.vel.H $(3DGreen)
	cp 3D.pars.par $(3DGreen)
	make greens_mpi
	./prepare_files.py 1 30 wav1-3D-wave.H 3DGreen
	./green.py 1 30 3DGreen
