# "make" will create all the .eps (or .ps) figures
# "make clean" will remove all the .v and .eps from Fig/ (but not .pdf)
# "make burn" will remove all files except .pdf in Fig/
# .pdf files were created with preview, so I didn't include them in the targets


# Common project definitions
S = ./Src
B = ./Bin
O = ./Obj
M = ./Mod
D = ./Data
R = ./Fig
P = ./Par
G = ./Grad

# SEPlib includes/libraries
SEPDIR = /opt/SEP/SEP7.0
SEPINC = $(SEPDIR)/include
SEPLIB = $(SEPDIR)/lib
SEPLIBFLAGS = -L$(SEPLIB) -lsepauxf90 -lsepgeef90 -lsepf90 -lsep2df90 -lsep3df90 -lsep3d -lsepf90 -lsep
 
# Common Libraries
CMNINCFLAGS = -I. -I$(SEPINC) -I$(S)
CMNLIBFLAGS = $(SEPLIBFLAGS)

# C Compiler Setup
CINCFLAGS = $(CMNINCFLAGS)
CLIBFLAGS = $(CMNLIBFLAGS)
CFLAGS    = -O3 -funroll-loops -fforce-addr -fopenmp -ffast-math -std=c99

# Cuda Runtime
CUDADIR = /opt/nvidia/cuda
CUDAINC = $(CUDADIR)/include
CUDALIB = $(CUDADIR)/lib64

# Cuda SDK
CUDASDKDIR = /opt/nvidia/cuda-2.3_sdk/C
CUDASDKINC = $(CUDASDKDIR)/common/inc
CUDASDKLIB = $(CUDASDKDIR)/lib

# CUDA-specific libraries from runtime and SDK
CUINCFLAGS = -I$(CUDAINC) -I$(CUDASDKINC)
CULIBFLAGS = -L$(CUDALIB) -L$(CUDASDKLIB) -lcuda -lcutil -lcublas -lcufft

CUFLAGS = -O3 --compiler-options $(CUINCFLAGS) $(CULIBFLAGS) $(CINCFLAGS) $(CLIBFLAGS)
#CUFLAGS = -maxrregcount=16 --compiler-options 

# Cuda Compiler
NVCC = /opt/nvidia/cuda/bin/nvcc


# Fortran Includes
F90INCFLAGS  = $(CMNINCFLAGS) $(CUINCFLAGS)
F90LIBFLAGS  = $(CMNLIBFLAGS) $(CULIBFLAGS)
F90FLAGS     = -O3 -module $(M)
CC   = gcc
CXX  = /lib/cpp
LINK = g++
F90  = ifort

default: results

################################################################################
# CUDA build rules for objects, library modules, and executables
$(O)/%.cu.o : $(S)/%.cu
	#Compiling CUDA $^
	$(NVCC) $< -c -o $@     $(CUFLAGS) $(CUINCPATH) $(CULIBPATH) $(CULIBS) #--ptxas-options=-v  --gpu-name sm_13

### C build rules ###
$(O)/%.c.o: $(S)/%.c
	#Compiling C $^
	$(CC) -c -o $@   $^  -DSOURCE="$^" $(CFLAGS) $(CINCFLAGS)

### Fortran rules ###
$(O)/%.f90.o: $(S)/%.f90
	#Preprocessing FORTRAN $^
	$(CXX) $(F90INCFLAGS) -DSOURCE="'$^'" < $(S)/$*.f90  > $(O)/$*.fix.f90
	#Compiling FORTRAN $@
	$(F90) $(F90INCFLAGS) $(F90FLAGS)  -c -o $(O)/$*.f90.o $(O)/$*.fix.f90

$(B)/%.x: $(O)/%.f90.o
	#Linking FORTRAN $@
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)


### Compile Programs #########################################################

$(B)/Forward.x: 	$(O)/gpu_model_wrap.c.o \
						$(O)/gpu_model.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/Forward.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/Forward2.x: 	$(O)/gpu_model_wrap.c.o \
						$(O)/gpu_model.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/Forward2.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/Forward_born.x: 	$(O)/gpu_model_wrap.c.o \
						$(O)/gpu_model.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/Forward_born.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/RTM.x: 	$(O)/gpu_prop_wrap.c.o \
						$(O)/gpu_prop.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/RTM.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/WEMVA.x: 	$(O)/gpu_prop_wrap.c.o \
						$(O)/gpu_prop.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/WEMVA.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/WEMVAA.x: 	$(O)/gpu_prop2_wrap.c.o \
						$(O)/gpu_prop2.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/WEMVAA.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

$(B)/WEMVAI.x: 	$(O)/gpu_prop2_wrap.c.o \
						$(O)/gpu_prop2.cu.o \
						$(O)/tomo_mod.f90.o \
						$(O)/angle_mod.f90.o \
						$(O)/WEMVAI.f90.o
	$(F90) -o $@ $^     $(F90FLAGS)  $(F90INCFLAGS) $(F90LIBFLAGS)

### Prepare input data #########################################################

wemva26_wavelet_15f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=15 tdelay=0.3 > wemva26_tmp.H
	Math file1=wemva26_tmp.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva26_tmp.H

wemva26_wavelet_5f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=5 tdelay=0.3 > wemva26_tmp2.H
	Math file1=wemva26_tmp2.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva26_tmp2.H

### Reflectivity #########################################################

wemva26_anom%_bg_pad.H: $(D)/wemva26_anom%.H
	Pad < $< beg1=300 end1=300 beg2=300 end2=300 extend=1 > $@

wemva26_anom0_rf_pad.H: $(D)/wemva26_anom0.H ${B}/Refl0.x
	${B}/Refl0.x < $< depth=110 | Pad beg1=300 end1=300 beg2=300 end2=300 extend=1 > $@

wemva26_anom2_drf_pad.H: $(D)/wemva26_anom2.H ${B}/Refl3.x
	${B}/Refl3.x < $< | Pad beg1=300 end1=300 beg2=300 end2=300 extend=1 > $@

### Run Modeling #########################################################

wemva26_data_anom2.H: wemva26_anom2_bg_pad.H wemva26_anom0_rf_pad.H wemva26_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model.p
	${B}/Forward_born.x < $< refl=wemva26_anom0_rf_pad.H wavfile=wemva26_wavelet_15f_3.0ms.H par=${P}/model.p > $@

wemva26_datad_anom2.H: wemva26_anom0_bg_pad.H wemva26_anom2_drf_pad.H wemva26_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model.p
	${B}/Forward_born.x < $< refl=wemva26_anom2_drf_pad.H wavfile=wemva26_wavelet_15f_3.0ms.H par=${P}/model.p > $@

wemva26_datad_lf_anom2.H: wemva26_anom0_bg_pad.H wemva26_anom2_drf_pad.H wemva26_wavelet_5f_3.0ms.H ${B}/Forward_born.x ${P}/model.p
	${B}/Forward_born.x < $< refl=wemva26_anom2_drf_pad.H wavfile=wemva26_wavelet_5f_3.0ms.H par=${P}/model.p > $@

### Run Imaging #########################################################

wemva26_imaged_bg_anom2.H: wemva26_anom0_bg_pad.H wemva26_datad_anom2.H wemva26_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/image.p
	${B}/RTM.x < $< data_obs=wemva26_datad_anom2.H wavfile=wemva26_wavelet_15f_3.0ms.H par=${P}/image.p > $@

wemva26_imaged_lf_bg_anom2.H: wemva26_anom0_bg_pad.H wemva26_datad_lf_anom2.H wemva26_wavelet_5f_3.0ms.H ${B}/RTM.x ${P}/image.p
	${B}/RTM.x < $< data_obs=wemva26_datad_lf_anom2.H wavfile=wemva26_wavelet_5f_3.0ms.H par=${P}/image.p > $@

wemva26_imaged_bg_anom2_scaled.H: wemva26_imaged_bg_anom2.H
	Scale < $< rscale=0.0007 > $@

wemva26_imaged_lf_bg_anom2_scaled.H: wemva26_imaged_lf_bg_anom2.H
	Scale < $< rscale=0.0007 > $@

wemva26_imaged_bg_anom2_fk.H: wemva26_imaged_bg_anom2_scaled.H
	Pad < $< n1out=2048 n2out=2048 | Rtoc | Ft3d sign1=1 sign2=1 center1=1 center2=1 maxsize=1000 | Cabs > $@

wemva26_imaged_lf_bg_anom2_fk.H: wemva26_imaged_lf_bg_anom2_scaled.H
	Pad < $< n1out=2048 n2out=2048 | Rtoc | Ft3d sign1=1 sign2=1 center1=1 center2=1 maxsize=1000 | Cabs > $@

wemva26_imaged_bg_anom2_stack.H: wemva26_imaged_bg_anom2_scaled.H
	Transp < $< | Stack3d > $@

wemva26_imaged_lf_bg_anom2_stack.H: wemva26_imaged_lf_bg_anom2_scaled.H
	Transp < $< | Stack3d > $@

wemva26_imaged_lf_bg_anom2_sstack.H: wemva26_imaged_lf_bg_anom2_scaled.H
	Transp < $< | Smooth rect1=20 | Stack3d > $@


### Run Forward WEMVA #########################################################

wemva26_deltaI_anom2.H: wemva26_anom0_bg_pad.H wemva26_anom2_bg_pad.H wemva26_data_anom2.H wemva26_wavelet_15f_3.0ms.H ${B}/WEMVA.x ${P}/parx.p
	${B}/WEMVA.x < $< data_obs=wemva26_data_anom2.H slow_true=wemva26_anom2_bg_pad.H wavfile=wemva26_wavelet_15f_3.0ms.H par=${P}/parx.p > $@

wemva26_deltaI_anom7.H: wemva26_deltaI_anom2.H
	Mute < $< tmute=1.7 vmute=0 > $@

