dip45.H:
	< /net/koko/claudio/Research/Kirmod1/4th-data.infill.H Window3d n2=255 max1=2.5 | Bandpass flo=2 fhi=36 | \
	Transf f_min=2 f_min1=4 f_max1=30 f_max=36 wei=y is2d=y maxsize=2000 > a1.H
	< a1.H Transp plane=12 reshape=2,4,5 > $@

dip30.H:
	< slowness.grad1.H Gpow gpow=-1|Vconvert n1out=1048 d1out=0.002 o1out=0 invelopt=1 outvelopt=1 intype=1 outtype=2 inpar=2 outpar=1 > vrms.grad1.H
	echo n3=9216 n4=1 >> /net/koko/claudio/Research/Kirmod1/4th-data.infill.H
	< /net/koko/claudio/Research/Kirmod1/4th-data.infill.H NMO velfile=vrms.grad1.H > a0.H
	<a0.H Interp o2out=-787.5 d2out=12.5 n2out=127 > a1.H
	< a1.H NMO velfile=vrms.grad1.H inv=1 > a2.H
	echo n3=96 n4=96 >> /net/koko/claudio/Research/Kirmod1/4th-data.infill.H
	echo n3=96 n4=96 >> a2.H
	< a2.H Pad beg1=0 end1=100 | Bandpass flo=2 fhi=36 | \
	Transf f_min=2 f_min1=4 f_max1=30 f_max=36 wei=y is2d=y maxsize=2000 > a1.H
	< a1.H Transp plane=12 reshape=2,4,5 > $@

dip45_CAM.vslow.H:
	Spike n1=151 d1=10 o1=0 n2=100 o2=0 d2=25 n3=100 o3=0 d3=25 nsp=1 k1=1 > j0.H
	Math file1=j0.H exp="file1*0+1" | Tpow tpow=1 > j1.H ; Math file1=j1.H exp="file1+1500" > j2.H
	< j2.H Transp plane=13 | Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 | Gpow gpow=-1 > slowness.grad1.H
	Gpow <  slowness.grad1.H gpow=-1 | Scale dscale=0.95 > 6velocity.H
	Lloyd_vel < 6velocity.H nref=1 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120 > 6slowref.H
	nohup Parallel command=${Prst_ExpRefl}/opteron/scr1/S_G_cam.x \
		files="dip45.H,6velocity.H,6slowref.H,$@" \
		tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		axis="5,4,4,5" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=86 device=eth0 \
		pverb=2 mach_mfile=mfile par=6sisfull45.par

dip30_CAM.vslow.H:
	Gpow < slowness.H gpow=-1 | Scale dscale=0.95 | Window3d n3=151 > 6velocity.H
	Lloyd_vel < 6velocity.H nref=1 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120 > 6slowref.H
	nohup Parallel command=${Prst_ExpRefl}/opteron/scr1/S_G_cam.x \
		files="dip30.H,6velocity.H,6slowref.H,$@" \
		tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		axis="5,4,4,5" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=79 device=eth0 \
		pverb=2 mach_mfile=mfile par=6sisfull.par

#dip45_CAM_vslow.rot.H:
#	echo o3=-3175 d3=25 >> dip45_CAM.vslow.H
#	< dip45_CAM.vslow.H Window3d min3=-800 n3=65 min2=1000 n2=1 max5=1500 squeeze=n | Transp plane=14 > a1.H
#	< a1.H Transp plane=12 reshape=2,4,5 > a2.H
#	< a2.H Transp plane=12 reshape=4,5 > a3.H
#	< a3.H ${Tese}/bin/OW_3dslant_2d.x nodes=14 nax=41 nahx=81 nay=1 nahy=1 oay=0 oahy=0 \
#	oax=0 oahx=-40 dax=1 dahx=1 nxw=41 rot=1 > j1.H;
#	echo d1=0.004 >> j1.H; < j1.H Bandpass flo=10 fhi=40 > $@; echo d1=10 d4=30 >> $@;

dip45_CAM_vslow.rot%:
	echo o3=-3175 d3=25 >> dip45_CAM.vslow.H
	< dip45_CAM.vslow.H Window3d min3=-250 n3=21 f2=$* n2=1 max5=1500 squeeze=n | Transp plane=14 > a1$*.H
	< a1$*.H Transp plane=12 reshape=2,4,5 > a2$*.H
	< a2$*.H Transp plane=12 reshape=4,5 > a3$*.H
	< a3$*.H Transp plane=23 > a4$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a4$*.H Cig_tilt=a4.rot$*.H par=Cig45-rot.P positive_rot=1

