################
#
# Importante !!! 04/04/10
#
# Falta rodar usando como input o window.scale
# elf.3D.bsch.window.scale.H
# n1=451        o1=0          d1=30          label1=
# n2=31         o2=-450       d2=30          label2=
# n3=133        o3=1600       d3=30          label3=
# n4=1          o4=0          d4=40          label4=
# n5=302        o5=0          d5=10          label5=
#
# Window3d n3=1 f3=$* < elf.3D.bsch.window.scale.H | Transp plane=12 reshape=2,3,5 | Rotation
#
# Feito em 04/05/10
#
#include Collect

Tese=/net/server/claudio/Research/Prst_ExpRefl/Tese
include 3Dprocfile_1

elf.L7f.vinit.grid30.H:
	< elf.L7f.vinit.grid20.H Interp d1out=30 d2out=30 o1out=-4500 o2out=-2000 n1out=800 n2out=420 type=1 > $@



topchalk02:
	< elf.L7f.vinit.grid30.H Window3d min1=-1200 n1=531 min2=700 n2=193|Transp plane=13|Gpow gpow=-1| \
	${Tese}/bin/OW_3deditvel_pik.x hor1=hor3dgt2500a.pik v=2500 > /dev/null



#########
#
# Improve elf.L7f.vinit.grid30.H with rmo picks
#

#	< j2.H Transp plane=13 | Mute vmute=0 tramp=0 tmute=100 | Clip chop=less to=-1 clip=0.001 > j2s.H

elf.L7f.vinit.rmo.grid30.H:
	< elf.L7f.vinit.grid30.H Window3d n1=448 min1=0 n2=133 min2=1930 | Gpow gpow=-1 | \
	Interp d1out=25 d2out=25 d3out=12.5 o2out=1950 n2out=133 n1out=540 maxsize=1000 type=1 > j0.H
	< j0.H Transp plane=13 | ${Tese}/bin/OW_vrms.x int2rms=1 | Transp plane=13 > j1.H
	< 3d.rmo.sedim.H Window3d max1=12000|Smooth rect1=31 rect2=31 rect3=31| Transp plane=13 | Mute vmute=0 tramp=0 tmute=100| \
	Clip chop=less to=-1 clip=0.001 | Interp d1out=12.5 d2out=25 d3out=25 o2out=1950 n2out=133 n3out=540 type=1 \
	n1out=401 maxsize=1000 |Transp plane=13 > j2.H 
	Math file1=j1.H file2=j2.H exp="file1*file2" | Transp plane=13 | ${Tese}/bin/OW_vrms.x int2rms=0 > j3.H
	< j0.H Transp plane=13 > j00.H
	< j3.H ../Tese/bin/OW_3deditvel_mut.x hor1=survey540 vvel=1 vsub=j00.H | Gpow gpow=-1 | Transp plane=13 > j4.H
	< j4.H Interp d3out=10 d2out=0 d1out=30 type=1 maxsize=1000 | Pad beg1=100 end1=100 beg2=100 end2=100 extend=1 > $@

elf.L7f.vinit.rmo2.grid30.H:
#	< elf.L7f.vinit.grid30.H Window3d n1=448 min1=0 n2=133 min2=1930 | Gpow gpow=-1 | \
#	Interp d1out=25 d2out=25 d3out=12.5 o2out=1950 n2out=133 n1out=540 maxsize=1000 type=1 > j0.H
#	< j0.H Transp plane=13 | ${Tese}/bin/OW_vrms.x int2rms=1 | Transp plane=13 > j1.H
#	< 3d.rmo.sedim2.H Window3d max1=12000 n3=115 |Pad beg3=0 end3=286 extend=1|Smooth rect1=31 rect2=31 rect3=31| Transp\
#	 plane=13 | Mute vmute=0 tramp=0 tmute=100 | \
#	Clip chop=less to=-1 clip=0.001 | Interp d1out=12.5 d2out=25 d3out=25 o2out=1950 n2out=133 n3out=540 type=1 \
#	n1out=401 maxsize=1000 |Transp plane=13 > j2.H 
#	Math file1=j1.H file2=j2.H exp="file1*file2" | Transp plane=13 | ${Tese}/bin/OW_vrms.x int2rms=0 > j3.H
#	< j0.H Transp plane=13 > j00.H
	< j3.H ../Tese/bin/OW_3deditvel_mut.x hor1=survey540 vvel=1 vsub=j00.H | Gpow gpow=-1 | Transp plane=13 > j4.H
	< j4.H Interp d3out=10 d2out=30 d1out=30 type=1 maxsize=1000 | Pad beg1=100 end1=100 beg2=100 end2=100 extend=1 > $@

elf.L7f.sinit.rmo.grid20.H:
	< elf.L7f.vinit.rmo.grid30.H Gpow gpow=-1 |Window3d max3=3015 | Interp d1out=20 d2out=20 d3out=10 maxsize=1000 type=1 > $@
	<$@ Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120> \
	elf.L7f.ref_vel.rmo.grid20.H

elf.L7f.sinit.rmo2.grid20.H:
	< elf.L7f.vinit.rmo2.grid30.H Gpow gpow=-1 |Window3d max3=3015 | Interp d1out=20 d2out=20 d3out=10 maxsize=1000 type=1 > $@
	<$@ Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120> \
	elf.L7f.ref_vel.rmo2.grid20.H

elf.3D.mig_adj.rmo2.H:
#	Window3d squeeze=n < elf.inv_freq.recip.H f5=1 > data.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data.H,elf.L7f.sinit.rmo2.grid20.H,elf.L7f.ref_vel.rmo2.grid20.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="5,4,4,6" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=359 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull360.par

elf.3D.mig_adj.agc.rmo2.H:
	Window3d squeeze=n < data.agc..H f5=1 > data.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data.H,elf.L7f.sinit.rmo2.grid20.H,elf.L7f.ref_vel.rmo2.grid20.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="5,4,4,6" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=359 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull360.par

elf.3D.mig.rmo2.int30.H:
	< elf.3D.mig_adj.rmo2.H Window3d min3=-450 n3=19 | Interp d1out=30 > temp00.H
	< temp00.H Transp | Interp d1out=30 > temp01.H
	< temp01.H Transp plane=31 | Interp d1out=30 > temp00.H
	< temp00.H Transp plane=21 reshape=1,3,5 > $@
	Rm temp*H

elf.3D.mig.rmo2.int30.new.H:
	echo o3=-1775 d3=25 >> elf.3D.mig_adj.rmo2.H
	< elf.3D.mig_adj.rmo2.H Window3d min3=-450 n3=37 | Interp d1out=30 maxsize=2000 > temp00.H
	< temp00.H Transp maxsize=2000 | Interp d1out=30 maxsize=2000 > temp01.H
	< temp01.H Transp plane=31 maxsize=2000 | Interp d1out=30 maxsize=2000 > temp00.H
	< temp00.H Transp plane=21 reshape=1,3,5 maxsize=2000 > $@
	Rm temp*H

# Base chalk left
elf.3D.bsch.window_left.H:
	< elf.3D.mig.rmo2.int30.H Transp plane=23 | ${Tese}/bin/OW3D_windowing.x nhorz=1 hor1=bsch_left_prstk wavz=5 | Window3d min1=510 min3=1990 max3=5050 | Pad beg1=17 beg3=13 end3=18> $@

elf.3D.bsch.window_left.mask%.H:
	Window3d < elf.3D.mig.rmo2.int30.H n3=1 f3=$* | Transp plane=13 > j0$*.H
	< elf.3D.bsch.window_left.H Window3d n2=1 f2=$* | Transp plane=13 | Agc window=200 | Envelope tc1=5 tc2=5 tc3=5 verb=1 > j1$*.H
	< j1$*.H Scale scale_to=1 | Clip to=1 clip=0.05 chop=greater | Clip to=0 clip=0.05 chop=less > j2$*.H
	Math file1=j0$*.H file2=j2$*.H exp="file1*file2" > $@
##	Rm j[0-2]$*.H 

elf.3D.bsch.window_left:
	make elf.3D.bsch.window_left.mask00.H;make elf.3D.bsch.window_left.mask01.H;make elf.3D.bsch.window_left.mask02.H;make elf.3D.bsch.window_left.mask03.H;make elf.3D.bsch.window_left.mask04.H;\
	make elf.3D.bsch.window_left.mask05.H;make elf.3D.bsch.window_left.mask06.H;make elf.3D.bsch.window_left.mask07.H;make elf.3D.bsch.window_left.mask08.H;make elf.3D.bsch.window_left.mask09.H;\
	make elf.3D.bsch.window_left.mask10.H;make elf.3D.bsch.window_left.mask11.H;make elf.3D.bsch.window_left.mask12.H;make elf.3D.bsch.window_left.mask13.H;make elf.3D.bsch.window_left.mask14.H;\
	make elf.3D.bsch.window_left.mask15.H;make elf.3D.bsch.window_left.mask16.H;make elf.3D.bsch.window_left.mask17.H;make elf.3D.bsch.window_left.mask18.H;make elf.3D.bsch.window_left.mask19.H;\
	make elf.3D.bsch.window_left.mask20.H;make elf.3D.bsch.window_left.mask21.H;make elf.3D.bsch.window_left.mask22.H;make elf.3D.bsch.window_left.mask23.H;make elf.3D.bsch.window_left.mask24.H;\
	make elf.3D.bsch.window_left.mask25.H;make elf.3D.bsch.window_left.mask26.H;make elf.3D.bsch.window_left.mask27.H;make elf.3D.bsch.window_left.mask28.H;make elf.3D.bsch.window_left.mask29.H;\
	make elf.3D.bsch.window_left.mask30.H

elf.3D.bsch.window_left.int30.H:
	Cat elf.3D.bsch.window_left.mask00.H elf.3D.bsch.window_left.mask01.H elf.3D.bsch.window_left.mask02.H elf.3D.bsch.window_left.mask03.H elf.3D.bsch.window_left.mask04.H \
	elf.3D.bsch.window_left.mask05.H elf.3D.bsch.window_left.mask06.H elf.3D.bsch.window_left.mask07.H elf.3D.bsch.window_left.mask08.H elf.3D.bsch.window_left.mask09.H \
	elf.3D.bsch.window_left.mask10.H elf.3D.bsch.window_left.mask11.H elf.3D.bsch.window_left.mask12.H elf.3D.bsch.window_left.mask13.H elf.3D.bsch.window_left.mask14.H \
	elf.3D.bsch.window_left.mask15.H elf.3D.bsch.window_left.mask16.H elf.3D.bsch.window_left.mask17.H elf.3D.bsch.window_left.mask18.H elf.3D.bsch.window_left.mask19.H \
	elf.3D.bsch.window_left.mask20.H elf.3D.bsch.window_left.mask21.H elf.3D.bsch.window_left.mask22.H elf.3D.bsch.window_left.mask23.H elf.3D.bsch.window_left.mask24.H \
	elf.3D.bsch.window_left.mask25.H elf.3D.bsch.window_left.mask26.H elf.3D.bsch.window_left.mask27.H elf.3D.bsch.window_left.mask28.H elf.3D.bsch.window_left.mask29.H \
	elf.3D.bsch.window_left.mask30.H axis=4 | Transp plane=23 reshape=1,2,4,5 > $@