wemva26_deltaS_anom7.H: wemva26_anom0_bg_pad.H wemva26_deltaI_anom7.H wemva26_data_anom2.H wemva26_wavelet_15f_3.0ms.H ${B}/WEMVAA.x ${P}/parx.p
	${B}/WEMVAA.x < $< data_obs=wemva26_data_anom2.H deltaI=wemva26_deltaI_anom7.H wavfile=wemva26_wavelet_15f_3.0ms.H par=${P}/parx.p > $@

wemva26_deltaS_anom7_scaled.H: wemva26_deltaS_anom7.H
	Scale < $< rscale=0.0000000003 > $@

wemva26_deltaS_anom7_fk.H: wemva26_deltaS_anom7_scaled.H
	Pad < $< n1out=2048 n2out=2048 | Rtoc | Ft3d sign1=1 sign2=1 center1=1 center2=1 maxsize=1000 | Cabs > $@

wemva26_deltaS_mfk0.H: wemva26_deltaS_anom7_fk.H wemva26_imaged_bg_anom2_fk.H 
	Math file1=$< file2=wemva26_imaged_bg_anom2_fk.H exp='file1*file2' > $@

wemva26_deltaS_mfk.H: wemva26_deltaS_anom7_fk.H wemva26_imaged_lf_bg_anom2_fk.H 
	Math file1=$< file2=wemva26_imaged_lf_bg_anom2_fk.H exp='file1*file2' > $@





### Prepare input data #########################################################

wemva30_wavelet_15f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=15 tdelay=0.3 > wemva30_tmp1.H
	Math file1=wemva30_tmp1.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva30_tmp1.H

wemva30_wavelet_5f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=5 tdelay=0.3 > wemva30_tmp2.H
	Math file1=wemva30_tmp2.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva30_tmp2.H

### Reflectivity #########################################################

wemva30_strue.H: $(D)/wemva30_anom1.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva30_sbg_iter0.H: $(D)/wemva30_anom0.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva30_rf.H: $(D)/wemva30_anom0.H ${B}/Refl.x
	${B}/Refl.x < $< depth=110 | Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

### Run Modeling #########################################################

wemva30_data_obs.H: wemva30_strue.H wemva30_rf.H wemva30_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model30.p
	${B}/Forward_born.x < $< refl=wemva30_rf.H wavfile=wemva30_wavelet_15f_3.0ms.H par=${P}/model30.p > $@

wemva30_data_lf_obs.H: wemva30_strue.H wemva30_wavelet_5f_3.0ms.H ${B}/Forward.x ${P}/model30.p
	${B}/Forward.x < $< wavfile=wemva30_wavelet_5f_3.0ms.H par=${P}/model30.p > $@

wemva30_data_iter%.H: wemva30_sbg_iter%.H wemva30_rf.H wemva30_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model30.p
	${B}/Forward_born.x < $< refl=wemva30_rf.H wavfile=wemva30_wavelet_15f_3.0ms.H par=${P}/model30.p > $@

wemva30_data_lf_iter%.H: wemva30_sbg_iter%.H wemva30_wavelet_5f_3.0ms.H ${B}/Forward.x ${P}/model30.p
	${B}/Forward.x < $< wavfile=wemva30_wavelet_5f_3.0ms.H par=${P}/model30.p > $@

wemva30_ddata_lf_iter%.H: wemva30_data_lf_iter%.H wemva30_data_lf_obs.H
	Add scale=-1.,1. $< wemva30_data_lf_obs.H > $@

### Run Imaging #########################################################

wemva30_image_flat.H: wemva30_sbg_iter0.H wemva30_data_iter0.H wemva30_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx30.p
	${B}/RTM.x < $< data_obs=wemva30_data_iter0.H wavfile=wemva30_wavelet_15f_3.0ms.H par=${P}/parx30.p > $@

wemva30_image_iter%.H: wemva30_sbg_iter%.H wemva30_data_obs.H wemva30_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx30.p
	${B}/RTM.x < $< data_obs=wemva30_data_obs.H wavfile=wemva30_wavelet_15f_3.0ms.H par=${P}/parx30.p > $@

wemva30_dimage_lf_iter%.H: wemva30_sbg_iter%.H wemva30_ddata_lf_iter%.H wemva30_wavelet_5f_3.0ms.H ${B}/RTM.x ${P}/parxlf30.p
	${B}/RTM.x < $< data_obs=wemva30_ddata_lf_iter$*.H wavfile=wemva30_wavelet_5f_3.0ms.H par=${P}/parxlf30.p > $@

### Run Forward WEMVA #########################################################

wemva30_deltaS_iter%.H: wemva30_sbg_iter%.H wemva30_image_iter%.H wemva30_image_flat.H wemva30_data_obs.H wemva30_wavelet_15f_3.0ms.H ${B}/WEMVAI.x ${P}/parx30.p
	${B}/WEMVAI.x < $< image_obs=wemva30_image_iter$*.H image_flat=wemva30_image_flat.H data_obs=wemva30_data_obs.H wavfile=wemva30_wavelet_15f_3.0ms.H par=${P}/parx30.p > $@

### wemva30 (forward- and back-scattered) iterate #########################################################

wemva30_sbg_iter1.H: ${D}/wemva30_anom1.H wemva30_sbg_iter0.H wemva30_deltaS_iter0.H wemva30_dimage_lf_iter0.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter0.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter0.H
	${B}/Comb2.x < $< iter=1 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter0.H dm1=wemva30_deltaS_iter0.H dm2=wemva30_dimage_lf_iter0.H > wemva30_tmp2_iter0.H
	Pad < wemva30_tmp2_iter0.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter0.H
	Add wemva30_tmp3_iter0.H wemva30_sbg_iter0.H > $@
	Rm wemva30_tmp_iter0.H wemva30_tmp2_iter0.H wemva30_tmp3_iter0.H

wemva30_sbg_iter2.H: ${D}/wemva30_anom1.H wemva30_sbg_iter1.H wemva30_deltaS_iter1.H wemva30_dimage_lf_iter1.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter1.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter1.H
	${B}/Comb2.x < $< iter=2 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter1.H dm1=wemva30_deltaS_iter1.H dm2=wemva30_dimage_lf_iter1.H > wemva30_tmp2_iter1.H
	Pad < wemva30_tmp2_iter1.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter1.H
	Add wemva30_tmp3_iter1.H wemva30_sbg_iter1.H > $@
	Rm wemva30_tmp_iter1.H wemva30_tmp2_iter1.H wemva30_tmp3_iter1.H

wemva30_sbg_iter3.H: ${D}/wemva30_anom1.H wemva30_sbg_iter2.H wemva30_deltaS_iter2.H wemva30_dimage_lf_iter2.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter2.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter2.H
	${B}/Comb2.x < $< iter=3 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter2.H dm1=wemva30_deltaS_iter2.H dm2=wemva30_dimage_lf_iter2.H > wemva30_tmp2_iter2.H
	Pad < wemva30_tmp2_iter2.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter2.H
	Add wemva30_tmp3_iter2.H wemva30_sbg_iter2.H > $@
	Rm wemva30_tmp_iter2.H wemva30_tmp2_iter2.H wemva30_tmp3_iter2.H

wemva30_sbg_iter4.H: ${D}/wemva30_anom1.H wemva30_sbg_iter3.H wemva30_deltaS_iter3.H wemva30_dimage_lf_iter3.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter3.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter3.H
	${B}/Comb2.x < $< iter=4 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter3.H dm1=wemva30_deltaS_iter3.H dm2=wemva30_dimage_lf_iter3.H > wemva30_tmp2_iter3.H
	Pad < wemva30_tmp2_iter3.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter3.H
	Add wemva30_tmp3_iter3.H wemva30_sbg_iter3.H > $@
	Rm wemva30_tmp_iter3.H wemva30_tmp2_iter3.H wemva30_tmp3_iter3.H

wemva30_sbg_iter5.H: ${D}/wemva30_anom1.H wemva30_sbg_iter4.H wemva30_deltaS_iter4.H wemva30_dimage_lf_iter4.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter4.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter4.H
	${B}/Comb2.x < $< iter=5 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter4.H dm1=wemva30_deltaS_iter4.H dm2=wemva30_dimage_lf_iter4.H > wemva30_tmp2_iter4.H
	Pad < wemva30_tmp2_iter4.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter4.H
	Add wemva30_tmp3_iter4.H wemva30_sbg_iter4.H > $@
	Rm wemva30_tmp_iter4.H wemva30_tmp2_iter4.H wemva30_tmp3_iter4.H

wemva30_sbg_iter6.H: ${D}/wemva30_anom1.H wemva30_sbg_iter5.H wemva30_deltaS_iter5.H wemva30_dimage_lf_iter5.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter5.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter5.H
	${B}/Comb2.x < $< iter=6 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter5.H dm1=wemva30_deltaS_iter5.H dm2=wemva30_dimage_lf_iter5.H > wemva30_tmp2_iter5.H
	Pad < wemva30_tmp2_iter5.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter5.H
	Add wemva30_tmp3_iter5.H wemva30_sbg_iter5.H > $@
	Rm wemva30_tmp_iter5.H wemva30_tmp2_iter5.H wemva30_tmp3_iter5.H

wemva30_sbg_iter7.H: ${D}/wemva30_anom1.H wemva30_sbg_iter6.H wemva30_deltaS_iter6.H wemva30_dimage_lf_iter6.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter6.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter6.H
	${B}/Comb2.x < $< iter=7 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter6.H dm1=wemva30_deltaS_iter6.H dm2=wemva30_dimage_lf_iter6.H > wemva30_tmp2_iter6.H
	Pad < wemva30_tmp2_iter6.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter6.H
	Add wemva30_tmp3_iter6.H wemva30_sbg_iter6.H > $@
	Rm wemva30_tmp_iter6.H wemva30_tmp2_iter6.H wemva30_tmp3_iter6.H

