norneG = $(PBS_O_WORKDIR)/norneH/
msp = cd ~/research.ln/Source/ShotMig/ ; make ; cd 
mgr = cd ~/research.ln/Source/green_new86_64/ ; make ; cd 
mgh = cd /home/goayeni/Research/Source/Hessian; make ; cd 
HH2 = $B/hessian_new.x
hess:
	$(mgh)

#compute
hess-norne-06-%: 
	$(mgh); 
	./hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grn-norne-06-$* hess-norne-06-$* hs.P 0 ; touch $@
hess-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grn-norne-06-$*  $@               hs-offd.P 0 ; touch $@
hesss-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grns-norne-06-$* $@               hs-offd-me.P      0 ; touch $@
hesss-hole-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grns-norne-06-$* $@               hs-offd-me-hole.P 0 ; touch $@
hessss-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grns-norne-06-$* $@               hs-offd-me.P      0 ; touch $@
hessss-hole-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grns-norne-06-$* $@               hss-offd-me-hole.P 0 ; touch $@
hessss-hole-sxgx-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grns-norne-06-$* $@               hss-offd-me-sxgx-hole.P 0 ; touch $@
#new--for dmult data
hessss-dmult-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 59 /data/sep/goayeni/workspace grns-norne-dmult-06-$* $@         hss-offd-dmult.P 0 ; touch $@
hessss-hole-sxgx-dmult-offd-norne-06-%: 
	$(mgr); 
	./hess2d.py $* 0 0 59 /data/sep/goayeni/workspace grns-norne-dmult-06-$* $@         hss-offd-dmult-sxgx-hole.P 0 ; touch $@
#add
add-hess-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grn-norne-06-$* hess-offd-norne-06-$*  	   	hs-offd.P         	0 ; touch $@
add-hesss-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grn-norne-06-$* hesss-offd-norne-06-$*       	hs-offd-me.P      	0 ; touch $@
add-hesss-hole-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 60 /data/sep/goayeni/workspace grn-norne-06-$* hesss-hole-offd-norne-06-$*  	hs-offd-me-hole.P 	0 ; touch $@
add-hessss-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grn-norne-06-$* hessss-offd-norne-06-$*  	   	hs-offd-mee.P     	0 ; touch $@
add-hessss-hole-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grn-norne-06-$* hessss-hole-offd-norne-06-$* 	hss-offd-me-hole.P 	0 ; touch $@
add-hessss-hole-sxgx-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grn-norne-06-$* hessss-hole-sxgx-offd-norne-06-$* hss-offd-me-sxgx-hole.P 0 ; touch $@
add-hessss-w4-offd-norne-06-%: 
	./add-w4-hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grn-norne-06-$* hessss-offd-norne-06-$* 	hs-offd-mee.P 	        0 ; touch $@
add-hessss-w4-hole-sxgx-offd-norne-06-%: 
	./add-w4-hess2d.py $* 0 0 80 /data/sep/goayeni/workspace grn-norne-06-$* hessss-hole-sxgx-offd-norne-06-$* hss-offd-me-sxgx-hole.P 0 ; touch $@
add-hessss-dmult-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 59 /data/sep/goayeni/workspace grn-norne-dmult-06-$* hessss-dmult-offd-norne-06-$*           hss-offd-dmult.P           0 ; touch $@
add-hessss-hole-sxgx-dmult-offd-norne-06-%: 
	./add-hess2d.py $* 0 0 59 /data/sep/goayeni/workspace grn-norne-dmult-06-$* hessss-hole-sxgx-dmult-offd-norne-06-$* hss-offd-dmult-sxgx-hole.P 0 ; touch $@


test_hess0: 
	$(mgh)
	@date
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grn-norne-06-2759/wave/gr1-441/w10-10/wav-gr1-441--w10-10.H > wav-$@.H 
	$(HH2) greens=/data/sep/goayeni/workspace/grn-norne-06-2759/green/gr1-441/w10-10/gr1-441--w10-10.H  hessian=$@.H 					\
	ngx=261      ogx=-3250 dgx=25 nsx=110  osx=500 dsx=25				\
	nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0 					\
	zmin_modl= 500 zmax_modl= 500							\
	nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0 				\
	marine=y									\
	nodes=7 spike=y verb=3 								\
	report=$@.log jw=1 wavelet=wav-$@.H

