Prst_ExpRefl=/net/koko/claudio/Research/Prst_ExpRefl

3000marm_vel_smooth2W:
	Math file1=${Prst_ExpRefl}/2Wpar/8000marmvel_24m.H exp="1/file1" > j1.H; < j1.H Window3d min2=0 max2=9000 > j2.H
	Math file1=j2.H exp="1/file1" | Transp plane=12 reshape=1,3 | Pad beg1=833 2000 end1=793 2000 extend=1 |\
	Window3d|Transp|Mute tmute=24 vmute=0 tramp=0| Clip clip=1500 chop=less to=1500 | Transp plane=12 reshape=1,3> j.H
	Math file1=j.H exp="1/@ABS(file1)" > 3000marm_vel_2W.H
	Math file1=${Prst_ExpRefl}/2Wpar/8000marmvel_24m.H exp="1/file1" > j1.H; < j1.H Window3d min2=0 max2=9000 | \
	Smooth rect1=48 rect2=16 > j2.H
	Math file1=j2.H exp="1/file1" | Transp plane=12 reshape=1,3 > 3000marm_vel_smooth2.H
	Math file1=3000marm_vel_2W.H exp="1/@ABS(file1)" | Window3d squeeze=n min1=0 max1=9000 |\
	${Prst_ExpRefl}/Tese/bin/OW_editvel_mut.x vvel=1 vsub=3000marm_vel_smooth2.H hor1=marmv0.orig > j1.H
	< j1.H /Tese/bin/OW_editvel_mut.x vper=1 vperc=.9 hor1=marmv0.orig | Pad beg1=833 2000 end1=793 2000 extend=1 |\
	Window3d|Transp|Mute tmute=24 vmute=0 tramp=0 | Clip clip=1500 chop=less to=1500 | Transp plane=12 reshape=1,3> j.H
	Math file1=j.H exp="1/@ABS(file1)" > 3000marm_vel_smooth.2W.H

########################################################################
#
# Perm w/2way data now w/FFD
#
8001marm_shot_fd%:
	Fdmod intag=${Prst_ExpRefl}/2Wpar/8000marmvel_4m.H  \
	no_stdout=1 oxs=$* dxs=24 nxs=25 tmax=4 ex=1650 bx=1650 jt=10 hsfile=${Prst_ExpRefl}/2Wpar/temp$*.H > /dev/null

8001marm_shots_fd.H:
	Cat ${Prst_ExpRefl}/2Wpar/temp0000.H \
	${Prst_ExpRefl}/2Wpar/temp0600.H \
	${Prst_ExpRefl}/2Wpar/temp1200.H \
	${Prst_ExpRefl}/2Wpar/temp1800.H \
	${Prst_ExpRefl}/2Wpar/temp2400.H \
	${Prst_ExpRefl}/2Wpar/temp3000.H \
	${Prst_ExpRefl}/2Wpar/temp3600.H \
	${Prst_ExpRefl}/2Wpar/temp4200.H \
	${Prst_ExpRefl}/2Wpar/temp4800.H \
	${Prst_ExpRefl}/2Wpar/temp5400.H \
	${Prst_ExpRefl}/2Wpar/temp6000.H \
	${Prst_ExpRefl}/2Wpar/temp6600.H \
	${Prst_ExpRefl}/2Wpar/temp7200.H \
	${Prst_ExpRefl}/2Wpar/temp7800.H \
	${Prst_ExpRefl}/2Wpar/temp8400.H axis=3 | \
	Bandpass flo=4 | Window3d j2=6 | Interp d1out=0.004 |Mute vmute=1500 tmute=0.05  > ${Prst_ExpRefl}/2Wpar/$@

3000marm_data_prep.H:
	< ${Prst_ExpRefl}/2Wpar/8001marm_shots_fd.H Window3d n1=951 | Transf shotmig=y is2d=y wei=y f_min=0 f_min1=4 f_max1=40 f_max=48 maxsize=1000 > $@

3000marmig2bckgffd:
	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=40 f_max=48 wei=y is2d=y shotmig=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 > 11111_2W.H
	Window3d < 3000marm2W.inv_vout11.perm.sd.H n4=1 f4=31 squeeze=n |\
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=11111_2W.H rec=3000marm_data_prep.H  \
	vel=3000marm_vel_smooth.2W.H bkgimg=3000marm.bkim2Wa.ffd.H \
	zmin=0 zmax=3000 nh=31 oh=-360 jh=1 migration=1 areal=0 \
	padtraces=50 traces_in_taper=50 node=14 \
	extrapolator=ffd

