include ${SEPINC}/SEP.top
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
BINDIR = ./Bin/${MTYPE}
OBJDIR = ./Obj/${MTYPE}
SRCDIR = ./Src8-31
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  


sigsbee2b_stratigraphy.H:
	Segy2sep tape=${SIGSMULT_DATA}/Tape/sigsbee2b_stratigraphy.segy par=${SIGSMULT_DATA}/loadstrat.p > sigsbee2b_stratigraphy.H

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$@ /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*$@

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 >> $@
#	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

vel_strati_pad.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" " >> $@ 

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

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

mig_Sis_refle_ang_f35_marine.H:
	< mig_Sis_refle_off_f35_marine.H Window3d f4=1 | Transp | Interp d2out=25 type=1 | Pad beg2=24 n2=97 n3=133> pppp.H
	Math file1=pppp.H exp="-1.0*file1" > pppp1.H
	< pppp1.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@

inv_Sis_az11ax21_ang_oldsolver_sym_mod.H:
	<inv_Sis_az11ax21_off_oldsolver_sym_mod.H Window3d f4=1 | Transp | Interp d2out=25 type=1 | Pad beg2=24 n2=97 n3=133> pppp.H
	Math file1=pppp.H exp="-1.0*file1" > pppp1.H
	< pppp1.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@


mig_Sis_refle_ang_f35_marine_test.H:
	< mig_Sis_refle_off_f35_marine_test.H Window3d f4=1 | Transp | Interp d2out=25 type=1 | Pad beg2=24 n2=97 n3=133> pppp.H
	Math file1=pppp.H exp="-1.0*file1" > pppp1.H
	< pppp1.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@
 
mig_Sis_refle_ang_f35_marine_strati.H:
	< mig_Sis_refle_off_f35_marine_strati.H Window3d f4=1 | Transp | Interp d2out=25 type=1 | Pad beg2=24 n2=97 n3=133> pppp.H
	Math file1=pppp.H exp="-1.0*file1" > pppp1.H
	< pppp1.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@

mig_Sis_pad_800_ang.H:
	< mig_Sis_pad_800.H Window3d min1=30000 n1=133 min5=7975 n5=381 |Transp plane=13 | Interp d2out=25 type=1> pppp.H
	< pppp.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@

mig_Sis_pad_2048_strati.H:  data_f.H wave.H vel_pad.H 
	@date
	nohup ../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=25 min_slow_dev=1. niter_lloyd=100 

mig_Sis_pad_2048_strati_ang.H:
	< mig_Sis_pad_2048_strati.H Window3d min1=30000 n1=133 min5=7975 n5=381|Transp plane=13 | Interp d2out=25 type=1> pppp.H
	< pppp.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@

inv_Sis_az11ax21_oldsolver_sym_strati_ang_reg.H:
	< inv_Sis_az11ax21_oldsolver_sym_strati_off_reg.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-10 dax=1. naax=71 style=1 adj=1 cosw=1> $@

mig_Sis_pad_800_strati_nointerp_ang.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> $@

mig_Sis_pad_800_strati_nointerp_ang_cow.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_pad_800_strati_nointerp_ang_inv.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 inv=1 niter=8 eps=0.01> $@

mig_Sis_pad_800_strati_nointerp_off_inv.H:
	< mig_Sis_pad_800_strati_nointerp_ang_inv.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x nhx=17 hxmin=-600 dhx=75 style=1 inv=1 niter=10 eps=0.01> $@ 

mig_Sis_pad_800_strati_nointerp_off.H:
	<mig_Sis_pad_800_strati_nointerp_ang.H  ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x adj=0 style=1 nhx=17 hxmin=-600 dhx=75> $@

mig_ang_unitary.H:
	Wavelet n1=100 wavelet=ricker2 tdelay=.05 fund=20> wave_test.H
	 Spike n1=17 d1=25 o1=-200 d2=25 n2=300 o2=0. n3=1. d3=1. o3=0. d4=1. o4=0.  d5=1. o5=0. n4=1 n5=1 k1=9 k2=150 k3=1 k4=1 k5=1| Smooth rect1=2 rect2=2 repeat=5 | Transp > pp6.H
	< pp6.H Filter filter=wave_test.H | Transp reshape=1,3| Scale scale_to=1. > pp5.H
	< pp5.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-80 dax=1 naax=161 style=1 adj=1> $@
mig_ang_unitary_cow.H:
	Wavelet n1=100 wavelet=ricker2 tdelay=.05 fund=20> wave_test.H
	 Spike n1=17 d1=25 o1=-200 d2=25 n2=300 o2=0. n3=1. d3=1. o3=0. d4=1. o4=0.  d5=1. o5=0. n4=1 n5=1 k1=9 k2=150 k3=1 k4=1 k5=1| Smooth rect1=2 rect2=2 repeat=5 | Transp > pp6.H
	< pp6.H Filter filter=wave_test.H | Transp reshape=1,3 | Scale scale_to=1.> pp5.H
	< pp5.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-80 dax=1 naax=161 style=1 adj=1 cosw=1> $@

mig_off_unitary.H:
	< mig_ang_unitary.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x nhx=17 hxmin=-200 dhx=25 style=1 adj=0> $@
mig_off_unitary_cow.H:
	< mig_ang_unitary_cow.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x nhx=17 hxmin=-200 dhx=25 style=1 adj=0 cosw=1> $@