test_hess: 
	$(mgh)
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-06-2759/wave/gr1-882/w11/wav-gr1-882--w11.H > wav-$@.H 
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-06-2759/green/gr1-882/w11/gr1-882--w11.H hessian=$@.H \
		par=hss-offd-me-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		nx_modl=   ox_modl=   nz_modl=  oz_modl=                                        \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H hole=n 
test_hess-sx: 
	$(mgh)
	mkdir -p /tmp/goayeni/$@/
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-06-2759/wave/gr1-882/w11/wav-gr1-882--w11.H > wav-$@.H datapath=/tmp/goayeni/$@/
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-06-2759/green/gr1-882/w11/gr1-882--w11.H hessian=$@.H \
		par=hss-offd-me-sx-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		nx_modl=   ox_modl=   nz_modl=  oz_modl=                                        \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H 
test_hess-sxgx: 
	mkdir -p /tmp/goayeni/$@/
	$(mgh)
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-06-2759/wave/gr1-882/w11/wav-gr1-882--w11.H > wav-$@.H datapath=/tmp/goayeni/$@/
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-06-2759/green/gr1-882/w11/gr1-882--w11.H hessian=$@.H \
		par=hss-offd-me-sxgx-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		nx_modl=   ox_modl=   nz_modl=  oz_modl=                                        \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H 


#try pre-offd
test_hess-hx_offd: 
	mkdir -p /tmp/goayeni/$@/
	$(mgh)
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-06-2759/wave/gr1-882/w11/wav-gr1-882--w11.H > wav-$@.H datapath=/tmp/goayeni/$@/
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-06-2759/green/gr1-882/w11/gr1-882--w11.H hessian=$@.H \
		par=hss-offd-me-sxgx-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=61 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H \


#single hole demult....
dmult_hess: 
	mkdir -p /data/sep/goayeni/workspace/norne/
	echo datapath=/data/sep/goayeni/workspace/norne/ >> .datapath
	rm -rf $@.H
	$(mgh)
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-dmult-06-2759/wave/gr1-882/w3/wav-gr1-882--w3.H > wav-$@.H datapath=/tmp/goayeni/
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-dmult-06-2759/green/gr1-882/w3/gr1-882--w3.H hessian=$@.H \
		par=hss-offd-dmult-sxgx-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=400    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		nx_modl=   ox_modl=   nz_modl=  oz_modl=                                        \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H  \
		hole=n

		
dmult_hess-sxgx: 
	mkdir -p /data/sep/goayeni/workspace/norne/
	echo datapath=/data/sep/goayeni/workspace/norne/ >> .datapath
	rm -rf $@.H
	$(mgh)
	Window3d n4=1 squeeze=n < /data/sep/goayeni/workspace/grns-norne-dmult-06-2759/wave/gr1-882/w3/wav-gr1-882--w3.H > wav-$@.H datapath=/tmp/goayeni/
	hessian_new_test.x greens=/data/sep/goayeni/workspace/grns-norne-dmult-06-2759/green/gr1-882/w3/gr1-882--w3.H hessian=$@.H \
		par=hss-offd-dmult-sxgx-hole.P \
		ngx=261    ogx=-3250 dgx=25 nsx=2  osx=2000 dsx=25 00                           \
		ngx=1      ogx=0     dgx=25 nsx=520    osx=500 dsx=25 00                      	\
		ngx=261    ogx=-3250 dgx=25 nsx=400    osx=500 dsx=25 00                      	\
		nx_offd=1  ox_offd=0  nz_offd=1 oz_offd=0                                       \
		nx_modl=   ox_modl=   nz_modl=  oz_modl=                                        \
		zmin_modl= 1500 zmax_modl= 500                                                  \
		nhx_offd=1  ohx_offd=0    nhx_modl=1 ohx_modl=0                                 \
		marine=y                                                                        \
	        nodes=15 spike= y verb=3                                                        \
	        report=$@.log jw=1 wavelet=wav-$@.H 

hess-all:
		make -n hessss-dmult-offd-norne-06-2759
		make -n hessss-hole-sxgx-dmult-offd-norne-06-2759
		make -n add-hessss-dmult-offd-norne-06-2759
		make -n add-hessss-dmult-offd-norne-08-2759