dip45.rot1:
	make dip45_CAM_vslow.rot00;	make dip45_CAM_vslow.rot01;	make dip45_CAM_vslow.rot02;	make dip45_CAM_vslow.rot03;
	make dip45_CAM_vslow.rot04;	make dip45_CAM_vslow.rot05;	make dip45_CAM_vslow.rot06;	make dip45_CAM_vslow.rot07;
	make dip45_CAM_vslow.rot08;	make dip45_CAM_vslow.rot09;	make dip45_CAM_vslow.rot10;	make dip45_CAM_vslow.rot11;
	make dip45_CAM_vslow.rot12;	make dip45_CAM_vslow.rot13;	make dip45_CAM_vslow.rot14;	make dip45_CAM_vslow.rot15;
dip45.rot2:
	make dip45_CAM_vslow.rot16;	make dip45_CAM_vslow.rot17;	make dip45_CAM_vslow.rot18;	make dip45_CAM_vslow.rot19;
	make dip45_CAM_vslow.rot20;	make dip45_CAM_vslow.rot21;	make dip45_CAM_vslow.rot22;	make dip45_CAM_vslow.rot23;
	make dip45_CAM_vslow.rot24;	make dip45_CAM_vslow.rot25;	make dip45_CAM_vslow.rot26;	make dip45_CAM_vslow.rot27;
	make dip45_CAM_vslow.rot28;	make dip45_CAM_vslow.rot29;	make dip45_CAM_vslow.rot30;	make dip45_CAM_vslow.rot31;
dip45.rot3:
	make dip45_CAM_vslow.rot32;	make dip45_CAM_vslow.rot33;	make dip45_CAM_vslow.rot34;	make dip45_CAM_vslow.rot35;
	make dip45_CAM_vslow.rot36;	make dip45_CAM_vslow.rot37;	make dip45_CAM_vslow.rot38;	make dip45_CAM_vslow.rot39;
	make dip45_CAM_vslow.rot40;	make dip45_CAM_vslow.rot41;	make dip45_CAM_vslow.rot42;	make dip45_CAM_vslow.rot43;
	make dip45_CAM_vslow.rot44;	make dip45_CAM_vslow.rot45;	make dip45_CAM_vslow.rot46;	make dip45_CAM_vslow.rot47;
dip45.rot4:
	make dip45_CAM_vslow.rot48;	make dip45_CAM_vslow.rot49;	make dip45_CAM_vslow.rot50;	make dip45_CAM_vslow.rot51;
	make dip45_CAM_vslow.rot52;	make dip45_CAM_vslow.rot53;	make dip45_CAM_vslow.rot54;	make dip45_CAM_vslow.rot55;
	make dip45_CAM_vslow.rot56;	make dip45_CAM_vslow.rot57;	make dip45_CAM_vslow.rot58;	make dip45_CAM_vslow.rot59;
	make dip45_CAM_vslow.rot60;	make dip45_CAM_vslow.rot61;	make dip45_CAM_vslow.rot62;	make dip45_CAM_vslow.rot63;
dip45.rot5:
	make dip45_CAM_vslow.rot64;	make dip45_CAM_vslow.rot65;	make dip45_CAM_vslow.rot66;	make dip45_CAM_vslow.rot67;
	make dip45_CAM_vslow.rot68;	make dip45_CAM_vslow.rot69;	make dip45_CAM_vslow.rot70;	make dip45_CAM_vslow.rot71;
	make dip45_CAM_vslow.rot72;	make dip45_CAM_vslow.rot73;	make dip45_CAM_vslow.rot74;	make dip45_CAM_vslow.rot75;
	make dip45_CAM_vslow.rot76;	make dip45_CAM_vslow.rot77;	make dip45_CAM_vslow.rot78;	make dip45_CAM_vslow.rot79;
dip45.rot6:
	make dip45_CAM_vslow.rot80;	make dip45_CAM_vslow.rot81;	make dip45_CAM_vslow.rot82;	make dip45_CAM_vslow.rot83;
	make dip45_CAM_vslow.rot84;	make dip45_CAM_vslow.rot85;	make dip45_CAM_vslow.rot86;	make dip45_CAM_vslow.rot87;
	make dip45_CAM_vslow.rot88;	make dip45_CAM_vslow.rot89;	make dip45_CAM_vslow.rot90;	make dip45_CAM_vslow.rot91;
	make dip45_CAM_vslow.rot92;	make dip45_CAM_vslow.rot93;	make dip45_CAM_vslow.rot94;	make dip45_CAM_vslow.rot95;

