# RTMVEL1akefile
# RTM imaging for velocity model 1
# for vel1.H file

# Figures
RTMVEL1RESULTSER = rtm1.vel1.mult.fullmuted rtm1.vel1.mult.fullmuted.lowcut \
		rtm1.vel1.mult.fullmuted.sousnap1 rtm1.vel1.mult.fullmuted.recsnap1 rtm1.vel1.mult.fullmuted.reflsnap1 \
		rtm1.vel1.mult.fullmuted.sousnap2 rtm1.vel1.mult.fullmuted.recsnap2 rtm1.vel1.mult.fullmuted.reflsnap2 \
		rtm1.vel1.mult.fullmuted.sousnap3 rtm1.vel1.mult.fullmuted.recsnap3 rtm1.vel1.mult.fullmuted.reflsnap3
RTMVEL1RESULTSCR = rtm0.vel1.mult.fullmuted rtm0.vel1.mult.fullmuted.lowcut
RTMVEL1RESULTSNR = 
RTMVEL1FIG_pdf = rtm1.vel1.mult.fullmuted.pdf rtm1.vel1.mult.fullmuted.lowcut.pdf \
	rtm1.vel1.mult.fullmuted.sousnap1.pdf rtm1.vel1.mult.fullmuted.recsnap1.pdf rtm1.vel1.mult.fullmuted.reflsnap1.pdf \
	rtm1.vel1.mult.fullmuted.sousnap2.pdf rtm1.vel1.mult.fullmuted.recsnap2.pdf rtm1.vel1.mult.fullmuted.reflsnap2.pdf \
	rtm1.vel1.mult.fullmuted.sousnap3.pdf rtm1.vel1.mult.fullmuted.recsnap3.pdf rtm1.vel1.mult.fullmuted.reflsnap3.pdf \
		 rtm0.vel1.mult.fullmuted.pdf rtm0.vel1.mult.fullmuted.lowcut.pdf

##########################
### generate data
## sh model2d.sh '$(m)' 'outputfile' 'modelprogram' 'velfile' 'soufile' 'parfile' '$(abstop)' 'jsx=4' 'sousnapoutfile' 'tmpdir=.'
##########################
## By default I will use $P/shot0.P file for the whole shots along 2D line.
##########
# with SRM
#wvfld%.mult.H:  vel%.H sou.HH ${P}/shot%.P ${B}/model2d_program.x 
#	<$< ${B}/model2d_program.x par=${P}/shot$*.P abstop=0 source=sou.HH $(dn) snap=$@
data%.vel1.mult.H souwvfld%.vel1.mult.H: vel1.H sou.H $P/shot%.vel1.P $B/model2d_program.x $(SH)/model2d.sh
	sh $(SH)/model2d.sh '$*.vel1.mult' 'data$*.vel1.mult.H' '$B/model2d_program.x' 'vel1.H' 'sou.H' '$P/shot$*.vel1.P' '0' '6' 'souwvfld$*.vel1.mult.H' './temp'
data%.vel1.mult: data%.vel1.mult.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='data$* with SRM' $(Axis) label1="Time (s)" | Tube &
souwvfld%.vel1.mult: souwvfld%.vel1.mult.H
	<$< Grey title='souwvfld$* with SRM' $(Axis) | Tube &

##########
# without SRM
#wvfld%.nomult.H:  vel%.H sou.HH ${P}/shot%.P ${B}/model2d_program.x 
#	<$< ${B}/model2d_program.x par=${P}/shot$*.P abstop=1 source=sou.HH $(dn) snap=$@
data%.vel1.nomult.H souwvfld%.vel1.nomult.H:  vel1.H sou.H $P/shot%.vel1.P $B/model2d_program.x $(SH)/model2d.sh
	sh $(SH)/model2d.sh '$*.vel1.nomult' 'data$*.vel1.nomult.H' '$B/model2d_program.x' 'vel1.H' 'sou.H' '$P/shot$*.vel1.P' '1' '6' 'souwvfld$*.vel1.nomult.H' './temp'
data%.vel1.nomult: data%.vel1.nomult.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='data$* without SRM' $(Axis) label1="Time (s)" | Tube &
souwvfld%.vel1.nomult: souwvfld%.vel1.nomult.H
	<$< Grey title='souwvfld$* without SRM' $(Axis) | Tube &

#################
## muting data
## sh mute2d.sh '$(m)' 'inputfile' 'outputfile' 'parfile' 'jsx=8' 'tmpdir=.'
#################
# with SRM
data%.vel1.mult.muted.H: data%.vel1.mult.H $P/mute.vel1.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > mute$*.vel1.mult.in.H
	sh $(SH)/mute2d.sh '$*.vel1.mult.muted' 'mute$*.vel1.mult.in.H' 'mute$*.vel1.mult.out.H' '$P/mute.vel1.P' '8' './temp'
	< mute$*.vel1.mult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm mute$*.vel1.mult.in.H mute$*.vel1.mult.out.H