wemva30_sbg_iter8.H: ${D}/wemva30_anom1.H wemva30_sbg_iter7.H wemva30_deltaS_iter7.H wemva30_dimage_lf_iter7.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter7.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter7.H
	${B}/Comb2.x < $< iter=8 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter7.H dm1=wemva30_deltaS_iter7.H dm2=wemva30_dimage_lf_iter7.H > wemva30_tmp2_iter7.H
	Pad < wemva30_tmp2_iter7.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter7.H
	Add wemva30_tmp3_iter7.H wemva30_sbg_iter7.H > $@
	Rm wemva30_tmp_iter7.H wemva30_tmp2_iter7.H wemva30_tmp3_iter7.H

wemva30_sbg_iter9.H: ${D}/wemva30_anom1.H wemva30_sbg_iter8.H wemva30_deltaS_iter8.H wemva30_dimage_lf_iter8.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter8.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter8.H
	${B}/Comb2.x < $< iter=9 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter8.H dm1=wemva30_deltaS_iter8.H dm2=wemva30_dimage_lf_iter8.H > wemva30_tmp2_iter8.H
	Pad < wemva30_tmp2_iter8.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter8.H
	Add wemva30_tmp3_iter8.H wemva30_sbg_iter8.H > $@
	Rm wemva30_tmp_iter8.H wemva30_tmp2_iter8.H wemva30_tmp3_iter8.H

wemva30_sbg_iter10.H: ${D}/wemva30_anom1.H wemva30_sbg_iter9.H wemva30_deltaS_iter9.H wemva30_dimage_lf_iter9.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter9.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter9.H
	${B}/Comb2.x < $< iter=10 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter9.H dm1=wemva30_deltaS_iter9.H dm2=wemva30_dimage_lf_iter9.H > wemva30_tmp2_iter9.H
	Pad < wemva30_tmp2_iter9.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter9.H
	Add wemva30_tmp3_iter9.H wemva30_sbg_iter9.H > $@
	Rm wemva30_tmp_iter9.H wemva30_tmp2_iter9.H wemva30_tmp3_iter9.H

wemva30_sbg_iter11.H: ${D}/wemva30_anom1.H wemva30_sbg_iter10.H wemva30_deltaS_iter10.H wemva30_dimage_lf_iter10.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter10.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter10.H
	${B}/Comb2.x < $< iter=11 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter10.H dm1=wemva30_deltaS_iter10.H dm2=wemva30_dimage_lf_iter10.H > wemva30_tmp2_iter10.H
	Pad < wemva30_tmp2_iter10.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter10.H
	Add wemva30_tmp3_iter10.H wemva30_sbg_iter10.H > $@
	Rm wemva30_tmp_iter10.H wemva30_tmp2_iter10.H wemva30_tmp3_iter10.H

wemva30_sbg_iter12.H: ${D}/wemva30_anom1.H wemva30_sbg_iter11.H wemva30_deltaS_iter11.H wemva30_dimage_lf_iter11.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter11.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter11.H
	${B}/Comb2.x < $< iter=12 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter11.H dm1=wemva30_deltaS_iter11.H dm2=wemva30_dimage_lf_iter11.H > wemva30_tmp2_iter11.H
	Pad < wemva30_tmp2_iter11.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter11.H
	Add wemva30_tmp3_iter11.H wemva30_sbg_iter11.H > $@
	Rm wemva30_tmp_iter11.H wemva30_tmp2_iter11.H wemva30_tmp3_iter11.H

wemva30_sbg_iter13.H: ${D}/wemva30_anom1.H wemva30_sbg_iter12.H wemva30_deltaS_iter12.H wemva30_dimage_lf_iter12.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter12.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter12.H
	${B}/Comb2.x < $< iter=13 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter12.H dm1=wemva30_deltaS_iter12.H dm2=wemva30_dimage_lf_iter12.H > wemva30_tmp2_iter12.H
	Pad < wemva30_tmp2_iter12.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter12.H
	Add wemva30_tmp3_iter12.H wemva30_sbg_iter12.H > $@
	Rm wemva30_tmp_iter12.H wemva30_tmp2_iter12.H wemva30_tmp3_iter12.H

wemva30_sbg_iter14.H: ${D}/wemva30_anom1.H wemva30_sbg_iter13.H wemva30_deltaS_iter13.H wemva30_dimage_lf_iter13.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter13.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter13.H
	${B}/Comb2.x < $< iter=14 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter13.H dm1=wemva30_deltaS_iter13.H dm2=wemva30_dimage_lf_iter13.H > wemva30_tmp2_iter13.H
	Pad < wemva30_tmp2_iter13.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter13.H
	Add wemva30_tmp3_iter13.H wemva30_sbg_iter13.H > $@
	Rm wemva30_tmp_iter13.H wemva30_tmp2_iter13.H wemva30_tmp3_iter13.H

wemva30_sbg_iter15.H: ${D}/wemva30_anom1.H wemva30_sbg_iter14.H wemva30_deltaS_iter14.H wemva30_dimage_lf_iter14.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter14.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter14.H
	${B}/Comb2.x < $< iter=15 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter14.H dm1=wemva30_deltaS_iter14.H dm2=wemva30_dimage_lf_iter14.H > wemva30_tmp2_iter14.H
	Pad < wemva30_tmp2_iter14.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter14.H
	Add wemva30_tmp3_iter14.H wemva30_sbg_iter14.H > $@
	Rm wemva30_tmp_iter14.H wemva30_tmp2_iter14.H wemva30_tmp3_iter14.H

wemva30_sbg_iter16.H: ${D}/wemva30_anom1.H wemva30_sbg_iter15.H wemva30_deltaS_iter15.H wemva30_dimage_lf_iter15.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter15.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter15.H
	${B}/Comb2.x < $< iter=16 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter15.H dm1=wemva30_deltaS_iter15.H dm2=wemva30_dimage_lf_iter15.H > wemva30_tmp2_iter15.H
	Pad < wemva30_tmp2_iter15.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter15.H
	Add wemva30_tmp3_iter15.H wemva30_sbg_iter15.H > $@
	Rm wemva30_tmp_iter15.H wemva30_tmp2_iter15.H wemva30_tmp3_iter15.H

wemva30_sbg_iter17.H: ${D}/wemva30_anom1.H wemva30_sbg_iter16.H wemva30_deltaS_iter16.H wemva30_dimage_lf_iter16.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter16.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter16.H
	${B}/Comb2.x < $< iter=17 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter16.H dm1=wemva30_deltaS_iter16.H dm2=wemva30_dimage_lf_iter16.H > wemva30_tmp2_iter16.H
	Pad < wemva30_tmp2_iter16.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter16.H
	Add wemva30_tmp3_iter16.H wemva30_sbg_iter16.H > $@
	Rm wemva30_tmp_iter16.H wemva30_tmp2_iter16.H wemva30_tmp3_iter16.H

wemva30_sbg_iter18.H: ${D}/wemva30_anom1.H wemva30_sbg_iter17.H wemva30_deltaS_iter17.H wemva30_dimage_lf_iter17.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter17.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter17.H
	${B}/Comb2.x < $< iter=18 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter17.H dm1=wemva30_deltaS_iter17.H dm2=wemva30_dimage_lf_iter17.H > wemva30_tmp2_iter17.H
	Pad < wemva30_tmp2_iter17.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter17.H
	Add wemva30_tmp3_iter17.H wemva30_sbg_iter17.H > $@
	Rm wemva30_tmp_iter17.H wemva30_tmp2_iter17.H wemva30_tmp3_iter17.H

wemva30_sbg_iter19.H: ${D}/wemva30_anom1.H wemva30_sbg_iter18.H wemva30_deltaS_iter18.H wemva30_dimage_lf_iter18.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter18.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter18.H
	${B}/Comb2.x < $< iter=19 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter18.H dm1=wemva30_deltaS_iter18.H dm2=wemva30_dimage_lf_iter18.H > wemva30_tmp2_iter18.H
	Pad < wemva30_tmp2_iter18.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter18.H
	Add wemva30_tmp3_iter18.H wemva30_sbg_iter18.H > $@
	Rm wemva30_tmp_iter18.H wemva30_tmp2_iter18.H wemva30_tmp3_iter18.H

wemva30_sbg_iter20.H: ${D}/wemva30_anom1.H wemva30_sbg_iter19.H wemva30_deltaS_iter19.H wemva30_dimage_lf_iter19.H ${D}/wemva30_a1.H ${D}/wemva30_a2.H ${B}/Comb2.x
	Window3d < wemva30_sbg_iter19.H f1=200 f2=200 n1=121 n2=201 > wemva30_tmp_iter19.H
	${B}/Comb2.x < $< iter=20 a1=${D}/wemva30_a1.H a2=${D}/wemva30_a2.H m0=wemva30_tmp_iter19.H dm1=wemva30_deltaS_iter19.H dm2=wemva30_dimage_lf_iter19.H > wemva30_tmp2_iter19.H
	Pad < wemva30_tmp2_iter19.H beg1=200 end1=200 beg2=200 end2=200 > wemva30_tmp3_iter19.H
	Add wemva30_tmp3_iter19.H wemva30_sbg_iter19.H > $@
	Rm wemva30_tmp_iter19.H wemva30_tmp2_iter19.H wemva30_tmp3_iter19.H






### Prepare input data #########################################################

wemva31_wavelet_15f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=15 tdelay=0.3 > wemva31_tmp1.H
	Math file1=wemva31_tmp1.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva31_tmp1.H

wemva31_wavelet_5f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=5 tdelay=0.3 > wemva31_tmp2.H
	Math file1=wemva31_tmp2.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva31_tmp2.H

### Reflectivity #########################################################