dip45.CAM_vslow.rot.H:
	Cat axis=5 a4.rot00.H a4.rot01.H a4.rot02.H a4.rot03.H a4.rot04.H a4.rot05.H a4.rot06.H a4.rot07.H a4.rot08.H a4.rot09.H a4.rot10.H a4.rot11.H a4.rot12.H a4.rot13.H a4.rot14.H a4.rot15.H a4.rot16.H a4.rot17.H a4.rot18.H a4.rot19.H a4.rot20.H a4.rot21.H a4.rot22.H a4.rot23.H a4.rot24.H a4.rot25.H a4.rot26.H a4.rot27.H a4.rot28.H a4.rot29.H a4.rot30.H a4.rot31.H a4.rot32.H a4.rot33.H a4.rot34.H a4.rot35.H a4.rot36.H a4.rot37.H a4.rot38.H a4.rot39.H a4.rot40.H a4.rot41.H a4.rot42.H a4.rot43.H a4.rot44.H a4.rot45.H a4.rot46.H a4.rot47.H a4.rot48.H a4.rot49.H a4.rot50.H a4.rot51.H a4.rot52.H a4.rot53.H a4.rot54.H a4.rot55.H a4.rot56.H a4.rot57.H a4.rot58.H a4.rot59.H a4.rot60.H a4.rot61.H a4.rot62.H a4.rot63.H a4.rot64.H a4.rot65.H a4.rot66.H a4.rot67.H a4.rot68.H a4.rot69.H a4.rot70.H a4.rot71.H a4.rot72.H a4.rot73.H a4.rot74.H a4.rot75.H a4.rot76.H a4.rot77.H a4.rot78.H a4.rot79.H a4.rot80.H a4.rot81.H a4.rot82.H a4.rot83.H a4.rot84.H a4.rot85.H a4.rot86.H a4.rot87.H a4.rot88.H a4.rot89.H a4.rot90.H a4.rot91.H a4.rot92.H a4.rot93.H a4.rot94.H a4.rot95.H > a10.H
	< a10.H Window3d | Transp plane=12 reshape=1,5 > $@
	
dip45.vslow.gx%.gy1.sort.H:
	echo d3=25 >> dip45.CAM_vslow.rot.H
	echo d4=25 >> dip45.CAM_vslow.rot.H
	< dip45.CAM_vslow.rot.H ${Tese}/bin/OW3d_sort.x ngapx=$* ngapy=1 > $@

dip45.gx%gy1:
	Gpow < slowness.H gpow=-1|Scale dscale=0.99|Gpow gpow=-1|Pad beg1=100 end1=100 beg2=100 end2=100 extend=1 > 6slowness.H
	Wavelet n1=1251 d1=0.002 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=1.25 domain=time > j0.H;echo o1=-1.25 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=30 f_max=36 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 > sou.H
	nohup Parallel command=${Tese}/bin/OW3d_arealsht.x \
	files="dip45.vslow.gx$*.gy1.sort.H,sou.H,6slowness.H,dip45.rec.gx$*gy1.H,dip45.sou.gx$*gy1.H" \
	tags="cig,sou,vel,rec_out,sou_out" usage="INPUT,INPUT,INPUT,OUTPUT,OUTPUT" \
	axis="7,4,4,4,4" file_type="DISTRIBUTE,COPY,COPY,DISTRIBUTE,DISTRIBUTE" \
	nblock=$* ngapx=$* ngapy=1 device=eth0 par=dip45_perm.P pverb=2 mach_mfile=mfile

dip45.48x1y: 
	echo d3=25 >> dip45.CAM_vslow.rot.H
	Wavelet n1=751 d1=0.002 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=0.75 domain=time > j0.H;echo o1=-0.75 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=42 f_max=48 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 > sou.H
	${Tese}/bin/OW3d_arealsht.x cig=dip45.CAM_vslow.rot.H sou=sou.H \
	vel=slowness.grad1.H rec_out=dip45.rec.gx48gy1f25.H dip45.rec.gx48gy1.H sou_out=dip45.sou.gx48gy1f25.H dip45.sou.gx48gy1.H \
	ngapx=48 ngapy=1 par=dip45_perm.P node=14 shape=1

migtest1:
	echo o6=0 >> 4th-plane/dip45.sou.gx48gy1.H
	echo o6=0 >> 4th-plane/dip45.rec.gx48gy1.H
	nohup /opt/MPI/mpich-1.2.7p1/bin/mpirun -v -np 17 -machinefile mfile ./Bin/steepest_ls3D.x \
	sou=4th-plane/dip45.sou.gx48gy1.H rec=4th-plane/dip45.rec.gx48gy1.H \
	vel=4th-plane/6slowness.H grad=4grad.H Dimage=4dimg.H bkgimg=4bimg.H Velout=4vout.H grad_=4grad_.H \
	zmin=0 zmax=1500 nhx=21 ohx=-250 jhx=1 nhy=1 ohy=0 jhy=1 niter=1 dso=1 inversion=1 areal=1 \
	xv_min=200 xv_max=1900 yv_min=200 yv_max=1900 \
	dxk= 1200 nxk= 6 oxk= -1200 dzk= 50 nzk= 39 ozk= 700 dxkmin= 1200 horv0= ${ELF}/hor3dgt2500a.pik \
	dyk= 1200 nyk= 4 oyk= 1090 dykmin= 1200 \
	extrapolator=pspi nrefvel=1 lloyd=1 node=4 bounds=1 cbounds=0 p_lo=60 p_up=60 nwblk=4 nzblk=10 \
	padtraces_x=0 traces_in_taper_x=30 padtraces_y=0 traces_in_taper_y=30

