include ${SEPINC}/SEP.top
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
BINDIR = ./Bin/${MTYPE}
OBJDIR = ./Obj/${MTYPE}
SRCDIR = ./Src
RESDIR = ./Fig
PARDIR = ./Par
TXTDIR = ./Txt
COLOR=y


RESULTSCR =  mig_Sis hesian_phase_Sis comp_Sis_full mig_inv mig_inv1 mig_off comp_Sis_full1 comp_Sis_full2 
RESULTSER = Sis_vel strati illu mig_Sis_pos inv_Sis_pos 
UF90FLAGS = -O3 -axP  -static-libcxa
UF90LIBS= -lsepgeef90 -lsep2df90  

data_over_f.H:
	< ${SIGSMULT_DATA}/sigsbee2b_nfs_shot.H Window3d max_memory=1000 max3=60000 |Mute vmute=4920 tmute=2 > t1$@
	< t1$@ Interp maxsize=2100  d3out=75 > t2$@
	< t2$@ 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=32 n1=400 n4=709> $@
	Rm t*$@

data_f.H:
	< ${SIGSMULT_DATA}/sigsbee2b_nfs_shot.H Window3d max_memory=1000 max3=60000 |Mute vmute=4920 tmute=2 > t1$@
	< t1$@ 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=132 n1=500 > $@
	Rm t*$@

wave.H:
	Wavelet n1=1024 o1=0 d1=0.011719 domain=time fund=20 wavelet=ricker2 > t1$@
	Add scale=-1 t1$@ > t2$@
	< t2$@ 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 >> $@
#	Rm t*$@

data_over_f_1.H:
	< data_over_f.H Window3d squeeze=n n4=1 f4=100 > t$@
	< t$@ Pad beg1=32 > $@

green_test_1.H:
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX/ShotMig  D_d=wave.H V=vel_pad.H W_d=$@  \
	nx=100 ny=1 ox=0 oy=0 dx=75 dy=75 nsx=1 nsy=1 osx=3000 osy=1  \
	nxt=50 nyt=1 oxt=2000 oyt=0 dxt=75 dyt=75 \
	nref=1  taper=y  taper_pcent=2 marine=n firstzW=400

vel_pad.H:
	< ${SIGSMULT_DATA}/sigsbee2b_migration_velocity.H Window> j.H
	echo "hff=-1 gff=-1 d1=25 o2=10025 d2=37.5"      >>j.H
	< j.H Window3d n1=701 j2=2|Pad beg2=2000 n2out=10000 extend=1|Transp> vel_pad.H
	/net/koko/alejandro/plib/pcs2/bin/Wedge.x vel_pad.H axis=2 nn=1 no=0. nd=75. 
	echo  "  label1="x" label2="y" label3="z" " >> vel_pad.H

mig_Sis_Parallel_16000.H:  data_f.H wave.H vel_pad.H ${PARDIR}/go_sis.P 
	@date
	nohup Parallel command=/net/koko/brad/bin/LINUX/ShotMig files="data_f.H,wave.H,vel_pad.H,$@"  \
	tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" nblock=361 axis=4,4,5,5 pverb=2  \
	file_type="DISTRIBUTE,COPY,COPY,COPY"  mach_mfile=machinefile_all nx_shots=327 ox_shots=10925. \
	par=${PARDIR}/go_sis.P  taper=y  taper_pcent=2 marine=y watercolumn=6000 #restart=y

mig_Sis_Parallel_pad.H:  data_f.H wave.H vel_pad.H ${PARDIR}/go_sis.P 
	@date
	nohup Parallel command=./ShotMigNew files="data_f_pad.H,wave.H,vel_pad.H,$@"  \
	tags="D,D_d,V,R" usage="INPUT,INPUT,INPUT,OUTPUT" nblock=37 axis=4,4,5,5 pverb=2  \
	file_type="DISTRIBUTE,COPY,COPY,COPY"  mach_mfile=machinefile_half nx_shots=327 ox_shots=10925. \
	par=${PARDIR}/go_sis.P  taper=y  taper_pcent=2 marine=y watercolumn=6000 minx=10925 maxx=60000   #restart=y

test.H:
	/net/koko/brad/bin/LINUX/ShotMigNew D=dat-b.H D_d=wav-b.H V=vel_pad.H R=test.H nref=1 taper=y  taper_pcent=2 marine=y watercolumn=6000 minx=10925 maxx=60040

mig_ang_97.H: 
	< mig_off64.H OFF2ANG  verb=y eps=1 > $@

green_Sis: wave_pad.H vel_pad.H ${PARDIR}/go_sis.P 
	@date
	nohup Parallel command=/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX/ShotMig files="wave_pad.H,vel_pad.H,$@.H"  \
	tags="D_d,V,W_d" usage="INPUT,INPUT,OUTPUT" nblock=22 axis=4,5,5 pverb=2  \
	file_type="DISTRIBUTE,COPY,DISTRIBUTE" mach_mfile=machinefile_half \
	nx=770 ny=1 ox=7925 oy=0 dx=75 dy=75 \
	par=${PARDIR}/go_sis.P  taper=y  taper_pcent=2 marine=n firstzW=320 nx_W=66 ny_W=1 ox_W=32000 oy_W=0 dx_W=75 dy_W=75 restart=y
	
green_Sis_OMP: wave_pad.H vel_pad.H 
	@date
	nohup /net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig.x D_d=wave_pad.H V=vel_pad.H W_d=pppp.H  \
	nx=1200 ny=1 ox=0 oy=0 dx=75 dy=75 \
	nref=8  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 restart=y
 

green_Sis_t.H:
	< green_Sis.H Transp reshape=2,3,6 maxmemory=1000 maxsise=1000 plane=23 > $@

green_Sis_test :  data_over_f.H wave_pad.H vel_pad.H ${PARDIR}/go_sis.P 
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX/ShotMig D_d=wave_pad.H V=vel_pad.H  W_d=$@.H \
	nx=750 ny=1 ox=7925 oy=0 dx=75 dy=75 \
	par=${PARDIR}/go_sis.P  taper=y  taper_pcent=2 marine=n firstzW=320 nx_W=66 ny_W=1 ox_W=32000 oy_W=0 dx_W=75 dy_W=75#restart=y 
#	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX/ShotMig D_d=wave_pad.H V=vel_pad.H  W_d=$@ \
  	/net/koko/brad/bin/LINUX/ShotMig D_d=wave_pad.H V=vel_pad.H  W_d=$@ \ 

LtL_Sis_az21_ax21_wave_old4.H: 
	date
	Window3d < green_Sis_OMP_f35_t.H n6=1 f6=6 squeeze=n > ppp.H
	< ppp.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_4_omp.x outtag=$@  \
	nax=21 naz=21 nay=1   marine=y noffx=5 \
	nry=1 nrx=348 orx=0. nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=64

LtL_Sis_az21_ax21_wave_old.H: 
	date
	Window3d < green_Sis_OMP_f35_t.H n6=1 f6=6 squeeze=n > ppp.H
	< ppp.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_2_omp.x outtag=$@  \
	nax=21 naz=21 nay=1   marine=y noffx=1 \
	nry=1 nrx=348 orx=0. nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=64
 