wemva31_strue.H: $(D)/wemva31_anom1.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva31_sbg_iter0.H: $(D)/wemva31_anom0.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva31_rf.H: $(D)/wemva31_anom0.H $(D)/wemva31_anom1.H ${B}/Refl2.x
	Add scale=1.,-1. $(D)/wemva31_anom1.H $(D)/wemva31_anom0.H | ${B}/Refl2.x depth=110 | Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

### Run Modeling #########################################################

wemva31_data_obs.H: wemva31_strue.H wemva31_rf.H wemva31_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model31.p
	${B}/Forward_born.x < $< refl=wemva31_rf.H wavfile=wemva31_wavelet_15f_3.0ms.H par=${P}/model31.p > $@

wemva31_data_lf_obs.H: wemva31_strue.H wemva31_rf.H wemva31_wavelet_5f_3.0ms.H ${B}/Forward_born.x ${P}/model31.p
	${B}/Forward_born.x < $< refl=wemva31_rf.H wavfile=wemva31_wavelet_5f_3.0ms.H par=${P}/model31.p > $@

wemva31_data_iter%.H: wemva31_sbg_iter%.H wemva31_rf.H wemva31_wavelet_15f_3.0ms.H ${B}/Forward_born.x ${P}/model31.p
	${B}/Forward_born.x < $< refl=wemva31_rf.H wavfile=wemva31_wavelet_15f_3.0ms.H par=${P}/model31.p > $@

wemva31_data_lf_iter%.H: wemva31_sbg_iter%.H wemva31_rf.H wemva31_wavelet_5f_3.0ms.H ${B}/Forward_born.x ${P}/model31.p
	${B}/Forward_born.x < $< refl=wemva31_rf.H wavfile=wemva31_wavelet_5f_3.0ms.H par=${P}/model31.p > $@

wemva31_ddata_lf_iter%.H: wemva31_data_lf_iter%.H wemva31_data_lf_obs.H
	Add scale=1.,-1. $< wemva31_data_lf_obs.H > $@

### Run Imaging #########################################################

wemva31_image_flat.H: wemva31_sbg_iter0.H wemva31_data_iter0.H wemva31_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx31.p
	${B}/RTM.x < $< data_obs=wemva31_data_iter0.H wavfile=wemva31_wavelet_15f_3.0ms.H par=${P}/parx31.p > $@

wemva31_image_iter%.H: wemva31_sbg_iter%.H wemva31_data_obs.H wemva31_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx31.p
	${B}/RTM.x < $< data_obs=wemva31_data_obs.H wavfile=wemva31_wavelet_15f_3.0ms.H par=${P}/parx31.p > $@

wemva31_dimage_lf_iter%.H: wemva31_sbg_iter%.H wemva31_ddata_lf_iter%.H wemva31_wavelet_5f_3.0ms.H ${B}/RTM.x ${P}/parxlf31.p
	${B}/RTM.x < $< data_obs=wemva31_ddata_lf_iter$*.H wavfile=wemva31_wavelet_5f_3.0ms.H par=${P}/parxlf31.p > $@

### Run Forward WEMVA #########################################################

wemva31_deltaS_iter%.H: wemva31_sbg_iter%.H wemva31_image_iter%.H wemva31_image_flat.H wemva31_data_obs.H wemva31_wavelet_15f_3.0ms.H ${B}/WEMVAI.x ${P}/parx31.p
	${B}/WEMVAI.x < $< image_obs=wemva31_image_iter$*.H image_flat=wemva31_image_flat.H data_obs=wemva31_data_obs.H wavfile=wemva31_wavelet_15f_3.0ms.H par=${P}/parx31.p > $@

### wemva31 (forward-scattered) iterate #########################################################

wemva31_sbg_iter1.H: ${D}/wemva31_anom1.H wemva31_sbg_iter0.H wemva31_deltaS_iter0.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter0.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter0.H
	${B}/Comb1.x < $< iter=1 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter0.H dm=wemva31_deltaS_iter0.H > wemva31_tmp2_iter0.H
	Pad < wemva31_tmp2_iter0.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter0.H
	Add wemva31_tmp3_iter0.H wemva31_sbg_iter0.H > $@
	Rm wemva31_tmp_iter0.H wemva31_tmp2_iter0.H wemva31_tmp3_iter0.H

wemva31_sbg_iter2.H: ${D}/wemva31_anom1.H wemva31_sbg_iter1.H wemva31_deltaS_iter1.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter1.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter1.H
	${B}/Comb1.x < $< iter=2 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter1.H dm=wemva31_deltaS_iter1.H > wemva31_tmp2_iter1.H
	Pad < wemva31_tmp2_iter1.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter1.H
	Add wemva31_tmp3_iter1.H wemva31_sbg_iter1.H > $@
	Rm wemva31_tmp_iter1.H wemva31_tmp2_iter1.H wemva31_tmp3_iter1.H

wemva31_sbg_iter3.H: ${D}/wemva31_anom1.H wemva31_sbg_iter2.H wemva31_deltaS_iter2.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter2.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter2.H
	${B}/Comb1.x < $< iter=3 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter2.H dm=wemva31_deltaS_iter2.H > wemva31_tmp2_iter2.H
	Pad < wemva31_tmp2_iter2.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter2.H
	Add wemva31_tmp3_iter2.H wemva31_sbg_iter2.H > $@
	Rm wemva31_tmp_iter2.H wemva31_tmp2_iter2.H wemva31_tmp3_iter2.H

wemva31_sbg_iter4.H: ${D}/wemva31_anom1.H wemva31_sbg_iter3.H wemva31_deltaS_iter3.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter3.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter3.H
	${B}/Comb1.x < $< iter=4 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter3.H dm=wemva31_deltaS_iter3.H > wemva31_tmp2_iter3.H
	Pad < wemva31_tmp2_iter3.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter3.H
	Add wemva31_tmp3_iter3.H wemva31_sbg_iter3.H > $@
	Rm wemva31_tmp_iter3.H wemva31_tmp2_iter3.H wemva31_tmp3_iter3.H

wemva31_sbg_iter5.H: ${D}/wemva31_anom1.H wemva31_sbg_iter4.H wemva31_deltaS_iter4.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter4.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter4.H
	${B}/Comb1.x < $< iter=5 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter4.H dm=wemva31_deltaS_iter4.H > wemva31_tmp2_iter4.H
	Pad < wemva31_tmp2_iter4.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter4.H
	Add wemva31_tmp3_iter4.H wemva31_sbg_iter4.H > $@
	Rm wemva31_tmp_iter4.H wemva31_tmp2_iter4.H wemva31_tmp3_iter4.H

wemva31_sbg_iter6.H: ${D}/wemva31_anom1.H wemva31_sbg_iter5.H wemva31_deltaS_iter5.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter5.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter5.H
	${B}/Comb1.x < $< iter=6 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter5.H dm=wemva31_deltaS_iter5.H > wemva31_tmp2_iter5.H
	Pad < wemva31_tmp2_iter5.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter5.H
	Add wemva31_tmp3_iter5.H wemva31_sbg_iter5.H > $@
	Rm wemva31_tmp_iter5.H wemva31_tmp2_iter5.H wemva31_tmp3_iter5.H

wemva31_sbg_iter7.H: ${D}/wemva31_anom1.H wemva31_sbg_iter6.H wemva31_deltaS_iter6.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter6.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter6.H
	${B}/Comb1.x < $< iter=7 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter6.H dm=wemva31_deltaS_iter6.H > wemva31_tmp2_iter6.H
	Pad < wemva31_tmp2_iter6.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter6.H
	Add wemva31_tmp3_iter6.H wemva31_sbg_iter6.H > $@
	Rm wemva31_tmp_iter6.H wemva31_tmp2_iter6.H wemva31_tmp3_iter6.H

wemva31_sbg_iter8.H: ${D}/wemva31_anom1.H wemva31_sbg_iter7.H wemva31_deltaS_iter7.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter7.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter7.H
	${B}/Comb1.x < $< iter=8 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter7.H dm=wemva31_deltaS_iter7.H > wemva31_tmp2_iter7.H
	Pad < wemva31_tmp2_iter7.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter7.H
	Add wemva31_tmp3_iter7.H wemva31_sbg_iter7.H > $@
	Rm wemva31_tmp_iter7.H wemva31_tmp2_iter7.H wemva31_tmp3_iter7.H

wemva31_sbg_iter9.H: ${D}/wemva31_anom1.H wemva31_sbg_iter8.H wemva31_deltaS_iter8.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter8.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter8.H
	${B}/Comb1.x < $< iter=9 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter8.H dm=wemva31_deltaS_iter8.H > wemva31_tmp2_iter8.H
	Pad < wemva31_tmp2_iter8.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter8.H
	Add wemva31_tmp3_iter8.H wemva31_sbg_iter8.H > $@
	Rm wemva31_tmp_iter8.H wemva31_tmp2_iter8.H wemva31_tmp3_iter8.H

wemva31_sbg_iter10.H: ${D}/wemva31_anom1.H wemva31_sbg_iter9.H wemva31_deltaS_iter9.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter9.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter9.H
	${B}/Comb1.x < $< iter=10 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter9.H dm=wemva31_deltaS_iter9.H > wemva31_tmp2_iter9.H
	Pad < wemva31_tmp2_iter9.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter9.H
	Add wemva31_tmp3_iter9.H wemva31_sbg_iter9.H > $@
	Rm wemva31_tmp_iter9.H wemva31_tmp2_iter9.H wemva31_tmp3_iter9.H

wemva31_sbg_iter11.H: ${D}/wemva31_anom1.H wemva31_sbg_iter10.H wemva31_deltaS_iter10.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter10.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter10.H
	${B}/Comb1.x < $< iter=11 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter10.H dm=wemva31_deltaS_iter10.H > wemva31_tmp2_iter10.H
	Pad < wemva31_tmp2_iter10.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter10.H
	Add wemva31_tmp3_iter10.H wemva31_sbg_iter10.H > $@
	Rm wemva31_tmp_iter10.H wemva31_tmp2_iter10.H wemva31_tmp3_iter10.H

