MAKE_DEPEND=yes
include ProcMakefile2

#include /usr/local/SEP/include/SEP.top
include /opt/SEP/SEP6.4.6/include/SEP.top
################# To run in gomorrah  ####################
#include /usr/local/share/make/SEP.top
##########################################################

# SEP report
LATOPTS  = style=html
# SEG abstract
#LATOPTS = option=cmfonts class=seg97
#LATOPTS = option=web 'style=jfc,makeidx,html,amsmath'
# Reviewer 
#LATOPTS=class=geophysics
SIGDIR = /mnt/sepdata/1a/sigsbee_multiple
OLEVEL=2
UF90FLAGS =-O3 -C 
MAINLIBS=-lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep -L${LIBDIR}  -lsepgeef90
UF77LIBS = ${BEILIB}

WHEREFFTW = /net/koko/linux_local/FFTW/FFTW2/lib/lib
UF90LIBS=${MPILIBS} ${MAINLIBS} -lsepparf90 -lseppar_mpi -lsepgeef90 -lsep2df90 -static-libcxa -lguide -lseppar  -lsuperset -lsupersetf90 -lsuperset -lsepauxf90 -lsepmathf90 -lsep3df90 -lsep3d -lsepf90 -lsep $(WHEREFFTW)fftw.a $(WHEREFFTW)rfftw.a #-Bstatic

COLOR=y

BIN=./Bin
RESDIR=./Fig
R = ${RESDIR}
FIG=./Fig
SRCDIR=./Src

RESULTSER = arm11phw arm11phz4 arm11phz armig11T armig11W armig41T armig41W arr11phw arr11phz arrec11 arrec41 ars11phw ars11phz arsou11 arsou41 mig08 \
arrec51 arsou51 armig51T armig51W arsource arrecev arrphz arrphw

ER_v = mig08.v \
arrec11.v arrec41.v arrec51.v \
arsou11.v arsou41.v arsou51.v \
arsource.v arrecev.v \
armig51W.v armig41W.v armig11W.v \
armig51T.v  armig41T.v armig11T.v \
ars11phw.v arr11phw.v arrphw.v \
ars11phz.v arr11phz.v arrphz.v \
arm11phw.v arm11phz.v arm11phz4.v \

ER_pdf = mig08.pdf \
arrec11.pdf arrec41.pdf arrec51.pdf \
arsou11.pdf arsou41.pdf arsou51.pdf \
arsource.pdf arrecev.pdf \
armig51W.pdf armig41W.pdf armig11W.pdf \
armig51T.pdf  armig41T.pdf armig11T.pdf \
ars11phw.pdf arr11phw.pdf arrphw.pdf \
ars11phz.pdf arr11phz.pdf arrphz.pdf \
arm11phw.pdf arm11phz.pdf arm11phz4.pdf \

er_burn: clean
er_build: ${ER_v}
er_pdf: ${ER_pdf}

#############Progs################
#
#
OW_arealsht.x:
	rm -f ${BIN}/OW_arealsht.x;rm -f OW_arealsht.x;make OW_arealsht.x
	mv -f OW_arealsht.x ${BIN}/OW_arealsht.x
	mv -f *.o ${OBJ}/.
	mv -f *.mod ${OBJ}/expld_refl/.

OWmig2d_areal.x:
	rm -f ${BIN}/OWmig2d_areal.x;rm -f OWmig2d_areal.x;make OWmig2d_areal.x
	mv -f OWmig2d_areal.x ${BIN}/OWmig2d_areal.x
	mv -f *.o ${OBJ}/owmig/.
	mv -f *.mod ${OBJ}/.

