sp = ~/research.ln/Source/Xtalk/.W/.G/Test/Data
norne-test-wave1:
	Wavelet n1=1251 d1=0.004 wavelet=ricker2 domain=time tdelay=0.04 fund=25 fhigh=50 | /opt/SEP/SEP6.4.6/bin/Transf f_min1=3 f_max=60 f_max1=45 | Transp plane=35 > t$@.H
	Pad end1=119 extend=1 < t$@.H > $@.H; echo n1=120  o1=0  d1=75 n2=1  o2=0 d2=1   label2=shot-y >> $@.H 
	touch $@

norne-test-vel: norne-test-wave
	Window3d min3=2916 n3=1 < ../Velocities/psdm-vel.H | Interp d1out=5 > 1; < 1 Transp plane=13 | Transp | Pad beg1=100 end1=100 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@

norne-test-dat: 
	#Window3d min3=0. max3=9000 j3=1 < 2916-shot.H > t$@.H; echo o3=0 o2=100. >> t$@.H 
	Window3d min3=2000. n3=400  < ../sail-line-06-2d-rmo.H | Pad beg2=178 end2=50 | Median_me.x wind2=1 wind3=1 nodes=15 node=15  > t$@.H ; echo o3=2000 >> t$@.H 
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24 > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@
norne-test-tdat: 
	Window3d min3=2000. n3=400  < ../sail-line-06-2d-rmo.H | Pad beg2=178 end2=50 | Median_me.x wind2=1 wind3=1 nodes=15 node=15 | Tpow tpow=1.2 pow=1.2 > t$@.H ; echo o3=2000 >> t$@.H 
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24 > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@

norne-test-wave: 
	Wavelet n1=1251 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 | Transp plane=24 | Pad end1=339 extend=1 > $@.H; echo o1=2000 d1=25 >> $@.H

%-mig: %-vel %-dat %-wave
	$B/joint_encode.x vel=$*-vel.H data=$*-dat.H wavelet=$*-wave.H report=$@.log adj=y \
		modl1=$*-mig.H  data1=$*-dat.H vel1=$*-vel.H \
		jw=1 lloyd=y nsurv=1 code_it=n small=y ltscale=n use_greens=n spike=y wfld_d_out=n wfld_u_out=n encode=n code_it=n spatio=n offset=y\
		taper_pct=0.3 taper_type=2 use_greens=n verb=6 nodes=15 lfreq=y  \
		jw=1 spike=n radian=  tmin=-0.5 tmax=0.5 spike_data=n lloyd=y wfld_d_out=n radian=n nsurv=1 code_it=n  small=y ltscale=y use_greens=n wfld_d_out=n wfld_u_out=n spatio=y \
		ogx=100 dgx=25 ngx=127 
%-tmig: %-vel %-tdat
	$B/joint_encode.x vel=$*-vel.H data=$*-tdat.H wavelet=$*-wave.H report=$@.log adj=y \
		modl1=$*-tmig.H  data1=$*-tdat.H vel1=$*-vel.H \
		jw=1 lloyd=y nsurv=1 code_it=n small=y ltscale=n use_greens=n spike=y wfld_d_out=n wfld_u_out=n encode=n code_it=n spatio=n offset=y\
		taper_pct=0.25 taper_type=2 use_greens=n verb=6 nodes=15 lfreq=y  \
		jw=1 spike=n radian=  tmin=-0.5 tmax=0.5 spike_data=n lloyd=y wfld_d_out=n radian=n nsurv=1 code_it=n  small=y ltscale=y use_greens=n wfld_d_out=y wfld_u_out=n spatio=y \
		ogx=100 dgx=25 ngx=127
norne-08-2d-wave:
	Wavelet n1=1251 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 | Transp plane=24 | Pad end1=339 extend=1 > $@.H; echo o1=2000 d1=25 >> $@.H

norne-08-2d-vel: norne-test-wave
	Window3d min3=2916 n3=1 < ../Velocities/psdm-vel.H | Interp d1out=5 > 1; < 1 Transp plane=13 | Transp | Pad beg1=200 end1=200 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@
norne-08-2d-dat: 
	Window3d min3=2000. n3=400  < ../sail-line-08-2d-rmo.H | Pad beg2=178 end2=50 > t$@.H ; echo o3=2000 >> t$@.H 
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24 > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@
Norne-%-2d-vel: norne-test-wave
	Window3d min3=2759 n3=1 < ../Velocities/psdm-vel.H | Interp d1out=5 | Median_me.x wind1=2 wind2=2 > 1; < 1 Transp plane=13 | Transp | Pad beg1=200 end1=400 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@