mig_inv_off_unitary.H:
	< mig_ang_unitary_inv.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x nhx=17 hxmin=-200 dhx=25 style=1 adj=0> $@
mig_ang_unitary_inv.H:
	Wavelet n1=100 wavelet=ricker2 tdelay=.05 fund=20> wave_test.H
	 Spike n1=17 d1=25 o1=-200 d2=25 n2=300 o2=0. n3=1. d3=1. o3=0. d4=1. o4=0.  d5=1. o5=0. n4=1 n5=1 k1=9 k2=150 k3=1 k4=1 k5=1| Smooth rect1=2 rect2=2 repeat=5 | Transp > pp6.H
	< pp6.H Filter filter=wave_test.H | Transp reshape=1,3 | Scale scale_to=1.> pp5.H
	< pp5.H  ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x amin=-80 dax=1 naax=161 style=1 inv=1 eps=1. niter=10 > $@

mig_off_unitary_inv.H:
	< mig_ang_unitary_inv.H ../NoStaDecon_save/Bin/LINUX86_64/off2ang_sd.x nhx=17 hxmin=-200 dhx=25 inv=1 eps=0.2 style=1 tan_max=5.76 adj=0> $@

#	Gauss n1=300 n2=17 var=3. > gauss.H
#	< gauss.H Transp reshape=1,3 > pp5.H
#	echo d1=25 o1=-200 d3=25 d4=1. o4=0.  d5=1. o5=0.>> pp5.H 




Imagertmax.H:
	<Imagertmhx.H Window3d n3=400 | ~claudio/Research/Off2ang/Yaxun/Bin/LINUX/off2ang_sd.x style=1 mode=1 niter=5> $@


mig_Sis_pad_800_strati.H:  data_f.H wave.H vel_pad.H 
	@date
	nohup ../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 

mig_Sis_pad_800_strati_ang.H:
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=7975 n5=381 |Transp plane=13 | Interp d2out=25 type=1> pppp.H
	< pppp.H ../Yaxun/Bin/LINUX86_64/off2ang_sd.x amin=-80 amax=80 style=1 > $@ 

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_strati_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_strati_pad.H W_d=green_Sis_strati_OMP.H  \
	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 
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

diagonal_marine_old_off17:
	Window3d <LtL_Sis_az21_ax21_wave_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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_marine_old_off17_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 \
n3d 	      n10.H n11.H n12.H n13.H n14.H n15.H n16.H  axis=3 > LtL_Sis_az21_ax21_wave_marine_old_off17_t_diag.H
	Rm n*.H

diagonal_17off:
	Window3d <LtL_Sis_az21_ax21_wave_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_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_marine_strati_off_t_diag.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_az1_ax1_wave_marine_old_off.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=1 naz=1 nay=1 marine=y noffx=33 \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=32000 nx=50
	date

LtL_Sis_az21_ax21_wave_marine_old_off: 
	date
	< green_Sis_OMP_f35.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$@
	nohup < t1$@ ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	nax=21 naz=11 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 #miniz=8225 nz=373
	date

LtL_Sis_az41_ax31_wave_marine_old_off: 
	date
	< green_Sis_OMP_f35.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$@ ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	nax=31 naz=21 nay=1 marine=y noffx=17 \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=30000 nx=133  miniz=7975 nz=381 #miniz=8225 nz=373
	date

LtL_Sis_az21_ax21_wave_marine_old_off_test: 
	date
#	< green_Sis_OMP_f35.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$@
	nohup < t1$@ ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	nax=21 naz=11 nay=1 marine=y noffx=17 \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=30000 nx=133 miniz=7975 nz=381 #miniz=8225 nz=373
	date

LtL_Sis_az21_ax21_wave_marine_old_off_test1: 
	date
	< t1LtL_Sis_az21_ax21_wave_marine_old_off_test Window3d n6=1 squeeze=n| ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	nax=1 naz=1 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=10000 nz=1 #miniz=8225 nz=373
	date
 

LtL_Sis_az21_ax21_wave_marine_strati_off: 
	date
	< green_Sis_strati_OMP_f35.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$@
	nohup < t1$@ ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	nax=21 naz=11 nay=1 marine=y noffx=17 \
	nry=1 nrx=348 orx=0 nsy=1 nsx=327 osx=10925 dsx=150 minix=30000 nx=133 miniz=7975 nz=381 #miniz=8225 nz=373
	date 


LtL_Sis_az41_ax41_wave_marine_strati: 
	date
	< green_Sis_strati_OMP_f35.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$@
	nohup < t1$@ ../Src_Hessian/Bin/LINUX86_64/Lt_L_old_3D_off_wave_6_omp.x outtag=$@.H  \
	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 #miniz=8225 nz=373
	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=1  \
	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_marine_old_off17.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_marine_old_off17_t.H
	echo d1=25 d2=25 o2=10000>> sigsbee2b_reflection_coefficients_77.H
	echo o8=8000 >> LtL_Sis_az21_ax21_wave_marine_old_off17_t.H
	< sigsbee2b_reflection_coefficients_77.H Window3d min1=8000 n1=383 min2=30000 n2=133 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Pad beg1=8 n1=17 | Window3d squeeze=n n3=381 > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=pp5.H >$@

mig_Sis_refle_off_f35_marine_test.H: 
#	< LtL_Sis_az21_ax21_wave_marine_old_off_test.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_marine_old_off_test_t.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 Pad beg1=8 n1=17 | Window3d squeeze=n > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_old_off_test_t.H ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=pp5.H >$@