LtL_Sis_az15_ax15_wave_old_marineno.H: 
	date
	Window3d < wave.H n3=1 f3=110 j3=200 squeeze=n >pp_wave.H
	Window3d < green_Sis_t.H n6=1 f6=110 j6=200 squeeze=n > ppp.H
#	Window3d < ppp.H min3=10000 n3=2 squeeze=n> pppp.H
	< ppp.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_diag_omp.x outtag=$@  \
	nax=1 naz=1 nay=1 nahx=1 nahy=1  marine=n noffx=33 noffy=1 \
	nry=1 nrx=709 orx=10925 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=64 

LtL_Sis_az21_ax21_wave_old_ang33_f35_border_sodom.H:
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_border_sodom.H  Transp reshape=1,2,4,10 plane=13 maxsize=10000> t$@
#	< t$@ ./OFF2ANG verb=y amax=45 maxmem=2000> t1$@
#	Rm t$@
	< t1$@ Transp maxsize=500 reshape=5,10 maxsize=2000> t2$@
	Rm t1$@
	< t2$@ Transp maxsize=500 reshape=4,5,10 maxsize=2000> t3$@
	Rm t2$@
	< t3$@ ./OFF2ANG verb=y amax=45 eps=10 maxmem=2000> $@
	Rm t3$@

diagonal_off:
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=0 n6=1 min6=0 n8=1 min8=0 n9=1 f9=0 | Transp > n0.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=1 n6=1 min6=0 n8=1 min8=0 n9=1 f9=1 | Transp > n1.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=2 n6=1 min6=0 n8=1 min8=0 n9=1 f9=2 | Transp > n2.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=3 n6=1 min6=0 n8=1 min8=0 n9=1 f9=3 | Transp > n3.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=4 n6=1 min6=0 n8=1 min8=0 n9=1 f9=4 | Transp > n4.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=5 n6=1 min6=0 n8=1 min8=0 n9=1 f9=5 | Transp > n5.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=6 n6=1 min6=0 n8=1 min8=0 n9=1 f9=6 | Transp > n6.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=7 n6=1 min6=0 n8=1 min8=0 n9=1 f9=7 | Transp > n7.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=8 n6=1 min6=0 n8=1 min8=0 n9=1 f9=8 | Transp > n8.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=9 n6=1 min6=0 n8=1 min8=0 n9=1 f9=9 | Transp > n9.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=10 n6=1 min6=0 n8=1 min8=0 n9=1 f9=10 | Transp > n10.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=11 n6=1 min6=0 n8=1 min8=0 n9=1 f9=11 | Transp > n11.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=12 n6=1 min6=0 n8=1 min8=0 n9=1 f9=12 | Transp > n12.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=13 n6=1 min6=0 n8=1 min8=0 n9=1 f9=13 | Transp > n13.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=14 n6=1 min6=0 n8=1 min8=0 n9=1 f9=14 | Transp > n14.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=15 n6=1 min6=0 n8=1 min8=0 n9=1 f9=15 | Transp > n15.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=16 n6=1 min6=0 n8=1 min8=0 n9=1 f9=16 | Transp > n16.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=17 n6=1 min6=0 n8=1 min8=0 n9=1 f9=17 | Transp > n17.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=18 n6=1 min6=0 n8=1 min8=0 n9=1 f9=18 | Transp > n18.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=19 n6=1 min6=0 n8=1 min8=0 n9=1 f9=19 | Transp > n19.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=20 n6=1 min6=0 n8=1 min8=0 n9=1 f9=20 | Transp > n20.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=21 n6=1 min6=0 n8=1 min8=0 n9=1 f9=21 | Transp > n21.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=22 n6=1 min6=0 n8=1 min8=0 n9=1 f9=22 | Transp > n22.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=23 n6=1 min6=0 n8=1 min8=0 n9=1 f9=23 | Transp > n23.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=24 n6=1 min6=0 n8=1 min8=0 n9=1 f9=24 | Transp > n24.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=25 n6=1 min6=0 n8=1 min8=0 n9=1 f9=25 | Transp > n25.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=26 n6=1 min6=0 n8=1 min8=0 n9=1 f9=26 | Transp > n26.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=27 n6=1 min6=0 n8=1 min8=0 n9=1 f9=27 | Transp > n27.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=28 n6=1 min6=0 n8=1 min8=0 n9=1 f9=28 | Transp > n28.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=29 n6=1 min6=0 n8=1 min8=0 n9=1 f9=29 | Transp > n29.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=30 n6=1 min6=0 n8=1 min8=0 n9=1 f9=30 | Transp > n30.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=31 n6=1 min6=0 n8=1 min8=0 n9=1 f9=31 | Transp > n31.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_border_t510.H n1=1 f1=32 n6=1 min6=0 n8=1 min8=0 n9=1 f9=32 | Transp > n32.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H n19.H \
          n20.H n21.H n22.H n23.H n24.H n25.H n26.H n27.H n28.H n29.H \
	      n30.H n31.H n32.H axis=3 >LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.H
	Rm n*.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.H Grey eout=1 gainpanel=a > LtL_Sis_az21_ax21_wave_old_off33_f35_border_diag.T

LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_newdiag_t.H: LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag.H Transp reshape=2,3,5,10 plane=13> pp1.H
	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H

diagonal_off_new:
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=0 n4=1 min4=0 n3=1 min3=0 n6=1 f6=0 | Transp > n0.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=1 n4=1 min4=0 n3=1 min3=0 n6=1 f6=1 | Transp > n1.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=2 n4=1 min4=0 n3=1 min3=0 n6=1 f6=2 | Transp > n2.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=3 n4=1 min4=0 n3=1 min3=0 n6=1 f6=3 | Transp > n3.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=4 n4=1 min4=0 n3=1 min3=0 n6=1 f6=4 | Transp > n4.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=5 n4=1 min4=0 n3=1 min3=0 n6=1 f6=5 | Transp > n5.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=6 n4=1 min4=0 n3=1 min3=0 n6=1 f6=6 | Transp > n6.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=7 n4=1 min4=0 n3=1 min3=0 n6=1 f6=7 | Transp > n7.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=8 n4=1 min4=0 n3=1 min3=0 n6=1 f6=8 | Transp > n8.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=9 n4=1 min4=0 n3=1 min3=0 n6=1 f6=9 | Transp > n9.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=10 n4=1 min4=0 n3=1 min3=0 n6=1 f6=10 | Transp > n10.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=11 n4=1 min4=0 n3=1 min3=0 n6=1 f6=11 | Transp > n11.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=12 n4=1 min4=0 n3=1 min3=0 n6=1 f6=12 | Transp > n12.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=13 n4=1 min4=0 n3=1 min3=0 n6=1 f6=13 | Transp > n13.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=14 n4=1 min4=0 n3=1 min3=0 n6=1 f6=14 | Transp > n14.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=15 n4=1 min4=0 n3=1 min3=0 n6=1 f6=15 | Transp > n15.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H n1=1 f1=16 n4=1 min4=0 n3=1 min3=0 n6=1 f6=16 | Transp > n16.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H axis=3 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_diag.H
	Rm n*.H
	<LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_diag.H  Grey eout=1 gainpanel=a > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.T

