MAKE_DEPEND=yes
include ${SEPINC}/SEP.top
include FigMakefile

UF90DEFINES=-DLINUX -DUSE_FFTW -DFFTW=1 -DSEP_OMP -DOMP # -DDBUG 
UF90INCLUDES=-I.
UF90FLAGS=-O3 -openmp -openmp_report2 -fast
FFTWLIBS=-L/opt/FFTW/lib -lfftw3f
MAINLIBS=-lsupersetf90 -lsuperset -lsepparf90 -lseppar -lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -lsepgeef90

VEC= -L/opt/intel/Compiler/11.1/059/mkl/lib/em64t -lmkl -lvml
UF90LIBS=-lsep2df90 -lsepparf90 -lseppar -lsupersetf90 -lsuperset ${FFTWLIBS} ${ACMLIB} ${MAINLIBS} -openmp -lpthread


BIN=./bin
PAR=./par
RESDIR=./Fig
R = ${RESDIR}
FIG=./Fig
DATA=./data
OBJ=./obj

RESULTER =
RESULTCR = intro01 intro02 intro03 intro04 intro05 intro06 intro07 intro08
RESULTNR =

ER_v = 
CR_v = intro01.v intro02.v intro03.v intro04.v intro05.v intro06.v intro07.v intro08.v

ER_pdf = 
CR_pdf = intro01.pdf intro02.pdf intro03.pdf intro04.pdf intro05.pdf intro06.pdf intro07.pdf intro08.pdf

NR_pdf = 

cr_burn: clean
er_burn: clean
er_build: ${ER_v}
cr_build: ${CR_v}
er_pdf: ${ER_pdf}
cr_pdf: ${CR_pdf}

rst_ExpRefl=/net/server/claudio/Research/Prst_ExpRefl
Tese=/net/server/claudio/Research/Prst_ExpRefl/Tese
PERM=/net/server/claudio/Research/Thesis/perm
##########################################
#
# 
#
marmplane:
	Wavelet n1=951 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1.9 domain=time > j2.H;echo o1=-1.9 >> j2.H
	<j2.H Transf f_min=2 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13|Window3d squeeze=n f3=1 > 11111.H
	< ${Prst_ExpRefl}/Tese/3000marmod_velsmoo.H Window3d squeeze=n > 2222.H;echo o3=2.103051 >> 2222.H
	${Prst_ExpRefl}/Tese/bin/OWplane_areal.x sou=11111.H rec=2222.H v0=1500 \
	rec_pl=3000marm.crec.plane.H sou_pl=3000marm.csou.plane.H ango=-40 angd=1. node=16 > /dev/null

marmmigplane00:
	< 3000marm.csou.plane.H Window3d squeeze=n min4=-0 n4=1 > s.H
	< 3000marm.crec.plane.H Window3d squeeze=n min4=-0 n4=1 > r.H
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=3000marm.bkim.plane000.H \
	zmin=0 zmax=3000 nh=1 oh=0 migration=1 areal=1 \
	padtraces=100 traces_in_taper=25 node=14 vsal=7000 nrefvel=4 > /dev/null

8100marm_areal_synthesis.H:
	Wavelet n1=951 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1.9 domain=time > j2.H;echo o1=-1.9 >> j2.H
	<j2.H Transf f_min=0 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13 |Transp plane=24 > 11111.H
	echo n4=1 d4=1 o4=0 >> 11111.H
	< ${Prst_ExpRefl}/Tese/3000marmod_velsmoo.H Window3d squeeze=n | Pad beg3=8 end3=0 > 2222.H;echo o3=0 >> 2222.H
	< 2222.H ${Prst_ExpRefl}/2Wpar/bin/OW_arealsht_s.x sou=11111.H vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.tr.H nrefvel=4 xmin=-10392 xmax=19008 ymin=0 \
	ymax=0 zmin=0 zmax=3000 nh=1 oh=0 dh=0 vsal=14700 traces_in_taper=0 depth_wav=2300 shot_in=2222.H \
	synth_op=8100marm_areal_synthesis.H synsh=8101marm_Synth_Areal_recip_da.H velout= velout.H nphi=1 dphi=1. ophi=-20 > $@