#	Rm elf.3D.bsch.window_left.mask*.H

elf.3D.bsch.left.rot%:
	< elf.3D.bsch.window_left.int30.H Window3d n4=1 f4=$* max2=5100 | ${Tese}/bin/OW_scale.x > a$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a$*.H Cig_tilt=b$*.H par=bsch_right-rot.P positive_rot=1
	< b$*.H Pad beg2=0 end2=278 > $@.H
	Rm a$*.H b$*.H

elf.3D.bschleft.rot:
	make elf.3D.bsch.left.rot000;make elf.3D.bsch.left.rot001;make elf.3D.bsch.left.rot002;make elf.3D.bsch.left.rot003;make elf.3D.bsch.left.rot004;\
	make elf.3D.bsch.left.rot005;make elf.3D.bsch.left.rot006;make elf.3D.bsch.left.rot007;make elf.3D.bsch.left.rot008;make elf.3D.bsch.left.rot009;\
	make elf.3D.bsch.left.rot010;make elf.3D.bsch.left.rot011;make elf.3D.bsch.left.rot012;make elf.3D.bsch.left.rot013;make elf.3D.bsch.left.rot014;\
	make elf.3D.bsch.left.rot015;make elf.3D.bsch.left.rot016;make elf.3D.bsch.left.rot017;make elf.3D.bsch.left.rot018;make elf.3D.bsch.left.rot019;\
	make elf.3D.bsch.left.rot020;make elf.3D.bsch.left.rot021;make elf.3D.bsch.left.rot022;make elf.3D.bsch.left.rot023;make elf.3D.bsch.left.rot024;\
	make elf.3D.bsch.left.rot025;make elf.3D.bsch.left.rot026;make elf.3D.bsch.left.rot027;make elf.3D.bsch.left.rot028;make elf.3D.bsch.left.rot029;\
	make elf.3D.bsch.left.rot030;make elf.3D.bsch.left.rot031;make elf.3D.bsch.left.rot032;make elf.3D.bsch.left.rot033;make elf.3D.bsch.left.rot034;\
	make elf.3D.bsch.left.rot035;make elf.3D.bsch.left.rot036;make elf.3D.bsch.left.rot037;make elf.3D.bsch.left.rot038;make elf.3D.bsch.left.rot039;\
	make elf.3D.bsch.left.rot040;make elf.3D.bsch.left.rot041;make elf.3D.bsch.left.rot042;make elf.3D.bsch.left.rot043;make elf.3D.bsch.left.rot044;\
	make elf.3D.bsch.left.rot045;make elf.3D.bsch.left.rot046;make elf.3D.bsch.left.rot047;make elf.3D.bsch.left.rot048;make elf.3D.bsch.left.rot049;\
	make elf.3D.bsch.left.rot050;make elf.3D.bsch.left.rot051;make elf.3D.bsch.left.rot052;make elf.3D.bsch.left.rot053;make elf.3D.bsch.left.rot054;\
	make elf.3D.bsch.left.rot055;make elf.3D.bsch.left.rot056;make elf.3D.bsch.left.rot057;make elf.3D.bsch.left.rot058;make elf.3D.bsch.left.rot059;\
	make elf.3D.bsch.left.rot060;make elf.3D.bsch.left.rot061;make elf.3D.bsch.left.rot062;make elf.3D.bsch.left.rot063;make elf.3D.bsch.left.rot064;\
	make elf.3D.bsch.left.rot065;make elf.3D.bsch.left.rot066;make elf.3D.bsch.left.rot067;make elf.3D.bsch.left.rot068;make elf.3D.bsch.left.rot069;\
	make elf.3D.bsch.left.rot070;make elf.3D.bsch.left.rot071;make elf.3D.bsch.left.rot072;make elf.3D.bsch.left.rot073;make elf.3D.bsch.left.rot074;\
	make elf.3D.bsch.left.rot075;make elf.3D.bsch.left.rot076;make elf.3D.bsch.left.rot077;make elf.3D.bsch.left.rot078;make elf.3D.bsch.left.rot079;\
	make elf.3D.bsch.left.rot080;make elf.3D.bsch.left.rot081;make elf.3D.bsch.left.rot082;make elf.3D.bsch.left.rot083;make elf.3D.bsch.left.rot084;\
	make elf.3D.bsch.left.rot085;make elf.3D.bsch.left.rot086;make elf.3D.bsch.left.rot087;make elf.3D.bsch.left.rot088;make elf.3D.bsch.left.rot089;\
	make elf.3D.bsch.left.rot090;make elf.3D.bsch.left.rot091;make elf.3D.bsch.left.rot092;make elf.3D.bsch.left.rot093;make elf.3D.bsch.left.rot094;\
	make elf.3D.bsch.left.rot095;make elf.3D.bsch.left.rot096;make elf.3D.bsch.left.rot097;make elf.3D.bsch.left.rot098;make elf.3D.bsch.left.rot099;\
	make elf.3D.bsch.left.rot100;make elf.3D.bsch.left.rot101;make elf.3D.bsch.left.rot102;make elf.3D.bsch.left.rot103;make elf.3D.bsch.left.rot104;\
	make elf.3D.bsch.left.rot105;make elf.3D.bsch.left.rot106;make elf.3D.bsch.left.rot107;make elf.3D.bsch.left.rot108;make elf.3D.bsch.left.rot109;\
	make elf.3D.bsch.left.rot110;make elf.3D.bsch.left.rot111;make elf.3D.bsch.left.rot112;make elf.3D.bsch.left.rot113;make elf.3D.bsch.left.rot114;\
	make elf.3D.bsch.left.rot115;make elf.3D.bsch.left.rot116;make elf.3D.bsch.left.rot117;make elf.3D.bsch.left.rot118;make elf.3D.bsch.left.rot119;\
	make elf.3D.bsch.left.rot120;make elf.3D.bsch.left.rot121;make elf.3D.bsch.left.rot122;make elf.3D.bsch.left.rot123;make elf.3D.bsch.left.rot124;\
	make elf.3D.bsch.left.rot125;make elf.3D.bsch.left.rot126;make elf.3D.bsch.left.rot127;make elf.3D.bsch.left.rot128;make elf.3D.bsch.left.rot129;\
	make elf.3D.bsch.left.rot130;make elf.3D.bsch.left.rot131;make elf.3D.bsch.left.rot132


elf.3D.bsch_left.rot:
	Cat axis=4 elf.3D.bsch.left.rot000.H elf.3D.bsch.left.rot001.H elf.3D.bsch.left.rot002.H elf.3D.bsch.left.rot003.H elf.3D.bsch.left.rot004.H elf.3D.bsch.left.rot005.H elf.3D.bsch.left.rot006.H elf.3D.bsch.left.rot007.H elf.3D.bsch.left.rot008.H elf.3D.bsch.left.rot009.H elf.3D.bsch.left.rot010.H elf.3D.bsch.left.rot011.H elf.3D.bsch.left.rot012.H elf.3D.bsch.left.rot013.H elf.3D.bsch.left.rot014.H elf.3D.bsch.left.rot015.H elf.3D.bsch.left.rot016.H elf.3D.bsch.left.rot017.H elf.3D.bsch.left.rot018.H elf.3D.bsch.left.rot019.H elf.3D.bsch.left.rot020.H elf.3D.bsch.left.rot021.H elf.3D.bsch.left.rot022.H elf.3D.bsch.left.rot023.H elf.3D.bsch.left.rot024.H elf.3D.bsch.left.rot025.H elf.3D.bsch.left.rot026.H elf.3D.bsch.left.rot027.H elf.3D.bsch.left.rot028.H elf.3D.bsch.left.rot029.H elf.3D.bsch.left.rot030.H elf.3D.bsch.left.rot031.H elf.3D.bsch.left.rot032.H elf.3D.bsch.left.rot033.H elf.3D.bsch.left.rot034.H elf.3D.bsch.left.rot035.H elf.3D.bsch.left.rot036.H elf.3D.bsch.left.rot037.H elf.3D.bsch.left.rot038.H elf.3D.bsch.left.rot039.H elf.3D.bsch.left.rot040.H elf.3D.bsch.left.rot041.H elf.3D.bsch.left.rot042.H elf.3D.bsch.left.rot043.H elf.3D.bsch.left.rot044.H elf.3D.bsch.left.rot045.H elf.3D.bsch.left.rot046.H elf.3D.bsch.left.rot047.H elf.3D.bsch.left.rot048.H > a.H
	Cat axis=4 elf.3D.bsch.left.rot049.H elf.3D.bsch.left.rot050.H elf.3D.bsch.left.rot051.H elf.3D.bsch.left.rot052.H elf.3D.bsch.left.rot053.H elf.3D.bsch.left.rot054.H elf.3D.bsch.left.rot055.H elf.3D.bsch.left.rot056.H elf.3D.bsch.left.rot057.H elf.3D.bsch.left.rot058.H elf.3D.bsch.left.rot059.H elf.3D.bsch.left.rot060.H elf.3D.bsch.left.rot061.H elf.3D.bsch.left.rot062.H elf.3D.bsch.left.rot063.H elf.3D.bsch.left.rot064.H elf.3D.bsch.left.rot065.H elf.3D.bsch.left.rot066.H elf.3D.bsch.left.rot067.H elf.3D.bsch.left.rot068.H elf.3D.bsch.left.rot069.H elf.3D.bsch.left.rot070.H elf.3D.bsch.left.rot071.H elf.3D.bsch.left.rot072.H elf.3D.bsch.left.rot073.H elf.3D.bsch.left.rot074.H elf.3D.bsch.left.rot075.H elf.3D.bsch.left.rot076.H elf.3D.bsch.left.rot077.H elf.3D.bsch.left.rot078.H elf.3D.bsch.left.rot079.H elf.3D.bsch.left.rot080.H elf.3D.bsch.left.rot081.H elf.3D.bsch.left.rot082.H elf.3D.bsch.left.rot083.H elf.3D.bsch.left.rot084.H elf.3D.bsch.left.rot085.H elf.3D.bsch.left.rot086.H elf.3D.bsch.left.rot087.H elf.3D.bsch.left.rot088.H elf.3D.bsch.left.rot089.H elf.3D.bsch.left.rot090.H elf.3D.bsch.left.rot091.H elf.3D.bsch.left.rot092.H elf.3D.bsch.left.rot093.H elf.3D.bsch.left.rot094.H elf.3D.bsch.left.rot095.H elf.3D.bsch.left.rot096.H elf.3D.bsch.left.rot097.H > b.H
	Cat axis=4 elf.3D.bsch.left.rot098.H elf.3D.bsch.left.rot099.H elf.3D.bsch.left.rot100.H elf.3D.bsch.left.rot101.H elf.3D.bsch.left.rot102.H elf.3D.bsch.left.rot103.H elf.3D.bsch.left.rot104.H elf.3D.bsch.left.rot105.H elf.3D.bsch.left.rot106.H elf.3D.bsch.left.rot107.H elf.3D.bsch.left.rot108.H elf.3D.bsch.left.rot109.H elf.3D.bsch.left.rot110.H elf.3D.bsch.left.rot111.H elf.3D.bsch.left.rot112.H elf.3D.bsch.left.rot113.H elf.3D.bsch.left.rot114.H elf.3D.bsch.left.rot115.H elf.3D.bsch.left.rot116.H elf.3D.bsch.left.rot117.H elf.3D.bsch.left.rot118.H elf.3D.bsch.left.rot119.H elf.3D.bsch.left.rot120.H elf.3D.bsch.left.rot121.H elf.3D.bsch.left.rot122.H elf.3D.bsch.left.rot123.H elf.3D.bsch.left.rot124.H elf.3D.bsch.left.rot125.H elf.3D.bsch.left.rot126.H elf.3D.bsch.left.rot127.H elf.3D.bsch.left.rot128.H elf.3D.bsch.left.rot129.H elf.3D.bsch.left.rot130.H elf.3D.bsch.left.rot131.H elf.3D.bsch.left.rot132.H > c.H
	Cat a.H b.H c.H axis=4 > $@.H