data%.vel1.mult.muted: data%.vel1.mult.muted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='muted data$* with SRM' $(Axis) label1="Time (s)" | Tube &
##full muted
data%.vel1.mult.fullmuted.H: data%.vel1.mult.H $P/fullmute.vel1.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > fullmute$*.vel1.mult.in.H
	sh $(SH)/mute2d.sh '$*.vel1.mult.fullmuted' 'fullmute$*.vel1.mult.in.H' 'fullmute$*.vel1.mult.out.H' '$P/fullmute.vel1.P' '8' './temp'
	< fullmute$*.vel1.mult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm fullmute$*.vel1.mult.in.H fullmute$*.vel1.mult.out.H
data%.vel1.mult.fullmuted: data%.vel1.mult.fullmuted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='fullmuted data$* with SRM' $(Axis) label1="Time (s)" | Tube &

# without SRM
data%.vel1.nomult.muted.H: data%.vel1.nomult.H $P/mute.vel1.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > mute$*.vel1.nomult.in.H
	sh $(SH)/mute2d.sh '$*.vel1.nomult.muted' 'mute$*.vel1.nomult.in.H' 'mute$*.vel1.nomult.out.H' '$P/mute.vel1.P' '8' './temp'
	< mute$*.vel1.nomult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm mute$*.vel1.nomult.in.H mute$*.vel1.nomult.out.H
data%.vel1.nomult.muted: data%.vel1.nomult.muted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='muted data$* without SRM' $(Axis) label1="Time (s)" | Tube &
##full muted
data%.vel1.nomult.fullmuted.H: data%.vel1.nomult.H $P/fullmute.vel1.P $(SH)/mute2d.sh
	<$< Transp plane=13 |Transp plane=34 > fullmute$*.vel1.nomult.in.H
	sh $(SH)/mute2d.sh '$*.vel1.nomult.fullmuted' 'fullmute$*.vel1.nomult.in.H' 'fullmute$*.vel1.nomult.out.H' '$P/fullmute.vel1.P' '8' './temp'
	< fullmute$*.vel1.nomult.out.H Transp plane=34 | Transp plane=13 > $@
	Rm fullmute$*.vel1.nomult.in.H fullmute$*.vel1.nomult.out.H
data%.vel1.nomult.fullmuted: data%.vel1.nomult.fullmuted.H
	<$< Transp plane=13 | Transp plane=34 | Grey title='fullmuted data$* with SRM' $(Axis) label1="Time (s)" | Tube &

##################################
### RTM images
##       sh rtm2d.sh '$(m)' 'in.H' 'out.H' 'rtmprogram' 'vel.mig.H' 'soufile' 'parfile' 'jsx=4' 'tmpdir=.' 'sousnapfile.out=sousnap${m}.H' 'recsnapfile.out=recsnap${m}.H' 'reflsnapfile.out=reflsnap${m}.H'
##################################
#########
## SRM
#########
## from unmuted data
rtm%.vel1.mult.H rtm%.vel1.mult.sousnap.H rtm%.vel1.mult.recsnap.H rtm%.vel1.mult.reflsnap.H: data%.vel1.mult.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.mult.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.mult' 'data$*.vel1.mult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.mult.sousnap.H' 'rtm$*.vel1.mult.recsnap.H' 'rtm$*.vel1.mult.reflsnap.H'
	Rm data$*.vel1.mult.interp.H
rtm%.vel1.mult.v: rtm%.vel1.mult.H
	<$< Grey title="An RTM image from data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.mult: rtm%.vel1.mult.H
	Sep_cube rtm$*.vel1.mult.sousnap.H rtm$*.vel1.mult.recsnap.H rtm$*.vel1.mult.reflsnap.H nviews=3 &

## from muted data
rtm%.vel1.mult.muted.H rtm%.vel1.mult.muted.sousnap.H rtm%.vel1.mult.muted.recsnap.H rtm%.vel1.mult.muted.reflsnap.H: data%.vel1.mult.muted.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.mult.muted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.mult.muted' 'data$*.vel1.mult.muted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.mult.muted.sousnap.H' 'rtm$*.vel1.mult.muted.recsnap.H' 'rtm$*.vel1.mult.muted.reflsnap.H'
	Rm data$*.vel1.mult.muted.interp.H