mig_Sis_refle_off_f35_marine_strati.H: 
	< LtL_Sis_az21_ax21_wave_marine_strati_off.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_marine_strati_off_t.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 Pad beg1=8 n1=17 | Window3d squeeze=n > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=pp5.H >$@

mig_Sis_refle_ang2off_f35_marine_strati.H: 
#	< LtL_Sis_az21_ax21_wave_marine_strati_off.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_marine_strati_off_t.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 Pad beg1=8 n1=17  extend=1 > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/test_ang2off.x m=pp5.H amax=60 amin=-60 adj=0>$@

mig_Sis_refle_off2ang_f35_marine_strati.H: # tested
#	< LtL_Sis_az21_ax21_wave_marine_strati_off.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_marine_strati_off_t.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 Pad beg1=8 n1=17 | Window3d squeeze=n > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/test_ang2off.x m=pp5.H amax=60 amin=-60 adj=1>$@

mig_Sis_refle_off2ang_f35_marine_strati_pad.H: # tested
#	< LtL_Sis_az21_ax21_wave_marine_strati_off.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_marine_strati_off_t.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 Pad beg1=12 n1=25 | Window3d squeeze=n > pp5.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Pad beg6=4 end6=4 extend=1> pppp_pad.H
	< pppp_pad.H ../NoStaDecon_save/Bin/LINUX86_64/test_ang2off.x m=pp5.H amax=60 amin=-60 adj=1>$@



mig_Sis_refle_off2ang_reg_f35_marine_strati.H: 
#	< LtL_Sis_az21_ax21_wave_marine_strati_off.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_marine_strati_off_t.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 Pad beg1=8 n1=17 extend=1| Window3d squeeze=n > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/test_reg_ang2off.x m=pp5.H amax=60 amin=-60 adj=0>$@
 
mig_Sis_refle_off_f35_marine_test_21.H: 
#	< LtL_Sis_az21_ax21_wave_marine_old_off_test.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_marine_old_off_test_t.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 Pad beg1=8 n1=17 | Window3d squeeze=n > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_old_off_test_t.H ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=pp5.H >$@


spike_ir_off_f35_marine_half.H: 
#	< LtL_Sis_az21_ax21_wave_marine_old_off17.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_marine_old_off17_t.H
#	Spike n1=17 n2=1 n3=381 n4=133 k1=8 k2=1 k3=160 k4=67 > spike.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=spike.H >$@


mig_Sis_refle_f35_marine.H: 
#	< LtL_Sis_az21_ax21_wave_marine_old_off17.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_marine_old_off17_t.H
	echo d1=25 d2=25 o2=10000>> sigsbee2b_reflection_coefficients_77.H
	< sigsbee2b_reflection_coefficients_77.H Window3d min1=7975 n1=383 min2=30000 n2=133 j2=3| Transp reshape=2,4> reflection_coefficients_target_conv.H 
	< reflection_coefficients_target_conv.H Window3d squeeze=n min3=9200 n3=332 > pp5.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d n1=1 min1=0 n6=1 min6=0 squeeze=n |../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=pp5.H >$@

LtL_Sis_az21_ax21_wave_marine_old_off17_t_prec.H:
	Math file1=LtL_Sis_az21_ax21_wave_marine_old_off17_t_diag.H exp="1./(file1+0.01)" | Transp reshape=2,4| Transp plane=34 > $@

LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H:
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_diag.H exp="1./(file1)" | Transp reshape=2,4| Transp plane=34 > $@ 
#	< pp.H Clip  pclip=90  | Scale >$@

inv_Sis_az11ax21_off_newsolver.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min1=7975 n1=381 j2=3 min2=-600 n2=17 squeeze=n  > ppp.H
	< ppp.H /opt/SEP/SEP6.4.6/bin/Transp reshape=1,3,4,5,6  > ppp1.H
	< ppp1.H  Scale scale_to=1. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_marine_old_off17_t_prec.H niter=5 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov.H wmov=wmov.H >$@

inv_Sis_az11ax21_off_oldsolver_sym.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min1=7975 n1=381 j2=3 min2=-600 n2=17 squeeze=n  > ppp.H
	< ppp.H /opt/SEP/SEP6.4.6/bin/Transp reshape=1,3,4,5,6  > ppp1.H
	< ppp1.H  Scale scale_to=1. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3.H   niter=500 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_old_sym_off.H wmov=wmov.H >$@

inv_Sis_az11ax21_off_oldsolver_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3.H   niter=10 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_old_sym_off_strati.H wmov=wmov.H >$@

inv_Sis_az11ax21_oldsolver_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. | Pad n3=400 > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2"| Pad n3=400 extend=1 | Window3d squeeze=n n1=1 min1=0 > res_pre.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_t.H Pad extend=1 n8=400 > LtL_Sis_az21_ax21_wave_marine_strati_t_pad.H 
	< LtL_Sis_az21_ax21_wave_marine_strati_t_pad.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3.H initial=0 resw=res_pre.H niter=100 eps=.0 alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_old_sym.H wmov=wmov.H >$@ #

inv_Sis_az11ax21_oldsolver_sym_strati_off.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. | Pad n3=400 > ppp3.H
	Math file1=ppp3.H exp="10.*file1" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2"| Pad n3=400 extend=1 | Window3d squeeze=n  > res_pre.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3.H initial=0  niter=200 eps=0. alpha=.001 na=1 reg_off=0. resw=res_pre.H remov=rmov.H mmov=mmov_old_sym.H wmov=wmov.H >$@ #