###################Data#######################
# Kirchhoff modeling
4001data_hor.H:
	susynlv ref="1:0,.5.;10.,.5" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=20 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk1.H
	susynlv ref="1:0,0.;10.,2.5" nt=1501 dt=0.002 dxs=0.02 nxs=201 fpeak=20 fxs=0 dxo=0.02 nxo=201 fxo=-2 v00=2 > junk1.su
	< junk1.su Su2sep > junk2.H
	Math file1=junk1.H file2=junk2.H exp="file1+file2" | Infill3d |Window3d j1=4 > $@; echo n2=201 o2=-2 d2=0.02 n3=201 o3=0 d3=0.02 >> $@

# Preparing data to migration
4003data_recip_prep.H: # 4001data_hor.H
	< 4001data_hor.H Window3d |  /net/koko/bob/SEP_BOB/bin/Transf shotmig=y \
	is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=48 maxsize=1000 > j2.H
	< j2.H Transp reshape=1,3,5 plane=23 | Pad beg1=100 end1=100 > $@

# Preparing source to migration
4004data_recip_source.H: # 4001data_hor.H
	Window3d < 4001data_hor.H n3=1 f3=200 | Mute vmute=2 tmute=0.7 | Window3d n2=1 min2=-1 > x.H
	< x.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=48 maxsize=1000 > y.H
	< y.H Transp plane=35 | Cabs | Stack |Window3d > 4004_spec.H;
	Math file1=4004_spec.H exp="file1*0" > 4004_phse.H; Cmplx 4004_spec.H 4004_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

# Preparing source areal modeling
4004data_recip_source_pad.H: # 4001data_hor.H
	Window3d < 4001data_hor.H n3=1 f3=200 | Mute vmute=2 tmute=0.7 | Window3d n2=1 min2=-1 | Pad end1=600 > x.H
	< x.H /net/koko/bob/SEP_BOB/bin/Transf shotmig=y is2d=y wei=y f_min=2 f_min1=6 f_max1=40 f_max=48 maxsize=1000 > y.H
	< y.H Transp plane=35 | Cabs | Stack |Window3d > 4004_spec.H;
	Math file1=4004_spec.H exp="file1*0" > 4004_phse.H; Cmplx 4004_spec.H 4004_phse.H |\
	Transp reshape=1,3,5 plane=23| Transp reshape=1,3,5 plane=12 > $@

# Velocity for areal stuff
4000vel.H:
	Spike n1=600 n2=4000 o1=0 o2=-8 d1=0.004 d2=0.02 > 1.H;Math file1=1.H exp="file1*0+2" > $@

# Velocity for shot profile migration
4000vel_brad_pad.H:
	Spike n1=4000 n3=500 o1=-8 o3=0 d1=0.02 d3=0.004 > 1.H;Math file1=1.H exp="file1*0+2" > $@

# Shot profile migration
4008_mig_brad.H: # 4003data_recip_prep.H 4004data_recip_source_pad.H 4000vel_brad_pad.H
	@date
	/net/koko/homes/sep/alejandro/research/SPWEI_SERIAL/bin/LINUX86_64/ShotMig \
	D=4003data_recip_prep.H D_d=4004data_recip_source_pad.H V=4000vel_brad_pad.H R=$@ \
	pverb=2  nx_shots=200 ox_shots=0 \
	par=difr_go_sis_plane.P taper=y taper_pcent=2 marine=y minx=-4 maxx=8.2 & #restart=y

# preparing data
4008_mig_brad_transp.H: # 4008_mig_brad.H 
	< 4008_mig_brad.H Window3d min1=0 max1=4 max5=1 | Transp plane=13 > $@

# correct velocity areal shots for every 51 CIG
3990areal_rec_comb1_enc0_vel.H 3990areal_sou_comb1_enc0_vel.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=3990areal_rec_comb1_enc0_vel.H sou_out=3990areal_sou_comb1_enc0_vel.H ngap=51 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 > /dev/null 

# correct velocity areal shots for every 41 CIG
4100areal_rec_comb1_enc0_vel.H 4100areal_sou_comb1_enc0_vel.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4100areal_rec_comb1_enc0_vel.H sou_out=4100areal_sou_comb1_enc0_vel.H ngap=41 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 > /dev/null 