wemva31_sbg_iter12.H: ${D}/wemva31_anom1.H wemva31_sbg_iter11.H wemva31_deltaS_iter11.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter11.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter11.H
	${B}/Comb1.x < $< iter=12 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter11.H dm=wemva31_deltaS_iter11.H > wemva31_tmp2_iter11.H
	Pad < wemva31_tmp2_iter11.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter11.H
	Add wemva31_tmp3_iter11.H wemva31_sbg_iter11.H > $@
	Rm wemva31_tmp_iter11.H wemva31_tmp2_iter11.H wemva31_tmp3_iter11.H

wemva31_sbg_iter13.H: ${D}/wemva31_anom1.H wemva31_sbg_iter12.H wemva31_deltaS_iter12.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter12.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter12.H
	${B}/Comb1.x < $< iter=13 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter12.H dm=wemva31_deltaS_iter12.H > wemva31_tmp2_iter12.H
	Pad < wemva31_tmp2_iter12.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter12.H
	Add wemva31_tmp3_iter12.H wemva31_sbg_iter12.H > $@
	Rm wemva31_tmp_iter12.H wemva31_tmp2_iter12.H wemva31_tmp3_iter12.H

wemva31_sbg_iter14.H: ${D}/wemva31_anom1.H wemva31_sbg_iter13.H wemva31_deltaS_iter13.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter13.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter13.H
	${B}/Comb1.x < $< iter=14 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter13.H dm=wemva31_deltaS_iter13.H > wemva31_tmp2_iter13.H
	Pad < wemva31_tmp2_iter13.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter13.H
	Add wemva31_tmp3_iter13.H wemva31_sbg_iter13.H > $@
	Rm wemva31_tmp_iter13.H wemva31_tmp2_iter13.H wemva31_tmp3_iter13.H

wemva31_sbg_iter15.H: ${D}/wemva31_anom1.H wemva31_sbg_iter14.H wemva31_deltaS_iter14.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter14.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter14.H
	${B}/Comb1.x < $< iter=15 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter14.H dm=wemva31_deltaS_iter14.H > wemva31_tmp2_iter14.H
	Pad < wemva31_tmp2_iter14.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter14.H
	Add wemva31_tmp3_iter14.H wemva31_sbg_iter14.H > $@
	Rm wemva31_tmp_iter14.H wemva31_tmp2_iter14.H wemva31_tmp3_iter14.H

wemva31_sbg_iter16.H: ${D}/wemva31_anom1.H wemva31_sbg_iter15.H wemva31_deltaS_iter15.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter15.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter15.H
	${B}/Comb1.x < $< iter=16 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter15.H dm=wemva31_deltaS_iter15.H > wemva31_tmp2_iter15.H
	Pad < wemva31_tmp2_iter15.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter15.H
	Add wemva31_tmp3_iter15.H wemva31_sbg_iter15.H > $@
	Rm wemva31_tmp_iter15.H wemva31_tmp2_iter15.H wemva31_tmp3_iter15.H

wemva31_sbg_iter17.H: ${D}/wemva31_anom1.H wemva31_sbg_iter16.H wemva31_deltaS_iter16.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter16.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter16.H
	${B}/Comb1.x < $< iter=17 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter16.H dm=wemva31_deltaS_iter16.H > wemva31_tmp2_iter16.H
	Pad < wemva31_tmp2_iter16.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter16.H
	Add wemva31_tmp3_iter16.H wemva31_sbg_iter16.H > $@
	Rm wemva31_tmp_iter16.H wemva31_tmp2_iter16.H wemva31_tmp3_iter16.H

wemva31_sbg_iter18.H: ${D}/wemva31_anom1.H wemva31_sbg_iter17.H wemva31_deltaS_iter17.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter17.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter17.H
	${B}/Comb1.x < $< iter=18 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter17.H dm=wemva31_deltaS_iter17.H > wemva31_tmp2_iter17.H
	Pad < wemva31_tmp2_iter17.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter17.H
	Add wemva31_tmp3_iter17.H wemva31_sbg_iter17.H > $@
	Rm wemva31_tmp_iter17.H wemva31_tmp2_iter17.H wemva31_tmp3_iter17.H

wemva31_sbg_iter19.H: ${D}/wemva31_anom1.H wemva31_sbg_iter18.H wemva31_deltaS_iter18.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter18.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter18.H
	${B}/Comb1.x < $< iter=19 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter18.H dm=wemva31_deltaS_iter18.H > wemva31_tmp2_iter18.H
	Pad < wemva31_tmp2_iter18.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter18.H
	Add wemva31_tmp3_iter18.H wemva31_sbg_iter18.H > $@
	Rm wemva31_tmp_iter18.H wemva31_tmp2_iter18.H wemva31_tmp3_iter18.H

wemva31_sbg_iter20.H: ${D}/wemva31_anom1.H wemva31_sbg_iter19.H wemva31_deltaS_iter19.H ${D}/wemva31_a1.H ${B}/Comb1.x
	Window3d < wemva31_sbg_iter19.H f1=200 f2=200 n1=121 n2=201 > wemva31_tmp_iter19.H
	${B}/Comb1.x < $< iter=20 a1=${D}/wemva31_a1.H m0=wemva31_tmp_iter19.H dm=wemva31_deltaS_iter19.H > wemva31_tmp2_iter19.H
	Pad < wemva31_tmp2_iter19.H beg1=200 end1=200 beg2=200 end2=200 > wemva31_tmp3_iter19.H
	Add wemva31_tmp3_iter19.H wemva31_sbg_iter19.H > $@
	Rm wemva31_tmp_iter19.H wemva31_tmp2_iter19.H wemva31_tmp3_iter19.H




### Prepare input data #########################################################

wemva32_wavelet_5f_3.0ms.H:
	Wavelet n1=3000 d1=0.003 wavelet=ricker2 phase=0 fund=5 tdelay=0.3 > wemva32_tmp2.H
	Math file1=wemva32_tmp2.H exp='file1*1000' > $@
	echo "label1='Time(s)'"  >> $@
	Rm wemva32_tmp2.H

### Reflectivity #########################################################

wemva32_strue.H: $(D)/wemva32_anom1.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva32_sbg_iter0.H: $(D)/wemva32_anom0.H
	Pad < $< beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

wemva32_rf.H: $(D)/wemva32_anom0.H $(D)/wemva32_anom1.H ${B}/Refl2.x
	Add scale=1.,-1. $(D)/wemva32_anom1.H $(D)/wemva32_anom0.H | ${B}/Refl2.x depth=110 | Pad beg1=200 end1=200 beg2=200 end2=200 extend=1 > $@

### Run Modeling #########################################################

wemva32_data_lf_obs.H: wemva32_strue.H wemva32_wavelet_5f_3.0ms.H ${B}/Forward2.x ${P}/model32.p
	${B}/Forward2.x < $< wavfile=wemva32_wavelet_5f_3.0ms.H par=${P}/model32.p > $@

wemva32_data_lf_iter%.H: wemva32_sbg_iter%.H wemva32_wavelet_5f_3.0ms.H ${B}/Forward2.x ${P}/model32.p
	${B}/Forward2.x < $< wavfile=wemva32_wavelet_5f_3.0ms.H par=${P}/model32.p > $@

wemva32_ddata_lf_iter%.H: wemva32_data_lf_iter%.H wemva32_data_lf_obs.H
	Add scale=-1.,1. $< wemva32_data_lf_obs.H > $@

### Run Imaging #########################################################

wemva32_image_flat.H: wemva32_sbg_iter0.H wemva32_data_iter0.H wemva32_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx32.p
	${B}/RTM.x < $< data_obs=wemva32_data_iter0.H wavfile=wemva32_wavelet_15f_3.0ms.H par=${P}/parx32.p > $@

wemva32_image_iter%.H: wemva32_sbg_iter%.H wemva32_data_obs.H wemva32_wavelet_15f_3.0ms.H ${B}/RTM.x ${P}/parx32.p
	${B}/RTM.x < $< data_obs=wemva32_data_obs.H wavfile=wemva32_wavelet_15f_3.0ms.H par=${P}/parx32.p > $@

wemva32_dimage_lf_iter%.H: wemva32_sbg_iter%.H wemva32_ddata_lf_iter%.H wemva32_wavelet_5f_3.0ms.H ${B}/RTM.x ${P}/parxlf32.p
	${B}/RTM.x < $< data_obs=wemva32_ddata_lf_iter$*.H wavfile=wemva32_wavelet_5f_3.0ms.H par=${P}/parxlf32.p > $@

### Run Forward WEMVA #########################################################

wemva32_deltaS_iter%.H: wemva32_sbg_iter%.H wemva32_image_iter%.H wemva32_image_flat.H wemva32_data_obs.H wemva32_wavelet_15f_3.0ms.H ${B}/WEMVAI.x ${P}/parx32.p
	${B}/WEMVAI.x < $< image_obs=wemva32_image_iter$*.H image_flat=wemva32_image_flat.H data_obs=wemva32_data_obs.H wavfile=wemva32_wavelet_15f_3.0ms.H par=${P}/parx32.p > $@

### wemva32 (back-scattered) iterate #########################################################

wemva32_sbg_iter1.H: ${D}/wemva32_anom1.H wemva32_sbg_iter0.H wemva32_dimage_lf_iter0.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter0.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter0.H
	${B}/Comb1.x < $< iter=1 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter0.H dm=wemva32_dimage_lf_iter0.H > wemva32_tmp2_iter0.H
	Pad < wemva32_tmp2_iter0.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter0.H
	Add wemva32_tmp3_iter0.H wemva32_sbg_iter0.H > $@
	Rm wemva32_tmp_iter0.H wemva32_tmp2_iter0.H wemva32_tmp3_iter0.H

