############################################
#
# Migrating Gauss anomaly data
#
00marmenc: # 
#	Wavelet n1=951 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1.9 domain=time > j2.H;echo o1=-1.9 >> j2.H
#	<j2.H Transf f_min=0 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
#	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13 |Transp plane=24 > 11111.H
#	echo n4=1 d4=1 o4=0 >> 11111.H
#	< ../Tese/3000marmod_velsmoo.H Window3d squeeze=n | Pad beg3=8 end3=0 > 2222.H;echo o3=0 >> 2222.H
	< 2222.H ./bin/OW_encode.x sou=11111.H vel=../Tese/3000marm_vel_smooth1.tr.H \
	encrec=encrec1.H encsou=encsou1.H xmin=-8000 xmax=20000 ngap=20 > /dev/null

00marmig:
	< encsou1.H Window3d n4=1 squeeze=n > s.H
	< encrec1.H Window3d n4=1 squeeze=n > r.H
	../Tese/bin/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=marm_random_mig.H \
	zmin=0 zmax=3000 nh=1 61 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 nrefvel=4 > /dev/null	


00marmigperm:
	< ../../Thesis/ispew/3000marm.csou35.perm.sd.H Window3d n4=1 squeeze=n > s.H
	< ../../Thesis/ispew/3000marm.crec35.perm.sd.H Window3d n4=1 squeeze=n > r.H
	../Tese/src/areal/wemvaop/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=marm_perm_mig.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 icwind=1 deltat=0.016 dt=0.008 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigpermf:
	../Tese/src/areal/wemvaop/OWmig2d_areal.x sou=../../Thesis/ispew/3000marm.csou35.perm.sd.H rec=../../Thesis/ispew/3000marm.crec35.perm.sd.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=marm_perm_migf.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 icwind=1 deltat=0.016 dt=0.008 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigispew:
	< ../../Thesis/ispew/3000marm.csou20.ispew.sd.H Window3d n4=1 squeeze=n > s.H
	< ../../Thesis/ispew/3000marm.crec20.ispew.sd.H Window3d n4=1 squeeze=n > r.H
	../Tese/src/areal/wemvaop/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=marm_ispew_mig.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigispewf:
	../Tese/src/areal/wemvaop/OWmig2d_areal.x sou=../../Thesis/ispew/3000marm.csou20.ispew.sd.H rec=../../Thesis/ispew/3000marm.crec20.ispew.sd.H  \
	vel=../Tese/3000marm_vel_smooth1.H bkgimg=marm_ispew_migf.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00gauss_noenc: # 
	< ../Tese/1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=../Tese/1005gauss_noenc.H encrec=encrec1.H encsou=encsou1.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=200 > /dev/null

00gauss_encall%: # 
	< ../Tese/1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=../Tese/1005gauss_encall$*.H noencode=0 \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 random=1 encall=1 nreal=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

00v2000_noenc: # 
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_noenc.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=200 > /dev/null

00v2000_encall%:
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_encall$*.H noencode=0 encrec=encrec.H encsou=encsou.H \
	vsal=3 nrefvel=1 xmin=-5.1 xmax=9.1 random=1 encall=1 nreal=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

00v2000_goldall%: # 
	< ../Tese/1004v2000_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004v2000_source.H vel=../Tese/1000_v2000mig.H \
	image=../Tese/1005v2000_encall$*.H noencode=0 encrec=encrec.H encsou=encsou.H \
	vsal=3 nrefvel=1 xmin=-5.1 xmax=9.1 gold=1 encall=1 nreal=$* shift1=$* zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

01sig_prep.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d max1=8 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y \
	f_min=3 f_min1=5 f_max1=36 f_max=42 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@

02sig_sou1.H:
	Window3d n4=1 < 01sig_prep.H | Transp | Cabs | Stack | Smooth rect1=9 > specsmoo.H
	Math file1=specsmoo.H exp="file1*0" > phsesmoo.H; Cmplx specsmoo.H phsesmoo.H | Transp reshape=1,3,5 plane=23| \
	Transp reshape=1,3,5 plane=12 > $@ 
	< $@ Window3d | Transp| Transf |Graph|Tube&