# correct velocity areal shots for every 11 CIG
4200areal_rec_comb1_enc0_vel.H 4200areal_sou_comb1_enc0_vel.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4200areal_rec_comb1_enc0_vel.H sou_out=4200areal_sou_comb1_enc0_vel.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=0 combine=1 > /dev/null 

# areal shots with encode=1 to let migration generate crosstalk between reflectors  for every 51 CIG
3991mig_comb1_enc0_vel.H: # 3990areal_rec_comb1_enc0_vel.H 3990areal_sou_comb1_enc0_vel.H  4000vel.H
	< 3990areal_rec_comb1_enc0_vel.H Transp plane=13 > 3990areal_rec_temp.H
	< 3990areal_sou_comb1_enc0_vel.H Transp plane=13 > 3990areal_sou_temp.H
	< 3990areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3990areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 > $@ 
	Rm 3990areal*temp.H

# areal shots with encode=1 to let migration generate crosstalk between reflectors  for every 41 CIG
4101mig_comb1_enc0_vel.H: # 4100areal_rec_comb1_enc0_vel.H 4100areal_sou_comb1_enc0_vel.H  4000vel.H
	< 4100areal_rec_comb1_enc0_vel.H Transp plane=13 > 4100areal_rec_temp.H
	< 4100areal_sou_comb1_enc0_vel.H Transp plane=13 > 4100areal_sou_temp.H
	< 4100areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4100areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 > $@ 
	Rm 4100areal*temp.H

# areal shots with encode=1 to let migration generate crosstalk between reflectors  for every 11 CIG
4201mig_comb1_enc0_vel.H: # 4200areal_rec_comb1_enc0_vel.H 4200areal_sou_comb1_enc0_vel.H  4000vel.H
	< 4200areal_rec_comb1_enc0_vel.H Transp plane=13 > 4200areal_rec_temp1.H
	< 4200areal_sou_comb1_enc0_vel.H Transp plane=13 > 4200areal_sou_temp1.H
	< 4200areal_rec_temp1.H ./bin/OWmig2d_areal.x sou=4200areal_sou_temp1.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=1 > $@ 
	Rm 4200areal*temp1.H

# areal shots with encode=0 time windowed imgaging condition for every 51 CIG
3990mig_comb1_enc0_vel.H: # 3990areal_rec_comb1_enc0_vel.H 3990areal_sou_comb1_enc0_vel.H  4000vel.H
	< 3990areal_rec_comb1_enc0_vel.H Transp plane=13 > 3990areal_rec_temp.H
	< 3990areal_sou_comb1_enc0_vel.H Transp plane=13 > 3990areal_sou_temp.H
	< 3990areal_rec_temp.H ./bin/OWmig2d_areal.x sou=3990areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 > $@ 
	Rm 3990areal*temp.H

# areal shots with encode=0 time windowed imgaging condition for every 41 CIG
4100mig_comb1_enc0_vel.H: # 4100areal_rec_comb1_enc0_vel.H 4100areal_sou_comb1_enc0_vel.H  4000vel.H
	< 4100areal_rec_comb1_enc0_vel.H Transp plane=13 > 4100areal_rec_temp.H
	< 4100areal_sou_comb1_enc0_vel.H Transp plane=13 > 4100areal_sou_temp.H
	< 4100areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4100areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 > $@ 
	Rm 4100areal*temp.H

# areal shots with encode=0 time windowed imgaging condition for every 11 CIG
4200mig_comb1_enc0_vel.H: # 4200areal_rec_comb1_enc0_vel.H 4200areal_sou_comb1_enc0_vel.H  4000vel.H
	< 4200areal_rec_comb1_enc0_vel.H Transp plane=13 > 4200areal_rec_temp.H
	< 4200areal_sou_comb1_enc0_vel.H Transp plane=13 > 4200areal_sou_temp.H
	< 4200areal_rec_temp.H ./bin/OWmig2d_areal.x sou=4200areal_sou_temp.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 deltat=0.02 dt=0.008 encode=0 > $@ 
	Rm 4200areal*temp.H