3000marm2Wffd.window.H:
	< 3000marm.bkim2Wa.ffd.H Window3d min1=0 max1=9000 min2=-192 n2=17 > 1.H
	< 1.H ./bin/OW_windowing.x nhorz=12 hor1=marm_ffd_h1 hor2=marm_ffd_h2 hor3=marm_ffd_h3 hor4=marm_ffd_h4 \
	hor5=marm_ffd_h5 hor6=marm_ffd_h6 hor7=marm_ffd_h7 hor8=marm_ffd_h8 hor9=marm_ffd_h9 hor10=marm_ffd_h10 \
	hor11=marm_ffd_h11 hor12=marm_ffd_h12 wavz=13 taperon=1 > $@

3000marm2W_rotsdffd:
	Window3d <3000marm2Wffd.window.H n1=70 f1=170 f3=200 n4=1 f4=0  |Transp plane=13|\
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=8 nxchunk=1 | Pad beg3=170 end3=135 beg1=200 > h1.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=70 f1=170 f3=200 n4=1 f4=0  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=8 nxchunk=1 | Pad beg3=170 end3=135 beg1=200 > h1.rec.sd.H
##
	Window3d <3000marm2Wffd.window.H n1=85 f1=290 f3=200 n4=1 f4=1  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=290 end3=0 beg1=200 > h2.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=85 f1=290 f3=200 n4=1 f4=1  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=290 end3=0 beg1=200 > h2.rec.sd.H
##
	Window3d <3000marm2Wffd.window.H n1=130 f1=170 f3=200 n4=1 f4=2  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=170 end3=75 beg1=200 > h3.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=130 f1=170 f3=200 n4=1 f4=2  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=170 end3=75 beg1=200 > h3.rec.sd.H
##
	Window3d <3000marm2Wffd.window.H n1=60 f1=120 f3=200 n4=1 f4=3  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=120 end3=195 beg1=200 > h4.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=60 f1=120 f3=200 n4=1 f4=3  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=120 end3=195 beg1=200 > h4.rec.sd.H
##
	Window3d <3000marm2Wffd.window.H n1=60 f1=90 f3=200 n4=1 f4=4  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=90 end3=225 beg1=200 > h5.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=60 f1=90 f3=200 n4=1 f4=4  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=90 end3=225 beg1=200 > h5.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=140 f1=150 f3=200 n4=1 f4=5  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=150 end3=85 beg1=200 > h6.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=140 f1=150 f3=200 n4=1 f4=5  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=150 end3=85 beg1=200 > h6.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=80 f1=200 f3=200 n4=1 f4=6  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=200 end3=95 beg1=200 > h7.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=80 f1=200 f3=200 n4=1 f4=6  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=200 end3=95 beg1=200 > h7.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=110 f1=220 f3=200 n4=1 f4=7  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=220 end3=45 beg1=200 > h8.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=110 f1=220 f3=200 n4=1 f4=7  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=220 end3=45 beg1=200 > h8.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=150 f1=180 f3=200 n4=1 f4=8  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=180 end3=45 beg1=200 > h9.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=150 f1=180 f3=200 n4=1 f4=8  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=180 end3=45 beg1=200 > h9.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=210 f1=0 f3=200 n4=1 f4=9  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=0 end3=165 beg1=200 > h10.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=210 f1=0 f3=200 n4=1 f4=9  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=0 end3=165 beg1=200 > h10.rec.sd.H
#
	Window3d <3000marm2Wffd.window.H n1=145 f1=230 f3=200 n4=1 f4=10  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=230 end3=0 beg1=200 > h11.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=145 f1=230 f3=200 n4=1 f4=10  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=230 end3=0 beg1=200 > h11.rec.sd.H

	Window3d <3000marm2Wffd.window.H n1=80 f1=150 f3=200 n4=1 f4=11  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 nodes=14 nxchunk=1 | Pad beg3=150 end3=145 beg1=200 > h12.sou.sd.H
	Window3d <3000marm2Wffd.window.H n1=80 f1=150 f3=200 n4=1 f4=11  |Transp plane=13 | \
	./bin/OW_3dslant.x nax=41 nah=31 oax=-40 oah=-30 dax=2 dah=2 \
	nxw=31 rot=1 upshft=0 rec=1 nodes=14 nxchunk=1 | Pad beg3=150 end3=145 beg1=200 > h12.rec.sd.H

	Cat h1.sou.sd.H h2.sou.sd.H h3.sou.sd.H h4.sou.sd.H h5.sou.sd.H h6.sou.sd.H h7.sou.sd.H \
	h8.sou.sd.H h9.sou.sd.H h10.sou.sd.H h11.sou.sd.H h12.sou.sd.H axis=4 >3000marm2Wffd.rot.sou.sd.H
	Cat h1.rec.sd.H h2.rec.sd.H h3.rec.sd.H h4.rec.sd.H h5.rec.sd.H h6.rec.sd.H h7.rec.sd.H \
	h8.rec.sd.H h9.rec.sd.H h10.rec.sd.H h11.rec.sd.H h12.rec.sd.H axis=4 >3000marm2Wffd.rot.rec.sd.H
	Rm h*.rec.sd.H h*.sou.sd.H 1.H