Norne-%-2d-wave:
	/opt/SEP/SEP7.0/bin/Wavelet n1=1251 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 | Transp plane=24 | Pad end1=519 extend=1 > $@.H; echo o1=500 d1=25 >> $@.H
Norne-%-2d-dat: 
#	Create3d < ../sail-line-$*-2d-rmo.H keyname1=offset_x keyname2=s_x > junk1.$@
#	Headermath < junk1.$@ > junk2.$@ key1=g_x eqn1="s_x+offset_x" key2=off_x eqn2="-offset_x"
#	Sort3d  < junk2.$@  > junk3.$@ nkeys=2 copy_data=1 grid_data=1 verb=1 key1=off_x og1=0.0 dg1=25 ng1=132 key2=g_x og2=-3498 dg2=25 ng2=953
#	Infill3d < junk3.$@ > junk4.$@
	Window3d min3=500. n3=520       < ../sail-line-$*-2d-rmo.H  | Pad beg2=50      > t1$@.H ; echo o3=500 >> t1$@.H 
	Window3d min3=500. n3=520   f2=1< junk4.$@ 	            | Pad end2=50 > t2$@.H ; echo o3=500 >> t2$@.H
	Cat3d t1$@.H t2$@.H axis=2 > t$@.H
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24  > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@

NNorne-%-2d-dat: 
	Cp Norne-$*-2d-dat.H $@.H
NNorne-%-2d-wave:
	Cp Norne-$*-2d-wave.H $@.H
NNorne-%-2d-vel: norne-test-wave
	Window3d min3=2759 n3=1 < ../Velocities/psdm-vel.H | Interp d1out=5 | Median_me.x wind1=2 wind2=2 > 1; < 1 Transp plane=13 | Transp | Pad beg1=200 end1=400 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@

NORNE-%-2d-vel: norne-test-wave
	Window3d min3=2759 n3=1 < ../Velocities/PSDM-vel.H | Interp d1out=5 o2out=-5000 n2out=1101 n1out=991 |  Median_me.x wind1=2 wind2=2 | Smooth rect1=2 rect2=2 > 1; < 1 Transp plane=13 | Transp  Pad beg1=200 end1=200 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@
NORNE-%-2d-wave:
	/opt/SEP/SEP7.0/bin/Wavelet n1=1251 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 | Transp plane=24 | Pad end1=519 extend=1 > $@.H; echo o1=500 d1=25 >> $@.H
NORNE-%-2d-dat: 
	Create3d < ../sail-line-$*-2d-rmo.H keyname1=offset_x keyname2=s_x > junk1.$@
	Headermath < junk1.$@ > junk2.$@ key1=g_x eqn1="s_x+offset_x" key2=off_x eqn2="-offset_x"
	Sort3d  < junk2.$@  > junk3.$@ nkeys=2 copy_data=1 grid_data=1 verb=1 key1=off_x og1=0.0 dg1=25 ng1=132 key2=g_x og2=-3498 dg2=25 ng2=953
	Infill3d < junk3.$@ > junk4.$@
	Window3d min3=500. n3=520       < ../sail-line-$*-2d-rmo.H  | Pad beg2=50      > t1$@.H ; echo o3=500 >> t1$@.H 
	Window3d min3=500. n3=520   f2=1< junk4.$@ 	            | Pad end2=50 > t2$@.H ; echo o3=500 >> t2$@.H
	Cat3d t1$@.H t2$@.H axis=2 > t$@.H
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24  > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@

NNorne-%-2d-dat: 
	Cp Norne-$*-2d-dat.H $@.H

Norne-%-shift:
	$B/joint_encode.x report=$@.log modl=Norne-08-2d-mig.H \
	modl1=Norne-08-2d-mig.H  rm1=Norne-08-2d-mig.H  smodl1=Norne-08-2d-mig-s.H ts1=Norne-08-2d-mig-ts1.H sd1=$@-tsd1.H cc1=$@-cc1.H rr1=$@-rr1.H \
	modl2=Norne-06-2d-mig.H  rm2=Norne-06-2d-mig.H  smodl2=Norne-06-2d-mig-s.H ts2=Norne-06-2d-mig-ts2.H sd2=$@-tsd2.H cc2=$@-cc2.H rr2=$@-rr2.H \
        par=shift.p
#half-vel-rmo
NORNE-%-2d-h-vel: #norne-test-wave
	Window3d min3=2759 n3=1 < ../Velocities/PSDM-vel.H | Interp d1out=5 o2out=-5000 n2out=1101 n1out=991 |  Median_me.x wind1=2 wind2=2 | Smooth rect1=2 rect2=2 > 1; < 1 Transp plane=13 | Transp  Pad beg1=200 end1=200 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@