diagonal_off_new1:LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_newdiag_t.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=0 n4=1 min4=0 n3=1 min3=0 n6=1 f6=0 | Transp > n0.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=1 n4=1 min4=0 n3=1 min3=0 n6=1 f6=1 | Transp > n1.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=2 n4=1 min4=0 n3=1 min3=0 n6=1 f6=2 | Transp > n2.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=3 n4=1 min4=0 n3=1 min3=0 n6=1 f6=3 | Transp > n3.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=4 n4=1 min4=0 n3=1 min3=0 n6=1 f6=4 | Transp > n4.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=5 n4=1 min4=0 n3=1 min3=0 n6=1 f6=5 | Transp > n5.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=6 n4=1 min4=0 n3=1 min3=0 n6=1 f6=6 | Transp > n6.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=7 n4=1 min4=0 n3=1 min3=0 n6=1 f6=7 | Transp > n7.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=8 n4=1 min4=0 n3=1 min3=0 n6=1 f6=8 | Transp > n8.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=9 n4=1 min4=0 n3=1 min3=0 n6=1 f6=9 | Transp > n9.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=10 n4=1 min4=0 n3=1 min3=0 n6=1 f6=10 | Transp > n10.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=11 n4=1 min4=0 n3=1 min3=0 n6=1 f6=11 | Transp > n11.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=12 n4=1 min4=0 n3=1 min3=0 n6=1 f6=12 | Transp > n12.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=13 n4=1 min4=0 n3=1 min3=0 n6=1 f6=13 | Transp > n13.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=14 n4=1 min4=0 n3=1 min3=0 n6=1 f6=14 | Transp > n14.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=15 n4=1 min4=0 n3=1 min3=0 n6=1 f6=15 | Transp > n15.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=16 n4=1 min4=0 n3=1 min3=0 n6=1 f6=16 | Transp > n16.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=17 n4=1 min4=0 n3=1 min3=0 n6=1 f6=17 | Transp > n17.H
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_t.H n1=1 f1=18 n4=1 min4=0 n3=1 min3=0 n6=1 f6=18 | Transp > n18.H
	Cat3d n0.H  n1.H  n2.H  n3.H  n4.H  n5.H  n6.H  n7.H  n8.H  n9.H \
	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H n17.H n18.H axis=3 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag_test.H
	Rm n*.H

LtL_Sis_az21_ax21_wave_old_off33_f35_border.H: 
	date
	nohup < green_Sis_OMP_f35_t.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_2_omp.x outtag=$@  \
	nax=21 naz=21 nay=1 marine=n noffx=33  \
	nry=1 nrx=709 orx=10925 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=70

LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_testdiag.H: 
	date
	< green_Sis_OMP_f35_t.H Window3d  n6=1 f6=5 squeeze=n |../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_4_omp.x outtag=$@  \
	nax=1 naz=1 nay=1 marine=y noffx=19  \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=30000 nx=33

LtL_Sis_az21_ax21_wave_off33_marine_test_reuse.H: 
	date
	< green_Sis_OMP_f35_t.H Window3d  n6=1 f6=5 n3=50 min3=11000 squeeze=n |../Src_Hessian/Bin/LINUX86_64/Lt_L_reuse_3D_off_wave_4_new.x outtag=$@  \
	nax=21 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

LtL_Sis_az21_ax21_wave_off33_marine_test_old.H: 
	date
	< green_Sis_OMP_f35.H  Transp reshape=2,3,6 plane=23 max_memory=1000 | Transp reshape=3,5,6 max_memory=1000> pp.H
	< pp.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_5_omp.x outtag=$@  \
	nax=21 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
	date

LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_alter_gomo.H: 
	date
	nohup < green_Sis_OMP_f35_alter_w_t.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_4_omp.x outtag=$@  \
	nax=17 naz=31 nay=1 marine=y noffx=17  \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=70 



LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_sodo.H: 
	date
	< green_Sis_OMP_f35_w_t_sodo.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_2_omp.x outtag=$@  \
	nax=17 naz=21 nay=1 marine=y noffx=17  \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=70
	date

LtL_Sis_az21_ax21_wave_old_off33_f35_diag.H: 
	date
	nohup < green_Sis_OMP_f35_t.H ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_diag_omp.x outtag=$@  \
	nax=21 naz=21 nay=1 nahx=11 marine=n noffx=33  \
	nry=1 nrx=709 orx=10925 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=70

LtL_Sis_az21_ax21_wave_old_fullf_test_sn.H: 
	date
	< green_Sis_OMP_f35_t.H Window3d n6=3 j6=3 f6=1 squeeze=n | ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_2_omp.x outtag=$@  \
	nax=1 naz=1 nay=1   marine=n noffx=33  \
	nry=1 nrx=709 orx=10925 nsy=1 nsx=327 osx=10925 dsx=150 minix=33000 nx=40

LtL_Sis_az21_ax21_wave_old_fullf_test_sp.H: 
	date
	< green_Sis_OMP_f35_t.H Window3d n6=3 j6=3 f6=1 squeeze=n | ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_2_sign_omp.x outtag=$@  \
	nax=1 naz=1 nay=1   marine=n noffx=33  \
	nry=1 nrx=709 orx=10925 nsy=1 nsx=327 osx=10925 dsx=150 minix=33000 nx=40

LtL_Sis_az21_ax21_f35_marine_offang17.H:
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine.H Transp reshape=2,3,10  maxsize=10000> t$@
	< t$@ ../sep124/OFF2ANG verb=y amax=45 maxmem=2000> t1$@
	< t1$@ Transp reshape=1,3,10 >  $@

LtL_Sis_az21_ax21_f35_marine_offang16_new.H:
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new.H   Window3d n1=16 n6=16 squeeze=n > t0$@
	< t0$@ Transp reshape=2,3,10  maxsize=10000> t$@
	< t$@ ../sep124/OFF2ANG verb=y amax=45 maxmem=2000> t1$@
	< t1$@ Transp reshape=1,3,10 >  $@


mig_Sis_refle_f5.H: 
#	< LtL_Sis_az21_ax21_wave_old_f5.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 > LtL_Sis_az21_ax21_wave_old_f5_t.H
	< LtL_Sis_az21_ax21_wave_old_f5_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_inv_off.x m=reflection_coefficients_target_conv.H >$@ 

mig_Sis_refle_f35.H: 
	< LtL_Sis_az21_ax21_wave_old_f35.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 > LtL_Sis_az21_ax21_wave_old_f35_t.H
	< LtL_Sis_az21_ax21_wave_old_f35_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_inv_off.x m=reflection_coefficients_target_conv.H >$@

mig_Sis_refle_f35_ba.H: 
#	< LtL_Sis_az21_ax21_wave_old_fullf_ba.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 > LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H
	< reflection_coefficients.H  Window3d min1=7975 n1=383 min2=31000 n2=80 j2=3 | Transp reshape=2,4> ppp.H
	< LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_inv_off.x m=ppp.H >$@

