#Figures

GREYD= Grey label1="z (ft) " label2="x (ft)" labelsz=4 titlesz=4 geophysics=1  crowd1=.75  crowd2=.75 title=" " min2=30600 max2=39600 > /dev/null 

diag-Sis.v: LtL_Sis_az41_ax41_wave_marine_strati_t.H 
	Window3d < LtL_Sis_az41_ax41_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 | ${GREYD} newclip=1 bclip=-.5 eclip=2.5 polarity=-1 allpos=y out=$@ 

refle-new.v: 
	Cp  ${SIGSMULT_DATA}/sigsbee2b_reflection_coefficients_77.H sigsbee2b_reflection_coefficients_77.H 
	echo d1=25 d2=25 o2=10000>> sigsbee2b_reflection_coefficients_77.H
	< sigsbee2b_reflection_coefficients_77.H Window3d min1=8000 n1=381 min2=30000 n2=133 j2=3| Transp reshape=2,4 > reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Window3d n3=381  > pp.H
	Math file1=pp.H exp="-1.*file1" > 1.H 
	< 1.H ${GREYD} title=" " pclip=95 o1num=9000  out=$@ 

mig-stack.v: mig_Sis_pad_800_strati_nointerp_ang_stack.H 
	Cp mig_Sis_pad_800_strati_nointerp_ang_stack.H 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=$@ 

mig-stack-filt.v: mig_Sis_strati_ang_claudio.H #filter migration in angle
	<  mig_Sis_strati_ang_claudio.H Window3d | Transp | Stack3d  > 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=$@ 
 
migoff1.v: mig_Sis_pad_800_strati.H 
	< mig_Sis_pad_800_strati.H Scale scale_to=1.| Window3d min1=30000 n1=133 min5=8000 n5=381 min3=0 n3=1 |Transp > 1.H 
	< 1.H Grey min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.03 polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< /net/koko/homes/sep/alejandro/research/sep124/mig_Sis_pad_800_strati.H Scale scale_to=1. | Window3d  n1=1 min1=36800 min5=8000 n5=381 min3=-600 n3=17 |Transp > 2.H
	< 2.H Grey label2="h (ft)" title=" " clip=0.03 crowd1=1. polarity=-1 screenratio=2. out=2.v geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=2.55,2.55 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> $@ 

migoff1-filt.v:  ~claudio/Research/Prst_ExpRefl/RTM/prim_01_subt_sort.H#filter migration offset needs rule
	< ~claudio/Research/Prst_ExpRefl/RTM/prim_01_subt_sort.H Scale scale_to=1. | Window3d min1=30000 n1=133 min5=8000 n5=381 min3=0 n3=1 |Transp > 1.H 
	< 1.H Grey   min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.03 polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< prim_01_subt_sort.H Scale scale_to=1. | Window3d  n1=1 min1=36800 min5=8000 n5=381 min3=-600 n3=17 |Transp > 2.H
	< 2.H Grey label2="h (ft)" title=" " clip=0.03 crowd1=1. polarity=-1 screenratio=2. out=2.v geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=2.55,2.55 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> $@ 

fragreya1=eout=1  gainpanel=a  label1="z (ft)" label2="x (ft)" label3="angle (degree)" title=" "  
frames1a= frame1=240 frame2=30 frame3=35 title=" "  crowd=.65 screenratio=1 point1=.75 point2=.75 polarity=-1 wantcoordlabel=-1 labelsz=3 titlesz=3

migang1.v: mig_Sis_pad_800_strati_nointerp_ang_cow.H 
	< mig_Sis_pad_800_strati_nointerp_ang_cow.H Scale scale_to=1.| Window3d f1=9 n1=1  > 1.H
	< 1.H Grey min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.05    polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< mig_Sis_pad_800_strati_nointerp_ang_cow.H Scale scale_to=1. | Window3d  n4=1 min4=36800 |Transp > 2.H
	< 2.H Grey label2="angle (degree)" title=" " clip=0.05 crowd1=1. polarity=-1 screenratio=2. out=2.v  geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=.8,.8 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> $@ 