wemva32_sbg_iter2.H: ${D}/wemva32_anom1.H wemva32_sbg_iter1.H wemva32_dimage_lf_iter1.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter1.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter1.H
	${B}/Comb1.x < $< iter=2 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter1.H dm=wemva32_dimage_lf_iter1.H > wemva32_tmp2_iter1.H
	Pad < wemva32_tmp2_iter1.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter1.H
	Add wemva32_tmp3_iter1.H wemva32_sbg_iter1.H > $@
	Rm wemva32_tmp_iter1.H wemva32_tmp2_iter1.H wemva32_tmp3_iter1.H

wemva32_sbg_iter3.H: ${D}/wemva32_anom1.H wemva32_sbg_iter2.H wemva32_dimage_lf_iter2.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter2.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter2.H
	${B}/Comb1.x < $< iter=3 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter2.H dm=wemva32_dimage_lf_iter2.H > wemva32_tmp2_iter2.H
	Pad < wemva32_tmp2_iter2.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter2.H
	Add wemva32_tmp3_iter2.H wemva32_sbg_iter2.H > $@
	Rm wemva32_tmp_iter2.H wemva32_tmp2_iter2.H wemva32_tmp3_iter2.H

wemva32_sbg_iter4.H: ${D}/wemva32_anom1.H wemva32_sbg_iter3.H wemva32_dimage_lf_iter3.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter3.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter3.H
	${B}/Comb1.x < $< iter=4 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter3.H dm=wemva32_dimage_lf_iter3.H > wemva32_tmp2_iter3.H
	Pad < wemva32_tmp2_iter3.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter3.H
	Add wemva32_tmp3_iter3.H wemva32_sbg_iter3.H > $@
	Rm wemva32_tmp_iter3.H wemva32_tmp2_iter3.H wemva32_tmp3_iter3.H

wemva32_sbg_iter5.H: ${D}/wemva32_anom1.H wemva32_sbg_iter4.H wemva32_dimage_lf_iter4.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter4.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter4.H
	${B}/Comb1.x < $< iter=5 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter4.H dm=wemva32_dimage_lf_iter4.H > wemva32_tmp2_iter4.H
	Pad < wemva32_tmp2_iter4.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter4.H
	Add wemva32_tmp3_iter4.H wemva32_sbg_iter4.H > $@
	Rm wemva32_tmp_iter4.H wemva32_tmp2_iter4.H wemva32_tmp3_iter4.H

wemva32_sbg_iter6.H: ${D}/wemva32_anom1.H wemva32_sbg_iter5.H wemva32_dimage_lf_iter5.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter5.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter5.H
	${B}/Comb1.x < $< iter=6 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter5.H dm=wemva32_dimage_lf_iter5.H > wemva32_tmp2_iter5.H
	Pad < wemva32_tmp2_iter5.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter5.H
	Add wemva32_tmp3_iter5.H wemva32_sbg_iter5.H > $@
	Rm wemva32_tmp_iter5.H wemva32_tmp2_iter5.H wemva32_tmp3_iter5.H

wemva32_sbg_iter7.H: ${D}/wemva32_anom1.H wemva32_sbg_iter6.H wemva32_dimage_lf_iter6.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter6.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter6.H
	${B}/Comb1.x < $< iter=7 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter6.H dm=wemva32_dimage_lf_iter6.H > wemva32_tmp2_iter6.H
	Pad < wemva32_tmp2_iter6.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter6.H
	Add wemva32_tmp3_iter6.H wemva32_sbg_iter6.H > $@
	Rm wemva32_tmp_iter6.H wemva32_tmp2_iter6.H wemva32_tmp3_iter6.H

wemva32_sbg_iter8.H: ${D}/wemva32_anom1.H wemva32_sbg_iter7.H wemva32_dimage_lf_iter7.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter7.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter7.H
	${B}/Comb1.x < $< iter=8 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter7.H dm=wemva32_dimage_lf_iter7.H > wemva32_tmp2_iter7.H
	Pad < wemva32_tmp2_iter7.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter7.H
	Add wemva32_tmp3_iter7.H wemva32_sbg_iter7.H > $@
	Rm wemva32_tmp_iter7.H wemva32_tmp2_iter7.H wemva32_tmp3_iter7.H

wemva32_sbg_iter9.H: ${D}/wemva32_anom1.H wemva32_sbg_iter8.H wemva32_dimage_lf_iter8.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter8.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter8.H
	${B}/Comb1.x < $< iter=9 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter8.H dm=wemva32_dimage_lf_iter8.H > wemva32_tmp2_iter8.H
	Pad < wemva32_tmp2_iter8.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter8.H
	Add wemva32_tmp3_iter8.H wemva32_sbg_iter8.H > $@
	Rm wemva32_tmp_iter8.H wemva32_tmp2_iter8.H wemva32_tmp3_iter8.H

wemva32_sbg_iter10.H: ${D}/wemva32_anom1.H wemva32_sbg_iter9.H wemva32_dimage_lf_iter9.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter9.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter9.H
	${B}/Comb1.x < $< iter=10 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter9.H dm=wemva32_dimage_lf_iter9.H > wemva32_tmp2_iter9.H
	Pad < wemva32_tmp2_iter9.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter9.H
	Add wemva32_tmp3_iter9.H wemva32_sbg_iter9.H > $@
	Rm wemva32_tmp_iter9.H wemva32_tmp2_iter9.H wemva32_tmp3_iter9.H

wemva32_sbg_iter11.H: ${D}/wemva32_anom1.H wemva32_sbg_iter10.H wemva32_dimage_lf_iter10.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter10.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter10.H
	${B}/Comb1.x < $< iter=11 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter10.H dm=wemva32_dimage_lf_iter10.H > wemva32_tmp2_iter10.H
	Pad < wemva32_tmp2_iter10.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter10.H
	Add wemva32_tmp3_iter10.H wemva32_sbg_iter10.H > $@
	Rm wemva32_tmp_iter10.H wemva32_tmp2_iter10.H wemva32_tmp3_iter10.H

wemva32_sbg_iter12.H: ${D}/wemva32_anom1.H wemva32_sbg_iter11.H wemva32_dimage_lf_iter11.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter11.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter11.H
	${B}/Comb1.x < $< iter=12 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter11.H dm=wemva32_dimage_lf_iter11.H > wemva32_tmp2_iter11.H
	Pad < wemva32_tmp2_iter11.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter11.H
	Add wemva32_tmp3_iter11.H wemva32_sbg_iter11.H > $@
	Rm wemva32_tmp_iter11.H wemva32_tmp2_iter11.H wemva32_tmp3_iter11.H

wemva32_sbg_iter13.H: ${D}/wemva32_anom1.H wemva32_sbg_iter12.H wemva32_dimage_lf_iter12.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter12.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter12.H
	${B}/Comb1.x < $< iter=13 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter12.H dm=wemva32_dimage_lf_iter12.H > wemva32_tmp2_iter12.H
	Pad < wemva32_tmp2_iter12.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter12.H
	Add wemva32_tmp3_iter12.H wemva32_sbg_iter12.H > $@
	Rm wemva32_tmp_iter12.H wemva32_tmp2_iter12.H wemva32_tmp3_iter12.H

wemva32_sbg_iter14.H: ${D}/wemva32_anom1.H wemva32_sbg_iter13.H wemva32_dimage_lf_iter13.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter13.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter13.H
	${B}/Comb1.x < $< iter=14 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter13.H dm=wemva32_dimage_lf_iter13.H > wemva32_tmp2_iter13.H
	Pad < wemva32_tmp2_iter13.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter13.H
	Add wemva32_tmp3_iter13.H wemva32_sbg_iter13.H > $@
	Rm wemva32_tmp_iter13.H wemva32_tmp2_iter13.H wemva32_tmp3_iter13.H

wemva32_sbg_iter15.H: ${D}/wemva32_anom1.H wemva32_sbg_iter14.H wemva32_dimage_lf_iter14.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter14.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter14.H
	${B}/Comb1.x < $< iter=15 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter14.H dm=wemva32_dimage_lf_iter14.H > wemva32_tmp2_iter14.H
	Pad < wemva32_tmp2_iter14.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter14.H
	Add wemva32_tmp3_iter14.H wemva32_sbg_iter14.H > $@
	Rm wemva32_tmp_iter14.H wemva32_tmp2_iter14.H wemva32_tmp3_iter14.H

wemva32_sbg_iter16.H: ${D}/wemva32_anom1.H wemva32_sbg_iter15.H wemva32_dimage_lf_iter15.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter15.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter15.H
	${B}/Comb1.x < $< iter=16 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter15.H dm=wemva32_dimage_lf_iter15.H > wemva32_tmp2_iter15.H
	Pad < wemva32_tmp2_iter15.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter15.H
	Add wemva32_tmp3_iter15.H wemva32_sbg_iter15.H > $@
	Rm wemva32_tmp_iter15.H wemva32_tmp2_iter15.H wemva32_tmp3_iter15.H

wemva32_sbg_iter17.H: ${D}/wemva32_anom1.H wemva32_sbg_iter16.H wemva32_dimage_lf_iter16.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter16.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter16.H
	${B}/Comb1.x < $< iter=17 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter16.H dm=wemva32_dimage_lf_iter16.H > wemva32_tmp2_iter16.H
	Pad < wemva32_tmp2_iter16.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter16.H
	Add wemva32_tmp3_iter16.H wemva32_sbg_iter16.H > $@
	Rm wemva32_tmp_iter16.H wemva32_tmp2_iter16.H wemva32_tmp3_iter16.H