inv_Sis_az11ax21_oldsolver_sym_strati_off_reg.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. | Pad n3=400 > ppp3.H
	Math file1=ppp3.H exp="10.*file1" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2"| Pad n3=400 extend=1 | Window3d squeeze=n  > res_pre.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3.H initial=0  niter=200 eps=.5 alpha=.001 na=1 reg_off=1. resw=res_pre.H remov=rmov.H mmov=mmov_old_sym_reg.H wmov=wmov.H >$@ #

inv_Sis_az11ax21_newsolver_sym_strati_off_reg.H: 
	< mig_Sis_pad_800_strati.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" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=121 min1=30600 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Window3d <  LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H  min4=30600 n4=121 squeeze=n > prec.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min9=30600 n9=121 squeeze=n > LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H 
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H /net/koko/homes/sep/alejandro/research/CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_off.x mig=ppp3.H initial=0  niter=40 eps=0. alpha=.001 na=1 reg_off=0 beta_max=.1 prec=prec.H  remov=rmov.H mmov=mmov_new_sym_reg1.H wmov=wmov.H >$@

inv_Sis_az11ax21_newsolver_sym_strati_off_reg_claudio.H: 
	< 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" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=121 min1=30600 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Window3d <  LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H  min4=30600 n4=121 squeeze=n > prec.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min9=30600 n9=121 squeeze=n > LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H 
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H /net/koko/homes/sep/alejandro/research/CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_off.x mig=ppp3.H initial=0  niter=10 eps=0. alpha=.001 na=1 reg_off=0 beta_max=10. prec=prec.H  remov=rmov_claudio.H mmov=mmov_new_sym_reg1_claudio.H wmov=wmov.H >$@

inv_Sis_az11ax21_newsolver_sym_strati_off_reg_test.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30600 n1=115 min5=13000 n5=180  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" > pp3.H
	< vel_strati_pad.H Window3d  min3=13000 n1=115 min1=30600 n3=180 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >pppppp1.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_diag.H exp="file1/file1" |Transp reshape=2,4| Transp plane=34 | Window3d  min3=13000 n3=180 min4=30600 n4=115 squeeze=n >  pppppp.H
	Math file1=pppppp1.H file2=pppppp.H exp="file1*file2" > resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_diag.H exp="1./(file1+0.2)" | Transp reshape=2,4| Transp plane=34 | Window3d  min3=13000 n3=180 min4=30600 n4=115 squeeze=n >  prec.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min9=30600 n9=115 min8=13000 n8=180 squeeze=n max_memory=500> LtL_Sis_az21_ax21_wave_marine_strati_off_t_w1.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_w1.H /net/koko/homes/sep/alejandro/research/CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_off.x mig=ppp3.H initial=0  niter=10  eps=0. alpha=.001 na=3 reg_off=0    prec=prec.H  remov=rmov.H mmov=mmov_new_sym_reg.H wmov=wmov.H >$@ #    prec=prec.H resw=resweigh.H value=1. beta_max=.1

   
inv_Sis_az11ax21_newsolver_sym_strati_off_reg_dp.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 extend=1>resweigh.H
	Cp LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H    prec.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H /net/koko/homes/sep/alejandro/research/CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_dp.x mig=ppp3.H initial=0  niter=10 eps=0. alpha=.001 na=1 reg_off=0  prec=prec.H remov=rmov_dp.H mmov=mmov_new_sym_reg_dp.H wmov=wmov.H >$@ # resw=resweigh.H   

inv_Sis_az11ax21_minres_sym_strati_off_reg.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=14000 n5=140  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" > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=14000 n8=140 squeeze=n > LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_w.H /net/koko/homes/sep/alejandro/research/MINRES/Bin/LINUX86_64/Sparse_Mxv_inv_minres_sym.x mig=ppp3.H  niter=100 eps=0.   >$@

##inv_Sis_az11ax21_newsolver_sym_strati_square.H:
#	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_square_sym.x mig=pp3.H niter=100 remov=rmov.H mmov=mmov_new_sym.H wmov=wmov.H >$@
#
#inv_Sis_az11ax21_newsolver_sym_strati_square_prec.H:
#	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
##   Preconditioner
##	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
#	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
#	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./file1" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
#	Cp LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv1.H
#	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H file2=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv1.H exp="file1*file2" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv_square.H
##   Weight
#	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
#	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
#	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 >resweigh.H
#	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv_square.H file2=resweigh.H exp="file1*file2" > pre_wei.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_square_sym.x mig=pp3.H eps=0.0 niter=50  remov=rmov.H mmov=mmov_old_sym_off.H w#mov=wmov.H prec=resweigh.H initial=0>$@ 

LtL_Sis_az21_ax21_wave_marine_strati_cr.H:
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=13000 n5=50  min3=0 n3=1 squeeze=n  > ppp.H
	< ppp.H Transp plane=13| Transp reshape=2,4,5,6 plane=23  > ppp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H Window3d squeeze=n n8=50 min8=13000 n9=133> pppp.H
	< pppp.H  ../NoStaDecon_save/Bin/LINUX86_64/test_symmetry.x m=ppp3.H >$@