01sig_prepw256.H:
	< ../OneWay/nfs-splitspread-shots.H Window3d max1=8 | Scale scale_to=1 | /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y \
	f_min=8 f_min1=12 f_max1=36 f_max=42 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Window3d squeeze=n n3=256 | Pad beg1=300 end1=300 > $@

02sig_sou1w256.H:
	Window3d n4=1 < 01sig_prepw256.H | Transp | Cabs | Stack | Smooth rect1=9 > specsmoo.H
	Math file1=specsmoo.H exp="file1*0" > phsesmoo.H; Cmplx specsmoo.H phsesmoo.H | Transp reshape=1,3,5 plane=23| \
	Transp reshape=1,3,5 plane=12 > $@ 
	< $@ Window3d | Transp| Transf |Graph|Tube&

00sig_velmig.H:
	<../2Wpar/1000velmig_my_areal.H Pad beg2=1000 end2=1000 extend=1 > $@

########################################
#
# Testing random phases
#
001sig_rand_%: 
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_rand_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

0015sig_rand_%: # new random function
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=0015sig_mig_rand_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

002sig_mseq_%: 
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=004sig_mig_mseq_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

003sig_golay_%: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=003sig_mig_golay_gap$* \
	nrefvel=4 xmin=-38000 xmax=70000 golay=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

004sig_gold_%: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=005sig_mig_gold_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=70000 gold=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* > /dev/null

05sig_mig_grid2: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< ../1Wpar/01sig_prep_grid2.H ./bin/OWmig2d_par.x sou=../1Wpar/02sig_sou1.H 02sig_sou1w256.H vel=../1Wpar/00sig_velmig_grid2.H 00sig_velmig_fix.H image=05sig_mig_grid2.H\
	encsou= 001rnd_encode_sou.H vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=194 400 > /dev/null

05sig_mig: # 
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=05sig_mig_noencode.H\
	vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 noencode=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=400 > /dev/null

05sig_mig_path: # 
	Window3d squeeze=n min4=32000 n4=50 n3=20 < 01sig_prepw256.H > z01sig_prep_select.H
	Window3d squeeze=n n3=20 < 02sig_sou1w256.H > z02sig_sou1w256.H
	< z01sig_prep_select.H ./bin/OWmig2d_par.x sou=z02sig_sou1w256.H vel=00sig_velmig_fix.H image=05sig_mig_noencode_path.H\
	vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 encode= 1 path=1 npath=3 5 wpath= 0.1 vrange= 0.05 zmin=0 zmax=24000 \
	nh=1 oh=0 dh=0 wbott=0 node= 1 > /dev/null

05sig_mig_nopath: # 
	Window3d squeeze=n min4=32000 n4=100 < 01sig_prepw256.H > z01sig_prep_select.H
#	Window3d squeeze=n min4=32000 n4=50 n3=20 < 01sig_prepw256.H > z01sig_prep_select.H
#	Window3d squeeze=n n3=20 < 02sig_sou1w256.H > z02sig_sou1w256.H
	< z01sig_prep_select.H ./bin/OWmig2d_par.x sou=z02sig_sou1w256.H vel=00sig_velmig_fix.H image=05sig_mig_noencode_nopath.H\
	vsal=14500 nrefvel=4 xmin=-38000 xmax=121000 encode= 1 path= 1 npath= 5 wpath= 0.1 vrange= 0.05 zmin=0 zmax=24000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

004tst: # one random code per shot
#	Window3d squeeze=n min4=12000 n4=400 < 01sig_prepw256.H > 01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=004tst.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=40 > /dev/null

z: 
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=z.H \
	nrefvel=4 xmin=-38000 xmax=121000 mseq=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=1 > /dev/null

#####################################
#
# Phase encoding all the shots 
#      in n realizations
000sig_rndall_tst_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=000sig_rndall_nre$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 encall=1 nreal=$* zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

000sig_goldall_tst_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=000sig_goldall_nre$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 gold=1 encall=1 nreal=$* shift=$* shift1=5 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

writephase_%:
#	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=wrt$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 random= 1 gold=1 encall=1 nreal=$* shift=$* shift1=5 zmin=0 zmax=24000 \
	nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 code=gold_all.H random_all.H > /dev/null