3000marm2Wffd_perm%:
	Window3d < 3000marm2Wffd.rot.sou.sd.H squeeze=n n4= 1 f4= 5 |Transp plane=13 > s1.H
	Window3d < 3000marm2Wffd.rot.rec.sd.H squeeze=n n4= 1 f4= 5 |Transp plane=13 > r1.H
	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=40 f_max=48 wei=y is2d=y shotmig=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 > 11111.H
	Math file1=3000marm_vel_smooth.2W.H exp="1/file1" > 3000marm_vel_smooth.2W.tr.H
	./bin/OW_arealsht.x rec_cig=r1.H sou_cig=s1.H sou=11111.H vel=3000marm_vel_smooth.2W.tr.H \
	rec_out=3000marm2Wffd.crec$*.perm.H sou_out=3000marm2Wffd.csou$*.perm.H extrapolator=ffd pspi nrefvel=4 vsal=7000 \
	padtraces=101 encode=1 encz=1 eps=1 random=1 injection=1 nreal=1 ngap=$* zmin=0 zmax=3000 node=14


3000marm2W02ffdperm%: # 
	Math file1=3000marm2Wffd.crec$*.perm.H exp="file1/100"|Pad beg1=267 end1=267>3000marm2Wffd.crec.perm.pad.H
	Math file1=3000marm2Wffd.csou$*.perm.H exp="file1/100"|Pad beg1=267 end1=267>3000marm2Wffd.csou.perm.pad.H
	${Prst_ExpRefl}/Tese/bin/OW_lbfgsc_areal.x sou=3000marm2Wffd.csou.perm.pad.H rec=3000marm2Wffd.crec.perm.pad.H \
	vel=3000marm_vel_smooth.2W.H grad=3000marm2W02ffd.inv_grad$*.perm.H \
	Dimage=3000marm2W02ffd.inv_dso$*.perm.H bkgimg=3000marm2W02ffd.inv_bkim$*.perm.H \
	Velout=3000marm2W02ffd.inv_vout$*.perm.H \
	zmin=0 zmax=3000 nh=17 oh=-192 jh=1 niter=15 dso=1 inversion=1 areal=1 mmax=10 xv_min=500 xv_max=8500 \
	dxk=240 nxk=109 oxk=-8400 dzk=16 nzk=188 horv0=marmv0.out0 \
	dipflt= 1 nang= 161 dang= 0.01 oang= -.8 indip= dip.smo9.H destroy= n \
	extrapolator=ffd padtraces=50 traces_in_taper=25 node=14 bounds=1 vmin=2000 vmax=5200 \
	log=log_bfgs_marm2W02ffd.perm$* list=listagem_bfgs_marm2W02ffd.perm$* 

3000marmig2W02vout11ffd:
	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=40 f_max=48 wei=y is2d=y shotmig=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 > 11111_2W.H
	Window3d < ${Prst_ExpRefl}/Tese/3000marm2W02ffd.inv_vout11.perm.H n4=1 f4=40 squeeze=n |\
	Pad extend=1 beg1=1000 end1=1000|Gpow gpow=-1 > 3000marm2W02ffd.inv_vout11_41.H
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=11111_2W.H rec=3000marm_data_prep.H  \
	vel=3000marm2W02ffd.inv_vout11_41.H bkgimg=3000marm.bkim2W02ffd.vout11_41.H \
	zmin=0 zmax=3000 nh=31 oh=-360 jh=1 migration=1 areal=0 padtraces=50 traces_in_taper=50 node=14 \
	extrapolator=ffd