NORNE-%-2d-h-wave:
	/opt/SEP/SEP7.0/bin/Wavelet n1=1251 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 | Transp plane=24 | Pad end1=519 extend=1 > $@.H; echo o1=500 d1=25 >> $@.H
NORNE-%-2d-h-dat: 
	#Create3d < ../sail-line-$*-2d-half-rmo.H keyname1=offset_x keyname2=s_x > junk1.$@
	#Headermath < junk1.$@ > junk2.$@ key1=g_x eqn1="s_x+offset_x" key2=off_x eqn2="-offset_x"
	#Sort3d  < junk2.$@  > junk3.$@ nkeys=2 copy_data=1 grid_data=1 verb=1 key1=off_x og1=0.0 dg1=25 ng1=132 key2=g_x og2=-3498 dg2=25 ng2=953
	#Infill3d < junk3.$@ > junk4.$@
	#Window3d min3=500. n3=520       < ../sail-line-$*-2d-half-rmo.H  | Pad beg2=50      > t1$@.H ; echo o3=500 >> t1$@.H 
	#Window3d min3=500. n3=520   f2=1< junk4.$@ 	                 | Pad end2=50 > t2$@.H ; echo o3=500 >> t2$@.H
	Cat3d t1$@.H t2$@.H axis=2 > t$@.H
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24  > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@

NORNE-%-shift:
	Transp plane=35 < $(sp)/norne-08-$*-mig-h1.H > mod1.H
	Transp plane=35 < $(sp)/norne-06-$*-mig-h1.H > mod2.H
	$B/joint_encode.x report=$@.log modl=mod1.H  \
	modl1=mod1.H rm1=mod1.H  smodl1=norne-06-$*-mig-h1-s.H ts1=norne-06-mig-h1-ts.H sd1=$@-tsd1.H cc1=$@-cc1.H rr1=$@-rr1.H \
	modl2=mod2.H rm2=mod2.H  smodl2=norne-08-$*-mig-h1-s.H ts2=norne-08-mig-h1-ts.H sd2=$@-tsd2.H cc2=$@-cc2.H rr2=$@-rr2.H \
        nmix=5 par=shift.p nodes=1 shift_nodes=1 


NORNE-%-2d-3225-h-vel: #norne-test-wave
	Window3d min3=3225 n3=1 < ../Velocities/PSDM-vel.H | Interp d1out=5 o2out=-5000 n2out=1101 n1out=991 |  Median_me.x wind1=2 wind2=2 | Smooth rect1=1 rect2=1 > 1; < 1 Transp plane=13 | Transp  Pad beg1=200 end1=200 extend=1 > $@.H; echo o2=0  >> $@.H
	touch $@
NORNE-%-2d-3225-h-wave:
	/opt/SEP/SEP7.0/bin/Wavelet n1=1251 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 | Transp plane=24 | Pad end1=519 extend=1 > $@.H; echo o1=500 d1=25 >> $@.H
NORNE-%-2d-3225-h-dat: 
	Create3d < ../sail-line-$*-2d-3225-half-rmo.H keyname1=offset_x keyname2=s_x > junk1.$@
	Headermath < junk1.$@ > junk2.$@ key1=g_x eqn1="s_x+offset_x" key2=off_x eqn2="-offset_x"
	Sort3d  < junk2.$@  > junk3.$@ nkeys=2 copy_data=1 grid_data=1 verb=1 key1=off_x og1=0.0 dg1=25 ng1=132 key2=g_x og2=-3498 dg2=25 ng2=953
	Infill3d < junk3.$@ > junk4.$@
	Window3d min3=500. n3=520       < ../sail-line-$*-2d-3225-half-rmo.H  | Pad beg2=50 > t1$@.H ; echo o3=500 >> t1$@.H 
	Window3d min3=500. n3=520   f2=1< junk4.$@ 	                      | Pad end2=50 > t2$@.H ; echo o3=500 >> t2$@.H
	Cat3d t1$@.H t2$@.H axis=2 > t$@.H
	/opt/SEP/SEP6.4.6/bin/Transf < t$@.H  maxsize=100 f_min1=3 f_max1=45 f_max=60 > tt$@.H
	Transp plane=35  < tt$@.H | Transp plane=24  > $@.H; 
	echo label1=offset-x  label2=offset-y label3=freq. label4=shot-x label5=shot-y o5=0 d5=1 n5=1 o2=0 n2=1 d2=1  >> $@.H;
	touch $@