migdip45.48x1y: 
	echo o6=0 d2=25 o2=0  >> dip45.sou.gx48gy1f25.H
	echo o6=0 d2=25 o2=0  >> dip45.rec.gx48gy1f25.H
	${Tese}/src/areal/3dwemva/OWmig3d_areal.x sou=dip45.sou.gx48gy1f25.H rec=dip45.rec.gx48gy1f25.H \
	vel=slowness.grad1.H bkgimg=4bimgf25.H \
	zmin=0 zmax=1500 nhx=21 ohx=-250 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=pspi nrefvel=1 lloyd=0 vsal=7000 nwblk= 1 > /dev/null

dip30_CAM_vslow.interp.H:
	echo o3=-787.5 d3=12.5 >> dip30_CAM.vslow.H
	< dip30_CAM.vslow.H Window3d min3=-250 n3=41 max5=1500 squeeze=n|Interp d3out=25 maxsize=2000 > $@

dip30_CAM_vslow.rot%:
	< dip30_CAM_vslow.interp.H Window3d f2=$* n2=1 squeeze=n | Transp plane=14 > a1$*.H
	< a1$*.H Transp plane=12 reshape=2,4,5 > a2$*.H
	< a2$*.H Transp plane=12 reshape=4,5 > a3$*.H
	< a3$*.H Transp plane=23 > a4$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a4$*.H Cig_tilt=a4.rot$*.H par=Cig45-rot.P positive_rot=1

dip30.rot1:
	make dip30_CAM_vslow.rot00;	make dip30_CAM_vslow.rot01;	make dip30_CAM_vslow.rot02;	make dip30_CAM_vslow.rot03;
	make dip30_CAM_vslow.rot04;	make dip30_CAM_vslow.rot05;	make dip30_CAM_vslow.rot06;	make dip30_CAM_vslow.rot07;
	make dip30_CAM_vslow.rot08;	make dip30_CAM_vslow.rot09;	make dip30_CAM_vslow.rot10;	make dip30_CAM_vslow.rot11;
	make dip30_CAM_vslow.rot12;	make dip30_CAM_vslow.rot13;	make dip30_CAM_vslow.rot14;	make dip30_CAM_vslow.rot15;
dip30.rot2:
	make dip30_CAM_vslow.rot16;	make dip30_CAM_vslow.rot17;	make dip30_CAM_vslow.rot18;	make dip30_CAM_vslow.rot19;
	make dip30_CAM_vslow.rot20;	make dip30_CAM_vslow.rot21;	make dip30_CAM_vslow.rot22;	make dip30_CAM_vslow.rot23;
	make dip30_CAM_vslow.rot24;	make dip30_CAM_vslow.rot25;	make dip30_CAM_vslow.rot26;	make dip30_CAM_vslow.rot27;
	make dip30_CAM_vslow.rot28;	make dip30_CAM_vslow.rot29;	make dip30_CAM_vslow.rot30;	make dip30_CAM_vslow.rot31;
dip30.rot3:
	make dip30_CAM_vslow.rot32;	make dip30_CAM_vslow.rot33;	make dip30_CAM_vslow.rot34;	make dip30_CAM_vslow.rot35;
	make dip30_CAM_vslow.rot36;	make dip30_CAM_vslow.rot37;	make dip30_CAM_vslow.rot38;	make dip30_CAM_vslow.rot39;
	make dip30_CAM_vslow.rot40;	make dip30_CAM_vslow.rot41;	make dip30_CAM_vslow.rot42;	make dip30_CAM_vslow.rot43;
	make dip30_CAM_vslow.rot44;	make dip30_CAM_vslow.rot45;	make dip30_CAM_vslow.rot46;	make dip30_CAM_vslow.rot47;