rtm%.vel1.mult.muted.v: rtm%.vel1.mult.muted.H
	<$< Grey title="An RTM image from muted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.mult.muted: rtm%.vel1.mult.muted.H
	Sep_cube rtm$*.vel1.mult.muted.sousnap.H rtm$*.vel1.mult.muted.recsnap.H rtm$*.vel1.mult.muted.reflsnap.H nviews=3 &
#full muted
rtm%.vel1.mult.fullmuted.H rtm%.vel1.mult.fullmuted.sousnap.H rtm%.vel1.mult.fullmuted.recsnap.H rtm%.vel1.mult.fullmuted.reflsnap.H: data%.vel1.mult.fullmuted.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.mult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.mult.fullmuted' 'data$*.vel1.mult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.mult.fullmuted.sousnap.H' 'rtm$*.vel1.mult.fullmuted.recsnap.H' 'rtm$*.vel1.mult.fullmuted.reflsnap.H'
	Rm data$*.vel1.mult.fullmuted.interp.H
rtm%.vel1.mult.fullmuted.v: rtm%.vel1.mult.fullmuted.H
	<$< Grey title="An RTM image from fullmuted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.mult.fullmuted: rtm%.vel1.mult.fullmuted.H
	Sep_cube rtm$*.vel1.mult.fullmuted.sousnap.H rtm$*.vel1.mult.fullmuted.recsnap.H rtm$*.vel1.mult.fullmuted.reflsnap.H nviews=3 &

rtm%.vel1.mult.fullmuted.spectrum.v: rtm%.vel1.mult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$@

#########
## NO SRM
#########
## from unmuted data
rtm%.vel1.nomult.H rtm%.vel1.nomult.sousnap.H rtm%.vel1.nomult.recsnap.H rtm%.vel1.nomult.reflsnap.H: data%.vel1.nomult.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.nomult.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.nomult' 'data$*.vel1.nomult.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.nomult.sousnap.H' 'rtm$*.vel1.nomult.recsnap.H' 'rtm$*.vel1.nomult.reflsnap.H'
	Rm data$*.vel1.nomult.interp.H
rtm%.vel1.nomult.v: rtm%.vel1.nomult.H
	<$< Grey title="An RTM image from data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.nomult: rtm%.vel1.nomult.H
	Sep_cube rtm$*.vel1.nomult.sousnap.H rtm$*.vel1.nomult.recsnap.H rtm$*.vel1.nomult.reflsnap.H nviews=3 &

## from muted data
rtm%.vel1.nomult.muted.H rtm%.vel1.nomult.muted.sousnap.H rtm%.vel1.nomult.muted.recsnap.H rtm%.vel1.nomult.muted.reflsnap.H: data%.vel1.nomult.muted.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.nomult.muted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.nomult.muted' 'data$*.vel1.nomult.muted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.nomult.muted.sousnap.H' 'rtm$*.vel1.nomult.muted.recsnap.H' 'rtm$*.vel1.nomult.muted.reflsnap.H'
	Rm data$*.vel1.nomult.muted.interp.H
rtm%.vel1.nomult.muted.v: rtm%.vel1.nomult.muted.H
	<$< Grey title="An RTM image from muted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.nomult.muted: rtm%.vel1.nomult.muted.H
	Sep_cube rtm$*.vel1.nomult.muted.sousnap.H rtm$*.vel1.nomult.muted.recsnap.H rtm$*.vel1.nomult.muted.reflsnap.H nviews=3 &
#full muted
rtm%.vel1.nomult.fullmuted.H rtm%.vel1.nomult.fullmuted.sousnap.H rtm%.vel1.nomult.fullmuted.recsnap.H rtm%.vel1.nomult.fullmuted.reflsnap.H: data%.vel1.nomult.fullmuted.H ${B}/crtm2d_program.x vel1.mig.H sou.H ${P}/rtm.vel1.P $(SH)/rtm2d.sh
	<$< Interp d2out=10 maxsize=1000 > data$*.vel1.nomult.fullmuted.interp.H
	sh ${SH}/rtm2d.sh '$*.vel1.nomult.fullmuted' 'data$*.vel1.nomult.fullmuted.interp.H' '$@' '${B}/crtm2d_program.x' 'vel1.mig.H' 'sou.H' '${P}/rtm.vel1.P' '6' './temp' 'rtm$*.vel1.nomult.fullmuted.sousnap.H' 'rtm$*.vel1.nomult.fullmuted.recsnap.H' 'rtm$*.vel1.nomult.fullmuted.reflsnap.H'
	Rm data$*.vel1.nomult.fullmuted.interp.H
rtm%.vel1.nomult.fullmuted.v: rtm%.vel1.nomult.fullmuted.H
	<$< Grey title="An RTM image from fullmuted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@
