norne3d = $(PBS_O_WORKDIR)/norne3d/
msp = cd ~/research.ln/Source/ShotMig/ ; make ; cd 
norne-3d-08-%: 
	$(msp)
	echo datapath=$(norne3d) > .datapath
	mkdir -p $(norne3d)
	./migrate3d.py $* 99 4 sail-line-08-2d-$*-half-rmo.H $@ 0 3d.p 1800; touch $@
norne-3d-06-%: 
	$(msp)
	echo datapath=$(norne3d) > .datapath
	mkdir -p $(norne3d)
	./migrate3d.py $* 99 4 sail-line-06-2d-$*-half-rmo.H $@ 0 3d.p 1800; touch $@
add-3d-08-%:
	./add3d.py $* 99 4 sail-line-08-2d-$*-half-rmo.H norne-3d-08-$* 0 ; touch $@
add-3d-06-%:
	./add3d.py $* 99 4 sail-line-06-2d-$*-half-rmo.H norne-3d-06-$* 0 ; touch $@
shift3d-%:
	Transp plane=12 reshape=1,3,6 < /data/sep/goayeni/workspace/norne3d-06-$*/modl/norne3d-06-$*.H > norne3d-06-$*-mig.H
	Transp plane=12 reshape=1,3,6 < /data/sep/goayeni/workspace/norne3d-08-$*/modl/norne3d-08-$*.H > norne3d-08-$*-mig.H
	cd ~/research.ln/Source/Oneway/; make $(bin)/joint_encode.x ; cd 
	$B/joint_encode.x report=$@.log modl=mod1.H  \
	modl1=norne3d-06-$*-mig.H rm1=norne3d-06-$*-mig.H  smodl1=norne3d-06-$*-mig-s.H ts1=norne3d-06-$*-mig-ts.H sd1=$@-tsd1.H cc1=$@-cc1.H rr1=$@-rr1.H \
	modl2=norne3d-08-$*-mig.H rm2=norne3d-08-$*-mig.H  smodl2=norne3d-08-$*-mig-s.H ts2=norne3d-08-$*-mig-ts.H sd2=$@-tsd2.H cc2=$@-cc2.H rr2=$@-rr2.H \
        nmix=3 par=shift.p nodes=1 shift_nodes=1
	touch $@

diff3d-%: shift3d-%
	Math file1=/data/sep/goayeni/workspace/norne3d-06-$*/modl/norne3d-06-$*.H file2=/data/sep/goayeni/workspace/norne3d-08-$*/modl/norne3d-08-$*.H exp=file2-file1 > b4-shift-$*.H 
	Math file1=norne3d-06-$*-mig.H file2=norne3d-08-$*-mig-s.H exp=file2-file1 | Window3d | Transp > af-shift-$*.H
	Cat3d b4-shift-$*.H af-shift-$*.H | Window3d min1=2000 max1=4500 min2=500 max2=12000 | Grey color=G newclip=1 wantscalebar=y title='xline: $*(m)'> /dev/null out=$@.v
	Window3d f5=0 f4=0 min3=2000 max3=4500 min1=500 max1=12000 n5=1 n4=1 < norne3d-08-$*-mig-ts.H | Window3d | Transp | Grey color=j newclip=1 wantscalebar=y title='xline: $*(m)'  >/dev/null out=tz-$*.v 
	tube $@.v tz-$*.v stretchy=y