##################
# Base chalk right
elf.3D.bsch.window_right.H:
	< elf.3D.mig.rmo2.int30.H Window3d min1=5130|Transp plane=23|${Tese}/bin/OW3D_windowing.x nhorz=1 hor1=bsch_right_prstk wavz=5|Window3d min1=5490 min3=1990 max3=5050|Pad beg1=183 beg3=13 end3=18 > $@

elf.3D.bsch.window_right.mask%.H:
	Window3d < elf.3D.mig.rmo2.int30.H n3=1 f3=$* | Transp plane=13 > h0$*.H
	< elf.3D.bsch.window_right.H Window3d n2=1 f2=$* | Transp plane=13 | Agc window=200 | Envelope tc1=5 tc2=5 tc3=5 verb=1 > h1$*.H
	< h1$*.H Scale scale_to=1 | Clip to=1 clip=0.05 chop=greater | Clip to=0 clip=0.05 chop=less > h2$*.H
	Math file1=h0$*.H file2=h2$*.H exp="file1*file2" > $@
#	Rm h[0-2]$*.H 

elf.3D.bsch.window_right:
	make elf.3D.bsch.window_right.mask00.H;make elf.3D.bsch.window_right.mask01.H;make elf.3D.bsch.window_right.mask02.H;make elf.3D.bsch.window_right.mask03.H;make elf.3D.bsch.window_right.mask04.H;\
	make elf.3D.bsch.window_right.mask05.H;make elf.3D.bsch.window_right.mask06.H;make elf.3D.bsch.window_right.mask07.H;make elf.3D.bsch.window_right.mask08.H;make elf.3D.bsch.window_right.mask09.H;\
	make elf.3D.bsch.window_right.mask10.H;make elf.3D.bsch.window_right.mask11.H;make elf.3D.bsch.window_right.mask12.H;make elf.3D.bsch.window_right.mask13.H;make elf.3D.bsch.window_right.mask14.H;\
	make elf.3D.bsch.window_right.mask15.H;make elf.3D.bsch.window_right.mask16.H;make elf.3D.bsch.window_right.mask17.H;make elf.3D.bsch.window_right.mask18.H;make elf.3D.bsch.window_right.mask19.H;\
	make elf.3D.bsch.window_right.mask20.H;make elf.3D.bsch.window_right.mask21.H;make elf.3D.bsch.window_right.mask22.H;make elf.3D.bsch.window_right.mask23.H;make elf.3D.bsch.window_right.mask24.H;\
	make elf.3D.bsch.window_right.mask25.H;make elf.3D.bsch.window_right.mask26.H;make elf.3D.bsch.window_right.mask27.H;make elf.3D.bsch.window_right.mask28.H;make elf.3D.bsch.window_right.mask29.H;\
	make elf.3D.bsch.window_right.mask30.H

elf.3D.bsch.window_right.int30.H:
	Cat elf.3D.bsch.window_right.mask00.H elf.3D.bsch.window_right.mask01.H elf.3D.bsch.window_right.mask02.H elf.3D.bsch.window_right.mask03.H elf.3D.bsch.window_right.mask04.H \
	elf.3D.bsch.window_right.mask05.H elf.3D.bsch.window_right.mask06.H elf.3D.bsch.window_right.mask07.H elf.3D.bsch.window_right.mask08.H elf.3D.bsch.window_right.mask09.H \
	elf.3D.bsch.window_right.mask10.H elf.3D.bsch.window_right.mask11.H elf.3D.bsch.window_right.mask12.H elf.3D.bsch.window_right.mask13.H elf.3D.bsch.window_right.mask14.H \
	elf.3D.bsch.window_right.mask15.H elf.3D.bsch.window_right.mask16.H elf.3D.bsch.window_right.mask17.H elf.3D.bsch.window_right.mask18.H elf.3D.bsch.window_right.mask19.H \
	elf.3D.bsch.window_right.mask20.H elf.3D.bsch.window_right.mask21.H elf.3D.bsch.window_right.mask22.H elf.3D.bsch.window_right.mask23.H elf.3D.bsch.window_right.mask24.H \
	elf.3D.bsch.window_right.mask25.H elf.3D.bsch.window_right.mask26.H elf.3D.bsch.window_right.mask27.H elf.3D.bsch.window_right.mask28.H elf.3D.bsch.window_right.mask29.H \
	elf.3D.bsch.window_right.mask30.H axis=4 | Transp plane=23 reshape=1,2,4,5 > $@
#	Rm elf.3D.bsch.window_right.mask*.H

elf.3D.bsch.right.rot%:
	< elf.3D.bsch.window_right.int30.H Window3d n4=1 f4=$* min2=5100 | ${Tese}/bin/OW_scale.x > c$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=c$*.H Cig_tilt=d$*.H par=bsch_right-rot.P positive_rot=1
	< d$*.H Pad beg2=170 end2=0  > $@.H
	Rm c$*.H d$*.H

elf.3D.bschright.rot:
	make elf.3D.bsch.right.rot000;make elf.3D.bsch.right.rot001;make elf.3D.bsch.right.rot002;make elf.3D.bsch.right.rot003;make elf.3D.bsch.right.rot004;\
	make elf.3D.bsch.right.rot005;make elf.3D.bsch.right.rot006;make elf.3D.bsch.right.rot007;make elf.3D.bsch.right.rot008;make elf.3D.bsch.right.rot009;\
	make elf.3D.bsch.right.rot010;make elf.3D.bsch.right.rot011;make elf.3D.bsch.right.rot012;make elf.3D.bsch.right.rot013;make elf.3D.bsch.right.rot014;\
	make elf.3D.bsch.right.rot015;make elf.3D.bsch.right.rot016;make elf.3D.bsch.right.rot017;make elf.3D.bsch.right.rot018;make elf.3D.bsch.right.rot019;\
	make elf.3D.bsch.right.rot020;make elf.3D.bsch.right.rot021;make elf.3D.bsch.right.rot022;make elf.3D.bsch.right.rot023;make elf.3D.bsch.right.rot024;\
	make elf.3D.bsch.right.rot025;make elf.3D.bsch.right.rot026;make elf.3D.bsch.right.rot027;make elf.3D.bsch.right.rot028;make elf.3D.bsch.right.rot029;\
	make elf.3D.bsch.right.rot030;make elf.3D.bsch.right.rot031;make elf.3D.bsch.right.rot032;make elf.3D.bsch.right.rot033;make elf.3D.bsch.right.rot034;\
	make elf.3D.bsch.right.rot035;make elf.3D.bsch.right.rot036;make elf.3D.bsch.right.rot037;make elf.3D.bsch.right.rot038;make elf.3D.bsch.right.rot039;\
	make elf.3D.bsch.right.rot040;make elf.3D.bsch.right.rot041;make elf.3D.bsch.right.rot042;make elf.3D.bsch.right.rot043;make elf.3D.bsch.right.rot044;\
	make elf.3D.bsch.right.rot045;make elf.3D.bsch.right.rot046;make elf.3D.bsch.right.rot047;make elf.3D.bsch.right.rot048;make elf.3D.bsch.right.rot049;\
	make elf.3D.bsch.right.rot050;make elf.3D.bsch.right.rot051;make elf.3D.bsch.right.rot052;make elf.3D.bsch.right.rot053;make elf.3D.bsch.right.rot054;\
	make elf.3D.bsch.right.rot055;make elf.3D.bsch.right.rot056;make elf.3D.bsch.right.rot057;make elf.3D.bsch.right.rot058;make elf.3D.bsch.right.rot059;\
	make elf.3D.bsch.right.rot060;make elf.3D.bsch.right.rot061;make elf.3D.bsch.right.rot062;make elf.3D.bsch.right.rot063;make elf.3D.bsch.right.rot064;\
	make elf.3D.bsch.right.rot065;make elf.3D.bsch.right.rot066;make elf.3D.bsch.right.rot067;make elf.3D.bsch.right.rot068;make elf.3D.bsch.right.rot069;\
	make elf.3D.bsch.right.rot070;make elf.3D.bsch.right.rot071;make elf.3D.bsch.right.rot072;make elf.3D.bsch.right.rot073;make elf.3D.bsch.right.rot074;\
	make elf.3D.bsch.right.rot075;make elf.3D.bsch.right.rot076;make elf.3D.bsch.right.rot077;make elf.3D.bsch.right.rot078;make elf.3D.bsch.right.rot079;\
	make elf.3D.bsch.right.rot080;make elf.3D.bsch.right.rot081;make elf.3D.bsch.right.rot082;make elf.3D.bsch.right.rot083;make elf.3D.bsch.right.rot084;\
	make elf.3D.bsch.right.rot085;make elf.3D.bsch.right.rot086;make elf.3D.bsch.right.rot087;make elf.3D.bsch.right.rot088;make elf.3D.bsch.right.rot089;\
	make elf.3D.bsch.right.rot090;make elf.3D.bsch.right.rot091;make elf.3D.bsch.right.rot092;make elf.3D.bsch.right.rot093;make elf.3D.bsch.right.rot094;\
	make elf.3D.bsch.right.rot095;make elf.3D.bsch.right.rot096;make elf.3D.bsch.right.rot097;make elf.3D.bsch.right.rot098;make elf.3D.bsch.right.rot099;\
	make elf.3D.bsch.right.rot100;make elf.3D.bsch.right.rot101;make elf.3D.bsch.right.rot102;make elf.3D.bsch.right.rot103;make elf.3D.bsch.right.rot104;\
	make elf.3D.bsch.right.rot105;make elf.3D.bsch.right.rot106;make elf.3D.bsch.right.rot107;make elf.3D.bsch.right.rot108;make elf.3D.bsch.right.rot109;\
	make elf.3D.bsch.right.rot110;make elf.3D.bsch.right.rot111;make elf.3D.bsch.right.rot112;make elf.3D.bsch.right.rot113;make elf.3D.bsch.right.rot114;\
	make elf.3D.bsch.right.rot115;make elf.3D.bsch.right.rot116;make elf.3D.bsch.right.rot117;make elf.3D.bsch.right.rot118;make elf.3D.bsch.right.rot119;\
	make elf.3D.bsch.right.rot120;make elf.3D.bsch.right.rot121;make elf.3D.bsch.right.rot122;make elf.3D.bsch.right.rot123;make elf.3D.bsch.right.rot124;\
	make elf.3D.bsch.right.rot125;make elf.3D.bsch.right.rot126;make elf.3D.bsch.right.rot127;make elf.3D.bsch.right.rot128;make elf.3D.bsch.right.rot129;\
	make elf.3D.bsch.right.rot130;make elf.3D.bsch.right.rot131;make elf.3D.bsch.right.rot132


