MAKE_DEPEND=yes
SEPINC = /opt/SEP/SEP6.4.6_INTEL/include
include $(SEPINC)/SEP.top
#include ../MAKERULES

SRCDIR = ./
BINDIR = /net/koko/claudio/Prst_ExpRefl/Tese/bin/
B = BINDIR
UPSTEXFLAGS=invras=y
LATOPTS  = style=html

LIBDIR = -L/opt/SEP/SEP6.4.6/lib 
#LIBDIR = -L/usr/local/SEP/SEP6.4.2_INTEL/lib 

FFT= -L/opt/FFTW/lib -lfftw3f 

UF90DEFINES = -DSEP_OMP
UF90FLAGS = -O3 -openmp #-CB #-axW -openmp #-static-libcxa -Bstatic # -vec_report3 -quiet
UF90LIBS =-lsepgeef90 -lsep2df90 -lvplotf90 -lvplot -lsep -lsep3df90 -lsupersetf90 -lsuperset -lsepdipf90 -lsepfilterf90 -lsep3d -lsepauxf90 -lsepvectorf90 -openmp
#UF90LIBS = -L/opt/intel_fce_80/lib/  -openmp -lsupersetf90 -lsuperset -lsepauxf90 -lsepmathf90 -lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -lsepgeef90  -lsepdipf90 -lsepfilterf90  -openmp


INCDIR = -I./ -I/opt/intel_fce_80/include/

default:${BINDIR}/Dip.x 

%.F90:  %.f90s
	ln -s $*.f90s $*.F90
Dip1=$(BINDIR)/Dip.x
Dip2=Dip/Dip.x
#Dip=dip/Dip
Dip3=Dip

dip = $(Dip1)

grey = Grey color=j wantscalebar=y newclip=1 eclip= 1.5 bclip= -1.5 > /dev/null 
median = Median.x wind1=3 wind2=3 

owdip: Dip.x
	rm -f ${B}/Dip.x; make Dip.x; mv -f Dip.x ${B}/Dip.x





#marmdip: Makefile $(dip)
#	< ../../Multiple/Marmousi/marm0mig.H  $(dip) method= puck_f  eps=0.00005 eps2=.00005 verb=y dave=y>1.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(dip) method= puck_f  eps=0.00005 eps2=.00005 verb=y dave=n>2.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip2) method= puck_f eps=0.000001 eps2=.000001 >3.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip3) method= puck_f eps=0.000001 eps2=.000001 >4.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip) method=sergey_s eps=0.002 eps2=.0002 | $(median) |$(grey) title='dips: marmousi2' out=$@2.v
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip) method=sergey_s eps=0.003 eps2=.0003 | $(median) |$(grey) title='dips: marmousi3' out=$@3.v
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip) method=sergey_s eps=0.004 eps2=.0004 | $(median) |$(grey) title='dips: marmousi4' out=$@4.v
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip) method=sergey_s eps=0.005 eps2=.0005 | $(median) |$(grey) title='dips: marmousi5' out=$@5.v
#	cat *.v  | tube
	Cat3d 1.H 2.H 3.H 4.H  ../../Multiple/Marmousi/marm1mig.H | Grey color=j gainpanel=e  wantscalebar=y| Tube	

marmdip: Makefile $(dip)
#	Window3d min3=0 n3=551 < ../../Multiple/Marmousi/marm0mig3.H | Transp |  $(dip) method= puck_f  eps=0.00005 eps2=.0001 verb=y dave=y>1.H
#	Window3d min2= 1000     < mod.H          		              |  $(dip) method= puck_f  eps=0.00005 eps2=.0001 verb=y dave=y>2.H
#	Sdip rect1=2 rect2=2 < mod2.H >1.H
#	Window3d min2= 1000     < mod.H		               	      	      |  $(dip) method=puck_f   eps=0.0001 eps2=.0001 nwind=3 verb=y dave= y>2.H
#	Sdip rect1=2 rect2=2 < tref.H >3.H
	Window3d min2= 1000     < tref.H 			              |  $(dip) method= puck_f  eps=0.00001 eps2=.00001 verb=y dave= y>3.H
	Window3d min2= 1000     < tref.H 			              |  $(dip) method =puck_f  eps=0.00001 eps2=.00001 verb=y dave=y>4.H
	Window3d min2= 1000     < tref.H 			              |  $(dip) method=puck_f   eps=0.00001 eps2=.00001 verb=y dave=y>5.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(dip) method= puck_f  eps=0.00005 eps2=.00005 verb=y dave=n>2.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip2) method= puck_f eps=0.000001 eps2=.000001 >3.H
#	< ../../Multiple/Marmousi/marm1mig.H  $(Dip3) method= puck_f eps=0.000001 eps2=.000001 >4.H
#	cat *.v  | tube
#	Cat3d 1.H | Grey color=j gainpanel=e  wantscalebar=y| Tube	&
#	Cat3d 2.H | Grey color=j gainpanel=e  wantscalebar=y| Tube	&
	Median.x < 1.H wind1=0 wind2=0 > test0.H 
	Median.x < 2.H wind1=1 wind2=1 > test1.H 
	Cp test0.H Fact/dip1.H
	Cp test1.H Fact/dip2.H
	Cp mod.H  Fact/in.H
	#Cat3d 4.H | Grey color=j gainpanel=e  wantscalebar=y polarity=1| Tube	&
	#Cat3d 5.H | Grey color=j gainpanel=e  wantscalebar=y polarity=1| Tube	&
	Cat3d [3-5].H | Grey color=j gainpanel=e  wantscalebar=y polarity=1| Tube	&
	#Cat3d [1-2].H | Grey color=j gainpanel=e  wantscalebar=y polarity=1| Tube	&
	#Cat3d test[0-1].H | Grey color=j gainpanel=e  wantscalebar=y polarity=1| Tube	&

tvel.H: 
	Vel vc=2000 n1=501 d1=5 o1=0 n2=501 d2=5 o2=0 n3=1 d3=5 o3=0      \
		xa1=0  xb1=2500   z1=100    alfa1=0    vr1=1200    const1=1 \
                xa2=0  xb2=2500   z2=200    alfa2=15   vr2=2400    const2=1 \
                xa4=0  xb4=2500   z4=300    alfa4=30   vr4=3500    const4=1 \
                xa5=0  xb5=2500   z5=400   alfa5=45   vr5=4800    const5=1 \
                xa6=0  xb6=2500   z6=600   alfa6=60   vr6=5000    const6=1 \
                xa7=0  xb7=2500   z7=800   alfa7=80   vr7=5500    const7=1 >  vel3.H
	mv vel3.H $@
	reflectivity.x vel=tvel.H | Smooth rect1=2 rect2=2 | Scale scale_to=1 > tref.H



clean:
	rm -f *.x *.fix.* *.mod *.o .make.dependencies.LINUX

include $(SEPINC)/SEP.bottom