3000marmig2Wtrueffd:
	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=40 f_max=48 wei=y is2d=y shotmig=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 > 11111_2W.H
	Window3d < ${Prst_ExpRefl}/2Wpar/8000marmvel_24m.H Transp plane=12 reshape=1,3|Pad extend=1 beg1=1000 end1=1000|Gpow gpow=-1 > 3000marm2W.true.H
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=11111_2W.H rec=3000marm_data_prep.H  \
	vel=3000marm2W.true.H bkgimg=3000marm.bkim2W02ffd.true.H \
	zmin=0 zmax=3000 nh=31 oh=-360 jh=1 migration=1 areal=0 padtraces=50 traces_in_taper=50 node=14 \
	extrapolator=ffd

3000marm.ang2W.vtrue.ffd.H:
	Window3d min1=0 max1=9000 min2=-192 n2=17 < 3000marm.bkim2W02ffd.true.H | Transp plane=13|Transp plane=34|Pad beg2=300 end2=300|OFF2ANG amax=40 eps=10 > x.H
	<x.H Window3d > $@ 

3000marm.ang2W.vsmoo2W.ffd.H:
	Window3d min1=0 max1=9000 min2=-192 n2=17 < 3000marm.bkim2W.vsmoo2W.ffd.H | Transp plane=13|Transp plane=34|Pad beg2=300 end2=300|OFF2ANG amax=40 eps=10 > x.H
	<x.H Window3d > $@ 

3000marm.ang2W02.ffd.H:
	Window3d min1=0 max1=9000 min2=-192 n2=17 <  3000marm.bkim2W02ffd.vout11_41.H| Transp plane=13|Transp plane=34|Pad beg2=300 end2=300|OFF2ANG amax=40 eps=10 > x.H
	<x.H Window3d > $@ 

# Figures

tslow.v:
	Window3d < ${Prst_ExpRefl}/2Wpar/8000marmvel_24m.H min2=1000 max2=8000 |Gpow gpow=-1 |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(a)" label1="z" barlabel= "slowness" \
	label2="x" screenratio= screenht= 4 screenwd= 5 wheretitle=b polarity= bartype=h \
	wantaxis1=y wantaxis2=y labelsz=6 titlesz=6 wantscalebar=1 barwidth=0.15 > /dev/null out=${R}/$@

tkim.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.bkim2W02ffd.true.H min2=0 n2=1 min1=1000 max1=8000 n4=1 min3= 1500 | Transp |\
	Grey pclip=95 title="(b)" label1="z" label2=" " screenratio= 4 screenht= 2 screenwd= 4 \
	wantaxis1=y wantaxis2=n labelsz=6 titlesz=6 wheretitle=b gainpanel=a > /dev/null out=${R}/$@

true.v: tslow.v tkim.v
	vppen < ${R}/tslow.v xcenter=0 ycenter=-9. vpstyle=n xsize=4 > junk1.v
	vppen < ${R}/tkim.v xcenter=0 ycenter=-0. vpstyle=n xsize=4 > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@
	
bslow.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm_vel_smooth.2W.H min1=1000 max1=8000 | Transp |\
	Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(a)" label1="z" barlabel= "slowness" \
	label2="x" screenratio= 4 screenht= 3.5 screenwd= 5 wheretitle=b polarity= bartype=h \
	wantaxis1=y wantaxis2=y labelsz=6 titlesz=6 wantscalebar=1 barwidth=0.15 > /dev/null out=${R}/$@

bkim.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.bkim2Wa.ffd.H min2=0 n2=1 min1=1000 max1=8000 n4=1 min3= 1500 | Transp |\
	Grey pclip=95 title="(b)" label1="z" label2="x" screenratio= 4 screenht= 3.5 screenwd= 4 \
	wantaxis1=y wantaxis2=n labelsz=6 titlesz=6 wheretitle=b gainpanel=a > /dev/null out=${R}/$@

bckg.v: bslow.v bkim.v
	vppen < ${R}/bslow.v xcenter=0 ycenter=-9. vpstyle=n xsize=4 > junk1.v
	vppen < ${R}/bkim.v xcenter=0 ycenter=-0. vpstyle=n xsize=4 > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