8101marm_areal_mig:
	< 8101marm_Synth_Areal_recip_da.H Transp plane=13 > receiver.H
	< 8100marm_areal_synthesis.H Window3d squeeze=n min3=-6600 n3=926 | Transp plane=13 > source.H
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=source.H rec=receiver.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=8101marm_areal_mig.H \
	zmin=0 zmax=3000 nh=1 61 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 nrefvel=4 > /dev/null	

00marmenc: # 
	Wavelet n1=951 d1=0.004 wavelet=ricker2 phase=0 fund=20 order=0.2 tdelay=1.9 domain=time > j2.H;echo o1=-1.9 >> j2.H
	<j2.H Transf f_min=0 f_min1=4 f_max1=42 f_max=48 wei=y is2d=y > j1.H; <j1.H Window3d | Cabs > j2.H
	Math file1=j2.H exp="file1*0">j1.H;Cmplx j2.H j1.H|Transp plane=13 |Transp plane=24 > 11111.H
	echo n4=1 d4=1 o4=0 >> 11111.H
	< ${Prst_ExpRefl}/Tese/3000marmod_velsmoo.H Window3d squeeze=n | Pad beg3=8 end3=0 > 2222.H;echo o3=0 >> 2222.H
	< 2222.H ${Prst_ExpRefl}/ShtPSPI/bin/OW_encode.x sou=11111.H vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.tr.H \
	encrec=encrec1.H encsou=encsou1.H xmin=-8000 xmax=20000 ngap=20 > /dev/null

00marmig:
	< encsou1.H Window3d n4=1 squeeze=n > s.H
	< encrec1.H Window3d n4=1 squeeze=n > r.H
	${Prst_ExpRefl}/Tese/bin/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=marm_random_mig.H \
	zmin=0 zmax=3000 nh=1 61 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 nrefvel=4 > /dev/null	


00marmigperm:
	< ../Thesis/ispew/3000marm.csou35.perm.sd.H Window3d n4=1 squeeze=n > s.H
	< ../Thesis/ispew/3000marm.crec35.perm.sd.H Window3d n4=1 squeeze=n > r.H
	${Prst_ExpRefl}/Tese/src/areal/wemvaop/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=marm_perm_mig.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 icwind=1 deltat=0.016 dt=0.008 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigpermf:
	${Prst_ExpRefl}/Tese/src/areal/wemvaop/OWmig2d_areal.x sou=../Thesis/ispew/3000marm.csou35.perm.sd.H rec=../Thesis/ispew/3000marm.crec35.perm.sd.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=marm_perm_migf.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 icwind=1 deltat=0.016 dt=0.008 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigispew:
	< ../Thesis/ispew/3000marm.csou20.ispew.sd.H Window3d n4=1 squeeze=n > s.H
	< ../Thesis/ispew/3000marm.crec20.ispew.sd.H Window3d n4=1 squeeze=n > r.H
	${Prst_ExpRefl}/Tese/src/areal/wemvaop/OWmig2d_areal.x sou=s.H rec=r.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=marm_ispew_mig.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	

00marmigispewf:
	${Prst_ExpRefl}/Tese/src/areal/wemvaop/OWmig2d_areal.x sou=../Thesis/ispew/3000marm.csou20.ispew.sd.H rec=../Thesis/ispew/3000marm.crec20.ispew.sd.H  \
	vel=${Prst_ExpRefl}/Tese/3000marm_vel_smooth1.H bkgimg=marm_ispew_migf.H \
	zmin=1500 zmax=3000 nh=1 oh=0 -300 jh= 1 migration=1 areal=1 \
	padtraces=50 traces_in_taper=25 node=14 vsal=7000 extrapolator=ffd  > /dev/null	


regclean_%:
	@-touch $*; Rm $*; touch $*; \rm $*

clean : texclean 
	@-touch j.o j.mod j.x *.fix.f90; \rm *.o *.mod *.x *.fix.f90

include ${SEPINC}/SEP.bottom