#define LINDIR ./../..
#define SLDDIR LINDIR/Slds/Resmig4
#define CHPDIR LINDIR/Resmig4
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define SLDLIST vprmig velsyn datspk imgspk datsyn imgsyn
#define USRLIBS VPLOTLIB
#define TLSZFAT  titlefat=5 titlesz=12 labelfat=3 labelsz=8 
#define THICK plotcol=0,0,0,0,0,0,0 plotfat=8
#define THIN plotcol=5,5,5,5,5,5,5 plotfat=4
#define FAST title=" " label1=" " label2=" " wantframe=y wantaxis=n

#include <SEP.defs>

SLDDIR/vprmig.v: vprmig.x
	vprmig.x out=junk.v > Junk
	vppen junk.v erase=n vpstyle=n > SLDDIR/vprmig.v
	rm junk.v Junk

SLDDIR/velsyn.v : DATDIR/strvel_syn.H DATDIR/str_syn.H
        Taplot < DATDIR/strvel_syn.H gpow=1 bias=2.7 clip=1. | Ta2vplot title="Velocity model" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 out=junk.v > Junk
	vp_Overlay back.v junk.v > junka.v
	Bar n2=200 type=h minval=1.89 maxval=3.01 > tempbar.H
	Taplot < tempbar.H gpow=1 bias=2.7 clip=1. | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.5 o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junkb.v > Junk
	vppen junka.v yshift=-0.3 > junk1.v
        vppen junkb.v yshift=-3.5 > junk2.v
	vppen junk1.v junk2.v erase=o vpstyle=n > SLDDIR/velsyn.v 
	rm back.v junk.v junka.v junkb.v junk1.v junk2.v Junk
	Rm tempbar.H

SLDDIR/datsyn.v : DATDIR/syncos.H DATDIR/synimg0.H
        Taplot < DATDIR/syncos.H | Ta2vplot title="Synthetic data" label1="Time (s)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=junk1.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THICK min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junka.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junkb.v > Junk
	vp_Overlay junka.v junkb.v > junk.v
        Taplot < DATDIR/synimg0.H | Ta2vplot title="Initially migrated image" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/datsyn.v 
	rm back.v junk.v junka.v junkb.v junk1.v junk2.v Junk

SLDDIR/imgsyn.v : DATDIR/synimg.H DATDIR/str_syn.H DATDIR/synrmg.H 
	Graph < DATDIR/str_syn.H FAST yreverse=y THICK min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junka.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junkb.v > Junk
	vp_Overlay junka.v junkb.v > junk.v
        Taplot < DATDIR/synrmg.H | Ta2vplot title="Residual depth migration" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junka.v
        Taplot < DATDIR/synimg.H | Ta2vplot title="Full depth migration" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/imgsyn.v 
	rm back.v junk.v junka.v junkb.v Junk

SLDDIR/datspk.v : DATDIR/spkcos.H DATDIR/spkimg0.H DATDIR/synoper.H
        Taplot < DATDIR/spkcos.H | Ta2vplot title="Synthetic data" label1="Time (s)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=junk1.v > Junk
        Taplot pclip=99.8 < DATDIR/spkimg0.H | Ta2vplot title="Initially migrated image" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	Graph < DATDIR/synoper.H FAST yreverse=y THICK min1=0. max1=4. min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.5,1.5,1.5,1.5,1.5,1.5,1.5 out=junka.v > Junk
	Graph < DATDIR/synoper.H FAST yreverse=y THIN min1=0. max1=4. min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.5,1.5,1.5,1.5,1.5,1.5,1.5 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/datspk.v 
	rm back.v junk1.v junk2.v junka.v junkb.v Junk
 
SLDDIR/imgspk.v : DATDIR/spkimg0.H DATDIR/str_syn.H DATDIR/spkrmg.H 
        Taplot pclip=99.8 < DATDIR/spkrmg.H | Ta2vplot title="Residual depth migraiton" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	vp_Overlay back.v > junka.v
        Taplot pclip=99.8 < DATDIR/spkimg.H | Ta2vplot title="Full depth migration" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i wheretitle=t xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=back.v > Junk
	vp_Overlay back.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/imgspk.v 
	rm back.v junka.v junkb.v Junk

clean&:
	-RM_CMD junk* *.o *.x a.out

#include "./../../Adm/LIN.sld.rules"

#include <SEP.obj.rules>

#include <SEP.prog.rules>
