include ${SEPINC}/SEP.top
include 3Dprocfile_old
include 3Dprocfile
include Recip50file

FFTWLIBS=-L/usr/local/FFTW/lib -lfftw3f
MKL= -L/opt/intel/mkl72/lib/em64t -lmkl -lvml

UF90LIBS=-lsupersetf90 -lsuperset -lsep2df90   ${FFTWLIBS} -lsepmathf90 -lsepmath -lsepauxf90 -lsepaux -openmp iipa #${MKL} -lsvml

#VECTOR= -L/usr/local/INTEL/mkl/8.0.1/lib/em64t -lvml 
#VECTOR=-ipa -openmp


UF90LIBS=-lsupersetf90 -lsuperset -lsep2df90   ${FFTWLIBS} -lsepmathf90 -lsepmath -lsepauxf90 -lsepaux ${VECTOR}  #-L/usr/local/ACML/pathscale64/lib -lacml  # -L. -lvml
UCINCLUDES=-I/usr/local/SU/include -O5
UF90FLAGS=-O3 -mx87-precision=32  -ipa -openmp
BINDIR=./

#elf.3D.mig_adj.H: 
#	cp mig_shell.sh mig_adj_shell.sh
#	echo "mig_off D=elf.inv_freq.recip.H R=$@ velocity=elf.L7f.sinit.grid20.cam.H nz=300 oz=0 dz=10 amx_o=0 amy_o=1600 amy_n=200 verb=6 program=./src/CAM.x pverb=2   device=ib0  mpich_hostfilecom="-hostfile" mpich_includelocalhost=cees-node-002" >>mig_adj_shell.sh
#	qsub -q 8GB.q mig_adj_shell.sh

aaa%:
	echo "$**3" | bc

elf.mig_adj.H:
#	Window3d squeeze=n < elf.inv_freq.recip.H n1=100 f1=100 n2=50 f2=50 n3=70 min3=0 n5=40 f5=100 max_memory=2000 > data.H 
	nohup Parallel command=./scr1/S_G_cam.x \
	               files="data.1f.H,elf.L7f.sinit.grid20.cam.H,elf.L7f.ref_vel.grid20.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="6,4,4,6" file_type="COPY,COPY,COPY,COPY" nblock=1 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull.par

elf.3D.mig_adj.%.H:
	cd ../$*
	rm -rf nohup.* stat
	runall "rm /scr1/claudio/*"
	Window3d squeeze=n < ../elf.inv_freq.recip.H n5=4 f5=$* > data.$*.H 
	nohup Parallel command=../scr1/S_G_cam.x \
	               files="data.$*.H,../elf.L7f.sinit.grid20.cam.H,../elf.L7f.ref_vel.grid20.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="6,4,4,6" file_type="COPY,COPY,COPY,COPY" nblock=4 device=eth0 \
		     pverb=2 mach_mfile=mfile par=../sisfull.par


elf.3D.mig_adj.rmo.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.rmo.grid20.H,elf.L7f.ref_vel.rmo.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.bob.H:
#	< ELF_vel.mig.bob.H Transp plane=12 reshape=1,3 | Interp d1out=20 d2out=2 type=1 maxsize=2000 >  vel.mig.bob20.H
#	< vel.mig.bob20.H Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120 > refvel.mig.bob20.H
#	Window3d squeeze=n < ../Elf3d/ELF_elf.inv_freq.H f5=1 > data.H
	Window3d n3=302 < vel.mig.bob20.H > vv.H
	Window3d n3=302 < refvel.mig.bob20.H > rr.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data.H,vv.H,rr.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.bob.agc.H:
	Window3d n3=302 < vel.mig.bob20.H > vv.H
	Window3d n3=302 < refvel.mig.bob20.H > rr.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data.agc.H,vv.H,rr.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

elf.3D.mig.velopt.agc.H:
	< elf.velopt_nl.gx45_20x10.H Window3d min1=-1000 max1=8000 min2=1300 max2=6000 | Gpow gpow=-1 > vv.H 
	< vv.H Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120 > rr.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data_agcbf.H,vv.H,rr.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="5,4,4,6" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=196 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull196.par

elf.3D.mig.velinit.agc.H:
	< elf.L7f.sinit.rmo2.grid20.H Window3d n3=281 min1=-1000 max1=8000 min2=1300 max2=6000 > vv1.H
	< vv1.H Lloyd_vel nref=4 vel_map=no min_region_pct=2 min_slow_dev=0.01 perc_start=no niter_lloyd=120 > rr1.H
	nohup Parallel command=scr1/S_G_cam.x \
	               files="data_agcbf.H,vv1.H,rr1.H,$@" \
	               tags="Data,S,ref_slow,R" usage="INPUT,INPUT,INPUT,OUTPUT" \
		     axis="5,4,4,6" file_type="DISTRIBUTE,COPY,COPY,COPY" nblock=196 device=eth0 \
		     pverb=2 mach_mfile=mfile par=sisfull196.par


#	< elf.L7f.vinit.grid20.cam.H Window3d min2=2420 n2=1 squeeze=n | Gpow gpow=-1 > elf.L7f.vinit.grid20.cam.2420.H;mv /tmp/elf.L7f.vinit.grid20.cam.2420.H@ /data/temp1/guerrac/.
#	echo sets next: in="/data/temp1/guerrac/elf.L7f.vinit.grid20.cam.2420.H@" >> elf.L7f.vinit.grid20.cam.2420.H
#	echo "mig_off D=elf.line2420.tr.H R=$@ velocity=elf.L7f.vinit.grid20.cam.2420.H nz=400 oz=0 dz=10 amx_o=0 amy_o=2420 amx_n=539 amy_n=1 verb=6 program=./src/CAM.x pverb=2   device=ib0  mpich_hostfilecom="-hostfile" mpich_includelocalhost=cees-node-002" >>mig_adj_shell.sh
#	qsub mig_adj_shell.sh
#sukdmig2d < z.su > mig2.su fzt=0 dzt=80 nzt=63 fxt=0 dxt=100 nxt=135 fs=0 ds=200 ns=67 ttfile=tfile.su fzo=0 dzo=10 nzo=400 doff=-50 noff=72)
#	echo o4=0 >> elf.inv_freq.25m.H
#	echo "mig_off D=elf.inv_freq.H  R=$@ velocity=elf.L7f.vinit.H nz=2 400 oz=0 dz=10 verb=6 program=/data/temp2/bob/wei/CAM.x pverb=2   device=ib0  mpich_hostfilecom="-hostfile" mpich_includelocalhost=cees-node-002" >>mig_adj_shell.sh
#	Window3d squeeze=n < elf.inv_freq.25m.H n5=1 f5=60 > data.H
#	Window3d squeeze=n < elf.inv_freq.25m.H j3=2 n5=4 > data.H
#	echo "mig_off D=data.H elf.inv_freq.25m.H  R=$@ velocity=elf.L7f.vinit.cam.H nz=2 400 oz=0 dz=12.5 amx_o=0 amy_o=1600 amx_n=539 amy_n=147 verb=6 program=./src/CAM.x pverb=2   device=ib0  mpich_hostfilecom="-hostfile" mpich_includelocalhost=cees-node-002" >>mig_adj_shell.sh
#include small.inc
#include fig.inc
include ${SEPINC}/SEP.bottom