LtL_Sis_az21_ax21_wave_marine_strati_cr_bin.H:
	< LtL_Sis_az21_ax21_wave_marine_strati_cr.H Transp max_memory=300 | Window3d n1=2 max_memory=300> headers.H
	< LtL_Sis_az21_ax21_wave_marine_strati_cr.H Transp max_memory=300 | Window3d n1=1 f1=2 max_memory=300 squeeze=n > 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=6650 dg1=1 og1=1 key2=row ng2=6650 dg2=1 og2=1 verb=1 max_memory=400 > t2$@
	< t2$@ Infill3d maxsize=200 verb=1 >$@

inv_Sis_az11ax21_chebi_sym_strati.H:
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=10000 n5=300   min3=0 n3=1| Transp > ppp3.H
	cp LtL_Sis_az21_ax21_wave_marine_strati_cr_bin.H LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H
	echo n1=39900 n2=39900 n3=1 >> LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H
	< LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H /net/koko/homes/sep/alejandro/research/CHEBI/Bin/LINUX86_64/driver_sep.x mig=ppp3.H niter=4 eps=100> $@
	echo n1=300 o1=10000      d1=25          label1=z n2=133        o2=29975      d2=75 >> $@

inv_Sis_az11ax21_cg_sym_strati.H:
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=10000 n5=300   min3=0 n3=1| Transp > ppp3.H
	cp LtL_Sis_az21_ax21_wave_marine_strati_cr_bin.H LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H
	echo n1=39900 n2=39900 n3=1 >> LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H
	< LtL_Sis_az21_ax21_wave_marine_strati_cr_bin_w.H /net/koko/homes/sep/alejandro/research/CG_MUSA/Bin/LINUX86_64/driver_cg.x mig=ppp3.H niter=10 eps=100> $@
	echo n1=300 o1=10000      d1=25          label1=z n2=133        o2=29975      d2=75  >> $@
#--------------------------------------------------------------------------------------------------------------------------------------------------------------
inv_Sis_az11ax21_newsolver_sym_strati_wei_prec.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./file1" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
#   Weight
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 >resweigh.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H resw=resweigh.H prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H  niter=15 eps=0. alpha=0.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_p.H wmov=wmov.H valuer=5. value=1 >$@  

inv_Sis_az11ax21_newsolver_sym_strati_prec.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./(file1)" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H  niter=20 eps=0.00 alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym.H wmov=wmov.H >$@ 

inv_Sis_az11ax21_newsolver_sym_strati_prec_claudio.H: 
	<prim_01_subt_sort.H  Window3d min1=30000 n1=133 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
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./(file1)" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H  niter=20 eps=0.00 alpha=.001 na=1 reg_off=0 remov=rmov_claudio0.H mmov=mmov_new_sym_claudio.H wmov=wmov.H beta_max=10.>$@

inv_Sis_az11ax21_newsolver_sym_strati_prec_test.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30600 n1=40 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
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 min9=30600 n9=40 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t_test.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 min9=30600 n9=40 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./file1" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t_test.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H    niter=20 eps=0.00 alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym.H wmov=wmov.H prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H >$@ #

inv_Sis_az11ax21_newsolver_sym_strati_prec_dp.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./file1" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_dp.x mig=p-p3.H  prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H  niter=5 eps=0.00 alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_dp.H wmov=wmov.H >$@ 

inv_Sis_az11ax21_minres_sym_strati_prec.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min1=0 n1=1 min6=0 n6=1 squeeze=n> LtL_Sis_az21_ax21_wave_marine_strati_t.H
	Window3d  < LtL_Sis_az21_ax21_wave_marine_strati_t.H n3=1 min3=0 n4=1 min4=0 squeeze=n | Transp reshape=5,10 > LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H exp="1./file1" > LtL_Sis_az21_ax21_wave_marine_strati_t_diag_inv.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../MINRES/Bin/LINUX86_64/Sparse_Mxv_inv_minres_sym.x mig=p-p3.H    niter=200 eps=0.  prec=LtL_Sis_az21_ax21_wave_marine_strati_t_diag.H > $@

inv_Sis_az11ax21_newsolver_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  niter=50 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_pppp.H wmov=wmov.H >$@

inv_Sis_az11ax21_minres_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../MINRES/Bin/LINUX86_64/Sparse_Mxv_inv_minres_sym.x mig=p-p3.H  niter=1000 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_minres.H wmov=wmov.H >$@

inv_Sis_az41ax21_newsolver_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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=1000. > p-p3.H
	< LtL_Sis_az41_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  niter=20 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_pppp_21.H wmov=wmov.H >$@


inv_Sis_az41ax41_newsolver_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
	< LtL_Sis_az41_ax41_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  niter=40 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_pppp_41.H wmov=wmov.H >$@

inv_Sis_az41ax41_minres_sym_strati.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
	< LtL_Sis_az41_ax41_wave_marine_strati_t.H ../MINRES/Bin/LINUX86_64/Sparse_Mxv_inv_minres_sym.x mig=p-p3.H  niter=40 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_pppp_41.H wmov=wmov.H >$@
 


inv_Sis_az11ax21_newsolver_sym_strati_w.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H Window3d squeeze=n min4=0 n4=1 j3=9  > LtL.H
	< LtL.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H  niter=100 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_pppp_w.H wmov=wmov.H >$@

inv_Sis_az11ax21_newsolver_sym_strati_wei.H: 
#   Weight
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 >resweigh.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=p-p3.H resw=resweigh.H  niter=50 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_ppppp.H wmov=wmov.H >$@
#--------------------------------------------------------------------------------------------------------------------------------------------------------------
inv_Sis_az11ax21_off_oldsolver_sym_mod.H: 
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=mig_Sis_refle_off_f35_marine.H niter=500 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_old_sym_off_mod.H wmov=wmov.H >$@