##################################
#
# Xiao's frequency independent
#       phase encoding
001sig_xiao_tst:
	Window3d squeeze=n min4=12000 n4=40 < 01sig_prepw256.H > z01sig_prep_select.H
	< z01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_xiao_tst.H \
	nrefvel=4 xmin=-38000 xmax=121000 xiao=1 random=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=4 node= 6 > /dev/null

001sig_xiao_%:
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=001sig_mig_xiao_gap$*.H \
	nrefvel=4 xmin=-38000 xmax=121000 xiao=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 ngap=$* node=6 > /dev/null


###################################
#
# Test of fftw + omp
#
tst_omp:
	Window3d squeeze=n n4=4 < ../Tese/1004gauss_prep.H > z1004gauss_prep.H
	< z1004gauss_prep.H ./bin/OWmig2d_par.x sou=../Tese/1004gauss_source.H vel=../Tese/1003gaussvel_.02pad.H \
	image=1.H \
	vsal=3 nrefvel=4 xmin=-5.1 xmax=9.1 noencode=1 zmin=0 zmax=1.15 nh=1 oh=0 dh=0 wbott=0 ngap=4 > /dev/null

tst_omp_sig:
	Window3d squeeze=n min4=12000 max4=70000 j4=80 < 01sig_prepw256.H > tst_omp_sig.H
	< tst_omp_sig.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=1.H \
	nrefvel=4 xmin=-38000 xmax=121000 noencode=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_rnd_omp_sig:
#	Window3d squeeze=n min4=12000 n4=1 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=2.H encrec=encrec.H encsou=encsou.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 encall=1 nreal=1 node= 1 eps=.25 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_rnd_omp_sig1:
#	Window3d squeeze=n min4=12000 n4=4 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=3.H encrec=encrec.H encsou=encsou.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 encall=0 nreal=1 node= 1 eps=1 nshots=10 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_rnd_omp_sig4:
#	Window3d squeeze=n min4=12000 n4=4 < 01sig_prepw256.H > z01sig_prep_select.H
	< 01sig_prep_select.H ./bin/OWmig2d_par.x sou=02sig_sou1w256.H vel=00sig_velmig_fix.H image=4.H encrec= encrec.H encsou= encsou.H \
	nrefvel=4 xmin=-38000 xmax=121000 random=1 encall=0 nreal=1 node= 1 eps=1 nshots=1 zmin=0 zmax=24000 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_marm:
#	Window3d squeeze=n min4=12000 n4=4 < 01sig_prepw256.H > z01sig_prep_select.H
	< ../2Wpar/9001marm_prep.H ./bin/OWmig2d_par.x sou=../2Wpar/9002marm_source.H vel=../2Wpar/9000marmvel_25m.H \
	image=1001marm_shtprof.H nrefvel=4 xmin=-8000 xmax=20000 noencode=0  node= 1 eps=1 zmin=0 zmax=3000 vsal=5500 \
	nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_marm4:
#	Window3d squeeze=n min4=12000 n4=4 < 01sig_prepw256.H > z01sig_prep_select.H
	< ../2Wpar/9001marm_prep.H ./bin/OWmig2d_par.x sou=../2Wpar/9002marm_source.H vel=../2Wpar/9000marmvel_25m.H image=5.H encrec= encrec.H encsou= encsou.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node= 1 eps=1 nshots=4 zmin=0 zmax=3000 vsal=5500 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

tst_marm5:
#	Window3d squeeze=n min4=12000 n4=4 < 01sig_prepw256.H > z01sig_prep_select.H
	< ../2Wpar/9001marm_prep.H ./bin/OWmig2d_par.x sou=../2Wpar/9002marm_source.H vel=../2Wpar/9000marmvel_25m.H image=6.H encrec= encrec.H encsou= encsou.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node= 1 eps=1 nshots=4 zmin=0 zmax=3000 vsal=5500 nh=1 oh=0 dh=0 wbott=0 6000 > /dev/null

###############################################
#
#
#
1000marmvel_25_nopad.H:
	Interp d1out=4 d2out=25 < /net/koko/claudio/Research/marmvel.H | Transp plane=12 reshape=1,3 > $@