wemva32_sbg_iter18.H: ${D}/wemva32_anom1.H wemva32_sbg_iter17.H wemva32_dimage_lf_iter17.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter17.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter17.H
	${B}/Comb1.x < $< iter=18 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter17.H dm=wemva32_dimage_lf_iter17.H > wemva32_tmp2_iter17.H
	Pad < wemva32_tmp2_iter17.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter17.H
	Add wemva32_tmp3_iter17.H wemva32_sbg_iter17.H > $@
	Rm wemva32_tmp_iter17.H wemva32_tmp2_iter17.H wemva32_tmp3_iter17.H

wemva32_sbg_iter19.H: ${D}/wemva32_anom1.H wemva32_sbg_iter18.H wemva32_dimage_lf_iter18.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter18.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter18.H
	${B}/Comb1.x < $< iter=19 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter18.H dm=wemva32_dimage_lf_iter18.H > wemva32_tmp2_iter18.H
	Pad < wemva32_tmp2_iter18.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter18.H
	Add wemva32_tmp3_iter18.H wemva32_sbg_iter18.H > $@
	Rm wemva32_tmp_iter18.H wemva32_tmp2_iter18.H wemva32_tmp3_iter18.H

wemva32_sbg_iter20.H: ${D}/wemva32_anom1.H wemva32_sbg_iter19.H wemva32_dimage_lf_iter19.H ${D}/wemva32_a1.H ${B}/Comb1.x
	Window3d < wemva32_sbg_iter19.H f1=200 f2=200 n1=121 n2=201 > wemva32_tmp_iter19.H
	${B}/Comb1.x < $< iter=20 a1=${D}/wemva32_a1.H m0=wemva32_tmp_iter19.H dm=wemva32_dimage_lf_iter19.H > wemva32_tmp2_iter19.H
	Pad < wemva32_tmp2_iter19.H beg1=200 end1=200 beg2=200 end2=200 > wemva32_tmp3_iter19.H
	Add wemva32_tmp3_iter19.H wemva32_sbg_iter19.H > $@
	Rm wemva32_tmp_iter19.H wemva32_tmp2_iter19.H wemva32_tmp3_iter19.H


### Output Results #########################################################

# Create ps, pdf and gif files
${R}/%.eps: ${R}/%.v
	pstexpen $< $@ color=y fat=1 fatmult=1.5 invras=n force=

${R}/%.pdf: ${R}/%.eps
	epstopdf $< --nocompress

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

wemva261= color=j newclip=1 bpclip=0.1 epclip=99 wantscalebar=y label1="Depth(km)" label2="Location(km)" label3="Slowness(s/km)" d1num=0.5 d2num=1 
wemva262= color=j newclip=1 bpclip=1 epclip=99.5 wantscalebar=y label1="kz(1/km)" label2="kx(1/km)" label3="Amplitude" min1=-10 min2=-10 max1=10 max2=10
wemva263= color=j newclip=1 bpclip=1 epclip=99.5 wantscalebar=y label1="kz(1/km)" label2="kx(1/km)" label3="Amplitude" min1=-10 min2=-10 max1=10 max2=10
wemva264= label1="Location(km)" label2="Slowness(s/km)"
wemva265= color=j newclip=1 bpclip=0 epclip=100 wantscalebar=y label1="Depth(km)" label2="Location(km)" label3="Slowness(s/km)" d1num=0.5 d2num=1 

${R}/wemva26deltaS.v: wemva26_deltaS_anom7_scaled.H 
	Grey < $< ${wemva261} title=' ' >/dev/null out=$@

${R}/wemva26deltaSfk.v: wemva26_deltaS_anom7_fk.H 
	Grey < $< ${wemva262} title=' ' >/dev/null out=$@

${R}/wemva27imaged.v: wemva26_imaged_bg_anom2_scaled.H 
	Grey < $< ${wemva261} title=' ' >/dev/null out=$@

${R}/wemva27imagedfk.v: wemva26_imaged_bg_anom2_fk.H 
	Grey < $< ${wemva262} title=' ' >/dev/null out=$@

${R}/wemva27imagedlf.v: wemva26_imaged_lf_bg_anom2_scaled.H 
	Grey < $< ${wemva261} title=' ' >/dev/null out=$@

${R}/wemva27imagedlffk.v: wemva26_imaged_lf_bg_anom2_fk.H 
	Grey < $< ${wemva262} title=' ' >/dev/null out=$@

${R}/wemva27deltaSmfk0.v: wemva26_deltaS_mfk0.H 
	Grey < $< ${wemva263} title=' ' >/dev/null out=$@

${R}/wemva27deltaSmfk.v: wemva26_deltaS_mfk.H 
	Grey < $< ${wemva263} title=' ' >/dev/null out=$@

${R}/wemva27stacks.v: wemva26_imaged_bg_anom2_stack.H wemva26_imaged_lf_bg_anom2_stack.H wemva26_imaged_lf_bg_anom2_sstack.H
	Graph < wemva26_imaged_bg_anom2_stack.H ${wemva264} title='(a)' >/dev/null out=${R}/1.v
	Graph < wemva26_imaged_lf_bg_anom2_stack.H ${wemva264} title='(b)' >/dev/null out=${R}/2.v
	Graph < wemva26_imaged_lf_bg_anom2_sstack.H ${wemva264} title='(c)' >/dev/null out=${R}/3.v
	vp_SideBySideIso ${R}/1.v ${R}/2.v ${R}/3.v > $@
	Rm ${R}/1.v ${R}/2.v ${R}/3.v

${R}/wemva30tv.v: ${D}/wemva30_anom1.H
	Window3d < $< min1=0 min2=0 n1=109 n2=201 | Grey ${wemva265} title=' ' >/dev/null out=$@

${R}/wemva30bp.v: wemva32_sbg_iter20.H
	Window3d < $< min1=0 min2=0 n1=109 n2=201 | Grey ${wemva265} title=' ' >/dev/null out=$@

${R}/wemva30fp.v: wemva31_sbg_iter20.H
	Window3d < $< min1=0 min2=0 n1=109 n2=201 | Grey ${wemva265} title=' ' >/dev/null out=$@

${R}/wemva30fbp.v: wemva30_sbg_iter20.H
	Window3d < $< min1=0 min2=0 n1=109 n2=201 | Grey ${wemva265} title=' ' >/dev/null out=$@

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