dip30.rot4:
	make dip30_CAM_vslow.rot48;	make dip30_CAM_vslow.rot49;	make dip30_CAM_vslow.rot50;	make dip30_CAM_vslow.rot51;
	make dip30_CAM_vslow.rot52;	make dip30_CAM_vslow.rot53;	make dip30_CAM_vslow.rot54;	make dip30_CAM_vslow.rot55;
	make dip30_CAM_vslow.rot56;	make dip30_CAM_vslow.rot57;	make dip30_CAM_vslow.rot58;	make dip30_CAM_vslow.rot59;
	make dip30_CAM_vslow.rot60;	make dip30_CAM_vslow.rot61;	make dip30_CAM_vslow.rot62;	make dip30_CAM_vslow.rot63;
dip30.rot5:
	make dip30_CAM_vslow.rot64;	make dip30_CAM_vslow.rot65;	make dip30_CAM_vslow.rot66;	make dip30_CAM_vslow.rot67;
	make dip30_CAM_vslow.rot68;	make dip30_CAM_vslow.rot69;	make dip30_CAM_vslow.rot70;	make dip30_CAM_vslow.rot71;
	make dip30_CAM_vslow.rot72;	make dip30_CAM_vslow.rot73;	make dip30_CAM_vslow.rot74;	make dip30_CAM_vslow.rot75;
	make dip30_CAM_vslow.rot76;	make dip30_CAM_vslow.rot77;	make dip30_CAM_vslow.rot78;	make dip30_CAM_vslow.rot79;
dip30.rot6:
	make dip30_CAM_vslow.rot80;	make dip30_CAM_vslow.rot81;	make dip30_CAM_vslow.rot82;	make dip30_CAM_vslow.rot83;
	make dip30_CAM_vslow.rot84;	make dip30_CAM_vslow.rot85;	make dip30_CAM_vslow.rot86;	make dip30_CAM_vslow.rot87;
	make dip30_CAM_vslow.rot88;	make dip30_CAM_vslow.rot89;	make dip30_CAM_vslow.rot90;	make dip30_CAM_vslow.rot91;
	make dip30_CAM_vslow.rot92;	make dip30_CAM_vslow.rot93;	make dip30_CAM_vslow.rot94;	make dip30_CAM_vslow.rot95;

dip30.CAM_vslow.rot.H:
	Cat axis=5 a4.rot00.H a4.rot01.H a4.rot02.H a4.rot03.H a4.rot04.H a4.rot05.H a4.rot06.H a4.rot07.H a4.rot08.H a4.rot09.H a4.rot10.H a4.rot11.H a4.rot12.H a4.rot13.H a4.rot14.H a4.rot15.H a4.rot16.H a4.rot17.H a4.rot18.H a4.rot19.H a4.rot20.H a4.rot21.H a4.rot22.H a4.rot23.H a4.rot24.H a4.rot25.H a4.rot26.H a4.rot27.H a4.rot28.H a4.rot29.H a4.rot30.H a4.rot31.H a4.rot32.H a4.rot33.H a4.rot34.H a4.rot35.H a4.rot36.H a4.rot37.H a4.rot38.H a4.rot39.H a4.rot40.H a4.rot41.H a4.rot42.H a4.rot43.H a4.rot44.H a4.rot45.H a4.rot46.H a4.rot47.H a4.rot48.H a4.rot49.H a4.rot50.H a4.rot51.H a4.rot52.H a4.rot53.H a4.rot54.H a4.rot55.H a4.rot56.H a4.rot57.H a4.rot58.H a4.rot59.H a4.rot60.H a4.rot61.H a4.rot62.H a4.rot63.H a4.rot64.H a4.rot65.H a4.rot66.H a4.rot67.H a4.rot68.H a4.rot69.H a4.rot70.H a4.rot71.H a4.rot72.H a4.rot73.H a4.rot74.H a4.rot75.H a4.rot76.H a4.rot77.H a4.rot78.H a4.rot79.H a4.rot80.H a4.rot81.H a4.rot82.H a4.rot83.H a4.rot84.H a4.rot85.H a4.rot86.H a4.rot87.H a4.rot88.H a4.rot89.H a4.rot90.H a4.rot91.H a4.rot92.H a4.rot93.H a4.rot94.H a4.rot95.H > a10.H
	< a10.H Window3d | Transp plane=12 reshape=1,5 | Pad end1=1 beg1=0  > $@


dip30.CAM_vslow.rot.sltstk.H:
#	< dip30_CAM.vslow.H Window3d min3=-250 n3=21 max5=1500 squeeze=n | Transp plane=14 > a.H
#	< a.H Transp plane=12 reshape=2,4,5 > x.H; < x.H Transp plane=12 reshape=4,5 > a.H
	< a.H Window3d squeeze=n n5=1 min5=800 | ${Tese}/bin/OW_3dslant_2d.x nodes=14 nax=41 oax=0 dax=1 nay=1 oay=0 day=1 \
	nahx=101 oahx=-50 dahx=1 nahy=1 oahy=0 dahy=1 nxw=31 rot=1 > j1.H;