1000marmvel_25.H:
	Interp d1out=4 o2out=-10000 d2out=25 n2out=1201 < ../2Wpar/marmvel.H > $@

1000marmvel_25brad.H:
	Transp plane=12 reshape=1,3 < 1000marmvel_25.H > 1000marmvel_25brad.H 

1000marmvel_25_smooth.H:
	Math file1=1000marmvel_25.H exp="1/file1" | Smooth rect1=5 rect2=5 > x1.H ; Math file1=x1.H exp="1/file1" >$@

1000marmvel_25_lloyd.H: # it is slowness !!!!!!
	Lloyd_vel < 1000marmvel_25_nopad.H nref=4 min_region_pct=2 niter_lloyd=120 min_slow_dev=.01 > $@

1001marm_prep.H:
	< ../2Wpar/marmdata.H Scale scale_to=1. | Pad end1=360 | \
	/net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=1 f_min1=4 f_max1=50 f_max=60 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=300 end1=300 > $@ ; echo d1=25 o1=-10075 >> $@

1002marm_source.H:
	< 1001marm_prep.H Window3d j4=10 | Transp plane=12 | Transp plane=13 > xxx.H
	 < xxx.H Cabs | Stack | Transp | Stack | Smooth rect1=15 > 9002marm_spec.H;
	Math file1=9002marm_spec.H exp="file1*0" > 9002marm_phse.H; Cmplx 9002marm_spec.H 9002marm_phse.H | Transp reshape=1,3,5 plane=23| \
	Transp reshape=1,3,5 plane=12 > $@

1003marm_shtpmig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_shtprof.H vref=1000marmvel_25_lloyd.H \
	nrefvel=4 xmin=-8000 xmax=20000 noencode=0  node= 1 eps=1 zmin=0 zmax=3000 vsal=5500 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_shtpmig_lloyd:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_shtprof_lloyd.H \
	vref=1000marmvel_25_lloyd.H nrefvel=4 xmin=-8000 xmax=20000 encode=0 node=16 eps=1 zmin=0 zmax=3000 vsal=5500 \
	nh=1 33 oh=0 -400 dh=25 wbott=0 > /dev/null

1003marm_shtpmig_lloyd_1:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_shtprof_lloyd_img1.H \
	vref=1000marmvel_25_lloyd.H nrefvel=4 xmin=-8000 xmax=20000 noencode=0  node=6 eps=1 zmin=0 zmax=3000 vsal=5500 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1003marm_renc4mig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_renc4mig.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node= 1 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_renc4mig1: # with a change in the way I define the reference velocities
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_renc4mig_vref2.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node=12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_renc4mig2: # with a change in the way I define the reference velocities and phase in +-pi/2
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_renc4mig_vref2_pi2.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node=12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_renc4mig16: # with a change in the way I define the reference velocities and phase in +-pi/2
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_renc16mig_vref2_pi2.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node=12 eps=1 nshots=16 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_renc4mig8: # with a change in the way I define the reference velocities and phase in +-pi/2
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_renc8mig_vref2_pi2.H \
	nrefvel=4 xmin=-8000 xmax=20000 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_genc4mig1: # with a change in the way I define the reference velocities
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_genc4mig_vref2.H \
	nrefvel=4 xmin=-8000 xmax=20000 gold=1 encall=0 nreal=1 node=12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_genc4mig16: # with a change in the way I define the reference velocities
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_genc16mig_vref2.H \
	nrefvel=4 xmin=-8000 xmax=20000 gold=1 encall=0 nreal=1 node=12 eps=1 nshots=16 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_genc4mig8: # with a change in the way I define the reference velocities
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1003marm_genc8mig_vref2.H \
	nrefvel=4 xmin=-8000 xmax=20000 gold=1 encall=0 nreal=1 node=12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1003marm_genc4mig2: # with a change in the way I define the reference velocities and smooth velocity !!! Bad !!!
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25_smooth.H image=1003marm_genc4mig_vref2_smoo.H \
	nrefvel=4 xmin=-8000 xmax=20000 gold=1 encall=0 nreal=1 node=12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1004marm_shtpmig_brad.H:
	@date
	nohup /net/koko/claudio/Research/Others/ShotMig/ShotMig \
	D=1001marm_prep.H D_d=1002marm_source.H V=1000marmvel_25brad.H R=$@ \
	pverb=2  nx_shots=240 ox_shots=3000 \
	par=marm_go_sis.P taper=y taper_pcent=2 marine=y minx=-8000 maxx=20000 & #restart=y