migang1-filt.v: mig_Sis_strati_ang_claudio.H 
	< mig_Sis_strati_ang_claudio.H Scale scale_to=1 | Window3d f1=9 n1=1  > 1.H
	< 1.H Grey  min2=30600 max2=39600 label1="z (ft)" label2="x (ft)" title=" " clip=0.05   polarity=-1 screenratio=.69 o1num=9000  out=1.v geophysics=1 > /dev/null
	< 1.H Overlay type=polyline x=10.,10. y=1.3,8.95 out=11.v >/dev/null
	vp_Overlay 1.v 11.v > 111.v
	< mig_Sis_strati_ang_claudio.H Scale scale_to=1 |Window3d n4=1 min4=36800 |Transp > 2.H
	< 2.H Grey label2="angle (degree)" title=" " clip=0.05 crowd1=1. polarity=-1 screenratio=2. out=2.v  geophysics=1> /dev/null
	< 2.H Overlay type=polyline x=.8,.8 y=1.3,8.95 out=22.v >/dev/null
	vp_Overlay 2.v 22.v > 222.v
	vp_SideBySideAniso 111.v 222.v | vppen vpstyle=n hshift=2 vshift=-.5> $@

GREYK=Grey clip=10 allpos=y label1="kx (1/ft)" label2="kh (1/ft)" title=" " geophysics=1 min1=-0.03 max1=0.03 min2=-0.03 max2=0.03> /dev/null
filtered.v: # need rules
	Window3d min1=14150 n1=1 < /net/koko/homes/sep/alejandro/research/sep124/filtered.H | Transp plane=12  >1.H
	< 1.H ${GREYK} out=$@

not-filtered.v: #need rules
	Window3d min1=14150 n1=1 < /net/koko/homes/sep/alejandro/research/sep124/not_filtered.H | Transp plane=12 >1.H
	< 1.H ${GREYK} out=$@

diff-filtered.v: # need rules
	Window3d min1=14150 n1=1 < /net/koko/homes/sep/alejandro/research/sep124/not_filtered.H | Transp plane=12 >1.H
	Window3d min1=14150 n1=1 < /net/koko/homes/sep/alejandro/research/sep124/filtered.H | Transp plane=12 >2.H
	Math file1=1.H file2=2.H exp="file1-file2" > 3.H
	< 3.H Grey ${GREYK} out=$@
	
invz.v: 
	Window3d </net/koko/homes/sep/alejandro/research/sep124/mmov_new_sym.H f6=6 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=85 o1num=9000  out=$@ 

invz-filt.v: 
	Window3d <inv_Sis_az41ax41_newsolver_sym_strati_prec_claudio2.0.H> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=87 o1num=9000  out=$@ 

resinvz.v: 
	Window3d </net/koko/homes/sep/alejandro/research/sep124/rmov0.H f6=9 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=100 o1num=9000  out=$@ 

resinvz-filt.v: 
	Window3d </net/koko/homes/sep/alejandro/research/sep124/rmov_claudio0.0.H f6=9 n6=1> 1s.H
	echo o1=8000 >> 1s.H
	< 1s.H ${GREYD} title=" " pclip=100 o1num=9000  out=$@
 

###

vel_strati_pad.H:
	Cp  ${SIGSMULT_DATA}/sigsbee2b_stratigraphy_f77.H sigsbee2b_stratigraphy_f77.H
	echo "d1=25 o2=10000 d2=25"      >> sigsbee2b_stratigraphy_f77.H
	< sigsbee2b_stratigraphy_f77.H Window3d n1=701 f2=1 j2=3|Pad beg2=2000 n2out=10000 extend=1|Transp> $@
	/net/koko/alejandro/plib/pcs2/bin/Wedge.x $@ axis=2 nn=1 no=0. nd=75. 
	echo  "  label1="x" label2="y" label3="z" " >> $@ 