# correct velocity areal shot random encode in z
4300areal_rec_comb1_enc1z_vel_rnd01_pad.H 4300areal_sou_comb1_enc1z_vel_rnd01_pad.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4300areal_rec_comb1_enc1z_vel_rnd01_pad.H sou_out=4300areal_sou_comb1_enc1z_vel_rnd01_pad.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=05 combine=1> /dev/null 

4300areal_rec_comb1_enc1z_vel_rnd1020_pad.H 4300areal_sou_comb1_enc1z_vel_rnd1020_pad.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4300areal_rec_comb1_enc1z_vel_rnd1020_pad.H sou_out=4300areal_sou_comb1_enc1z_vel_rnd1020_pad.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=1020 combine=1> /dev/null 

4300areal_rec_comb1_enc1z_vel_rnd0100_pad.H 4300areal_sou_comb1_enc1z_vel_rnd0100_pad.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4300areal_rec_comb1_enc1z_vel_rnd0100_pad.H sou_out=4300areal_sou_comb1_enc1z_vel_rnd0100_pad.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=0100 combine=1> /dev/null 

4300areal_rec_comb1_enc1z_vel_rnd020_pad.H 4300areal_sou_comb1_enc1z_vel_rnd020_pad.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H 
	./bin/OW_arealsht.x rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4300areal_rec_comb1_enc1z_vel_rnd020_pad.H sou_out=4300areal_sou_comb1_enc1z_vel_rnd020_pad.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=0 encz=1 eps=0.15 rndcte=020 combine=1> /dev/null 

# correct velocity areal shot random encode in w
#4301areal04:
4301areal_rec_comb1_enc1w_vel_04pad.H 4301areal_sou_comb1_enc1w_vel_04pad.H: # 4008_mig_brad_transp.H 4004data_recip_source_pad.H 4000vel.H
	./bin/ rec_cig=4008_mig_brad_transp.H sou_cig=4008_mig_brad_transp.H sou=4004data_recip_source_pad.H vel=4000vel.H \
	rec_out=4301areal_rec_comb1_enc1w_vel_04pad.H sou_out=4301areal_sou_comb1_enc1w_vel_04pad.H ngap=11 nrefvel=1 vsal=6 \
	padtraces=41 encode=1 encw=1 encz=0 combine=1> /dev/null 

# migration with correct velocity areal shot random encode in z
4300mig_comb1_enc1z_vel_rnd01_pad.H: # 4300areal_sou_comb1_enc1z_vel_rnd01_pad.H 4300areal_sou_comb1_enc1z_vel_rnd01_pad.H 4000vel.H 
	< 4300areal_rec_comb1_enc1z_vel_rnd01_pad.H Transp plane=13 > 4300areal_rec_temp_rnd01.H
	< 4300areal_sou_comb1_enc1z_vel_rnd01_pad.H Transp plane=13 > 4300areal_sou_temp_rnd01.H
	< 4300areal_rec_temp_rnd01.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp_rnd01.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 > $@ 
	Rm 4300*temp_rnd01.H

# correct velocity areal shot random encode in z
4300mig_comb1_enc1z_vel_rnd1020_pad.H: # 4300areal_sou_comb1_enc1z_vel_rnd1020_pad.H 4300areal_sou_comb1_enc1z_vel_rnd1020_pad.H 4000vel.H 
	< 4300areal_rec_comb1_enc1z_vel_rnd1020_pad.H Transp plane=13 > 4300areal_rec_temp_rnd1020.H
	< 4300areal_sou_comb1_enc1z_vel_rnd1020_pad.H Transp plane=13 > 4300areal_sou_temp_rnd1020.H
	< 4300areal_rec_temp_rnd1020.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp_rnd1020.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 > $@ 
	Rm 4300*temp_rnd1020.H