dslow.v:
	Window3d < ${Prst_ExpRefl}/2Wpar/8000marmvel_24m.H min2=1000 max2=8000 min1= 1500 |Gpow gpow=-1 > x1.H
	Window3d < ${Prst_ExpRefl}/Tese/3000marm_vel_smooth.2W.H min1=1000 max1=8000 min3= 1500 |Transp > x2.H
	Math file1=x1.H file2=x2.H exp="file1/file2" |\
	Grey newclip=1 color=j bclip=0.7 eclip=1.3 bias=1. allpos=1 title=" " label1="z" barlabel=" " \
	label2="x" screenratio= 4 screenht= 3.5 screenwd= 5 wheretitle=b bartype=h \
	wantaxis1=y wantaxis2=y labelsz=6 titlesz=6 wantscalebar=1 barwidth=0.15 > /dev/null out=${R}/$@
	
bwkim.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm2Wffd.window.H min2=0 n2=1 min1=1000 max1=8000 min3= 1500 | \
	Transp plane=23|Stack|Transp| \
	Grey gainpanel=a title=" " label1="z" label2="x" screenratio= 4 screenht= 3.5 screenwd= 4 pclip=100 \
	wantaxis1=y wantaxis2=y labelsz=6 titlesz=6 wheretitle=b > /dev/null out=${R}/$@

oslow.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm2W02ffd.inv_vout11.perm.H n4=1 f4=40 min1=1000 max1=8000 | Transp |\
	Gpow gpow=-1|Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title="(a)" label1="z" barlabel=" " \
	label2="x" screenratio= 4 screenht= 3.5 screenwd= 5 wheretitle=b polarity= -1 bartype=h \
	wantaxis1=y wantaxis2=y labelsz=6 titlesz=6 wantscalebar=1 barwidth=0.15 > /dev/null out=${R}/$@

okim.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.bkim2W02ffd.vout11_41.H min2=0 n2=1 min1=1000 max1=8000 n4=1 min3= | Transp |\
	Grey pclip=95 title="(b)" label1="z" label2="x" screenratio= 4 screenht= 3.5 screenwd= 4 \
	wantaxis1=y wantaxis2=n labelsz=6 titlesz=6 wheretitle=b gainpanel=a > /dev/null out=${R}/$@

opti.v: oslow.v okim.v
	vppen < ${R}/oslow.v xcenter=0 ycenter=-9. vpstyle=n xsize=4 > junk1.v
	vppen < ${R}/okim.v xcenter=0 ycenter=-0. vpstyle=n xsize=4 > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

islow.v: tslow.v bslow.v
	vppen < ${R}/tslow.v xcenter=0 ycenter=-3. vpstyle=n > junk1.v
	vppen < ${R}/bslow.v xcenter=0 ycenter=-0. vpstyle=n > junk2.v
	vp_Overlay junk1.v junk2.v > ${R}/$@

sfperm11.v:
	< ${Prst_ExpRefl}/Tese/3000marm2W02ffd.inv_vout11.perm.H Window3d min1=1000 max1=8000 min3=1500 n4=1 f4=40 | Transp |\
	Gpow gpow=-1 | Grey newclip=1 color=j bclip=0.0001 eclip=0.000666667 title=" " label1="z" barlabel="  " \
	label2="x" screenratio=4 screenht=2.5 screenwd=5 wheretitle=b \
	wantaxis1=y wantaxis2=n labelsz=3 titlesz=3 wantscalebar=1 barwidth=0.1 > /dev/null out=${R}/$@

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

fangtrue.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=1500 min1=1500|\
	Grey title="1500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=y wantaxis2=y labelsz=2 > /dev/null out=0.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=2500 min1=1500|\
	Grey title="2500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=1.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=3500 min1=1500|\
	Grey title="3500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=2.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=4500 min1=1500|\
	Grey title="4500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=3.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=5500 min1=1500|\
	Grey title="5500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=4.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=6500 min1=1500|\
	Grey title="6500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=5.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vtrue.ffd.H n3=1 min3=7500 min1=1500|\
	Grey title="7500" pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2="angle" screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=y labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-1.0 ycenter=-6.5 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1.5 ycenter=-6.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2.0 ycenter=-6.5 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-2.5 ycenter=-6.5 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-3.0 ycenter=-6.5 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-3.5 ycenter=-6.5 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-4.0 ycenter=-6.5 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk2.v junk4.v junk6.v junk1.v junk3.v junk5.v> ${R}/$@	
	rm -f [0-6].v junk[0-6].v