inv_Sis_az11ax21_off_newsolver_w.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min1=7975 n1=381 j2=3 min2=-600 n2=17 squeeze=n  > ppp.H
	< ppp.H /opt/SEP/SEP6.4.6/bin/Transp reshape=1,3,4,5,6  > ppp1.H
	< ppp1.H  Scale scale_to=1.|Window3d  min4=32000 n4=80 squeeze=n  > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t_prec.H Window3d  min4=32000 n4=80 squeeze=n > pp4.H
#	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d  min9=32000 n9=80 squeeze=n > pp6.H
	< pp6.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  prec=pp4.H niter=20 eps=0. 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 squeeze=n |../NoStaDecon_save/Bin/LINUX86_64/test_symmetry.x m=pp3.H >$@





test_sym_w.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 f8=50 n8=200 f9=30 n9=70 n3=11 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 | 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=20 >$@

inv_Sis_az21ax21_eps10_50_newsolver_test_reg.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | 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=.1 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 | 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_weig1.H  eps=.0 niter=20>$@

inv_Sis_az21ax21_eps10_50_newsolver_test_prec_sym.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1.> pp3.H
#	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 > pp4.H
	< pp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H  eps=.0  mmov=m_mov.H remov=r_mov.H niter=100>$@

inv_Sis_az21ax21_eps10_50_newsolver_test_prec_sym_cd.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1.> pp3.H
#	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 > pp4.H
	< pp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H  eps=.05  mmov=m_mov0.H remov=r_mov0.H niter=25 cd=1 cg=0>$@
#---------------------------------------------------------------------------------------------------------------------------------------------------
inv_Sis_az21ax21_eps10_50_newsolver_test_prec_sym_cg.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1.> pp30.H
#	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 > pp40.H
	< pp40.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp30.H  prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H  eps=.0  mmov=m_mov0.H remov=r_mov0_prec.H niter=50 cd=0 cg=1>$@

inv_Sis_az21ax21_eps10_50_oldsolver.H:
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 n1=381 | Transp reshape=2,4| Scale scale_to=1.| Pad n3=400> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H Window3d n8=381 squeeze=n | Pad n8=400> ppp4.H
	<  ppp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H  mmov=m_mov_old.H remov=r_mov.H eps=.0 niter=500>$@


# Same operator for the data
inv_Sis_az11ax21_off_newsolver_sym_mod_strati.H: 
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_diag.H exp="1./file1" | Transp reshape=2,4| Transp plane=34 | Window3d n1=1 min1=0 squeeze=n > prec0.H
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../CG_new/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x  mig=mig_Sis_refle_f35_marine_strati.H   niter=20 eps=0. alpha=.001 na=1 reg_off=0 remov=rmov.H mmov=mmov_new_sym_off_mod.H wmov=wmov.H >$@ #prec=prec0.H

inv_Sis_az21ax21_eps10_50_oldsolver_mod_strati.H: 
	< LtL_Sis_az21_ax21_wave_marine_strati_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=mig_Sis_refle_f35_marine_strati.H  mmov=m_mov_mod_old.H remov=r_mov.H eps=.0 niter=100>$@

inv_Sis_az21ax21_eps10_50_newsolver_mod_sym_cg.H: 
	< mig_Sis_refle_off_f35_marine.H Window3d min1=0 n1=1 squeeze=n| Scale scale_to=1.| Pad n3=400> pp3_mod.H
#	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d squeeze=n n1=1 min1=0 n6=1 min6=0 max_memory=2000| Pad n8=400> pp4_mod.H
	< pp4_mod.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3_mod.H  eps=.0  mmov=m_mov_mod.H remov=r_mov_mod.H niter=500 cd=0 cg=1>$@

inv_Sis_az41ax21_eps10_50_newsolver_mod_sym_cg.H: 
	< mig_Sis_refle_off_f35_marine.H Window3d min1=0 n1=1 squeeze=n| Scale scale_to=1.| Pad n3=410> pp3_mod.H
#	< LtL_Sis_az41_ax21_wave_marine_old_off1.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_az41_ax21_wave_marine_old_off1_t.H
	< LtL_Sis_az41_ax21_wave_marine_old_off1_t.H Window3d squeeze=n n1=1 min1=0 n6=1 min6=0 max_memory=2000 f3=8| Pad n8=410> pp41_mod.H
	< pp41_mod.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3_mod.H  eps=.0  mmov=m_mov_mod41.H remov=r_mov_mod.H niter=50 cd=0 cg=1>$@ 

inv_Sis_az21ax21_eps10_50_newsolver_mod_sym_cg_taper.H: 
	< mig_Sis_refle_off_f35_marine.H Window3d min1=0 n1=1 squeeze=n| Scale scale_to=1.| Pad n3=400> pp3_mod.H
#	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d squeeze=n n1=1 min1=0 n6=1 min6=0 max_memory=2000| Pad n8=400> pp4_mod.H
	< pp4_mod.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3_mod.H  eps=.0  mmov=m_mov_mod_taper.H remov=r_mov_mod.H niter=50 cd=0 cg=1>$@