#	< a.H Window3d squeeze=n n5=1 min5=800 | ${Tese}/bin/OW_3dslant.x nax=31 nah=41 oax=-30 oah=-40 dax=2 dah=2 nxw=25 rot=1 upshft=0 nodes=14 nxchunk=1 > j1.H

#	< dip.bslo.H Window3d min1=-1000 max1=1010 | Transp plane=13 | ${Prst_ExpRefl}/Tese/bin/OW_3dslant.x nax=61 nah=81 oax=-30 oah=-40 dax=1 dah=1 nxw=41 rot=1 upshft=0 nodes=14 nxchunk=1 > j2.sou.H


dip30.48x1y: 
	echo d3=25 >> dip30.CAM_vslow.rot.H
	Gpow < slowness.H gpow=-1|Scale dscale=0.95|Gpow gpow=-1|Pad beg1=100 end1=100 beg2=100 end2=100 extend=1 > 6slowness.H
	Wavelet n1=751 d1=0.002 wavelet=ricker2 phase=0 fund=20 order=8 tdelay=0.75 domain=time > j0.H;echo o1=-0.75 >> 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 > 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=dip30.CAM_vslow.rot.H sou=sou.H \
	vel=6slowness.H rec_out=dip30.rec.gx48gy1f20.H sou_out=dip30.sou.gx48gy1f20.H \
	ngapx=48 ngapy=1 par=dip30_perm.P node=14 shape=1

migdip30.48x1y: 
	echo o6=0 d2=25 o2=0 >> dip30.sou.gx48gy1f20.H
	echo o6=0 d2=25 o2=0 >> dip30.rec.gx48gy1f20.H
	Gpow < slowness.H gpow=-1|Scale dscale=0.95|Gpow gpow=-1|Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 > 6slowness.H
	${Tese}/src/areal/3dwemva/OWmig3d_areal.x sou=dip30.sou.gx48gy1f20.H rec=dip30.rec.gx48gy1f20.H \
	vel=6slowness.H bkgimg=dip30.bimgf25.H \
	zmin=0 zmax=1500 nhx=21 ohx=-250 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=pspi nrefvel=1 lloyd=0 vsal=7000 nwblk= 1 > /dev/null

dip30.48x1y25: 
	echo d3=25 d4=25 >> dip30.CAM_vslow.rot.H
	Gpow < slowness.H gpow=-1|Scale dscale=0.95|Gpow gpow=-1|Pad beg1=100 end1=100 beg2=100 end2=100 extend=1 > 6slowness.H
	Wavelet n1=751 d1=0.002 wavelet=ricker2 phase=0 fund=25 order=8 tdelay=0.75 domain=time > j0.H;echo o1=-0.75 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=50 f_max=56 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 > sou.H
	${Tese}/bin/OW3d_arealsht.x cig=dip30.CAM_vslow.rot.H sou=sou.H \
	vel=6slowness.H rec_out=dip30.rec.gx48gy1f25.H sou_out=dip30.sou.gx48gy1f25.H \
	ngapx=48 ngapy=1 par=dip30_perm.P node=14 shape=1

migdip30.48x1y25: 
	echo o6=0 d2=25 o2=0 >> dip30.sou.gx48gy1f25.H
	echo o6=0 d2=25 o2=0 >> dip30.rec.gx48gy1f25.H
	Gpow < slowness.H gpow=-1|Scale dscale=0.95|Gpow gpow=-1|Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 > 6slowness.H
	${Tese}/src/areal/3dwemva/OWmig3d_areal.x sou=dip30.sou.gx48gy1f25.H rec=dip30.rec.gx48gy1f25.H \
	vel=6slowness.H bkgimg=dip30.bimgf25.H \
	zmin=0 zmax=1500 nhx=21 ohx=-250 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=pspi nrefvel=1 lloyd=0 vsal=7000 nwblk= 1 > /dev/null

######################
#
# Stretch by cos ay'
#
dip30_CAM.vslow.stretch.H:
	Window3d < dip30_CAM_vslow.interp.H min3=-250 n3=21 squeeze=n | Transp plane=23 maxsize=2000 > 1.H
	echo o2=-288.67 d2=28.867 >> 1.H
	Interp < 1.H o2out=-300 d2out=25 n2out=25 maxsize=2000 | Transp plane=23 maxsize=2000 > $@

dip30_CAM_vslow.strch.rot%:
	< dip30_CAM.vslow.stretch.H Window3d f2=$* n2=1 squeeze=n | Transp plane=14 > a1$*.H
	< a1$*.H Transp plane=12 reshape=2,4,5 > a2$*.H
	< a2$*.H Transp plane=12 reshape=4,5 > a3$*.H
	< a3$*.H Transp plane=23 > a4$*.H
	~claudio/research.ln/Others/Biondo/rotation/v3/RotCig_v3.x Cig_hx=a4$*.H Cig_tilt=a4.rot$*.H par=Cig30-rot.str.P positive_rot=1