elf.3D.bsch_right.rot:
	Cat axis=4 elf.3D.bsch.right.rot000.H elf.3D.bsch.right.rot001.H elf.3D.bsch.right.rot002.H elf.3D.bsch.right.rot003.H elf.3D.bsch.right.rot004.H elf.3D.bsch.right.rot005.H elf.3D.bsch.right.rot006.H elf.3D.bsch.right.rot007.H elf.3D.bsch.right.rot008.H elf.3D.bsch.right.rot009.H elf.3D.bsch.right.rot010.H elf.3D.bsch.right.rot011.H elf.3D.bsch.right.rot012.H elf.3D.bsch.right.rot013.H elf.3D.bsch.right.rot014.H elf.3D.bsch.right.rot015.H elf.3D.bsch.right.rot016.H elf.3D.bsch.right.rot017.H elf.3D.bsch.right.rot018.H elf.3D.bsch.right.rot019.H elf.3D.bsch.right.rot020.H elf.3D.bsch.right.rot021.H elf.3D.bsch.right.rot022.H elf.3D.bsch.right.rot023.H elf.3D.bsch.right.rot024.H elf.3D.bsch.right.rot025.H elf.3D.bsch.right.rot026.H elf.3D.bsch.right.rot027.H elf.3D.bsch.right.rot028.H elf.3D.bsch.right.rot029.H elf.3D.bsch.right.rot030.H elf.3D.bsch.right.rot031.H elf.3D.bsch.right.rot032.H elf.3D.bsch.right.rot033.H elf.3D.bsch.right.rot034.H elf.3D.bsch.right.rot035.H elf.3D.bsch.right.rot036.H elf.3D.bsch.right.rot037.H elf.3D.bsch.right.rot038.H elf.3D.bsch.right.rot039.H elf.3D.bsch.right.rot040.H elf.3D.bsch.right.rot041.H elf.3D.bsch.right.rot042.H elf.3D.bsch.right.rot043.H elf.3D.bsch.right.rot044.H elf.3D.bsch.right.rot045.H elf.3D.bsch.right.rot046.H elf.3D.bsch.right.rot047.H elf.3D.bsch.right.rot048.H > a.H
	Cat axis=4 elf.3D.bsch.right.rot049.H elf.3D.bsch.right.rot050.H elf.3D.bsch.right.rot051.H elf.3D.bsch.right.rot052.H elf.3D.bsch.right.rot053.H elf.3D.bsch.right.rot054.H elf.3D.bsch.right.rot055.H elf.3D.bsch.right.rot056.H elf.3D.bsch.right.rot057.H elf.3D.bsch.right.rot058.H elf.3D.bsch.right.rot059.H elf.3D.bsch.right.rot060.H elf.3D.bsch.right.rot061.H elf.3D.bsch.right.rot062.H elf.3D.bsch.right.rot063.H elf.3D.bsch.right.rot064.H elf.3D.bsch.right.rot065.H elf.3D.bsch.right.rot066.H elf.3D.bsch.right.rot067.H elf.3D.bsch.right.rot068.H elf.3D.bsch.right.rot069.H elf.3D.bsch.right.rot070.H elf.3D.bsch.right.rot071.H elf.3D.bsch.right.rot072.H elf.3D.bsch.right.rot073.H elf.3D.bsch.right.rot074.H elf.3D.bsch.right.rot075.H elf.3D.bsch.right.rot076.H elf.3D.bsch.right.rot077.H elf.3D.bsch.right.rot078.H elf.3D.bsch.right.rot079.H elf.3D.bsch.right.rot080.H elf.3D.bsch.right.rot081.H elf.3D.bsch.right.rot082.H elf.3D.bsch.right.rot083.H elf.3D.bsch.right.rot084.H elf.3D.bsch.right.rot085.H elf.3D.bsch.right.rot086.H elf.3D.bsch.right.rot087.H elf.3D.bsch.right.rot088.H elf.3D.bsch.right.rot089.H elf.3D.bsch.right.rot090.H elf.3D.bsch.right.rot091.H elf.3D.bsch.right.rot092.H elf.3D.bsch.right.rot093.H elf.3D.bsch.right.rot094.H elf.3D.bsch.right.rot095.H elf.3D.bsch.right.rot096.H elf.3D.bsch.right.rot097.H > b.H
	Cat axis=4 elf.3D.bsch.right.rot098.H elf.3D.bsch.right.rot099.H elf.3D.bsch.right.rot100.H elf.3D.bsch.right.rot101.H elf.3D.bsch.right.rot102.H elf.3D.bsch.right.rot103.H elf.3D.bsch.right.rot104.H elf.3D.bsch.right.rot105.H elf.3D.bsch.right.rot106.H elf.3D.bsch.right.rot107.H elf.3D.bsch.right.rot108.H elf.3D.bsch.right.rot109.H elf.3D.bsch.right.rot110.H elf.3D.bsch.right.rot111.H elf.3D.bsch.right.rot112.H elf.3D.bsch.right.rot113.H elf.3D.bsch.right.rot114.H elf.3D.bsch.right.rot115.H elf.3D.bsch.right.rot116.H elf.3D.bsch.right.rot117.H elf.3D.bsch.right.rot118.H elf.3D.bsch.right.rot119.H elf.3D.bsch.right.rot120.H elf.3D.bsch.right.rot121.H elf.3D.bsch.right.rot122.H elf.3D.bsch.right.rot123.H elf.3D.bsch.right.rot124.H elf.3D.bsch.right.rot125.H elf.3D.bsch.right.rot126.H elf.3D.bsch.right.rot127.H elf.3D.bsch.right.rot128.H elf.3D.bsch.right.rot129.H elf.3D.bsch.right.rot130.H elf.3D.bsch.right.rot131.H elf.3D.bsch.right.rot132.H > c.H
	Cat a.H b.H c.H axis=4 > $@.H
	Rm a.H b.H c.H 

#
# Base chalk
elf.3D.bsch.rot.H:
	echo d4=30 o4=1600 >> elf.3D.bsch_right.rot.H
	echo d4=30 o4=1600 >> elf.3D.bsch_left.rot.H
	Math file1=elf.3D.bsch_right.rot.H file2=elf.3D.bsch_left.rot.H exp="file1+file2"|Transp plane=12 reshape=1,5 > j0.H
	< j0.H Pad beg1=0 end1=3 > $@
	echo d4=30 >> $@

elf.3D.bsch.rot.scale.H:
	< elf.3D.bsch.rot.H ${Tese}/bin/OW_3dscale.x > $@

elf.3D.bsch.window.H:
	Math file1=elf.3D.bsch.window_right.int30.H file2=elf.3D.bsch.window_left.int30.H exp=file1+file2 | Transp plane=12 reshape=1,5 > a.H
	<a.H Pad beg1=0 end1=3 > $@

elf.3D.bsch.window.scale.H:
	Math file1=elf.3D.bsch.window_right.int30.H file2=elf.3D.bsch.window_left.int30.H exp=file1+file2 | ${Tese}/bin/OW_3dscale.x | Transp plane=12 reshape=1,5 > a.H
	<a.H Pad beg1=0 end1=3 > $@

#elf.3D.bsch.rot%:
#	Window3d n3=1 f3=$* < elf.3D.bsch.window.scale.H | Transp plane=12 reshape=2,3,5 > a$*.H
#	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a$*.H Cig_tilt=elf.3D.bsch.rot$*.H par=bsch_right-rot.P positive_rot=1
#	Rm a$*.H

elf.3D.bschrot:
	make elf.3D.bsch.rot000;make elf.3D.bsch.rot001;make elf.3D.bsch.rot002;make elf.3D.bsch.rot003;make elf.3D.bsch.rot004;\
	make elf.3D.bsch.rot005;make elf.3D.bsch.rot006;make elf.3D.bsch.rot007;make elf.3D.bsch.rot008;make elf.3D.bsch.rot009;\
	make elf.3D.bsch.rot010;make elf.3D.bsch.rot011;make elf.3D.bsch.rot012;make elf.3D.bsch.rot013;make elf.3D.bsch.rot014;\
	make elf.3D.bsch.rot015;make elf.3D.bsch.rot016;make elf.3D.bsch.rot017;make elf.3D.bsch.rot018;make elf.3D.bsch.rot019;\
	make elf.3D.bsch.rot020;make elf.3D.bsch.rot021;make elf.3D.bsch.rot022;make elf.3D.bsch.rot023;make elf.3D.bsch.rot024;\
	make elf.3D.bsch.rot025;make elf.3D.bsch.rot026;make elf.3D.bsch.rot027;make elf.3D.bsch.rot028;make elf.3D.bsch.rot029;\
	make elf.3D.bsch.rot030;make elf.3D.bsch.rot031;make elf.3D.bsch.rot032;make elf.3D.bsch.rot033;make elf.3D.bsch.rot034;\
	make elf.3D.bsch.rot035;make elf.3D.bsch.rot036;make elf.3D.bsch.rot037;make elf.3D.bsch.rot038;make elf.3D.bsch.rot039;\
	make elf.3D.bsch.rot040;make elf.3D.bsch.rot041;make elf.3D.bsch.rot042;make elf.3D.bsch.rot043;make elf.3D.bsch.rot044;\
	make elf.3D.bsch.rot045;make elf.3D.bsch.rot046;make elf.3D.bsch.rot047;make elf.3D.bsch.rot048;make elf.3D.bsch.rot049;\
	make elf.3D.bsch.rot050;make elf.3D.bsch.rot051;make elf.3D.bsch.rot052;make elf.3D.bsch.rot053;make elf.3D.bsch.rot054;\
	make elf.3D.bsch.rot055;make elf.3D.bsch.rot056;make elf.3D.bsch.rot057;make elf.3D.bsch.rot058;make elf.3D.bsch.rot059;\
	make elf.3D.bsch.rot060;make elf.3D.bsch.rot061;make elf.3D.bsch.rot062;make elf.3D.bsch.rot063;make elf.3D.bsch.rot064;\
	make elf.3D.bsch.rot065;make elf.3D.bsch.rot066;make elf.3D.bsch.rot067;make elf.3D.bsch.rot068;make elf.3D.bsch.rot069;\
	make elf.3D.bsch.rot070;make elf.3D.bsch.rot071;make elf.3D.bsch.rot072;make elf.3D.bsch.rot073;make elf.3D.bsch.rot074;\
	make elf.3D.bsch.rot075;make elf.3D.bsch.rot076;make elf.3D.bsch.rot077;make elf.3D.bsch.rot078;make elf.3D.bsch.rot079;\
	make elf.3D.bsch.rot080;make elf.3D.bsch.rot081;make elf.3D.bsch.rot082;make elf.3D.bsch.rot083;make elf.3D.bsch.rot084;\
	make elf.3D.bsch.rot085;make elf.3D.bsch.rot086;make elf.3D.bsch.rot087;make elf.3D.bsch.rot088;make elf.3D.bsch.rot089;\
	make elf.3D.bsch.rot090;make elf.3D.bsch.rot091;make elf.3D.bsch.rot092;make elf.3D.bsch.rot093;make elf.3D.bsch.rot094;\
	make elf.3D.bsch.rot095;make elf.3D.bsch.rot096;make elf.3D.bsch.rot097;make elf.3D.bsch.rot098;make elf.3D.bsch.rot099;\
	make elf.3D.bsch.rot100;make elf.3D.bsch.rot101;make elf.3D.bsch.rot102;make elf.3D.bsch.rot103;make elf.3D.bsch.rot104;\
	make elf.3D.bsch.rot105;make elf.3D.bsch.rot106;make elf.3D.bsch.rot107;make elf.3D.bsch.rot108;make elf.3D.bsch.rot109;\
	make elf.3D.bsch.rot110;make elf.3D.bsch.rot111;make elf.3D.bsch.rot112;make elf.3D.bsch.rot113;make elf.3D.bsch.rot114;\
	make elf.3D.bsch.rot115;make elf.3D.bsch.rot116;make elf.3D.bsch.rot117;make elf.3D.bsch.rot118;make elf.3D.bsch.rot119;\
	make elf.3D.bsch.rot120;make elf.3D.bsch.rot121;make elf.3D.bsch.rot122;make elf.3D.bsch.rot123;make elf.3D.bsch.rot124;\
	make elf.3D.bsch.rot125;make elf.3D.bsch.rot126;make elf.3D.bsch.rot127;make elf.3D.bsch.rot128;make elf.3D.bsch.rot129;\
	make elf.3D.bsch.rot130;make elf.3D.bsch.rot131;make elf.3D.bsch.rot132