wave.H:
	Wavelet n1=1024 o1=0 d1=0.011719 domain=time fund=20 wavelet=ricker2 > t1$@
	Add scale=-1 t1$@ > t2$@
	< t2$@ /opt/SEP/SEP6.4.6/bin/Transf shotmig=y is2d=y wei=y f_min=5 f_min1=10 f_max1=30 f_max=35 maxsize=1000  >$@
	echo o3=4.91667    d3=0.0833333 >> $@
	Rm t*$@
 
wave_pad.H:
	Wavelet n1=1024 o1=0 d1=0.011719 domain=time fund=20 wavelet=ricker2 > t1$@
	Add scale=-1 t1$@ > t2$@
	< t2$@ /opt/SEP/SEP6.4.6/bin/Transf shotmig=y is2d=y wei=y f_min=5 f_min1=10 f_max1=30 f_max=35 maxsize=1000  >t$@
	<t$@ Pad n4out=1005 extend=1>> $@
	echo o3=4.91667    d3=0.0833333 o4=10925      d4=75 >> $@

data_f.H:
	< ${SIGSMULT_DATA}/sigsbee2b_nfs_shot.H Window3d max_memory=1000 max3=60000 |Mute vmute=4920 tmute=2 > t1$@
	< t1$@ /opt/SEP/SEP6.4.6/bin/Transf shotmig=y is2d=y wei=y f_min=5 f_min1=10 f_max1=30 f_max=35 maxsize=1000 >t3$@
	< t3$@ Pad beg1=400 n1=800 > $@
	Rm t*$@
 
green_Sis_strati_OMP.H: wave_pad.H vel_strati_pad.H 
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig.x D_d=wave_pad.H V=vel_strati_pad.H W_d=$@  \
	nx=1200 ny=1 ox=0 oy=0 dx=75 dy=75 min_slow_dev=1. niter_lloyd=100 watercolumn=6000\
	nref=15  taper=y  taper_pcent=2 marine=n firstzW=320 nx_W=133 ny_W=1 ox_W=30000 oy_W=0 dx_W=75 dy_W=75 

mig_Sis_pad_800_strati.H:  data_f.H wave.H vel_strati_pad.H 
	../SPWEI_SERIAL/bin/LINUX86_64/ShotMig.x D=data_f.H D_d=wave.H V=vel_strati_pad.H R=$@  \
	taper=y  taper_pcent=2 marine=y watercolumn=6000  nhx=33 nref=15 min_slow_dev=1. niter_lloyd=100 

LtL_Sis_az41_ax41_wave_marine_strati.H: green_Sis_strati_OMP.H /Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x 
	date
	< green_Sis_strati_OMP.H Window3d squeeze=n f3=2 n3=7 > t$@ 
	< t$@ Transp reshape=2,3,6 plane=23 max_memory=1000 | Transp reshape=3,5,6 max_memory=1000> t1$@
	< t1$@ ~valencia/research/Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@  \
	nax=41 naz=21 nay=1 marine=y noffx=1 \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=30000 nx=133 miniz=7975 nz=381 
	date 

LtL_Sis_az41_ax41_wave_marine_strati_t.H: LtL_Sis_az41_ax41_wave_marine_strati.H 
	< LtL_Sis_az41_ax41_wave_marine_strati.H Transp reshape=3,6,10 > pp1.H
	< pp1.H Transp reshape=2,3,5,6,10 plane=24 > pp2.H
	< pp2.H Transp reshape=7,9,10 plane=23 > $@

mig_Sis_pad_800_strati_nointerp_ang_stack.H: mig_Sis_pad_800_strati.H 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=7975 n5=381 min3=-600 n3=17 squeeze=n | Transp reshape=2,4,5 | Transp reshape=2,4,5 plane=23> pppp.H
	< pppp.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-10 dax=1. naax=71 style=1 adj=1> 1.H 
	< 1.H Window3d | Transp | Stack3d > $@ 

mig_Sis_pad_800_strati_nointerp_ang_cow.H: mig_Sis_pad_800_strati.H 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=7975 n5=381 min3=-600 n3=17 squeeze=n | Transp reshape=2,4,5 | Transp reshape=2,4,5 plane=23> pppp.H
	< pppp.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-10 dax=1. naax=71 style=1 adj=1 cosw=1> $@
 