mig_Sis_refle_off_f35_border.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_border.H Window3d squeeze=n min10=11000 n10=200 >pp0.H
#	< pp0.H Transp reshape=2,3,5,10 plane=13> pp1.H
	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_f35_border_t.H
	< reflection_coefficients.H Window3d min1=7975 n1=383 min2=32000 n2=70 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Pad beg1=16 n1=33 | Window3d squeeze=n min3=11000 n3=200 > pp5.H
	< LtL_Sis_az21_ax21_wave_old_f35_border_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=pp5.H >$@

mig_Sis_refle_off_f35_flipoff.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff.H Window3d squeeze=n min10=11000 n10=200 >pp0.H
#	< pp0.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_f35_flipoff_t.H
#	< reflection_coefficients.H Window3d min1=7975 n1=383 min2=32000 n2=70 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
#	< reflection_coefficients_target_conv.H Pad beg1=16 n1=33 | Window3d squeeze=n min3=11000 n3=200 > pp5.H
	< LtL_Sis_az21_ax21_wave_old_f35_flipoff_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=pp5.H >$@ 

mig_Sis_refle_off_f35_marine.H: 
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine.H Transp reshape=2,3,5,10 plane=13> pp1.H
	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_t.H
	< reflection_coefficients.H Window3d min1=7975 n1=383 min2=32000 n2=70 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Pad beg1=8 n1=17 | Window3d squeeze=n min3=11000 n3=200 > pp5.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=pp5.H >$@

mig_Sis_refle_off_f35_marine_new.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H
#	< reflection_coefficients.H Window3d min1=7975 n1=383 min2=32000 n2=70 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
#	< reflection_coefficients_target_conv.H Pad beg1=8 n1=17 | Window3d squeeze=n min3=11000 n3=200 > pp5.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=pp5.H >$@

mig_Sis_refle_offang_f35_marine_new.H: 
#	< LtL_Sis_az21_ax21_f35_marine_offang16_new.H Transp reshape=7,9,10 plane=23 >LtL_Sis_az21_ax21_f35_marine_offang16_t_new.H
	< reflection_coefficients.H Window3d min1=7975 n1=383 min2=32000 n2=70 j2=3 | Transp reshape=2,4> reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Pad beg1=8 n1=16 | Window3d squeeze=n min3=11000 n3=200 > pp5.H
	< LtL_Sis_az21_ax21_f35_marine_offang16_t_new.H ../NoStaDecon/Bin/LINUX86_64/convolve_Mxv_fly_off.x m=pp5.H >$@



inv_Sis_az21ax21_eps10_50_off.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff.H Window3d squeeze=n min10=11000 n10=200 n4=1 min4=0 n6=1 min6=0 >pp0.H
#	< pp0.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > pp6.H 
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 squeeze=n n3=1 min3=0 | Transp reshape=2,5,6  | Scale scale_to=1. > pp3.H
	<  pp6.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  niter=70 eps=.1 alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov.H wmov=wmov.H >$@

inv_Sis_az21ax21_off_marine_new_reg.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 squeeze=n min3=-600 n3=17 | Transp reshape=2,5,6  | Scale scale_to=1. > pp3s.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3s.H  niter=20 initial=1 eps=25. alpha=.001 na=1 reg_off=1 remov=rmov3.H mmov=mmov3.H wmov=wmov.H >$@

inv_Sis_az21ax21_off_marine_new.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 squeeze=n min3=-600 n3=17 | Transp reshape=2,5,6  | Scale scale_to=1. > pp3s.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3s.H  niter=20 initial=1 eps=8. alpha=.001 na=1 reg_damp=1 remov=rmov3.H mmov=mmov3.H wmov=wmov.H >$@
 
inv_Sis_az21ax21_eps10_50_off_marine.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 squeeze=n min3=-600 n3=17 | Transp reshape=2,5,6   > pp3.H
	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  niter=20 initial=0 eps=1. alpha=.001 na=1 reg_off=0 remov=rmov1.H mmov=mmov1.H wmov=wmov.H >$@

inv_Sis_az21ax21_eps10_50_off_reg_sign.H: 
#	< LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff.H Window3d squeeze=n min10=11000 n10=200 >pp0.H
#	< pp0.H Transp reshape=2,3,5,10 plane=13> pp1.H
#	< pp1.H Transp reshape=7,9,10 plane=23 > LtL_Sis_az21_ax21_wave_old_f35_flipoff_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 squeeze=n  | Transp reshape=2,5,6  | Scale scale_to=1. > pp3.H
	Math file1=pp3.H exp="-1.*file1" > pp4.H
	< LtL_Sis_az21_ax21_wave_old_f35_flipoff_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp4.H  niter=80 eps=.1 alpha=.001 na=1 reg_off=1 remov=rmov2.H mmov=mmov2.H wmov=wmov.H >$@


inv_Sis_az21ax21_eps10_50_ba.H: 
#	< LtL_Sis_az21_ax21_wave_old_fullf_ba.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 > pp4.H
#	< pp4.H Window3d f9=1 squeeze=n >LtL_Sis_az21_ax21_wave_old_f35_ba_t.H
#	< mig_Sis_Parallel_pad_800.H Window3d min1=31075 n1=79 min3=0 n3=1 min5=7975 | Pad n2=383 |Window3d f1=1| Transp | Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=31075 n4=79 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_f35_ba_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  niter=40 eps=0.01 alpha=.001 na=1 reg_off=0 remov=rmovsis.H mmov=mmovsis.H wmov=wmov.H>$@

test_sym.H: 
	< LtL_Sis_az21_ax21_wave_off33_marine_test_old.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 > LtL_Sis_az21_ax21_wave_old_t.H
	< LtL_Sis_az21_ax21_wave_old_t.H Window3d n4=1 min4=0  squeeze=n |../NoStaDecon_save/Bin/LINUX86_64/test_symmetry.x m=pp3.H >$@


inv_Sis_az21ax21_eps10_50_newsolver_test.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver.x mig=pp3.H  eps=.0 niter=3>$@

inv_Sis_az21ax21_eps10_50_newsolver_test_prec.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_old_t_weigh.H  eps=.0 niter=4>$@

inv_Sis_az21ax21_eps10_50_newsolver_diag.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H Window3d n4=1 min4=0 n3=1 min3=0 squeeze=n |../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_old_t_weigh.H eps=.0 niter=1>$@

inv_Sis_az21ax21_eps10_50_test.H: 
	< LtL_Sis_az21_ax21_wave_off33_marine_test_old.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 > LtL_Sis_az21_ax21_wave_old_t.H
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  eps=.0 niter=3>$@

test_sym_bin.H:
#	< test_sym.H Window3d n2=2 | Transp > headers.H
#	< test_sym.H Window3d n2=1 f2=2 | Transp > values.H
#	< values.H Create3d headers=headers.H keyname1='row' keyname2='col' verb=0> t$@
#	< t$@ Headermath key1=off eqn1=(row-col) verb=1 > t1$@
#	< t1$@ Sort3d nkeys=2 key1=col ng1=50939 dg1=1 og1=1 key2=row ng2=50939 dg2=1 og2=1 verb=1 max_memory=400 > t2$@
	< t2$@ Infill3d maxsize=1000 verb=1 >$@
 