elf.3Dbschrot:
	Cat axis=4 elf.3D.bsch.rot000.H elf.3D.bsch.rot001.H elf.3D.bsch.rot002.H elf.3D.bsch.rot003.H elf.3D.bsch.rot004.H elf.3D.bsch.rot005.H elf.3D.bsch.rot006.H elf.3D.bsch.rot007.H elf.3D.bsch.rot008.H elf.3D.bsch.rot009.H elf.3D.bsch.rot010.H elf.3D.bsch.rot011.H elf.3D.bsch.rot012.H elf.3D.bsch.rot013.H elf.3D.bsch.rot014.H elf.3D.bsch.rot015.H elf.3D.bsch.rot016.H elf.3D.bsch.rot017.H elf.3D.bsch.rot018.H elf.3D.bsch.rot019.H elf.3D.bsch.rot020.H elf.3D.bsch.rot021.H elf.3D.bsch.rot022.H elf.3D.bsch.rot023.H elf.3D.bsch.rot024.H elf.3D.bsch.rot025.H elf.3D.bsch.rot026.H elf.3D.bsch.rot027.H elf.3D.bsch.rot028.H elf.3D.bsch.rot029.H elf.3D.bsch.rot030.H elf.3D.bsch.rot031.H elf.3D.bsch.rot032.H elf.3D.bsch.rot033.H elf.3D.bsch.rot034.H elf.3D.bsch.rot035.H elf.3D.bsch.rot036.H elf.3D.bsch.rot037.H elf.3D.bsch.rot038.H elf.3D.bsch.rot039.H elf.3D.bsch.rot040.H elf.3D.bsch.rot041.H elf.3D.bsch.rot042.H elf.3D.bsch.rot043.H elf.3D.bsch.rot044.H elf.3D.bsch.rot045.H elf.3D.bsch.rot046.H elf.3D.bsch.rot047.H elf.3D.bsch.rot048.H > a.H
	Cat axis=4 elf.3D.bsch.rot049.H elf.3D.bsch.rot050.H elf.3D.bsch.rot051.H elf.3D.bsch.rot052.H elf.3D.bsch.rot053.H elf.3D.bsch.rot054.H elf.3D.bsch.rot055.H elf.3D.bsch.rot056.H elf.3D.bsch.rot057.H elf.3D.bsch.rot058.H elf.3D.bsch.rot059.H elf.3D.bsch.rot060.H elf.3D.bsch.rot061.H elf.3D.bsch.rot062.H elf.3D.bsch.rot063.H elf.3D.bsch.rot064.H elf.3D.bsch.rot065.H elf.3D.bsch.rot066.H elf.3D.bsch.rot067.H elf.3D.bsch.rot068.H elf.3D.bsch.rot069.H elf.3D.bsch.rot070.H elf.3D.bsch.rot071.H elf.3D.bsch.rot072.H elf.3D.bsch.rot073.H elf.3D.bsch.rot074.H elf.3D.bsch.rot075.H elf.3D.bsch.rot076.H elf.3D.bsch.rot077.H elf.3D.bsch.rot078.H elf.3D.bsch.rot079.H elf.3D.bsch.rot080.H elf.3D.bsch.rot081.H elf.3D.bsch.rot082.H elf.3D.bsch.rot083.H elf.3D.bsch.rot084.H elf.3D.bsch.rot085.H elf.3D.bsch.rot086.H elf.3D.bsch.rot087.H elf.3D.bsch.rot088.H elf.3D.bsch.rot089.H elf.3D.bsch.rot090.H elf.3D.bsch.rot091.H elf.3D.bsch.rot092.H elf.3D.bsch.rot093.H elf.3D.bsch.rot094.H elf.3D.bsch.rot095.H elf.3D.bsch.rot096.H elf.3D.bsch.rot097.H > b.H
	Cat axis=4 elf.3D.bsch.rot098.H elf.3D.bsch.rot099.H elf.3D.bsch.rot100.H elf.3D.bsch.rot101.H elf.3D.bsch.rot102.H elf.3D.bsch.rot103.H elf.3D.bsch.rot104.H elf.3D.bsch.rot105.H elf.3D.bsch.rot106.H elf.3D.bsch.rot107.H elf.3D.bsch.rot108.H elf.3D.bsch.rot109.H elf.3D.bsch.rot110.H elf.3D.bsch.rot111.H elf.3D.bsch.rot112.H elf.3D.bsch.rot113.H elf.3D.bsch.rot114.H elf.3D.bsch.rot115.H elf.3D.bsch.rot116.H elf.3D.bsch.rot117.H elf.3D.bsch.rot118.H elf.3D.bsch.rot119.H elf.3D.bsch.rot120.H elf.3D.bsch.rot121.H elf.3D.bsch.rot122.H elf.3D.bsch.rot123.H elf.3D.bsch.rot124.H elf.3D.bsch.rot125.H elf.3D.bsch.rot126.H elf.3D.bsch.rot127.H elf.3D.bsch.rot128.H elf.3D.bsch.rot129.H elf.3D.bsch.rot130.H elf.3D.bsch.rot131.H elf.3D.bsch.rot132.H > c.H
	Cat a.H b.H c.H axis=4 > $@.H

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

bschlk00.permr%:
	Wavelet n1=1001 d1=0.004 wavelet=ricker2 phase=0 fund=15 order=8 tdelay=2 domain=time > j0.H;echo o1=-2 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=32 f_max=34 wei=y is2d=y | Window3d | Cabs > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > bsch.sou.H
	${Tese}/bin/OW3d_arealsht.x cig=elf.3D.bsch.rot.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.perm.rec.gx$*gy1.rot.H sou_out=elf.perm.sou.gx$*gy1.rot.H \
	ngapx=$* ngapy=1 par=bsch_perm.P node=14 shape=1

bschlk_eps0.permr.%:
	Wavelet n1=1001 d1=0.004 wavelet=ricker2 phase=0 fund=15 order=8 tdelay=2 domain=time > j0.H;echo o1=-2 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=32 f_max=34 wei=y is2d=y | Window3d | Cabs > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > bsch.sou.H
	Window3d < elf.3D.bsch.rot.H min2=-300 n2=21 squeeze=n > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H elf.3D.bsch.rot.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.perm.rec.gx$*gy1.rot.H sou_out=elf.perm.sou.gx$*gy1.rot.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

#bschlk_eps0a.permr.%:
#	Wavelet n1=1001 d1=0.004 wavelet=ricker2 phase=0 fund=25 order=8 tdelay=2 domain=time > j0.H;echo o1=-2 >> j0.H
#	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=42 f_max=46 wei=y is2d=y | Window3d | Cabs | Smooth rect1=11 > j1.H
#	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > bsch.sou.H
#	Window3d min3=-300 n3=21 squeeze=n < elf.3Dbschrot.H | ${Tese}/bin/OW_3dscale.x | Transp plane=12 reshape=1,5 > cig.H
#	echo d3=30 d4=30 >> cig.H
#	${Tese}/bin/OW3d_arealsht.x cig=cig.H elf.3D.bsch.rot.H sou=bsch.sou.H \
#	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.perm.rec.gx$*gy1.f25.rot.H sou_out=elf.perm.sou.gx$*gy1.f25.rot.H \
#	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

bschlk_eps0.perm.%:
	Wavelet n1=1001 d1=0.004 wavelet=ricker2 phase=0 fund=15 order=8 tdelay=2 domain=time > j0.H;echo o1=-2 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=32 f_max=34 wei=y is2d=y | Window3d | Cabs > j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > bsch.sou.H
	Window3d < elf.3D.bsch.window.scale.H min2=-300 n2=21 squeeze=n > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H elf.3D.bsch.rot.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.perm.rec.gx$*gy1.H sou_out=elf.perm.sou.gx$*gy1.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

elf.3D.mig_adj.bsch01.H:
#	Window3d squeeze=n < elf.inv_freq.recip.H f5=1 > data.H
#	< elf.velopt01.CAM.H Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120> \
#	elf.refvelopt01.CAM.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data.H,elf.velopt01.CAM.H,elf.refvelopt01.CAM.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="5,4,4,6" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=355 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull360.par

#bschlk00.perm:
##	Scale < elf.L7f.vinit.grid30.H dscale=2 > vel.H
#	Wavelet n1=501 d1=0.008 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=2. domain=time > j0.H;echo o1=-2 >> j0.H
#	< j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=32 f_max=38 wei=y is2d=y | Window3d | Cabs | Smooth rect1=11 > j1.H
#	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > sou.H
#	${Tese}/bin/OW3d_arealsht.x cig=elf.3D.bschlk01.flt.int30.sort.H sou=sou.H vel=vel.H \
#	rec_out=elf.3D.bsck01.rec.perm.H sout_out=elf.3D.bsck01.sou.perm.H par=3dperm10_par.P node=10 > /dev/null