inv_Sis_az21ax21_eps10_50_oldsolver_mod.H: 
	< mig_Sis_refle_off_f35_marine.H Window3d min1=0 n1=1 squeeze=n| Scale scale_to=1. | Pad n3=400> pp3_mod.H;
	< pp4_mod.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly_sym.x mig=pp3_mod.H  mmov=m_mov_mod_old.H remov=r_mov.H eps=.0 niter=100>$@

inv_Sis_az21ax21_eps10_50_newsolver_mod_prec_sym_cg.H: 
	< mig_Sis_refle_off_f35_marine.H Window3d min1=0 n1=1 squeeze=n| Scale scale_to=1.| Pad n3=400> pp3_mod.H
#	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d squeeze=n n1=1 min1=0 n6=1 min6=0 max_memory=2000 | Pad n8=400> pp4_mod.H
	< LtL_Sis_az21_ax21_wave_old_t_weig2.H Pad n3=400 > pp5_mod.H
	< pp4_mod.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3_mod.H  eps=.0 prec=pp5_mod.H mmov=m_mov_mod_prec.H remov=r_mov_mod.H wmov=w_mov_mod.H niter=50 cd=0 cg=1 na=1 alpha=0.001>$@

inv_Sis_az21ax21_eps10_50_newsolver_off_mod_sym_cg.H: 
	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=mig_Sis_refle_off_f35_marine.H  eps=.0  mmov=m_mov_mod_off.H remov=r_mov_mod_off.H niter=5 cd=0 cg=1>$@


inv_Sis_az21ax21_oldsolver_ang2off_mod_sym_strati_%.H: 
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_sym_ang2off.x mig=mig_Sis_refle_off_f35_marine_strati.H  eps=.0  mmov=m_mov_mod_ang2off_strati.H remov=r_mov_mod_off.H niter=$* >$@

inv_Sis_az21ax21_oldsolver_ang2off_mod_reg_strati_%.H: 
	< mig_Sis_pad_800_strati_nointerp_ang.H Transp reshape=1,3 | Transp plane=24 > pp6.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_sym_ang2off.x mig=mig_Sis_refle_off_f35_marine_strati.H  eps=0.01  mmov=m_mov_mod_ang2off_strati_reg.H remov=r_mov_mod_off.H niter=$* >$@

inv_Sis_az21ax21_oldsolver_ang2off_sym_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_sym_ang2off.x mig=pp3.H  eps=.0  mmov=m_mov_ang2off_strati.H remov=r_mov_mod_off.H niter=$* reg_off=0. amax=10.>$@

#inv_Sis_az21ax21_newsolver_ang2off_sym_strati_%.H: 
#	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym_ang2off.x mig=pp3.H  eps=.0  mmov=m_mov_ang2off_strati_newsolv.H remov=r_mov_mod_off_newsolv.H niter=$* reg_off=0. amax=10.>$@
#
inv_Sis_az21ax21_oldsolver_ang2off_reg_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=1000.  mmov=m_mov_ang2off_strati_reg.H remov=r_mov_mod_off_reg.H niter=$* reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@

inv_Sis_az21ax21_oldsolver_ang2off_reg_prec_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp3.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=1000.  mmov=m_mov_ang2off_strati_prec_reg.H remov=r_mov_mod_off_prec_reg.H niter=$*  resw=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H  reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@

inv_Sis_az21ax21_oldsolver_ang2off_reg_prec_mask_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. | Pad n3=400 > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2"| Pad n3=400 extend=1 > res_pre.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Pad n8=400 extend=1 > LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=1000  mmov=m_mov_ang2off_strati_prec_reg_mask.H remov=r_mov_mod_ang2off_prec_reg_mask.H niter=$*  resw=res_pre.H reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1 >$@ 

inv_Sis_az21ax21_oldsolver_ang2off_prec_mask_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. | Pad n3=400 > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2"| Pad n3=400 extend=1 > res_pre.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Pad n8=400 extend=1 > LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t_pad.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=0.  mmov=m_mov_ang2off_strati_prec_mask.H remov=r_mov_mod_ang2off_prec_mask.H niter=$*  resw=res_pre.H reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1 >$@  

inv_Sis_az21ax21_oldsolver_ang2off_mask_strati_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
	< vel_strati_pad.H Window3d  min3=8000 n1=133 min1=30000 n3=381 | Transp reshape=2,4 > ppp1.H
	Math file1=ppp1.H exp="@ABS(file1-14800)" | Scale scale_to=1. | Clip chop=greater clip=0.7 | Scale scale_to=1.| Window3d | Transp | Smooth rect1=30 rect2=30|Clip chop=less clip=.6 > ppp.H
	Math file1=ppp.H exp="file1-0.6"| Scale scale_to=1.| Transp reshape=2,4 | Pad extend=1 n1out=17 >resweigh.H
	Math file1=LtL_Sis_az21_ax21_wave_marine_strati_off_t_prec.H file2=resweigh.H exp="file1*file2" > res_pre.H
	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=0.  mmov=m_mov_ang2off_strati_mask.H remov=r_mov_mod_ang2off_mask.H niter=$*  resw=resweigh.H reg_off=0  amax=10.>$@  
 
inv_Sis_az21ax21_oldsolver_ang2off_strati_window_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=67 min5=10000 n5=301  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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=9970 n8=301 n9=67 squeeze=n > pppppw.H
	< pppppw.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=0.  mmov=m_mov_ang2off_strati_window.H remov=r_mov_mod_ang2off_window.H niter=$*  reg_off=0  amax=10. amin=-50  naax=61>$@