# correct velocity areal shot random encode in z
4300mig_comb1_enc1z_vel_rnd0100_pad.H: # 4300areal_sou_comb1_enc1z_vel_rnd0100_pad.H 4300areal_sou_comb1_enc1z_vel_rnd0100_pad.H 4000vel.H 
	< 4300areal_rec_comb1_enc1z_vel_rnd0100_pad.H Transp plane=13 > 4300areal_rec_temp_rnd0100.H
	< 4300areal_sou_comb1_enc1z_vel_rnd0100_pad.H Transp plane=13 > 4300areal_sou_temp_rnd0100.H
	< 4300areal_rec_temp_rnd0100.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp_rnd0100.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 > $@ 
	Rm 4300*temp_rnd0100.H

# correct velocity areal shot random encode in z
4300mig_comb1_enc1z_vel_rnd020_pad.H: # 4300areal_sou_comb1_enc1z_vel_rnd020_pad.H 4300areal_sou_comb1_enc1z_vel_rnd020_pad.H 4000vel.H 
	< 4300areal_rec_comb1_enc1z_vel_rnd020_pad.H Transp plane=13 > 4300areal_rec_temp_rnd020.H
	< 4300areal_sou_comb1_enc1z_vel_rnd020_pad.H Transp plane=13 > 4300areal_sou_temp_rnd020.H
	< 4300areal_rec_temp_rnd020.H ./bin/OWmig2d_areal.x sou=4300areal_sou_temp_rnd020.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 > $@ 
	Rm 4300*temp_rnd020.H

# migration with correct velocity areal shot random encode in w
4301mig_comb1_enc1w_vel_04.H: # 4301areal_rec_comb1_enc1w_vel_04.H 4301areal_sou_comb1_enc1w_vel_04.H  4000vel.H
	< 4301areal_rec_comb1_enc1w_vel_04.H Transp plane=13 > 4301areal_rec_temp_04.H
	< 4301areal_sou_comb1_enc1w_vel_04.H Transp plane=13 > 4301areal_sou_temp_04.H
	< 4301areal_rec_temp_04.H ./bin/OWmig2d_areal.x sou=4301areal_sou_temp_04.H \
	vel=4000vel.H nrefvel=1 xmin=-1.22 xmax=5.2 ymin=0 ymax=0 zmin=0 zmax=1. nh=41 oh=-.4 dh=0.02 wbott=0 vsal=5 \
	traces_in_taper=0 every=0 encode=1 > $@ 
	Rm 4301*temp_04.H

# migration with correct velocity areal shot random encode in z - summing 4 different realizations
4300mig_comb1_enc1z_vel_rnd4real.H: # 4300mig_comb1_enc1z_vel_rnd01_pad.H 4300mig_comb1_enc1z_vel_rnd1020_pad.H 4300mig_comb1_enc1z_vel_rnd0100_pad.H 4300mig_comb1_enc1z_vel_rnd020_pad.H
	Math file1=4300mig_comb1_enc1z_vel_rnd01_pad.H file2=4300mig_comb1_enc1z_vel_rnd1020_pad.H file3=4300mig_comb1_enc1z_vel_rnd0100_pad.H \
	file4=4300mig_comb1_enc1z_vel_rnd020_pad.H exp=file1+file2+file3+file4 > $@

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

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

clean:
	rm -f *.o *.mod
	rm -f ${RESDIR}/*.pdf ${RESDIR}/*.ps ${RESDIR}/*.v
#	rm -f ${BINDIR}/*.x
#	rm -f *.H

include /opt/SEP/SEP6.4.6/include/SEP.bottom
#include /usr/local/share/make/SEP.bottom
#include /usr/local/SEP/include/SEP.bottom