#	files="elf.3D.bschlk01.flt.int30.sort.H,sou.H,vel.H,elf.3D.bsck00.rec.perm.H,elf.3D.bsck00.sou.perm.H" \


###################################################
#
# Residual moveout analysis
#
#0250elf.rho.H:
#	Window3d min1=0 max1=13000 < 0200elf.bkim.4000.plane.H|Transp plane=13|Transp plane=23 reshape=2,4,5|\
#	Interp d2out=3.375 Pad beg1=100 end1=100 beg2=250 end2=250|OFF2ANG amax=60 eps=10 verb=y maxmem=1000 > x1.H
#	Window3d < x1.H min1=0 max1=4000 min2=-40 max2=40 j4=10 | Agc window=50 | ../Tese/bin/OW_rhopanel.x method=gcg \
#	exponent=.5 eps=1 niter=50 orho=0.8 drho=0.02 nrho=201 > m.H

3Delf.sediment.scan.rho.H:
	< elf.3D.000to356.flt.int30.input2window.H Window3d n4=181 min1= 5000 max1= 7000 min3= 2000 max3= 3000|Transp plane=34 > xxx.H
	${Tese}/bin/OW_3dresmig.x bkgimg=xxx.H Dimage=$@ node=1 2 6 nrho=101 orho=0.8 drho=0.004 > $@


############################################
#
# Amplitude analysis
#
amp_expl:
	Window3d < data.H n3=1 min3=0 squeeze=n | Transf maxsize=2000 | Window3d | Transp plane=23 > elf.off0000.H
	< elf.off0000.H Bandpass flo=4 fhi=10 | Envelope rect1=101 rect2=1 rect3=1 | Grey eout=1 pclip=90 gainpanel=a > elf.off0000.bp.env.T
	< elf.off0000.H Bandpass flo=4 fhi=10|Agc window=400|Envelope rect1=101 rect2=1 rect3=1|Grey eout=1 pclip= 100 gainpanel=a>elf.off0000.bp.agc400.env.T
	cube elf.off0000.bp.env.T elf.off0000.bp.agc400.env.T nviews=2 &
	Window3d < data.H n3=1 min3=1000 squeeze=n | Transf maxsize=2000 | Window3d | Transp plane=23 > elf.off1000.H
	< elf.off1000.H Bandpass flo=4 fhi=10 | Envelope rect1=101 rect2=1 rect3=1 | Grey eout=1 pclip= 100 gainpanel=a > elf.off1000.bp.env.T
	< elf.off1000.H Bandpass flo=4 fhi=10|Agc window=400|Envelope rect1=101 rect2=1 rect3=1|Grey eout=1 pclip= 100 gainpanel=a>elf.off1000.bp.agc400.env.T
	cube elf.off2000.bp.env.T elf.off2000.bp.agc400.env.T nviews=2 &

analysis%:
	Window3d < data.H n3=1 f3=$* squeeze=n|Transf maxsize=2000|Window3d|Transp plane=23|Bandpass flo=4 fhi=10|Envelope rect1=101 rect2=1 rect3=1|Transp|\
	Stack|Window3d|Scale scale_to=1|Clip clip=0.05 to=1|Gpow gpow=10> o$*.H

elf.mask.agc.H:
	Window3d < o70.H f2=250 max2=12000 | Pad beg2=250 end2=0 extend=1 | Pad beg1=0 end2=72 > o70a.H
	Window3d < o00.H f2=100 | Pad beg2=100 end2=0 extend=1 > o00a.H
	Add o00a.H o70a.H mode=prod |Transp> $@

data.mask.H:
	< data.H ${Tese}/bin/OW_3dmask.x mask=elf.mask.agc.H > $@

data.agc.H:
#	< data.mask.H Transf | Agc window=400 > j.H
	< j.H Transf  f_min=1 f_min1=4 f_max1=52 f_max=60 wei=1 > $@

data.agcbf%.H:
	< data.mask.H Window3d squeeze=n n3=3 f3=$* | Transf maxsize=2000 > a$*.H
	< a$*.H Agc window=400 | Window3d max1=5 | Bandpass flo=1 fhi=40 > j$*.H
	< j$*.H Transf f_min=1 f_min1=4 f_max1=32 f_max=40 wei=1 maxsize=2000 > a$*.H
	< a$*.H Transp plane=35 maxsize=2000 > $@
	Rm a$*.H j$*.H

ddd0:
	make data.agcbf000.H; make data.agcbf003.H;make data.agcbf006.H; make data.agcbf009.H;make data.agcbf012.H; make data.agcbf015.H
	
ddd1:
	make data.agcbf018.H; make data.agcbf021.H;make data.agcbf024.H; make data.agcbf027.H;make data.agcbf030.H; make data.agcbf033.H

ddd2:
	make data.agcbf036.H; make data.agcbf039.H;make data.agcbf042.H; make data.agcbf045.H;make data.agcbf048.H; make data.agcbf051.H

ddd3:
	make data.agcbf054.H; make data.agcbf057.H;make data.agcbf060.H; make data.agcbf063.H;make data.agcbf066.H; make data.agcbf069.H

	
data_agcbf.H:
	Cat data.agcbf000.H data.agcbf003.H data.agcbf006.H data.agcbf009.H data.agcbf012.H data.agcbf015.H \
	data.agcbf018.H data.agcbf021.H data.agcbf024.H data.agcbf027.H data.agcbf030.H data.agcbf033.H \
	data.agcbf036.H data.agcbf039.H data.agcbf042.H data.agcbf045.H data.agcbf048.H data.agcbf051.H \
	data.agcbf054.H data.agcbf057.H data.agcbf060.H data.agcbf063.H data.agcbf066.H data.agcbf069.H \
	axis=5 max_memory=2000 > $@

data_agcbf.recip.H:
#	echo n2=42336 n3=72 n4=1 n5=1 o3=0 o4=0.800364 o5=0 d3=25 d4=0.200091 d5=25 >> data_agcbf.H
#	Reverse < data_agcbf.H which=4 opt=no | Window3d max3=0 > w.H
	< data_agcbf.H ${Tese}/bin/merge.x > z.H
	Transp < z.H plane=35 maxsize=2000 > $@

######################
#
# Abstract SEG2010
#
# Migration run on cees-cluster
#
elf.3D.mig.velinit.agc.ang.H:
#	Transp < elf.3D.mig.velinit.agc.H reshape=2,4,5 plane=13 maxsize=2000 > a00.H
	< a00.H Pad beg2=50 end2=50 | OFF2ANG amax=40 verb=1 eps=10 maxmem=2000 > x00.H
	< x00.H Window3d min2=0 > $@

elf.3D.mig.velopt01.bsch.agc.ang.H:
	Transp < elf.3D.mig.velopt01.bsch.agc.H reshape=2,4,5 plane=13 maxsize=2000 > a00.H
	< a00.H Pad beg2=50 end2=50 | OFF2ANG amax=40 verb=1 eps=10 > x00.H
	< x00.H Window3d min2=0 > $@

#########################
#
# After rms scaling
#
elf.3Dbschrot.scale:
	echo o4=1600 >> elf.3Dbschrot.H
	Window3d min3=-300 n3=21 < elf.3Dbschrot.H|Transp plane=12 reshape=1,5 maxsize=2000 > a.H
	< a.H Scale scale_to=1 | ${Tese}/bin/OW_3dscale.x wx=5 wy=5 scale=elf.scale.bsch.H > $@.H ;Rm a.H
	echo o4=0 d4=30 >> $@.H

bsch.sou.H:
	Wavelet n1=1001 d1=0.002 wavelet=ricker2 phase=0 fund=25 order=8 tdelay=2 domain=time > j0.H;
	echo o1=-2 >> j0.H
	<j0.H Scale dclip=-1|Transf f_min=0 f_min1=4 f_max1=52 f_max=60 wei=y|\
	Window3d|Cabs|Smooth rect1=11>j1.H
	Math file1=j1.H exp="file1*0" > j2.H ; Cmplx j1.H j2.H | Transp plane=13 > $@

bschlk_eps0a.permr.%:
	${Tese}/bin/OW3d_arealsht.x cig=elf.3Dbschrot.scale.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.perm.rec.gx$*gy1.f25.rot.scale.H \
	sou_out=elf.perm.sou.gx$*gy1.f25.rot.scale.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

bschlk.ispew.%:
	${Tese}/bin/OW3d_arealsht.x cig=elf.3Dbschrot.scale.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf.ispew.rec.gx$*gy1.f25.rot.H \
	sou_out=elf.ispew.sou.gx$*gy1.f25.rot.H \
	ngapx=$* ngapy=1 par=bsch_perm.P node=14 shape=1

bsch_agc%:
	<elf.perm.sou.gx$*gy1.f25.rot.scale.H Transp plane=35 maxsize=2000 > 0.H
	< 0.H Transf maxsize=2000 > a.H
	< a.H Agc window=254 > b.H
	< b.H Transf f_min=0 f_min1=4 f_max1=42 f_max=46 wei=y maxsize=2000 > c.H
	< c.H Transp plane=35 maxsize=2000 > elf.perm.sou.gx$*gy1.f25.rot.scale.rms.H
	< elf.perm.rec.gx$*gy1.f25.rot.scale.H Transp plane=35 maxsize=2000 > 0.H
	< 0.H Transf maxsize=2000 > a.H
	< a.H Agc window=254 > b.H
	< b.H Transf f_min=0 f_min1=4 f_max1=42 f_max=46 wei=y maxsize=2000 > c.H
	< c.H Transp plane=35 maxsize=2000 > elf.perm.rec.gx$*gy1.f25.rot.scale.rms.H

##########################################
#
# New windowing. The previous was too big
# and generated reflector crosstalk
#
elf.3D.mig.rmo2.int30a.H:
	 < elf.3D.mig.rmo2.int30.H ../Tese/bin/OW_square.x > $@

elf.3D.mig.rmo2.int30a.new.H:
	 < elf.3D.mig.rmo2.int30.new.H ../Tese/bin/OW_square.x > $@

elf000.bsch.window_right.mask%.H:
	Window3d < elf.3D.mig.rmo2.int30a.H n3=1 f3=$* | Transp plane=13 > h0$*.H
	< elf.3D.bsch.window_right.H Window3d n2=1 f2=$* | Transp plane=13 | Agc window=200 | Envelope tc1=1 tc2=2 tc3=2 verb=1 > h1$*.H
	< h1$*.H Scale scale_to=1|Clip to=1 clip=0.05 chop=greater | Clip to=0 clip=0.05 chop=less|Smooth rect1=1 rect2=2 rect3=2 > h2$*.H
	Math file1=h0$*.H file2=h2$*.H exp="file1*file2" > $@
#	Rm h[0-2]$*.H 