dip30.rot1s:
	make dip30_CAM_vslow.strch.rot00;	make dip30_CAM_vslow.strch.rot01;	make dip30_CAM_vslow.strch.rot02;	make dip30_CAM_vslow.strch.rot03;
	make dip30_CAM_vslow.strch.rot04;	make dip30_CAM_vslow.strch.rot05;	make dip30_CAM_vslow.strch.rot06;	make dip30_CAM_vslow.strch.rot07;
	make dip30_CAM_vslow.strch.rot08;	make dip30_CAM_vslow.strch.rot09;	make dip30_CAM_vslow.strch.rot10;	make dip30_CAM_vslow.strch.rot11;
	make dip30_CAM_vslow.strch.rot12;	make dip30_CAM_vslow.strch.rot13;	make dip30_CAM_vslow.strch.rot14;	make dip30_CAM_vslow.strch.rot15;
dip30.rot2s:
	make dip30_CAM_vslow.strch.rot16;	make dip30_CAM_vslow.strch.rot17;	make dip30_CAM_vslow.strch.rot18;	make dip30_CAM_vslow.strch.rot19;
	make dip30_CAM_vslow.strch.rot20;	make dip30_CAM_vslow.strch.rot21;	make dip30_CAM_vslow.strch.rot22;	make dip30_CAM_vslow.strch.rot23;
	make dip30_CAM_vslow.strch.rot24;	make dip30_CAM_vslow.strch.rot25;	make dip30_CAM_vslow.strch.rot26;	make dip30_CAM_vslow.strch.rot27;
	make dip30_CAM_vslow.strch.rot28;	make dip30_CAM_vslow.strch.rot29;	make dip30_CAM_vslow.strch.rot30;	make dip30_CAM_vslow.strch.rot31;
dip30.rot3s:
	make dip30_CAM_vslow.strch.rot32;	make dip30_CAM_vslow.strch.rot33;	make dip30_CAM_vslow.strch.rot34;	make dip30_CAM_vslow.strch.rot35;
	make dip30_CAM_vslow.strch.rot36;	make dip30_CAM_vslow.strch.rot37;	make dip30_CAM_vslow.strch.rot38;	make dip30_CAM_vslow.strch.rot39;
	make dip30_CAM_vslow.strch.rot40;	make dip30_CAM_vslow.strch.rot41;	make dip30_CAM_vslow.strch.rot42;	make dip30_CAM_vslow.strch.rot43;
	make dip30_CAM_vslow.strch.rot44;	make dip30_CAM_vslow.strch.rot45;	make dip30_CAM_vslow.strch.rot46;	make dip30_CAM_vslow.strch.rot47;
dip30.rot4s:
	make dip30_CAM_vslow.strch.rot48;	make dip30_CAM_vslow.strch.rot49;	make dip30_CAM_vslow.strch.rot50;	make dip30_CAM_vslow.strch.rot51;
	make dip30_CAM_vslow.strch.rot52;	make dip30_CAM_vslow.strch.rot53;	make dip30_CAM_vslow.strch.rot54;	make dip30_CAM_vslow.strch.rot55;
	make dip30_CAM_vslow.strch.rot56;	make dip30_CAM_vslow.strch.rot57;	make dip30_CAM_vslow.strch.rot58;	make dip30_CAM_vslow.strch.rot59;
	make dip30_CAM_vslow.strch.rot60;	make dip30_CAM_vslow.strch.rot61;	make dip30_CAM_vslow.strch.rot62;	make dip30_CAM_vslow.strch.rot63;
dip30.rot5s:
	make dip30_CAM_vslow.strch.rot64;	make dip30_CAM_vslow.strch.rot65;	make dip30_CAM_vslow.strch.rot66;	make dip30_CAM_vslow.strch.rot67;
	make dip30_CAM_vslow.strch.rot68;	make dip30_CAM_vslow.strch.rot69;	make dip30_CAM_vslow.strch.rot70;	make dip30_CAM_vslow.strch.rot71;
	make dip30_CAM_vslow.strch.rot72;	make dip30_CAM_vslow.strch.rot73;	make dip30_CAM_vslow.strch.rot74;	make dip30_CAM_vslow.strch.rot75;
	make dip30_CAM_vslow.strch.rot76;	make dip30_CAM_vslow.strch.rot77;	make dip30_CAM_vslow.strch.rot78;	make dip30_CAM_vslow.strch.rot79;