#	< t1$@ Sort3d nkeys=2 key1=col ng1=1900 dg1=1 og1=1 key2=off ng2=21 dg2=1 og2=-10 verb=1> t2$@
#	< t1$@ Sort3d nkeys=2 key1=col ng1=4000 dg1=1 og1=1 key2=off ng2=2020 dg2=1 og2=-1010 verb=1> t2$@

inv_Sis_az21ax21_eps10_50_newsolver.H: 
#	< LtL_Sis_az21_ax21_wave_old_fullf_ba.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 > pp4.H
#	< pp4.H Window3d f9=1 squeeze=n >LtL_Sis_az21_ax21_wave_old_f35_ba_t.H
#	< mig_Sis_Parallel_pad_800.H Window3d min1=31075 n1=79 min3=0 n3=1 min5=7975 | Pad n2=383 |Window3d f1=1| Transp | Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=31075 n4=79 min2=0 n2=1 min1=7975 | Pad n1=383 |Window3d f2=1| Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_f35_ba_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver.x mig=pp3.H  eps=.0 niter=20>$@


inv_Sis_az21ax21_eps10_50_ba_weigh.H: 
#	< LtL_Sis_az21_ax21_wave_old_fullf_ba.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 > pp4.H
#	< pp4.H Window3d f9=1 squeeze=n >LtL_Sis_az21_ax21_wave_old_f35_ba_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=31075 n1=79 min3=0 n3=1 min5=7975 | Pad n2=383 |Window3d f1=1| Transp | Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< vel_pad.H Window3d  min3=7975 n1=78 min1=31175 | Pad n2out=383 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.4 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=20 rect2=20|Clip chop=less clip=.5 > ppp.H
	Math file1=ppp.H exp="file1-0.5"| Scale scale_to=1.| Transp reshape=2,4 >resweigh.H
	< LtL_Sis_az21_ax21_wave_old_f35_ba_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  niter=30 eps=.5 alpha=.001 na=1 reg_off=0 resw=resweigh.H remov=rmovsis.H mmov=mmovsis.H wmov=wmov.H>$@

inv_Sis_az21ax21_eps10_50_ba1.H: 
#	< LtL_Sis_az21_ax21_wave_old_fullf_ba.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 > pp4.H
#	< pp4.H Window3d f9=1 squeeze=n >LtL_Sis_az21_ax21_wave_old_f35_ba_t.H
	< mig_Sis_Parallel_pad_800.H Window3d min1=31075 n1=79 min3=0 n3=1 min5=7975 | Pad n2=383 |Window3d f1=1| Transp | Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< LtL_Sis_az21_ax21_wave_old_f35_ba_t.H ../NoStaDecon/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  niter=70 eps=3. alpha=.001 na= reg_off=0 remov=rmov.H mmov=mmov.H wmov=wmov.H>$@ 

#mig_Sis_off17_pad_800_ang.H:
#	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=7975 j3=2 squeeze=n |Pad n5out=383 | Transp reshape=2,5,6  | Scale scale_to=1. | Window3d |Transp reshape=1,3,4 plane=23 > pp3r.H
#	Math file1=pp3r.H exp="-1.*file1" > pp4.H
#	< pp4.H Transp reshape=2,3,4 |  Interp d2out=25 type=1> pp4_off.H
#	< pp4_off.H ./OFF2ANG verb=1 amax=45 nzpad=1024> $@

division.H:
	< LtL_Sis_az21_ax21_wave_old_f35_ba_t.H Window3d n4=1 min4=0 n3=1 min3=0 f9=1 > illu0.H
	Window3d < pp3.H > mig0.H
	Math file1=mig0.H file2=illu0.H exp="file1/(file2)" >$@

mig_Sis_off17_pad_800_ang_all.H:
	< mig_Sis_Parallel_pad_800.H Window3d min3=-600 n3=17 squeeze=n | Transp reshape=2,5,6  > pp4.H
	< pp4.H Transp reshape=2,3,4 |  Interp d2out=50 type=1 | Pad n2out=100 > pp4_off_all.H
	< pp4_off_all.H ./OFF2ANG verb=1 amax=90  > $@

mig_Sis_off17_pad_800_ang.H:
	< mig_Sis_Parallel_pad_800.H Window3d min1=32000 n1=70 min5=11000 n5=200 min3=-600 n3=17 squeeze=n | Transp reshape=2,5,6  > pp4.H
	< pp4.H Transp reshape=2,3,4 |  Interp d2out=50 type=1 | Pad n2out=100 > pp4_off.H
	< pp4_off.H ./OFF2ANG verb=1 amax=45 nzpad=512 > $@

inv_Sis_az21ax21_eps10_50_ang_reg.H:
	< inv_Sis_az21ax21_eps10_50_off_reg.H Transp reshape=2,3,4 |  Interp d2out=25 > inv_Sis_az21ax21_eps10_50_off_reg_1.H
	< inv_Sis_az21ax21_eps10_50_off_reg_1.H ./OFF2ANG amax=45 verb=1 nzpad=1024> $@

inv_Sis_az21ax21_ang_marine_new_reg.H:
	< inv_Sis_az21ax21_off_marine_new_reg.H Window3d  squeeze=n | Transp reshape=2,3,4 |  Interp d2out=50 | Pad n2out=100> ppreg.H
	< ppreg.H ./OFF2ANG verb=1 amax=45 nzpad=512 > $@

inv_Sis_az21ax21_ang_marine_new.H:
	< inv_Sis_az21ax21_off_marine_new.H Window3d  squeeze=n | Transp reshape=2,3,4 |  Interp d2out=50 | Pad n2out=100> ppnoreg.H
	< ppnoreg.H ./OFF2ANG verb=1 amax=45 nzpad=512 > $@ 

inv_Sis_az21ax21_eps10_50_ang.H:
	< inv_Sis_az21ax21_eps10_50_off.H Reverse  which=1 opt=ignore |Transp reshape=2,3,4 |  Interp d2out=25 type=1> inv_Sis_az21ax21_eps10_50_off1.H
	<  inv_Sis_az21ax21_eps10_50_off1.H ./OFF2ANG verb=1 amax=45 nzpad=512> $@

compare_ang_w:
	Window3d < mig_Sis_off17_pad_800_ang.H | Transp plane=23 | Grey eout=1 gainpanel=a pclip=98> mig_Sis_off17_pad_800_ang_w.T
	Window3d < inv_Sis_az21ax21_ang_marine_new_reg.H | Transp plane=23 | Grey eout=1 gainpanel=a pclip=98> inv_Sis_az21ax21_ang_marine_new_reg.T
	Window3d < inv_Sis_az21ax21_ang_marine_new.H | Transp plane=23 | Grey eout=1 gainpanel=a pclip=98> inv_Sis_az21ax21_ang_marine_new.T
	Rickmulti mig_Sis_off17_pad_800_ang_w.T inv_Sis_az21ax21_ang_marine_new_reg.T inv_Sis_az21ax21_ang_marine_new.T 