elf000.bsch.window_right1:
	make elf000.bsch.window_right.mask00.H;make elf000.bsch.window_right.mask01.H;make elf000.bsch.window_right.mask02.H;make elf000.bsch.window_right.mask03.H;make elf000.bsch.window_right.mask04.H;\
	make elf000.bsch.window_right.mask05.H;make elf000.bsch.window_right.mask06.H;make elf000.bsch.window_right.mask07.H;make elf000.bsch.window_right.mask08.H;make elf000.bsch.window_right.mask09.H;

elf000.bsch.window_right2:
	make elf000.bsch.window_right.mask10.H;make elf000.bsch.window_right.mask11.H;make elf000.bsch.window_right.mask12.H;make elf000.bsch.window_right.mask13.H;make elf000.bsch.window_right.mask14.H;\
	make elf000.bsch.window_right.mask15.H;make elf000.bsch.window_right.mask16.H;make elf000.bsch.window_right.mask17.H;make elf000.bsch.window_right.mask18.H;make elf000.bsch.window_right.mask19.H;

elf000.bsch.window_right3:
	make elf000.bsch.window_right.mask20.H;make elf000.bsch.window_right.mask21.H;make elf000.bsch.window_right.mask22.H;make elf000.bsch.window_right.mask23.H;make elf000.bsch.window_right.mask24.H;\
	make elf000.bsch.window_right.mask25.H;make elf000.bsch.window_right.mask26.H;make elf000.bsch.window_right.mask27.H;make elf000.bsch.window_right.mask28.H;make elf000.bsch.window_right.mask29.H;\
	make elf000.bsch.window_right.mask30.H

elf000.bsch.window_right.int30.H:
	Cat elf000.bsch.window_right.mask00.H elf000.bsch.window_right.mask01.H elf000.bsch.window_right.mask02.H elf000.bsch.window_right.mask03.H elf000.bsch.window_right.mask04.H \
	elf000.bsch.window_right.mask05.H elf000.bsch.window_right.mask06.H elf000.bsch.window_right.mask07.H elf000.bsch.window_right.mask08.H elf000.bsch.window_right.mask09.H \
	elf000.bsch.window_right.mask10.H elf000.bsch.window_right.mask11.H elf000.bsch.window_right.mask12.H elf000.bsch.window_right.mask13.H elf000.bsch.window_right.mask14.H \
	elf000.bsch.window_right.mask15.H elf000.bsch.window_right.mask16.H elf000.bsch.window_right.mask17.H elf000.bsch.window_right.mask18.H elf000.bsch.window_right.mask19.H \
	elf000.bsch.window_right.mask20.H elf000.bsch.window_right.mask21.H elf000.bsch.window_right.mask22.H elf000.bsch.window_right.mask23.H elf000.bsch.window_right.mask24.H \
	elf000.bsch.window_right.mask25.H elf000.bsch.window_right.mask26.H elf000.bsch.window_right.mask27.H elf000.bsch.window_right.mask28.H elf000.bsch.window_right.mask29.H \
	elf000.bsch.window_right.mask30.H axis=4 | Transp plane=23 reshape=1,2,4,5 > $@
	Rm elf000.bsch.window_right.mask*.H

elf000.bsch.window_left.mask%.H:
	Window3d < elf.3D.mig.rmo2.int30a.H n3=1 f3=$* | Transp plane=13 > j0$*.H
	< elf.3D.bsch.window_left.H Window3d n2=1 f2=$* | Transp plane=13 | Agc window=200 | Envelope tc1=1 tc2=2 tc3=2 verb=1 > j1$*.H
	< j1$*.H Scale scale_to=1|Clip to=1 clip=0.05 chop=greater | Clip to=0 clip=0.05 chop=less|Smooth rect1=1 rect2=2 rect3=2 > j2$*.H
	Math file1=j0$*.H file2=j2$*.H exp="file1*file2" > $@
#	Rm j[0-2]$*.H 

elf000.bsch.window_left1:
	make elf000.bsch.window_left.mask00.H;make elf000.bsch.window_left.mask01.H;make elf000.bsch.window_left.mask02.H;make elf000.bsch.window_left.mask03.H;make elf000.bsch.window_left.mask04.H;\
	make elf000.bsch.window_left.mask05.H;make elf000.bsch.window_left.mask06.H;make elf000.bsch.window_left.mask07.H;make elf000.bsch.window_left.mask08.H;make elf000.bsch.window_left.mask09.H;\

elf000.bsch.window_left2:
	make elf000.bsch.window_left.mask10.H;make elf000.bsch.window_left.mask11.H;make elf000.bsch.window_left.mask12.H;make elf000.bsch.window_left.mask13.H;make elf000.bsch.window_left.mask14.H;\
	make elf000.bsch.window_left.mask15.H;make elf000.bsch.window_left.mask16.H;make elf000.bsch.window_left.mask17.H;make elf000.bsch.window_left.mask18.H;make elf000.bsch.window_left.mask19.H;\
	make elf000.bsch.window_left.mask20.H;make elf000.bsch.window_left.mask21.H;

elf000.bsch.window_left3:
	make elf000.bsch.window_left.mask22.H;make elf000.bsch.window_left.mask23.H;make elf000.bsch.window_left.mask24.H;\
	make elf000.bsch.window_left.mask25.H;make elf000.bsch.window_left.mask26.H;make elf000.bsch.window_left.mask27.H;make elf000.bsch.window_left.mask28.H;make elf000.bsch.window_left.mask29.H;\
	make elf000.bsch.window_left.mask30.H

elf000.bsch.window_left.int30.H:
	Cat elf000.bsch.window_left.mask00.H elf000.bsch.window_left.mask01.H elf000.bsch.window_left.mask02.H elf000.bsch.window_left.mask03.H elf000.bsch.window_left.mask04.H \
	elf000.bsch.window_left.mask05.H elf000.bsch.window_left.mask06.H elf000.bsch.window_left.mask07.H elf000.bsch.window_left.mask08.H elf000.bsch.window_left.mask09.H \
	elf000.bsch.window_left.mask10.H elf000.bsch.window_left.mask11.H elf000.bsch.window_left.mask12.H elf000.bsch.window_left.mask13.H elf000.bsch.window_left.mask14.H \
	elf000.bsch.window_left.mask15.H elf000.bsch.window_left.mask16.H elf000.bsch.window_left.mask17.H elf000.bsch.window_left.mask18.H elf000.bsch.window_left.mask19.H \
	elf000.bsch.window_left.mask20.H elf000.bsch.window_left.mask21.H elf000.bsch.window_left.mask22.H elf000.bsch.window_left.mask23.H elf000.bsch.window_left.mask24.H \
	elf000.bsch.window_left.mask25.H elf000.bsch.window_left.mask26.H elf000.bsch.window_left.mask27.H elf000.bsch.window_left.mask28.H elf000.bsch.window_left.mask29.H \
	elf000.bsch.window_left.mask30.H axis=4 | Transp plane=23 reshape=1,2,4,5 > $@
	Rm elf000.bsch.window_left.mask*.H

elf000.bsch.window.H:
	Math file1=elf000.bsch.window_right.int30.H file2=elf000.bsch.window_left.int30.H exp=file1+file2|Transp plane=12 reshape=1,5 > a.H
	<a.H Pad beg1=0 end1=3 > $@

elf000.bsch.window.scale.H:
	Math file1=elf000.bsch.window_right.int30.H file2=elf000.bsch.window_left.int30.H exp=file1+file2|\
	Transp plane=12 reshape=1,5 maxsize=1000 > b.H;< b.H ${Tese}/bin/OW_3dscale.x wx=1 wy=1 > a.H
	<a.H Pad beg1=0 end1=3 > $@

bschlk000_eps0a.perm.%:
	echo d4=30 >> elf000.bsch.window.scale.H
	Window3d squeeze=n < elf000.bsch.window.scale.H n5=281 > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf000.perm.rec.gx$*.scale.H \
	sou_out=elf000.perm.sou.gx$*.scale.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

##############################################
elf000.bsch.rot%:
	Window3d n3=1 f3=$* < elf000.bsch.window.scale.H | Transp plane=12 reshape=2,3,5 > a$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a$*.H Cig_tilt=elf000.bsch.rot$*.H par=bsch_right-rot.P positive_rot=1
	Rm a$*.H

elf000.bschrot1:
	make elf000.bsch.rot000;make elf000.bsch.rot001;make elf000.bsch.rot002;make elf000.bsch.rot003;make elf000.bsch.rot004;\
	make elf000.bsch.rot005;make elf000.bsch.rot006;make elf000.bsch.rot007;make elf000.bsch.rot008;make elf000.bsch.rot009;\
	make elf000.bsch.rot010;make elf000.bsch.rot011;make elf000.bsch.rot012;make elf000.bsch.rot013;make elf000.bsch.rot014;

elf000.bschrot2:
	make elf000.bsch.rot015;make elf000.bsch.rot016;make elf000.bsch.rot017;make elf000.bsch.rot018;make elf000.bsch.rot019;\
	make elf000.bsch.rot020;make elf000.bsch.rot021;make elf000.bsch.rot022;make elf000.bsch.rot023;make elf000.bsch.rot024;

elf000.bschrot3:
	make elf000.bsch.rot025;make elf000.bsch.rot026;make elf000.bsch.rot027;make elf000.bsch.rot028;make elf000.bsch.rot029;\
	make elf000.bsch.rot030;make elf000.bsch.rot031;make elf000.bsch.rot032;make elf000.bsch.rot033;make elf000.bsch.rot034;

elf000.bschrot4:
	make elf000.bsch.rot035;make elf000.bsch.rot036;make elf000.bsch.rot037;make elf000.bsch.rot038;make elf000.bsch.rot039;\
	make elf000.bsch.rot040;make elf000.bsch.rot041;make elf000.bsch.rot042;make elf000.bsch.rot043;make elf000.bsch.rot044;\
	make elf000.bsch.rot045;make elf000.bsch.rot046;make elf000.bsch.rot047;make elf000.bsch.rot048;make elf000.bsch.rot049;

elf000.bschrot5:
	make elf000.bsch.rot050;make elf000.bsch.rot051;make elf000.bsch.rot052;make elf000.bsch.rot053;make elf000.bsch.rot054;\
	make elf000.bsch.rot055;make elf000.bsch.rot056;make elf000.bsch.rot057;make elf000.bsch.rot058;make elf000.bsch.rot059;\
	make elf000.bsch.rot060;make elf000.bsch.rot061;make elf000.bsch.rot062;make elf000.bsch.rot063;make elf000.bsch.rot064;

elf000.bschrot6:
	make elf000.bsch.rot065;make elf000.bsch.rot066;make elf000.bsch.rot067;make elf000.bsch.rot068;make elf000.bsch.rot069;\
	make elf000.bsch.rot070;make elf000.bsch.rot071;make elf000.bsch.rot072;make elf000.bsch.rot073;make elf000.bsch.rot074;