mig_Sis_strati_ang_claudio.H: ~claudio/Research/Prst_ExpRefl/RTM/prim_01_subt_sort.H#needs rule filter
	< ~claudio/Research/Prst_ExpRefl/RTM/prim_01_subt_sort.H Window3d min1=30600 n1=121 min5=8000 n5=381  min3=-600 n3=17 squeeze=n  > ppp.H
	< ppp.H Transp plane=13| Transp reshape=2,4,5,6 plane=23  > ppp1.H
	< ppp1.H  Scale scale_to=1.  > ppp3.H
	Math file1=ppp3.H exp="10.*file1" > pppp.H
	< pppp.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-10 dax=1. naax=71 style=1 adj=1> $@

# Inversion
 
inv_Sis_az41ax41_newsolver_sym_strati_prec.H rmov0.H mmov_new_sym.H: mig_Sis_pad_800_strati.H LtL_Sis_az41_ax41_wave_marine_strati_t.H 
	< mig_Sis_pad_800_strati.H Window3d min1=30600 n1=121 min5=8000 n5=381  min3=0 n3=1 squeeze=n  > ppp.H
	< ppp.H Transp plane=13| Transp reshape=2,4,5,6 plane=23  > ppp1.H
	< ppp1.H  Scale scale_to=1. > p-p3.H
#   Preconditioner
	Window3d  < LtL_Sis_az41_ax41_wave_marine_strati_t.H min9=30600 n9=121 squeeze=n  > LtL_Sis_az41_ax41_wave_marine_strati_t_w.H
	Window3d  < LtL_Sis_az41_ax41_wave_marine_strati_t_w.H n3=1 min3=0 n4=1 min4=0 min9=30600 n9=121 squeeze=n | Transp reshape=5,10 > LtL_Sis_az41_ax41_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az41_ax41_wave_marine_strati_t_diag.H exp="1./(file1)" > LtL_Sis_az41_ax41_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az41_ax41_wave_marine_strati_t_w.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_off.x mig=p-p3.H  prec=LtL_Sis_az41_ax41_wave_marine_strati_t_diag_inv.H  niter=200 eps=0.00 alpha=.001 na=1 reg_off=0 reg=0 remov=rmov0.H mmov=mmov_new_sym.H wmov=wmov.H beta_max=.001>$@
 
inv_Sis_az41ax41_newsolver_sym_strati_prec_claudio%.H rmovclaudio0.H: LtL_Sis_az41_ax41_wave_marine_strati_t.H  prim_01_subt_sort.H     
	< prim_01_subt_sort.H  Window3d min1=30600 n1=121 min5=8000 n5=381  min3=0 n3=1 squeeze=n  > ppp.H
	< ppp.H Transp plane=13| Transp reshape=2,4,5,6 plane=23  > ppp1.H
	< ppp1.H  Scale scale_to=1. > p-p3.H
#   Preconditioner
	Window3d  < LtL_Sis_az41_ax41_wave_marine_strati_t.H min9=30600 n9=121 squeeze=n  > LtL_Sis_az41_ax41_wave_marine_strati_t_w.H
	Window3d  < LtL_Sis_az41_ax41_wave_marine_strati_t_w.H n3=1 min3=0 n4=1 min4=0 min9=30600 n9=121 squeeze=n | Transp reshape=5,10 > LtL_Sis_az41_ax41_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az41_ax41_wave_marine_strati_t_diag.H exp="1./(file1)" > LtL_Sis_az41_ax41_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az41_ax41_wave_marine_strati_t_w.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_off.x mig=p-p3.H  prec=LtL_Sis_az41_ax41_wave_marine_strati_t_diag_inv.H  niter=200 eps=0.00 alpha=.001 na=1 reg_off=0 reg=0 remov=rmovclaudio0.H mmov=mmov_new_sym_claudio.H wmov=wmov.H beta_max=.001 >$@