compare_ang:
	< mig_Sis_off17_pad_800_ang.H Window3d | Transp plane=23 | Grey eout=1 gainpanel=a pclip=90 > pp4_ang.T
	< inv_Sis_az21ax21_eps10_50_ang.H Window3d | Transp plane=23 | Grey eout=1 gainpanel=a pclip=90 > inv_Sis_az21ax21_eps10_50_ang.T
	< inv_Sis_az21ax21_eps10_50_ang_reg.H Window3d | Transp plane=23 | Grey eout=1 gainpanel=a pclip=90 > inv_Sis_az21ax21_eps10_50_ang_reg.T
	Rickmulti pp4_ang.T inv_Sis_az21ax21_eps10_50_ang.T inv_Sis_az21ax21_eps10_50_ang_reg.T

compare_off:
	< pp4_off.H Window3d | Transp plane=23 | Grey eout=1 gainpanel=a > pp4_off.T
	< ppreg.H  Window3d| Transp plane=23 | Grey eout=1 gainpanel=a  > ppreg.T
	< ppnoreg.H  Window3d| Transp plane=23 | Grey eout=1 gainpanel=a  > pp.T
	Rickmulti pp4_off.T ppreg.T pp.T

fragreyo1=eout=1  gainpanel=a  label1="z (ft)" label2="x (ft)" label3="h (ft)" title=" "     gainpanel=a
fragreyo2=eout=1  gainpanel=a  label1="z (ft)" label2="x (ft)" label3="angle (degree)" title=" "   gainpanel=a
frames1= frame1=160 frame2=8 frame3=12 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=1.8 point1=.85 point2=.75 
frames11= frame1=160 frame2=8 frame3=12 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=.6 point1=1 point2=.75 
frames1a= frame1=160 frame2=68 frame3=16 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=1. point1=1 point2=.75 
frames1illu= frame1=160 frame2=8 frame3=8 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=.6 point1=1 point2=.75 
frames2= frame1=160 frame2=8 frame3=1 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=.6 point1=1 point2=.75 
frames2illu= frame1=160 frame2=8 frame3=1 title=" "  labelfat=5 labelsz=6 crowd=.65 screenratio=.6 point1=1 point2=.75

${RESDIR}/compare_off1.v:
	Window3d <pp4_off.H n2=25 | Transp plane=23 |   Grey ${fragreyo1} |Cubeplot ${frames11}  out=$@ > /dev/null

${RESDIR}/compare_off_art.v:
	 Window3d < mig_Sis_Parallel_pad_800.H  min5=7975 min1=31075 n1=79 | Transp reshape=2,3 |   Grey ${fragreyo1} pclip=95  |Cubeplot ${frames1a}  out=$@ > /dev/null

${RESDIR}/compare_ang_art.v:
	< mig_Sis_off17_pad_800_ang_all.H Window3d min4=31075 n4=79 min1=7975  | Transp plane=23 |   Grey ${fragreyo1} pclip=95 |Cubeplot ${frames1a}  frame3=50 out=$@ > /dev/null

${RESDIR}/compare_off2.v:
	Window3d < ppreg.H  n2=25 | Transp plane=23 |   Grey ${fragreyo1} |Cubeplot ${frames11}  out=$@ > /dev/null

${RESDIR}/compare_off3.v:
	Window3d < ppnoreg.H n2=25 | Transp plane=23 |   Grey ${fragreyo1} |Cubeplot ${frames11}  out=$@ > /dev/null

${RESDIR}/illu_off.v:
	Window3d < LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_diag_off.H n2=17| Transp plane=23 |   Grey ${fragreyo1} allpos=y|Cubeplot ${frames1illu} color=j  out=$@ > /dev/null

fragreyfo=eout=1  gainpanel=a  label1="az (ft)" label2="ax (ft)" label3="h (ft)" title=" " labelfat=4 color=j pclip=100
framesof=frame1=15 frame2=8 frame3=8 flat=1 title=" " polarity=-1 color=j labelfat=5 titlefat=5 labelsize=12  titlesz=12

${RESDIR}/filt_off17.v:
	Window3d < LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H min6=0 n6=1 min9=33000 n8=1 min8=12000 n9=1 | Transp reshape=1,3 | Grey ${fragreyfo} |Cubeplot ${framesof} out=$@1 > /dev/null
	Window3d < LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H min6=0 n6=1 min9=34000 n8=1 min8=12000 n9=1 | Transp reshape=1,3 | Grey ${fragreyfo} |Cubeplot ${framesof} out=$@2 > /dev/null
	Window3d < LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H min6=0 n6=1 min9=35000 n8=1 min8=12000 n9=1 | Transp reshape=1,3 | Grey ${fragreyfo} |Cubeplot ${framesof} out=$@3 > /dev/null
	Window3d < LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_t.H min6=0 n6=1 min9=36000 n8=1 min8=12000 n9=1 | Transp reshape=1,3 | Grey ${fragreyfo} |Cubeplot ${framesof} out=$@4 > /dev/null
	vp_Movie $@1 $@2 $@3 $@4 > $@_movie
	< $@1 vppen vpstyle=n yscale=2  > $@11
	< $@2 vppen vpstyle=n yscale=2  > $@22
	< $@3 vppen vpstyle=n yscale=2  > $@33
	< $@4 vppen vpstyle=n yscale=2  > $@44 
	vp_SideBySideIso $@11 $@22  > t.v 
	vp_SideBySideIso $@33 $@44  > b.v
	vp_OverUnderAniso t.v b.v | vppen vpstyle=n txscale=3 tyscale=3 > $@ 


${RESDIR}/compare_ang1.v:
	Window3d < mig_Sis_off17_pad_800_ang.H min2=0| Transp plane=23 |   Grey ${fragreyo2} |Cubeplot ${frames2}  out=$@ > /dev/null

${RESDIR}/compare_ang2.v:
	Window3d < inv_Sis_az21ax21_ang_marine_new_reg.H min2=0| Transp plane=23 |   Grey ${fragreyo2} |Cubeplot ${frames2}  out=$@ > /dev/null

${RESDIR}/compare_ang3.v:
	Window3d < inv_Sis_az21ax21_ang_marine_new.H min2=0| Transp plane=23 |   Grey ${fragreyo2} |Cubeplot ${frames2}  out=$@ > /dev/null

${RESDIR}/illu_ang.v:
	Window3d <LtL_Sis_az21_ax21_wave_old_off33_f35_flipoff_marine_new_diag_ang.H min2=0| Transp plane=23 |   Grey ${fragreyo2} allpos=y|Cubeplot ${frames2illu} color=j  out=$@ > /dev/null
#----------------------------------------------------------

#----------------------------------------------------------


GREYP= Grey  label1="az (ft)" label2="ax (ft)" label3=" " labelsz=12 titlesz=12 crowd2=.7 color=j  \
	    grid=y > /dev/null 
GREYP1= Grey  label1="az (ft)" label2="ax (ft)" label3=" " labelsz=12 titlesz=12 crowd2=.7 color=j  \
	    grid=y gridfat=10 screenratio=.33 labelfat=5 titlefat=5 framenumfat=5 title=" "> /dev/null
GREYS= Grey  label1="z (ft) " label2="x (ft)" label3=" "  title=" "  \
	  title=" " labelsz=12 titlesz=12   > /dev/null \