fangback.v:
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=1500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=y wantaxis2=n labelsz=2 > /dev/null out=0.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=2500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=1.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=3500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=2.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=4500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=3.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=5500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=4.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=6500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=5.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W.vsmoo2W.ffd.H n3=1 min3=7500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-1.0 ycenter=-5 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1.5 ycenter=-5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2.0 ycenter=-5 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-2.5 ycenter=-5 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-3.0 ycenter=-5 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-3.5 ycenter=-5 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-4.0 ycenter=-5 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk2.v junk4.v junk6.v junk1.v junk3.v junk5.v> ${R}/$@	
	rm -f [0-6].v junk[0-6].v

fang11.v:
	Window3d < ${Prst_ExpRefl}/Tese/ n3=1 min3=1500 min1=1500|\
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=1500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=y wantaxis2=n labelsz=2 > /dev/null out=0.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=2500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=1.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=3500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=2.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=4500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=3.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=5500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=4.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=6500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=5.v
	Window3d < ${Prst_ExpRefl}/Tese/3000marm.ang2W02.ffd.H n3=1 min3=7500 min1=1500|\
	Grey title=" " pclip=95 label1="z" titlesz=2 labelsz=2 o2num=-40 d2num=40 n2tic=2 \
	label2=" " screenratio=4 screenht=1.75 screenwd=.5 wheretitle=t \
	wantaxis1=n wantaxis2=n labelsz=2 > /dev/null out=6.v
	vppen < 0.v xcenter=-1.0 ycenter=-3.5 vpstyle=n > junk0.v
	vppen < 1.v xcenter=-1.5 ycenter=-3.5 vpstyle=n > junk1.v
	vppen < 2.v xcenter=-2.0 ycenter=-3.5 vpstyle=n > junk2.v
	vppen < 3.v xcenter=-2.5 ycenter=-3.5 vpstyle=n > junk3.v
	vppen < 4.v xcenter=-3.0 ycenter=-3.5 vpstyle=n > junk4.v
	vppen < 5.v xcenter=-3.5 ycenter=-3.5 vpstyle=n > junk5.v
	vppen < 6.v xcenter=-4.0 ycenter=-3.5 vpstyle=n > junk6.v
	vp_Overlay junk0.v junk2.v junk4.v junk6.v junk1.v junk3.v junk5.v> ${R}/$@	
	rm -f [0-6].v junk[0-6].v


fang1.v: fangback.v fangtrue.v fang11.v
	vp_Overlay ${R}/fangback.v ${R}/fangtrue.v ${R}/fang11.v > ${R}/$@

###pdf
bkim.pdf:
	pstexpen ${R}/bkim.v ${R}/bkim.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/bkim.ps

okim.pdf:
	pstexpen ${R}/okim.v ${R}/okim.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/okim.ps

tkim.pdf:
	pstexpen ${R}/tkim.v ${R}/tkim.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/tkim.ps

bwkim.pdf:
	pstexpen ${R}/bwkim.v ${R}/bwkim.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/bwkim.ps

dimg1.pdf:
	pstexpen ${R}/dimg1.v ${R}/dimg1.ps color=n fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/dimg1.ps

islow.pdf:
	pstexpen ${R}/islow.v ${R}/islow.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/islow.ps

tslow.pdf:
	pstexpen ${R}/tslow.v ${R}/tslow.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/tslow.ps

bslow.pdf:
	pstexpen ${R}/bslow.v ${R}/bslow.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/bslow.ps

oslow.pdf:
	pstexpen ${R}/oslow.v ${R}/oslow.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/oslow.ps

dslow.pdf:
	pstexpen ${R}/dslow.v ${R}/dslow.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/dslow.ps

sfperm11.pdf:
	pstexpen ${R}/sfperm11.v ${R}/sfperm11.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/sfperm11.ps

fang1.pdf:
	pstexpen ${R}/fang1.v ${R}/fang1.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/fang1.ps

fangback.pdf:
	pstexpen ${R}/fangback.v ${R}/fangback.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/fangback.ps

fangtrue.pdf:
	pstexpen ${R}/fangtrue.v ${R}/fangtrue.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/fangtrue.ps

fang11.pdf:
	pstexpen ${R}/fang11.v ${R}/fang11.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/fang11.ps

true.pdf:
	pstexpen ${R}/true.v ${R}/true.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/true.ps

bckg.pdf:
	pstexpen ${R}/bckg.v ${R}/bckg.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/bckg.ps

opti.pdf:
	pstexpen ${R}/opti.v ${R}/opti.ps color=y fat=1 fatmult=1.5 invras=n force=  invras=y
	epstopdf ${R}/opti.ps