########################################
#
# Phase encoding-- now using Lloyd
#
1005marm_renc4mig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_renc4mig_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_renc4mig1:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_renc4mig1_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_renc8mig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_renc8mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1005marm_renc8mig2:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_renc8mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_renc16mig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_renc16mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node= 12 eps=1 nshots=16 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1005marm_rencallmig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_rencallmig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=1 nreal=1 60 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 node=10 > /dev/null

1005marm_genc4mig: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc4mig1: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig1_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc4mig1a: # after changing the computation of preferred pairs
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig1a_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node=12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc4mig2: # combining with no encoding 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig2_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=0 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc4migpi: # ver como fica com +-pi
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig_ll_pi.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc8migpi: # ver como fica com +-pi
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc8mig_ll_pi.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc8mig: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc8mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1005marm_genc8mig2: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc8mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc8miglin: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc8miglin_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 16 eps=1 nshots=8 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 linint=1 > /dev/null

1005marm_genc16mig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc16mig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=16 zmin=0 zmax=3000 vsal=4000 \
	nh=1 oh=0 dh=0 wbott=0 > /dev/null

1005marm_gencallmig:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_gencallmig_ll.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=1 nreal=1 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_gencallmigperf:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_gencallmig_ll_perf.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 perfect=1 encall=1 nreal=1 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_genc4migperf_new:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig_ll_perf_new.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 perfect=1 encall=0 nreal=1 nshots=4 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 jshift=1 > /dev/null

1005marm_gencallmigperf_new: # Cp 1005marm_gencallmig_ll_perf_new.H 1005marm_genc4mig_ll_perf_new.H
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_gencallmig_ll_perf_new.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 perfect=1 encall=1 nreal=1 10 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 node=10 > /dev/null

#1005marm_gencallmig:
#	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_gencallmig_ll.H \
#	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=1 nreal=60 eps=1 zmin=0 zmax=3000 \
#	nh=1 oh=0 dh=0 wbott=0 > /dev/null
################################
################################
#
# Biondo's simulated annealing random codes
#
1005marm_anneal4mig: 
	<~biondo/Code/Anneal/GenAnneal/Test/Claudio/Anneal_Corr_indfreq_wide.H Window3d | Transp > Anneal_Corr_indfreq_wide.H
	Window3d n2=40 <Anneal_Corr_indfreq_wide.H > Anneal_Corr_indfreq_wide1.H;echo o2=201 >> Anneal_Corr_indfreq_wide1.H
	Cat Anneal_Corr_indfreq_wide.H Anneal_Corr_indfreq_wide1.H axis=2 > Anneal_Corr_indfreq_wide240.H
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_anneal4mig_ll.H InCode=Anneal_Corr_indfreq_wide240.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 random=1 encall=0 nreal=1 node=8 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_anneal4migall:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_optmig_ll_all.H InCode=Anneal_Corr_indfreq_wide240.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=1 nreal=1 eps=1 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_noenc_migall:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_noenc_mig_ll_all.H InCode=Anneal_Corr_indfreq_wide240.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 random=1 encall=1 nreal=1 eps=0 zmin=0 zmax=3000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

#
#
#
################################
################################
#
# Gold codes
#
1005marm_genc4mig_perf: 
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_genc4mig_perf_ll.H \
	vref=1000marmvel_25_lloyd.H  xmin=-8000 xmax=20000 encode=1 gold=1 perfect=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

#
#
#
################################


################################
#
# Writing codes
#
codegold:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=tstcode.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nshots=8 eps=1 zmin=0 zmax=3000 \
	code=goldcode.H nh=1 oh=0 dh=0 wbott=0 > /dev/null


coderand:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=tstcode.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nshots=8 eps=1 zmin=0 zmax=3000 \
	code=randcode.H nh=1 oh=0 dh=0 wbott=0 > /dev/null