GREYS1= Grey  label1="z (ft) " label2="x (ft)" label3=" "  labelsz=8 title=" "  \
		 title=" "  > /dev/null \

#		labelsz=14 titlesz=14  crowd1=.7  crowd2=.6
GREYC2= Grey   pclip=100 label1="z (ft) " label2="x (ft)" label3=" "  \
	    title=" " > /dev/null \
		labelsz=8
# titlesz=14   
GREYV2= Grey allpos=y  pclip=100 label1="z (ft) " label2="x (ft)" label3="velociy ft/s" labelsz=8 titlesz=8 \
	    wantscalebar=y  crowd1=.75  crowd2=.75 bias=4920 title=" "  > /dev/null
GREYSS= Grey label1="z (ft) " label2="h (ft)"   title=" " labelsz=16 titlesz=16 d1num=300 o1num=16000 crowd=.65 #xll=3 yll=1.2 xur=10 yur=9   crowd=.85   
frames= frame1=0 frame2=100 frame3=0 flat=1 point1=0.668539 point2=0.935094 titlesz=8 wantcoordlabel=-1 crowd=.65
 
${RESDIR}/Sis_vel.v:
	< ${SIGSMULT_DATA}/sigsbee2b_migration_velocity.H Window> j.H
	echo "hff=-1 gff=-1 d1=25 o2=10025 d2=37.5"      >>j.H
	Window3d < j.H  min2=10925 max2=65000 min1=6000 max1=17525| ${GREYV2}  out=$@
#	< pp.v vp_annotate batch=y text=${TXTDIR}/text_Sis  | vppen vpstyle=n   > $@
#	Rm j.H; \rm pp.v 

${RESDIR}/hesian_phase_Sis.v: 
	cp 	LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H LtL_1.H
	Window3d n8=1 min8=12000 n9=1 min9=32000 n1=1 n6=1  < LtL_1.H | ${GREYP} title="(a)"  out=1.v
	Window3d n8=1 min8=12000 n9=1 min9=33500 n1=1 n6=1  < LtL_1.H | ${GREYP} title="(b)"  out=2.v
	Window3d n8=1 min8=12000 n9=1 min9=35000 n1=1 n6=1  < LtL_1.H | ${GREYP} title="(c)"  out=3.v
	Window3d n8=1 min8=12000 n9=1 min9=36500 n1=1 n6=1  < LtL_1.H | ${GREYP} title="(d)"  out=4.v
	vp_SideBySideAniso 1.v 2.v | vppen vpstyle=n txscale=1.8 tyscale=1.8 > top.v
	vp_SideBySideAniso 3.v 4.v | vppen vpstyle=n txscale=1.8 tyscale=1.8 > bot.v
	vp_OverUnderAniso top.v bot.v > $@
	\rm 1.v 2.v 3.v 4.v  top.v bot.v

${RESDIR}/hesian_phase_Sis_movie.v: 
	cp 	LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H LtL_1.H
	Window3d n8=1 min8=11000 n9=1 min9=31750 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=0.v
	Window3d n8=1 min8=11000 n9=1 min9=32250 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=1.v
	Window3d n8=1 min8=11000 n9=1 min9=32750 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=2.v
	Window3d n8=1 min8=11000 n9=1 min9=33250 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=3.v
	Window3d n8=1 min8=11000 n9=1 min9=33750 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=4.v
	Window3d n8=1 min8=11000 n9=1 min9=34250 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=5.v
	Window3d n8=1 min8=11000 n9=1 min9=34750 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=6.v
	Window3d n8=1 min8=11000 n9=1 min9=35250 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=7.v
	Window3d n8=1 min8=11000 n9=1 min9=35750 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=8.v
	Window3d n8=1 min8=11000 n9=1 min9=36250 n1=1 n6=1  < LtL_1.H | ${GREYP1}   out=9.v
	vp_Movie 0.v 1.v 2.v 3.v 4.v 5.v 6.v 7.v 8.v 9.v > $@
	\rm 0.v 1.v 2.v 3.v 4.v 5.v 6.v 7.v 8.v 9.v

${RESDIR}/hesian_phase_Sis_movie.gif:
	vplot2gif ${RESDIR}/hesian_phase_Sis_movie.v $@

${RESDIR}/mig_Sis.v: mig_Sis_Parallel_pad_800.H
	Window3d <mig_Sis_Parallel_pad_800.H min3=0 n3=1| Transp | ${GREYS1} pclip=97 out=$@ > /dev/null

${RESDIR}/comp_Sis_full.v: 
	< reflection_coefficients.H  Window3d min1=7975 n1=383 min2=31000 n2=80 j2=3 | ${GREYS}  out=0.v pclip=97 title="(a)"
	Window3d n3=1 min3=0 n4=1 min4=0  < LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H | ${GREYC2} out=00.v  allpos=y title="(b)"
	Window3d < mig_Sis_Parallel_pad_800.H min3=0 n3=1 min5=7975 min1=31075 n1=79| Pad n2out=383 |Transp| ${GREYS} pclip=90   out=1.v title="(c)"
	Window3d <   inv_Sis_az21ax21_eps10_50_ba_weigh.H | ${GREYS} pclip=90  out=2.v title="(d)"
	vp_SideBySideAniso 0.v 00.v |  vppen vpstyle=n> 000.v
	vp_SideBySideAniso 1.v  2.v |  vppen vpstyle=n> 12.v
	vp_OverUnderAniso 000.v 12.v | vppen vpstyle=n > $@

${RESDIR}/comp_Sis_full1.v: 
	< reflection_coefficients.H  Window3d min1=7975 n1=383 min2=31000 n2=80 j2=3 | ${GREYS3} out=0.v pclip=97 title="(a)"
	Window3d n3=1 min3=0 n4=1 min4=0  < LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H | ${GREYC3} out=00.v  allpos=y title="(b)"
	vp_OverUnderAniso 0.v 00.v |  vppen vpstyle=n xscale=.5 >  $@

${RESDIR}/strati.v: 
	< reflection_coefficients.H  Window3d min1=7975 n1=383 min2=31000 n2=79 j2=3 > pp.H
	echo  o2=31075 >>pp.H
	< pp.H ${GREYS1} pclip=97 polarity=-1 o1num=9000 out=pp.v
	<  pp.v vppen vpstyle=n xscale=.5 > $@

${RESDIR}/illu.v: 
	Window3d n3=1 min3=0 n4=1 min4=0  min9=31075< LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H | ${GREYC2} out=pp.v  allpos=y color=j o1num=9000
	<  pp.v vppen vpstyle=n xscale=.5 > $@

${RESDIR}/illu_pre.v: 
	Window3d n3=1 min3=0 n4=1 min4=0  min9=31075< LtL_Sis_az21_ax21_wave_old_fullf_ba_t.H | ${GREYC2} out=pp.v  allpos=y color=j o1num=9000
	<  pp.v vppen vpstyle=n xscale=.5 > $@