dip30.rot6s:
	make dip30_CAM_vslow.strch.rot80;	make dip30_CAM_vslow.strch.rot81;	make dip30_CAM_vslow.strch.rot82;	make dip30_CAM_vslow.strch.rot83;
	make dip30_CAM_vslow.strch.rot84;	make dip30_CAM_vslow.strch.rot85;	make dip30_CAM_vslow.strch.rot86;	make dip30_CAM_vslow.strch.rot87;
	make dip30_CAM_vslow.strch.rot88;	make dip30_CAM_vslow.strch.rot89;	make dip30_CAM_vslow.strch.rot90;	make dip30_CAM_vslow.strch.rot91;
	make dip30_CAM_vslow.strch.rot92;	make dip30_CAM_vslow.strch.rot93;	make dip30_CAM_vslow.strch.rot94;	make dip30_CAM_vslow.strch.rot95;

dip30.CAM_vslow.rot.stretch.H:
	Cat axis=5 a4.rot00.H a4.rot01.H a4.rot02.H a4.rot03.H a4.rot04.H a4.rot05.H a4.rot06.H a4.rot07.H a4.rot08.H a4.rot09.H a4.rot10.H a4.rot11.H a4.rot12.H a4.rot13.H a4.rot14.H a4.rot15.H a4.rot16.H a4.rot17.H a4.rot18.H a4.rot19.H a4.rot20.H a4.rot21.H a4.rot22.H a4.rot23.H a4.rot24.H a4.rot25.H a4.rot26.H a4.rot27.H a4.rot28.H a4.rot29.H a4.rot30.H a4.rot31.H a4.rot32.H a4.rot33.H a4.rot34.H a4.rot35.H a4.rot36.H a4.rot37.H a4.rot38.H a4.rot39.H a4.rot40.H a4.rot41.H a4.rot42.H a4.rot43.H a4.rot44.H a4.rot45.H a4.rot46.H a4.rot47.H a4.rot48.H a4.rot49.H a4.rot50.H a4.rot51.H a4.rot52.H a4.rot53.H a4.rot54.H a4.rot55.H a4.rot56.H a4.rot57.H a4.rot58.H a4.rot59.H a4.rot60.H a4.rot61.H a4.rot62.H a4.rot63.H a4.rot64.H a4.rot65.H a4.rot66.H a4.rot67.H a4.rot68.H a4.rot69.H a4.rot70.H a4.rot71.H a4.rot72.H a4.rot73.H a4.rot74.H a4.rot75.H a4.rot76.H a4.rot77.H a4.rot78.H a4.rot79.H a4.rot80.H a4.rot81.H a4.rot82.H a4.rot83.H a4.rot84.H a4.rot85.H a4.rot86.H a4.rot87.H a4.rot88.H a4.rot89.H a4.rot90.H a4.rot91.H a4.rot92.H a4.rot93.H a4.rot94.H a4.rot95.H > a10.H
	< a10.H Window3d | Transp plane=12 reshape=1,5 | Pad end1=1 beg1=0 > $@

dip30.48x1y.strch:
	echo d4=25 d3=25 n6=1 >> dip30.CAM_vslow.rot.stretch.H
	Wavelet n1=751 d1=0.002 wavelet=ricker2 phase=0 fund=25 order=8 tdelay=0.75 domain=time > j0.H;echo o1=-0.75 >> j0.H
	<j0.H Scale dclip=-1 | Transf f_min=0 f_min1=4 f_max1=50 f_max=56 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 > sou.H
	${Tese}/bin/OW3d_arealsht.x cig=dip30.CAM_vslow.rot.stretch.H sou=sou.H \
	vel=6slowness.H rec_out=dip30.rec.gx48gy1f25.strch.H sou_out=dip30.sou.gx48gy1f25.strch.H \
	ngapx=48 ngapy=1 par=dip30_perm.P node=14 shape=1

migdip30.48x1y25.strch: 
	echo o6=0 d2=25 >> dip30.sou.gx48gy1f25.strch.H
	echo o6=0 d2=25 >> dip30.rec.gx48gy1f25.strch.H
	${Tese}/src/areal/3dwemva/OWmig3d_areal.x sou=dip30.sou.gx48gy1f25.strch.H rec=dip30.rec.gx48gy1f25.strch.H \
	vel=6slowness.H bkgimg=dip30.bimgf25.strch.H \
	zmin=0 zmax=1500 nhx=21 ohx=-250 jhx=1 nhy=1 ohy=0 jhy=1 areal=1 migration=1 \
	padtraces_x=25 traces_in_taper_x=25 padtraces_y=25 traces_in_taper_y=25 node=14 \
	extrapolator=pspi nrefvel=1 lloyd=0 vsal=7000 nwblk= 1 > /dev/null