inv_Sis_az21ax21_oldsolver_ang2off_strati_reg_window_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=67 min5=10000 n5=301  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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=9970 n8=301 n9=67 squeeze=n > pppppw.H
	< pppppw.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=100.  mmov=m_mov_ang2off_strati_reg_window.H remov=r_mov_mod_ang2off_reg_window.H niter=$*  reg_off=1 amin=-10 dax=1.  naax=71 style=1 >$@

inv_Sis_az21ax21_oldsolver_ang2off_strati_coswreg_window_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=67 min5=10000 n5=301  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. > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=9970 n8=301 n9=67 squeeze=n > pppppw.H
	< pppppw.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=1000.  mmov=m_mov_ang2off_strati_coswreg_window.H remov=r_mov_mod_ang2off_coswreg_window.H niter=$*  reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@

inv_Sis_az21ax21_oldsolver_ang2off_strati_coswreg_window1_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=13000 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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=12975 n8=181  squeeze=n > ppppp.H
	< ppppp.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=1000.  mmov=m_mov_ang2off_strati_coswreg_window.H remov=r_mov_mod_ang2off_coswreg_window.H niter=$*  reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@

inv_Sis_az21ax21_oldsolver_ang2off_strati_window1_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=13000 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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=12975 n8=181  squeeze=n > ppppp.H
	< ppppp.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3.H  eps=0.  mmov=m_mov_ang2off_strati_coswreg_window.H remov=r_mov_mod_ang2off_coswreg_window.H niter=$*  reg_off=0 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@

inv_Sis_az21ax21_oldsolver_ang2off_strati_coswreg_window2_%.H: 
	< mig_Sis_pad_800_strati.H Window3d min1=30000 n1=133 min5=14700 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. > pp33.H
	Math file1=pp33.H exp="3.0*file1" > pp3w2.H
#	< LtL_Sis_az21_ax21_wave_marine_strati_off_t.H Window3d min8=14675 n8=113  squeeze=n > pppppw2.H
	< pppppw2.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_oldsolver_prec_ang2off.x mig=pp3w2.H  eps=1000.  mmov=m_mov_ang2off_strati_coswreg_window2.H remov=r_mov_mod_ang2off_coswreg_window2.H niter=$*  reg_off=1 amin=-10 dax=1.  naax=71  cosw=1 style=1>$@
#---------------------------------------------------------------------------------------------------------------------------------------------------
ir_taper_sym.H: 
	Spike n1=1 n2=1 n3=400 n4=133 k1=1 k2=1 k3=100 k4=10 > spike2d.H
#	< LtL_Sis_az21_ax21_wave_marine_old_off17_t.H Window3d squeeze=n n1=1 min1=0 n6=1 min6=0 max_memory=2000 | Pad n8=400> pp4_mod.H
	< pp4_mod.H  ../NoStaDecon_save/Bin/LINUX86_64/convolve_Mxv_sym.x m=spike2d.H >$@


inv_Sis_az21ax21_eps10_50_newsolver_test_sym.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1.> pp3.H
#	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 > pp4.H
	< pp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  eps=.0 mmov=m_mov.H remov=r_mov.H  niter=20>$@

inv_Sis_az21ax21_eps10_50_newsolver_test_prec_sym_wind.H: 
	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1. | Window3d min3=13000 squeeze=n> pp3.H
	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 min8=13000 > pp4.H
	< LtL_Sis_az21_ax21_wave_old_t_weig1.H Window3d min3=13000 squeeze=n> pp5.H
	< pp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H eps=.0 mmov=m_mov.H remov=r_mov.H  niter=50 cd=1 cg=0 prec=pp5.H >$@

inv_Sis_az21ax21_eps10_50_newsolver_test_prec_sym_reg.H: 
#	< mig_Sis_Parallel_pad_2048_off.H Window3d min4=30000 n4=133 min2=0 n2=1 min1=7975 | Pad n1=383 | Transp reshape=2,4| Scale scale_to=1.> pp3.H
#	< LtL_Sis_az21_ax21_wave_old_t.H Window3d squeeze=n n3=11 > pp4.H
	< pp4.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver_sym.x mig=pp3.H  prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H  eps=.01 niter=10>$@ 

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 | 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 > pp5.H
	< pp5.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_newsolver.x mig=pp3.H  eps=.0 niter=100>$@

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 | 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  mmov=m_mov.H remov=r_mov.H eps=.0 niter=100>$@

inv_Sis_az21ax21_eps10_50_test_prec_mask.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 | Transp reshape=2,4| Scale scale_to=1.> pp3.H
	< vel_pad.H Window3d  min3=7975 n1=133 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_t.H ../NoStaDecon_save/Bin/LINUX86_64/Sparse_Mxv_inv_off_fly.x mig=pp3.H prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H resw=resweigh.H mmov=m_mov.H remov=r_mov.H eps=.0 niter=30>$@

inv_Sis_az21ax21_eps10_50_test_prec.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 | 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 prec=LtL_Sis_az21_ax21_wave_old_t_weig1.H  mmov=m_mov.H remov=r_mov.H eps=.0 niter=50>$@

test_sym_bin.H:
#	< test_sym.H Transp max_memory=300 | Window3d n1=2 max_memory=300> headers.H
	< test_sym.H Transp max_memory=300 | Window3d n1=1 f1=2 max_memory=300 squeeze=n > 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 | 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