${RESDIR}/mig_Sis_pos.v: 
	Window3d < mig_Sis_Parallel_pad_800.H min3=0 n3=1 min5=7975 min1=31175 n1=79| Pad n2out=383 |Transp> pp.H
	echo o2=31075 >> pp.H
	<pp.H ${GREYS1} pclip=95   out=pp.v o1num=9000
	<  pp.v vppen vpstyle=n xscale=.5 > $@

${RESDIR}/inv_Sis_pos.v: 
	Window3d <  inv_Sis_az21ax21_eps10_50_ba_weigh.H  | ${GREYS1} pclip=95  out=pp.v o1num=9000
	<  pp.v vppen vpstyle=n xscale=.5 > $@
  


${RESDIR}/inv_Sis_pos_movie.v: 
	Transp <  mmovsis.H reshape=2,4,6| Window3d j6=2| ${GREYS1} screenratio=1.7 out=$@ o1num=9000 labelfat=5 titlefat=5 framenumfat=5

${RESDIR}/inv_Sis_pos_movie.gif:
	vplot2gif ${RESDIR}/inv_Sis_pos_movie.v  $@

${RESDIR}/compare_inv_pp.v: 
	Window3d < mig_Sis_Parallel_pad_800.H min3=0 n3=1 min5=11000 n5=200  min1=32075 n1=67 | Transp > mig_Sis_Parallel_pad_800_comp.H
	< mig_Sis_Parallel_pad_800_comp.H ${GREYS1} gainpanel=e out=${RESDIR}/compare_inv_pp1.v > /dev/null
	< inv_Sis_az21ax21_off_marine_new_reg.H Window3d n1=1 min1=0 f4=1 n4=67 > inv_Sis_az21ax21_off0_marine_new_reg.H
	< inv_Sis_az21ax21_off0_marine_new_reg.H ${GREYS1}  gainpanel=e out=${RESDIR}/compare_inv_pp2.v > /dev/null
	Window3d < inv_Sis_az21ax21_eps10_50_ba.H min3=11000 max3=16000 min4=32000 > inv_Sis_az21ax21_eps10_50_ba_comp.H
	echo o2=32075 >> inv_Sis_az21ax21_eps10_50_ba_comp.H
	<  inv_Sis_az21ax21_eps10_50_ba_comp.H ${GREYS1}  gainpanel=e out=${RESDIR}/compare_inv_pp3.v > /dev/null
	< reflection_coefficients.H  Window3d min1=11025 n1=200 min2=32150 n2=67 j2=3 > pp.H
	Math file1=pp.H exp="-1*file1"  >pp1.H
	Cat3d pp1.H mig_Sis_Parallel_pad_800_comp.H inv_Sis_az21ax21_eps10_50_ba_comp.H inv_Sis_az21ax21_off0_marine_new_reg.H  axis=3 | ${GREYS1} gainpanel=e out=$@ > /dev/null
 

${RESDIR}/comp_Sis_full2.v: 
	Window3d < mig_Sis_Parallel_pad_800.H min3=0 n3=1 min5=7975 min1=31075 n1=79| Pad n2out=383 |Transp| ${GREYS3} pclip=90   out=1.v title="(a)"
	Window3d <  inv_Sis_az21ax21_eps10_50_ba.H  | ${GREYS3} pclip=90  out=2.v title="(b)"
	vp_OverUnderAniso 1.v  2.v |  vppen vpstyle=n xscale=.5> $@
 
${RESDIR}/comp_off_mig.v: mig_Sis_Parallel_16000_w_stack_t.H
    Window3d < mig_Sis_Parallel_16000_w_stack_t.H n1=5 j1=2 min1=0. min3=16000 n3=61 min4=38775 n4=230 squeeze=n| Scale scale_to=1.> ppp.H
	Window3d < ppp.H min4=45000 n4=1| Transp | ${GREYSS} title="(a)"  out=1.v >/dev/null
	Window3d < ppp.H min4=50000 n4=1| Transp | ${GREYSS} title="(b)"  out=2.v >/dev/null
	Window3d < ppp.H min4=55000 n4=1| Transp | ${GREYSS} title="(c)"  out=3.v >/dev/null
	vp_SideBySideAniso 1.v 2.v 3.v| vppen  vpstyle=n  > $@


${RESDIR}/comp_off_inv_damp.v: inv_Sis_az21ax7zo_off_eps10_50.H
	Window3d < inv_Sis_az21ax7zo_off_eps10_50.H min4=45000 n4=1| Transp | ${GREYSS} title="(a)"  out=1.v >/dev/null
	Window3d < inv_Sis_az21ax7zo_off_eps10_50.H min4=50000 n4=1| Transp | ${GREYSS} title="(b)"  out=2.v >/dev/null
	Window3d < inv_Sis_az21ax7zo_off_eps10_50.H min4=55000 n4=1| Transp | ${GREYSS} title="(c)"  out=3.v >/dev/null
	vp_SideBySideAniso 1.v 2.v 3.v| vppen  vpstyle=n  > $@

${RESDIR}/comp_off_inv_semb.v: inv_Sis_az21ax7Diff_off_eps10_50.H
	Window3d < inv_Sis_az21ax7Diff_off_eps10_50.H min4=45000 n4=1| Transp | ${GREYSS} title="(a)"  out=1.v >/dev/null
	Window3d < inv_Sis_az21ax7Diff_off_eps10_50.H min4=50000 n4=1| Transp | ${GREYSS} title="(b)"  out=2.v >/dev/null
	Window3d < inv_Sis_az21ax7Diff_off_eps10_50.H min4=55000 n4=1| Transp | ${GREYSS} title="(c)"  out=3.v >/dev/null
	vp_SideBySideAniso 1.v 2.v 3.v| vppen  vpstyle=n  > $@


fragreyo=eout=1  gainpanel=a  label1="z (feet)" label2="x (feet)" label3="h (feet)" title=" " labelfat=4 pclip=99
frames= flat=1  title=" " 

${RESDIR}/mig_inv1.v:
	Window3d < mmov3.H n6=1 f6=9 | Transp reshape=1,3 |Grey ${fragreyo} |Cubeplot flat=1 point1=.8 point2=.8 frame1=157 frame2=22 frame3=8 title=" " out=$@ > /dev/null

${RESDIR}/mig_off.v:
	Window3d < pp3.H | Transp reshape=1,3 |Grey ${fragreyo} |Cubeplot flat=1 point1=.8 point2=.8 frame1=157 frame2=22 frame3=8 title=" " out=$@ > /dev/null

  gs= gs -dBATCH -sPAPERSIZE=letter -r200x200 -dEPSCrop -sDEVICE=jpeg
${RESDIR}/%.jpg:${RESDIR}/%.ps
        ${gs} -sOutputFile=$@ ${RESDIR}/$*.ps
                                                                                 
 
clean: jclean
	rm -f ./*.stb ./*fix*
	rm -f ./Obj/LINUX/*.o ./Obj/SGI64/*.o
	rm -f ./Bin/LINUX/*.x ./Bin/SGI64/*.x
	rm -f ./Lib/LINUX/* ./Lib/SGI64/*
	rm -f ./Src/*.o ./Src/*.stb ./Src/*.mod ./Src/*fix*

include ${SEPINC}/SEP.bottom