elf000.bschrot7:
	make elf000.bsch.rot075;make elf000.bsch.rot076;make elf000.bsch.rot077;make elf000.bsch.rot078;make elf000.bsch.rot079;\
	make elf000.bsch.rot080;make elf000.bsch.rot081;make elf000.bsch.rot082;make elf000.bsch.rot083;make elf000.bsch.rot084;\
	make elf000.bsch.rot085;make elf000.bsch.rot086;make elf000.bsch.rot087;make elf000.bsch.rot088;make elf000.bsch.rot089;

elf000.bschrot8:
	make elf000.bsch.rot090;make elf000.bsch.rot091;make elf000.bsch.rot092;make elf000.bsch.rot093;make elf000.bsch.rot094;\
	make elf000.bsch.rot095;make elf000.bsch.rot096;make elf000.bsch.rot097;make elf000.bsch.rot098;make elf000.bsch.rot099;

elf000.bschrot9:
	make elf000.bsch.rot100;make elf000.bsch.rot101;make elf000.bsch.rot102;make elf000.bsch.rot103;make elf000.bsch.rot104;\
	make elf000.bsch.rot105;make elf000.bsch.rot106;make elf000.bsch.rot107;make elf000.bsch.rot108;make elf000.bsch.rot109;\
	make elf000.bsch.rot110;make elf000.bsch.rot111;make elf000.bsch.rot112;make elf000.bsch.rot113;make elf000.bsch.rot114;

elf000.bschrot10:
	make elf000.bsch.rot115;make elf000.bsch.rot116;make elf000.bsch.rot117;make elf000.bsch.rot118;make elf000.bsch.rot119;\
	make elf000.bsch.rot120;make elf000.bsch.rot121;make elf000.bsch.rot122;make elf000.bsch.rot123;make elf000.bsch.rot124;

elf000.bschrot11:
	make elf000.bsch.rot125;make elf000.bsch.rot126;make elf000.bsch.rot127;make elf000.bsch.rot128;make elf000.bsch.rot129;\
	make elf000.bsch.rot130;make elf000.bsch.rot131;make elf000.bsch.rot132


elf.3Dbschrot00:
	Cat axis=4 elf000.bsch.rot000.H elf000.bsch.rot001.H elf000.bsch.rot002.H elf000.bsch.rot003.H elf000.bsch.rot004.H elf000.bsch.rot005.H elf000.bsch.rot006.H elf000.bsch.rot007.H elf000.bsch.rot008.H elf000.bsch.rot009.H elf000.bsch.rot010.H elf000.bsch.rot011.H elf000.bsch.rot012.H elf000.bsch.rot013.H elf000.bsch.rot014.H elf000.bsch.rot015.H elf000.bsch.rot016.H elf000.bsch.rot017.H elf000.bsch.rot018.H elf000.bsch.rot019.H elf000.bsch.rot020.H elf000.bsch.rot021.H elf000.bsch.rot022.H elf000.bsch.rot023.H elf000.bsch.rot024.H elf000.bsch.rot025.H elf000.bsch.rot026.H elf000.bsch.rot027.H elf000.bsch.rot028.H elf000.bsch.rot029.H elf000.bsch.rot030.H elf000.bsch.rot031.H elf000.bsch.rot032.H elf000.bsch.rot033.H elf000.bsch.rot034.H elf000.bsch.rot035.H elf000.bsch.rot036.H elf000.bsch.rot037.H elf000.bsch.rot038.H elf000.bsch.rot039.H elf000.bsch.rot040.H elf000.bsch.rot041.H elf000.bsch.rot042.H elf000.bsch.rot043.H elf000.bsch.rot044.H elf000.bsch.rot045.H elf000.bsch.rot046.H elf000.bsch.rot047.H elf000.bsch.rot048.H > a.H
	Cat axis=4 elf000.bsch.rot049.H elf000.bsch.rot050.H elf000.bsch.rot051.H elf000.bsch.rot052.H elf000.bsch.rot053.H elf000.bsch.rot054.H elf000.bsch.rot055.H elf000.bsch.rot056.H elf000.bsch.rot057.H elf000.bsch.rot058.H elf000.bsch.rot059.H elf000.bsch.rot060.H elf000.bsch.rot061.H elf000.bsch.rot062.H elf000.bsch.rot063.H elf000.bsch.rot064.H elf000.bsch.rot065.H elf000.bsch.rot066.H elf000.bsch.rot067.H elf000.bsch.rot068.H elf000.bsch.rot069.H elf000.bsch.rot070.H elf000.bsch.rot071.H elf000.bsch.rot072.H elf000.bsch.rot073.H elf000.bsch.rot074.H elf000.bsch.rot075.H elf000.bsch.rot076.H elf000.bsch.rot077.H elf000.bsch.rot078.H elf000.bsch.rot079.H elf000.bsch.rot080.H elf000.bsch.rot081.H elf000.bsch.rot082.H elf000.bsch.rot083.H elf000.bsch.rot084.H elf000.bsch.rot085.H elf000.bsch.rot086.H elf000.bsch.rot087.H elf000.bsch.rot088.H elf000.bsch.rot089.H elf000.bsch.rot090.H elf000.bsch.rot091.H elf000.bsch.rot092.H elf000.bsch.rot093.H elf000.bsch.rot094.H elf000.bsch.rot095.H elf000.bsch.rot096.H elf000.bsch.rot097.H > b.H
	Cat axis=4 elf000.bsch.rot098.H elf000.bsch.rot099.H elf000.bsch.rot100.H elf000.bsch.rot101.H elf000.bsch.rot102.H elf000.bsch.rot103.H elf000.bsch.rot104.H elf000.bsch.rot105.H elf000.bsch.rot106.H elf000.bsch.rot107.H elf000.bsch.rot108.H elf000.bsch.rot109.H elf000.bsch.rot110.H elf000.bsch.rot111.H elf000.bsch.rot112.H elf000.bsch.rot113.H elf000.bsch.rot114.H elf000.bsch.rot115.H elf000.bsch.rot116.H elf000.bsch.rot117.H elf000.bsch.rot118.H elf000.bsch.rot119.H elf000.bsch.rot120.H elf000.bsch.rot121.H elf000.bsch.rot122.H elf000.bsch.rot123.H elf000.bsch.rot124.H elf000.bsch.rot125.H elf000.bsch.rot126.H elf000.bsch.rot127.H elf000.bsch.rot128.H elf000.bsch.rot129.H elf000.bsch.rot130.H elf000.bsch.rot131.H elf000.bsch.rot132.H > c.H
	Cat a.H b.H c.H axis=4 > d.H
	echo d4=30 o4=1600 d5=30 o5=0 >> d.H
	< d.H Transp plane=12 reshape=1,5 maxsize=2000 > $@.H

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

bschlk000_eps0a.rot.perm.%:
	echo o4=1600 d4=30 d5=30 >> elf.3Dbschrot00.H
	Window3d squeeze=n < elf.3Dbschrot00.H n1=281 min3=-300 n3=21 | Transp plane=12 reshape=1,5 \
	maxsize=1000 > cig.H
#	Window3d squeeze=n < elf.3Dbschrot00.H n1=281 min3=-300 n3=21 | Transp plane=12 reshape=1,5 \
#	maxsize=1000 > b.H;< b.H ${Tese}/bin/OW_3dscale.x > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf000.perm.rec.gx$*.scale.rot.H \
	sou_out=elf000.perm.sou.gx$*.scale.rot.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1


bschlk000.rms:
#	Gpow gpow=2 < elf.3Dbschrot00.H|Transp maxsize=2000|Stack|Window3d|Stack|Window3d|Gpow gpow=0.5|Scale scale_to=1 > c.H
	< elf.3Dbschrot00.H ${Tese}/bin/OW_3dCscale.x w=c.H > elf.3Dbschrot00.rms.H

bschlk000_eps0a.rot.perms.%:
	echo o4=1600 d4=30 d5=30 >> elf.3Dbschrot00.rms.H
	Window3d squeeze=n < elf.3Dbschrot00.rms.H n1=281 | Transp plane=12 reshape=1,5 maxsize=1000 > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf000.perm.rec.gx$*.scale.rot.H \
	sou_out=elf000.perm.sou.gx$*.scale.rot.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1

bschlk000.rna%:
	Window3d < elf.3Dbschrot00.rms.H n3=1 f3=$* | Transp plane=23 | Fx2d lambda=10 > j$*.H

bschlk000.rna.exe1:
	make bschlk000.rna000;make bschlk000.rna001;make bschlk000.rna002;make bschlk000.rna003;make bschlk000.rna004;\
	make bschlk000.rna005;make bschlk000.rna006;make bschlk000.rna007;make bschlk000.rna008;make bschlk000.rna009;\

bschlk000.rna.exe2:
	make bschlk000.rna010;make bschlk000.rna011;make bschlk000.rna012;make bschlk000.rna013;make bschlk000.rna014;\
	make bschlk000.rna015;make bschlk000.rna016;make bschlk000.rna017;make bschlk000.rna018;make bschlk000.rna019;\

bschlk000.rna.exe3:
	make bschlk000.rna020;make bschlk000.rna021;make bschlk000.rna022;make bschlk000.rna023;make bschlk000.rna024;\
	make bschlk000.rna025;make bschlk000.rna026;make bschlk000.rna027;make bschlk000.rna028;make bschlk000.rna029;\
	make bschlk000.rna030

bschlk000.rna.compose:
	Cat axis=4 j000.H j001.H j002.H j003.H j004.H j005.H j006.H j007.H j008.H j009.H j010.H j011.H j012.H j013.H j014.H j015.H j016.H j017.H j018.H j019.H j020.H j021.H j022.H j023.H j024.H j025.H j026.H j027.H j028.H j029.H j030.H > $@.H

bschlk000_eps0a.rot.permrna.%:
	Window3d squeeze=n < bschlk000.rna.compose.H n1=281 | Transp plane=32 reshape=1,2,4,5 maxsize=1000 > a.H
	< a.H Transp plane=12 reshape=1,5 maxsize=1000 > cig.H
	${Tese}/bin/OW3d_arealsht.x cig=cig.H sou=bsch.sou.H \
	vel=elf.L7f.vinit.rmo2.grid30.H rec_out=elf000.perm.rec.gx$*.scale.rot.rna.H \
	sou_out=elf000.perm.sou.gx$*.scale.rot.rna.H \
	ngapx=$* ngapy=1 par=bsch_permeps0.P node=14 shape=1


##############################################
elf001.bsch.rot%:
	Window3d n3=1 f3=$* n1=100 f1=220 < elf000.bsch.window.scale.H squeeze=n| Transp maxsize=2000 > a1$*.H
	< a1$*.H Transp plane=34 maxsize=2000 > a2$*.H
	< a2$*.H Transp plane=12 reshape=4,5 maxsize=2000 > a3$*.H
	< a3$*.H ${Tese}/bin/OW_3dslant_2d.x nodes=14 nax=81 nahx=81 nay=1 nahy=1 oay=0 oahy=0 \
	oax=-40 oahx=-40 dax=1 dahx=1 nxw=11 rot=1 rec=1 > $@.H




	