results: $(B)/Forward_born.x $(B)/RTM.x $(B)/WEMVA.x $(B)/WEMVAA.x ${B}/Refl0.x ${B}/Refl3.x wemva26_anom0_bg_pad.H wemva26_anom2_bg_pad.H wemva26_anom0_rf_pad.H wemva26_anom2_drf_pad.H wemva26_data_anom2.H wemva26_datad_anom2.H wemva26_datad_lf_anom2.H wemva26_imaged_bg_anom2.H wemva26_imaged_lf_bg_anom2.H wemva26_imaged_bg_anom2_scaled.H wemva26_imaged_lf_bg_anom2_scaled.H wemva26_imaged_bg_anom2_fk.H wemva26_imaged_lf_bg_anom2_fk.H wemva26_imaged_bg_anom2_stack.H wemva26_imaged_lf_bg_anom2_stack.H wemva26_imaged_lf_bg_anom2_sstack.H wemva26_deltaI_anom2.H wemva26_deltaI_anom7.H wemva26_deltaS_anom7.H wemva26_deltaS_anom7_scaled.H wemva26_deltaS_anom7_fk.H wemva26_deltaS_mfk0.H wemva26_deltaS_mfk.H ${R}/wemva26deltaS.eps ${R}/wemva26deltaSfk.eps ${R}/wemva27imaged.eps ${R}/wemva27imagedfk.eps ${R}/wemva27imagedlf.eps ${R}/wemva27imagedlffk.eps ${R}/wemva27deltaSmfk0.eps ${R}/wemva27deltaSmfk.eps ${R}/wemva27stacks.eps $(B)/Forward.x $(B)/Forward2.x $(B)/Refl.x $(B)/Refl2.x $(B)/WEMVAI.x $(B)/Comb1.x $(B)/Comb2.x wemva30_wavelet_15f_3.0ms.H wemva30_wavelet_5f_3.0ms.H wemva30_strue.H wemva30_rf.H wemva30_data_obs.H wemva30_image_flat.H wemva30_data_lf_obs.H wemva30_sbg_iter0.H wemva30_data_iter0.H wemva30_image_iter0.H wemva30_deltaS_iter0.H wemva30_data_lf_iter0.H wemva30_ddata_lf_iter0.H wemva30_dimage_lf_iter0.H wemva30_sbg_iter1.H wemva30_data_iter1.H wemva30_image_iter1.H wemva30_deltaS_iter1.H wemva30_data_lf_iter1.H wemva30_ddata_lf_iter1.H wemva30_dimage_lf_iter1.H wemva30_sbg_iter2.H wemva30_data_iter2.H wemva30_image_iter2.H wemva30_deltaS_iter2.H wemva30_data_lf_iter2.H wemva30_ddata_lf_iter2.H wemva30_dimage_lf_iter2.H wemva30_sbg_iter3.H wemva30_data_iter3.H wemva30_image_iter3.H wemva30_deltaS_iter3.H wemva30_data_lf_iter3.H wemva30_ddata_lf_iter3.H wemva30_dimage_lf_iter3.H wemva30_sbg_iter4.H wemva30_data_iter4.H wemva30_image_iter4.H wemva30_deltaS_iter4.H wemva30_data_lf_iter4.H wemva30_ddata_lf_iter4.H wemva30_dimage_lf_iter4.H wemva30_sbg_iter5.H wemva30_data_iter5.H wemva30_image_iter5.H wemva30_deltaS_iter5.H wemva30_data_lf_iter5.H wemva30_ddata_lf_iter5.H wemva30_dimage_lf_iter5.H wemva30_sbg_iter6.H wemva30_data_iter6.H wemva30_image_iter6.H wemva30_deltaS_iter6.H wemva30_data_lf_iter6.H wemva30_ddata_lf_iter6.H wemva30_dimage_lf_iter6.H wemva30_sbg_iter7.H wemva30_data_iter7.H wemva30_image_iter7.H wemva30_deltaS_iter7.H wemva30_data_lf_iter7.H wemva30_ddata_lf_iter7.H wemva30_dimage_lf_iter7.H wemva30_sbg_iter8.H wemva30_data_iter8.H wemva30_image_iter8.H wemva30_deltaS_iter8.H wemva30_data_lf_iter8.H wemva30_ddata_lf_iter8.H wemva30_dimage_lf_iter8.H wemva30_sbg_iter9.H wemva30_data_iter9.H wemva30_image_iter9.H wemva30_deltaS_iter9.H wemva30_data_lf_iter9.H wemva30_ddata_lf_iter9.H wemva30_dimage_lf_iter9.H wemva30_sbg_iter10.H wemva30_data_iter10.H wemva30_image_iter10.H wemva30_deltaS_iter10.H wemva30_data_lf_iter10.H wemva30_ddata_lf_iter10.H wemva30_dimage_lf_iter10.H wemva30_sbg_iter11.H wemva30_data_iter11.H wemva30_image_iter11.H wemva30_deltaS_iter11.H wemva30_data_lf_iter11.H wemva30_ddata_lf_iter11.H wemva30_dimage_lf_iter11.H wemva30_sbg_iter12.H wemva30_data_iter12.H wemva30_image_iter12.H wemva30_deltaS_iter12.H wemva30_data_lf_iter12.H wemva30_ddata_lf_iter12.H wemva30_dimage_lf_iter12.H wemva30_sbg_iter13.H wemva30_data_iter13.H wemva30_image_iter13.H wemva30_deltaS_iter13.H wemva30_data_lf_iter13.H wemva30_ddata_lf_iter13.H wemva30_dimage_lf_iter13.H wemva30_sbg_iter14.H wemva30_data_iter14.H wemva30_image_iter14.H wemva30_deltaS_iter14.H wemva30_data_lf_iter14.H wemva30_ddata_lf_iter14.H wemva30_dimage_lf_iter14.H wemva30_sbg_iter15.H wemva30_data_iter15.H wemva30_image_iter15.H wemva30_deltaS_iter15.H wemva30_data_lf_iter15.H wemva30_ddata_lf_iter15.H wemva30_dimage_lf_iter15.H wemva30_sbg_iter16.H wemva30_data_iter16.H wemva30_image_iter16.H wemva30_deltaS_iter16.H wemva30_data_lf_iter16.H wemva30_ddata_lf_iter16.H wemva30_dimage_lf_iter16.H wemva30_sbg_iter17.H wemva30_data_iter17.H wemva30_image_iter17.H wemva30_deltaS_iter17.H wemva30_data_lf_iter17.H wemva30_ddata_lf_iter17.H wemva30_dimage_lf_iter17.H wemva30_sbg_iter18.H wemva30_data_iter18.H wemva30_image_iter18.H wemva30_deltaS_iter18.H wemva30_data_lf_iter18.H wemva30_ddata_lf_iter18.H wemva30_dimage_lf_iter18.H wemva30_sbg_iter19.H wemva30_data_iter19.H wemva30_image_iter19.H wemva30_deltaS_iter19.H wemva30_data_lf_iter19.H wemva30_ddata_lf_iter19.H wemva30_dimage_lf_iter19.H wemva30_sbg_iter20.H ${R}/wemva30tv.eps ${R}/wemva30fbp.eps wemva31_wavelet_15f_3.0ms.H wemva31_strue.H wemva31_rf.H wemva31_data_obs.H wemva31_image_flat.H wemva31_sbg_iter0.H wemva31_data_iter0.H wemva31_image_iter0.H wemva31_deltaS_iter0.H wemva31_sbg_iter1.H wemva31_data_iter1.H wemva31_image_iter1.H wemva31_deltaS_iter1.H wemva31_sbg_iter2.H wemva31_data_iter2.H wemva31_image_iter2.H wemva31_deltaS_iter2.H wemva31_sbg_iter3.H wemva31_data_iter3.H wemva31_image_iter3.H wemva31_deltaS_iter3.H wemva31_sbg_iter4.H wemva31_data_iter4.H wemva31_image_iter4.H wemva31_deltaS_iter4.H wemva31_sbg_iter5.H wemva31_data_iter5.H wemva31_image_iter5.H wemva31_deltaS_iter5.H wemva31_sbg_iter6.H wemva31_data_iter6.H wemva31_image_iter6.H wemva31_deltaS_iter6.H wemva31_sbg_iter7.H wemva31_data_iter7.H wemva31_image_iter7.H wemva31_deltaS_iter7.H wemva31_sbg_iter8.H wemva31_data_iter8.H wemva31_image_iter8.H wemva31_deltaS_iter8.H wemva31_sbg_iter9.H wemva31_data_iter9.H wemva31_image_iter9.H wemva31_deltaS_iter9.H wemva31_sbg_iter10.H wemva31_data_iter10.H wemva31_image_iter10.H wemva31_deltaS_iter10.H wemva31_sbg_iter11.H wemva31_data_iter11.H wemva31_image_iter11.H wemva31_deltaS_iter11.H wemva31_sbg_iter12.H wemva31_data_iter12.H wemva31_image_iter12.H wemva31_deltaS_iter12.H wemva31_sbg_iter13.H wemva31_data_iter13.H wemva31_image_iter13.H wemva31_deltaS_iter13.H wemva31_sbg_iter14.H wemva31_data_iter14.H wemva31_image_iter14.H wemva31_deltaS_iter14.H wemva31_sbg_iter15.H wemva31_data_iter15.H wemva31_image_iter15.H wemva31_deltaS_iter15.H wemva31_sbg_iter16.H wemva31_data_iter16.H wemva31_image_iter16.H wemva31_deltaS_iter16.H wemva31_sbg_iter17.H wemva31_data_iter17.H wemva31_image_iter17.H wemva31_deltaS_iter17.H wemva31_sbg_iter18.H wemva31_data_iter18.H wemva31_image_iter18.H wemva31_deltaS_iter18.H wemva31_sbg_iter19.H wemva31_data_iter19.H wemva31_image_iter19.H wemva31_deltaS_iter19.H wemva31_sbg_iter20.H ${R}/wemva30fp.eps wemva32_wavelet_5f_3.0ms.H wemva32_strue.H wemva32_data_lf_obs.H wemva32_sbg_iter0.H wemva32_data_lf_iter0.H wemva32_ddata_lf_iter0.H wemva32_dimage_lf_iter0.H wemva32_sbg_iter1.H wemva32_data_lf_iter1.H wemva32_ddata_lf_iter1.H wemva32_dimage_lf_iter1.H wemva32_sbg_iter2.H wemva32_data_lf_iter2.H wemva32_ddata_lf_iter2.H wemva32_dimage_lf_iter2.H wemva32_sbg_iter3.H wemva32_data_lf_iter3.H wemva32_ddata_lf_iter3.H wemva32_dimage_lf_iter3.H wemva32_sbg_iter4.H wemva32_data_lf_iter4.H wemva32_ddata_lf_iter4.H wemva32_dimage_lf_iter4.H wemva32_sbg_iter5.H wemva32_data_lf_iter5.H wemva32_ddata_lf_iter5.H wemva32_dimage_lf_iter5.H wemva32_sbg_iter6.H wemva32_data_lf_iter6.H wemva32_ddata_lf_iter6.H wemva32_dimage_lf_iter6.H wemva32_sbg_iter7.H wemva32_data_lf_iter7.H wemva32_ddata_lf_iter7.H wemva32_dimage_lf_iter7.H wemva32_sbg_iter8.H wemva32_data_lf_iter8.H wemva32_ddata_lf_iter8.H wemva32_dimage_lf_iter8.H wemva32_sbg_iter9.H wemva32_data_lf_iter9.H wemva32_ddata_lf_iter9.H wemva32_dimage_lf_iter9.H wemva32_sbg_iter10.H wemva32_data_lf_iter10.H wemva32_ddata_lf_iter10.H wemva32_dimage_lf_iter10.H wemva32_sbg_iter11.H wemva32_data_lf_iter11.H wemva32_ddata_lf_iter11.H wemva32_dimage_lf_iter11.H wemva32_sbg_iter12.H wemva32_data_lf_iter12.H wemva32_ddata_lf_iter12.H wemva32_dimage_lf_iter12.H wemva32_sbg_iter13.H wemva32_data_lf_iter13.H wemva32_ddata_lf_iter13.H wemva32_dimage_lf_iter13.H wemva32_sbg_iter14.H wemva32_data_lf_iter14.H wemva32_ddata_lf_iter14.H wemva32_dimage_lf_iter14.H wemva32_sbg_iter15.H wemva32_data_lf_iter15.H wemva32_ddata_lf_iter15.H wemva32_dimage_lf_iter15.H wemva32_sbg_iter16.H wemva32_data_lf_iter16.H wemva32_ddata_lf_iter16.H wemva32_dimage_lf_iter16.H wemva32_sbg_iter17.H wemva32_data_lf_iter17.H wemva32_ddata_lf_iter17.H wemva32_dimage_lf_iter17.H wemva32_sbg_iter18.H wemva32_data_lf_iter18.H wemva32_ddata_lf_iter18.H wemva32_dimage_lf_iter18.H wemva32_sbg_iter19.H wemva32_data_lf_iter19.H wemva32_ddata_lf_iter19.H wemva32_dimage_lf_iter19.H wemva32_sbg_iter20.H ${R}/wemva30bp.eps

clean:
	rm -f $(R)/*.v
	rm -f $(R)/*.eps

burn:
	rm -f *.H
	rm -f $(R)/*.v
	rm -f $(R)/*.eps
	rm -f $(M)/*
	rm -f $(O)/*
	rm -f $(B)/*