rtm%.vel1.nomult.fullmuted: rtm%.vel1.nomult.fullmuted.H
	Sep_cube rtm$*.vel1.nomult.fullmuted.sousnap.H rtm$*.vel1.nomult.fullmuted.recsnap.H rtm$*.vel1.nomult.fullmuted.reflsnap.H nviews=3 &

rtm%.vel1.nomult.fullmuted.spectrum.v: rtm%.vel1.nomult.fullmuted.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of an RTM image from fullmuted data without SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$

#################################
### RTM after lowcut filtered
#################################
## from unmuted data
rtm%.vel1.mult.lowcut.H: rtm%.vel1.mult.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel1.mult.lowcut.v: rtm%.vel1.mult.lowcut.H
	<$< Grey title="An lowcut-filtered RTM image from data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@

## from muted data
rtm%.vel1.mult.muted.lowcut.H: rtm%.vel1.mult.muted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel1.mult.muted.lowcut.v: rtm%.vel1.mult.muted.lowcut.H
	<$< Grey title="An lowcut-filtered RTM image from muted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@

## from fullmuted data
rtm%.vel1.mult.fullmuted.lowcut.H: rtm%.vel1.mult.fullmuted.H $B/Light1.x
	<$< $B/Light1.x >$@
rtm%.vel1.mult.fullmuted.lowcut.v: rtm%.vel1.mult.fullmuted.lowcut.H
	<$< Window n3=1 f3=0 | Grey title="An lowcut-filtered RTM image from muted data with SRM" $(Axis) screenratio=$(vel1screenratio) | Vppen vpstyle=n scale=$(vel1screenratio) $(dn) out=$R/$@

rtm%.vel1.mult.fullmuted.lowcut.spectrum.v: rtm%.vel1.mult.fullmuted.lowcut.H $B/Fftw2D_program.x
	<$< Rtoc | $B/Fftw2D_program.x flag=1 dim=2 axis1=1 axis2=1 | Cabs | Grey title="The spectrum of a lowcut-filtered RTM image from fullmuted data with SRM" $(Axis) label1="Kz (/m)" label2="Kx (/m)" $(dn) out=$R/$


################################
## get the image showing origins of RTM noises from rtm1.vel1 where sx=3500
################################
rtm1.vel1.mult.fullmuted.imagingsnaps: rtm1.vel1.mult.fullmuted.sousnap1.pdf rtm1.vel1.mult.fullmuted.recsnap1.pdf rtm1.vel1.mult.fullmuted.reflsnap1.pdf \
	rtm1.vel1.mult.fullmuted.sousnap2.pdf rtm1.vel1.mult.fullmuted.recsnap2.pdf rtm1.vel1.mult.fullmuted.reflsnap2.pdf \
	rtm1.vel1.mult.fullmuted.sousnap3.pdf rtm1.vel1.mult.fullmuted.recsnap3.pdf rtm1.vel1.mult.fullmuted.reflsnap3.pdf

## % must be sou, rec or refl
rtm1.vel1.mult.fullmuted.%snap1.v: rtm1.vel1.mult.fullmuted.H ${TXTDIR}/text-imagingsnap1-vel1-mult
	< rtm1.vel1.mult.fullmuted.$*snap.H Window n3=1 f3=12 | Grey $(Axis) screenratio=$(vel1screenratio) $(dn) out=$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap1.v
	vp_annotate <$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap1.v batch=y vpstyle=n scale=$(vel1screenratio) text=${TXTDIR}/text-imagingsnap1-vel1-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap1.v
rtm1.vel1.mult.fullmuted.%snap2.v: rtm1.vel1.mult.fullmuted.H ${TXTDIR}/text-imagingsnap2-vel1-mult
	< rtm1.vel1.mult.fullmuted.$*snap.H Window n3=1 f3=8 | Grey $(Axis) screenratio=$(vel1screenratio) $(dn) out=$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap2.v
	vp_annotate <$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap2.v batch=y vpstyle=n scale=$(vel1screenratio) text=${TXTDIR}/text-imagingsnap2-vel1-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap2.v
rtm1.vel1.mult.fullmuted.%snap3.v: rtm1.vel1.mult.fullmuted.H ${TXTDIR}/text-imagingsnap3-vel1-mult
	< rtm1.vel1.mult.fullmuted.$*snap.H Window n3=1 f3=4 | Grey $(Axis) screenratio=$(vel1screenratio) $(dn) out=$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap3.v
	vp_annotate <$R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap3.v batch=y vpstyle=n scale=$(vel1screenratio) text=${TXTDIR}/text-imagingsnap3-vel1-mult >$R/$@
	rm -f $R/tmp1x.rtm1.vel1.mult.fullmuted.$*snap3.v