randcovar.H:
	Window3d <randcode.H squeeze=n n3=1 | Transp plane=12 reshape=1,3 > xxx.H
	/net/koko/claudio/Research/Others/Biondo/PhaseCorr_v2.x < xxx.H par=corr2.P Covar_mat=randcovar_mat.H > $@

goldcovar.H:
	Window3d < goldcode.H squeeze=n n3=1 | Transp plane=12 reshape=1,3 > xxx.H
	/net/koko/claudio/Research/Others/Biondo/PhaseCorr_v2.x < xxx.H par=corr2.P Covar_mat=goldcovar_mat.H > $@

#
# Finished writing codes
#
################################
################################
#
# Angle gathers
#
1006marm_shtprofang_ll.H:
	Window3d < 1003marm_shtprof_lloyd.H min1=3000 max1=8600 | Transp plane=13 | Transp plane=34 | Scale scale_to=1 |\
	Pad beg2=44 end2=44 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=5 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=3000 > $@

1006marm_renc4ang_ll.H:
	Window3d < 1005marm_renc4mig_ll.H min1=3000 max1=8600 | Transp plane=13 | Transp plane=34 | Scale scale_to=1 |\
	Pad beg2=44 end2=44 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=5 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=3000 > $@

1006marm_genc4ang_ll.H:
	Window3d < 1005marm_genc4mig_ll_pi.H min1=3000 max1=8600 | Transp plane=13 | Transp plane=34 | Scale scale_to=1 |\
	Pad beg2=44 end2=44 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=5 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=3000 > $@

1006marm_renc8ang_ll.H:
	Window3d < 1005marm_renc8mig_ll.H min1=3000 max1=8600 | Transp plane=13 | Transp plane=34 | Scale scale_to=1 |\
	Pad beg2=44 end2=44 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=5 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=3000 > $@

1006marm_genc8ang_ll.H:
	Window3d < 1005marm_genc8mig_ll_pi.H min1=3000 max1=8600 | Transp plane=13 | Transp plane=34 | Scale scale_to=1 |\
	Pad beg2=44 end2=44 end1=100 > xx.H
	< xx.H /net/koko/bob/SEP_BOB/bin/OFF2ANG amax=40 eps=5 maxmem=1000 verb=y > yy.H; <yy.H Window3d max1=3000 > $@

###################
#
# tests
1005marm_tst:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_tst.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_tst1:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_tst1.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=1 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null

1005marm_tst2:
	< 1001marm_prep.H ./bin/OWmig2d_par.x sou=1002marm_source.H vel=1000marmvel_25.H image=1005marm_tst2.H \
	vref=1000marmvel_25_lloyd.H xmin=-8000 xmax=20000 encode=1 gold=1 encall=0 nreal=1 node= 12 eps=.75 nshots=4 zmin=0 zmax=3000 vsal=4000 \
	nh=33 oh=-400 dh=25 wbott=0 > /dev/null
#######################
#
# Figures
shtprf.v: # 1003marm_shtprof_lloyd.H
	Window3d < 1003marm_shtprof_lloyd.H min1=3000 max1=8600 | Transp plane=13 | Scale scale_to=1 | \
	Grey gainpanel=a pclip=95 eout=1 label1="depth(m)" label2="offset(m)" label3="distance(m)" title="shot profile" | Transp plane=23 > xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=750 n2=224 n3=33 n4=1 n5=1  >a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=228 frame2=180 frame3=16 title=" " > /dev/null out=./Fig/$@

renc4.v: # 1003marm_shtprof_lloyd.H
	Window3d < 1003marm_shtprof_lloyd.H min1=3000 max1=8600 | Transp plane=13 | Scale scale_to=1 | \
	Grey gainpanel=a pclip=95 eout=1 label1="depth(m)" label2="offset(m)" label3="distance(m)" title="shot profile" | Transp plane=23 > xxx.T
	< xxx.T Window3d  f1=0 f2=0 f3=0 f4=0 f5=0 n1=750 n2=224 n3=33 n4=1 n5=1  >a0.T
	< a0.T Cubeplot wantcoordlabel=-1 point1=.7 point2=0.7 flat=1 frame1=228 frame2=180 frame3=16 title=" " > /dev/null out=./Fig/$@