#define LINDIR ./../..
#define SLDDIR LINDIR/Slds/Fiedat5
#define CHPDIR LINDIR/Fiedat5
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB
#define SLDLIST fiewlt fienf fiecmp fiecom fievel \
		fieimgavg fieimg fieimgv0 fiermg
#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/datacube.v : DATDIR/fiecos.A
$	Cubeplot < DATDIR/fiecos.A frame1=350 frame2=200 frame3=0 point1=0.8 point2=0.8 title="Data cube" label1="Time (s)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.4 d2num=1. o1num=0.8 o2num=0 grid=n polarity=1 xinch=10.24 yinch=7.5 out=SLDDIR/datacube.v > Junk
$	rm Junk

$SLDDIR/imagecube.v : DATDIR/fiemig.A
$	Cubeplot < DATDIR/fiemig.A frame1=186 frame2=200 frame3=0 point1=0.8 point2=0.8 title="Image cube" label1="Time (s)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.4 d2num=1. o1num=0.8 o2num=0 grid=n polarity=1 xinch=10.24 yinch=7.5 out=SLDDIR/imagecube.v > Junk
$	rm Junk

SLDDIR/fiewlt.v : DATDIR/fiewlt.H
	Graph < DATDIR/fiewlt.H TLSZFAT title="Estimated source wavelet" label1="Time (s)" label2="Amplitude" plotfat=4 xinch=10.24 yinch=6.827 out=junk.v > Junk
	vppen junk.v vpstyle=n scale=1. > SLDDIR/fiewlt.v
	rm junk.v Junk

SLDDIR/fienf.v : DATDIR/fienear.H DATDIR/fiefar.H
	Taplot tpow=1. < DATDIR/fienear.H | Ta2vplot title="Near offset section" label1="Time (s)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.4 d2num=1. o1num=0.8 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junk1.v > Junk
	Taplot tpow=1. < DATDIR/fiefar.H | Ta2vplot title="Far offset section" label1="Time (s)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.4 d2num=1. o1num=0.8 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junk2.v > Junk
	vp_Movie junk1.v junk2.v > SLDDIR/fienf.v
	rm junk1.v junk2.v Junk

SLDDIR/fiecmp.v : DATDIR/fiecmp1.H DATDIR/fiecmp2.H
	Taplot tpow=1. < DATDIR/fiecmp1.H | Ta2vplot title=" " label1="Time (s)" label2="Offset (km)" wheretitle=t TLSZFAT color=i d1num=0.4 d2num=0.4 o1num=0.8 o2num=0.1 grid=n polarity=1 xinch=3. yinch=7.5 out=junka.v > Junk
	Box x0=6.4 y0=3.1 size=0.1 label="Non-hyperbolic\\\\n    moveout" lab_fat=3 lab_color=4 xt=0.000000 yt=-1.00000 out=point.v > Junk
	vp_Overlay junka.v point.v > junk1.v
	Taplot tpow=1. < DATDIR/fiecmp2.H | Ta2vplot title=" " label1=" " label2="Offset (km)" wantaxis1=n wheretitle=t TLSZFAT color=i d1num=0.4 d2num=0.4 o1num=0.8 o2num=0.1 grid=n polarity=1 xinch=3. yinch=7.5 out=junkb.v > Junk
	Box x0=6.85 y0=2.4 size=0.1 label="Water-bottom\\\\n    multiple" lab_fat=3 lab_color=4 xt=0.000000 yt=1.000 out=point.v > Junk
	vp_Overlay junkb.v point.v > junk2.v
	vppen junk1.v xshift=-1 > junka.v
	vppen junk2.v xshift=2.5 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=1. > SLDDIR/fiecmp.v
	rm point.v junka.v junkb.v junk1.v junk2.v Junk
	
SLDDIR/fiecom.v : DATDIR/fierawsamp.H DATDIR/fiesamp.H 
	< DATDIR/fierawsamp.H Tpow tpow=2 | Wiggle zplot=1.2 TLSZFAT title="Before" label1="Time (s)" label2="Midpoint (km)" wherexlabel=t poly=y transp=1 yreverse=1 xinch=4.5 yinch=7.5 out=junka.v > Junk
	< DATDIR/fiesamp.H Tpow tpow=2 | Wiggle zplot=1.2 TLSZFAT title="After" label1="Time (s)" label2="Midpoint (km)" wherexlabel=t poly=y transp=1 yreverse=1 xinch=4.5 yinch=7.5 out=junkb.v > Junk
	Scale axis=12 < DATDIR/fierawsamp.H | Spectrum | Graph TLSZFAT title="Spectrum" label1="Frequency (Hz)" label2="Amplitude" plotfat=4 xinch=4.5 yinch=3 out=junkc.v > Junk
	Scale axis=12 < DATDIR/fiesamp.H | Spectrum | Graph TLSZFAT title="Spectrum" label1="Frequency (Hz)" label2="Amplitude" plotfat=4 xinch=4.5 yinch=3 out=junkd.v > Junk
	vppen junka.v xshift=-2 yshift=0 > junk1.v
	vppen junkb.v xshift=4 yshift=0 > junk2.v
	vp_Overlay junk1.v junk2.v > junka.v
	vppen junkc.v xshift=-2 yshift=0 > junk3.v
	vppen junkd.v xshift=4 yshift=0 > junk4.v
	vp_Overlay junk3.v junk4.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/fiecom.v
	rm junka.v junkb.v junkc.v junkd.v junk1.v junk2.v junk3.v junk4.v Junk

SLDDIR/fievel.v : DATDIR/fievelf.H DATDIR/fiestrf.H
	Graph < DATDIR/fiestrf.H TLSZFAT title="Picked reflectors" label2="Depth (km)" label1="Surface (km)" wherexlabel=t min2=0.6 max2=1.3 yreverse=y pad=n THIN xinch=10.24 yinch=6.12 out=interface.v > Junk
	vp_Overlay interface.v > junk1.v
	Window n1=167 f1=150 < DATDIR/fievelf.H > fievelf.temp
        Taplot < fievelf.temp bias=1.94 clip=0.46 | Ta2vplot title="Estimated interval velocity model" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=6.12 out=junk.v > Junk
	Contour nc=10 c0=1.482 dc=0.1 < fievelf.temp THICK FAST xinch=10.24 yinch=6.12 out=junka.v > Junk
	Contour nc=10 c0=1.482 dc=0.1 < fievelf.temp THIN FAST xinch=10.24 yinch=6.12 out=junkb.v > Junk
	vp_Overlay junk.v junka.v junkb.v > back.v
        Bar n2=100 type=h minval=1.479 maxval=2.41 | Taplot bias=1.94 clip=0.46 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.3 o2num=1.48 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3.5 > bar.v
	vp_Overlay back.v bar.v > junk2.v
        vp_Movie junk1.v junk2.v > SLDDIR/fievel.v
        rm back.v bar.v junk.v junka.v junkb.v junk1.v junk2.v interface.v Junk
	Rm fievelf.temp

SLDDIR/fieimgavg.v : DATDIR/fiemigavg.H  DATDIR/fieimgavg.H
	Taplot tpow=0. < DATDIR/fieimgavg.H | Ta2vplot title="Stacked image after migration" label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junka.v > Junk
	Window n2=1 f2=37 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis2=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-4.6 > junk1.v
	Window n2=1 f2=46 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-3.285 > junk2.v
	Window n2=1 f2=124 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-1.971 > junk3.v
	Window n2=1 f2=169 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title="\\\\k350 Common surface location gathers" label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-0.657 > junk4.v
	Window n2=1 f2=197 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=0.657 > junk5.v
	Window n2=1 f2=269 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=1.971 > junk6.v
	Window n2=1 f2=368 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=3.285 > junk7.v
	Window n2=1 f2=386 < DATDIR/fiemigavg.H > fiemigavg.temp
	Taplot tpow=0. < fiemigavg.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=4.6 > junk8.v
	vp_Overlay junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v > junkb.v
        vp_Movie junka.v junkb.v > SLDDIR/fieimgavg.v
	rm junka.v junkb.v junk.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v Junk
	Rm fiemigavg.temp

SLDDIR/fieimg.v : DATDIR/fiemig.H  DATDIR/fieimg.H
	Taplot tpow=0. < DATDIR/fieimg.H | Ta2vplot title="Stacked image after migration" label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junka.v > Junk
	Window n2=1 f2=37 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis2=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-4.6 > junk1.v
	Window n2=1 f2=46 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-3.285 > junk2.v
	Window n2=1 f2=124 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-1.971 > junk3.v
	Window n2=1 f2=169 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title="\\\\k350 Common surface location gathers" label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=-0.657 > junk4.v
	Window n2=1 f2=197 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=0.657 > junk5.v
	Window n2=1 f2=269 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=1.971 > junk6.v
	Window n2=1 f2=368 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=3.285 > junk7.v
	Window n2=1 f2=386 < DATDIR/fiemig.H > fiemig.temp
	Taplot tpow=0. < fiemig.temp | Ta2vplot title=" " label1="Depth (km)" label2="(km)" wheretitle=t TLSZFAT color=i wantaxis=n d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=1.2 yinch=7.5 out=junk.v > Junk
	vppen junk.v xshift=4.6 > junk8.v
	vp_Overlay junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v > junkb.v
        vp_Movie junka.v junkb.v > SLDDIR/fieimg.v
	rm junka.v junkb.v junk.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v Junk
	Rm fiemig.temp

SLDDIR/fieimgv0.v : DATDIR/fieimgv0.H  DATDIR/fiermg.H
	Taplot tpow=0. < DATDIR/fieimgv0.H | Ta2vplot title="Constant velocity migration" label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junka.v > Junk
	Taplot tpow=0.4 < DATDIR/fiermg.H | Ta2vplot title="Residual depth migration" label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junkb.v > Junk
        vp_Movie junka.v junkb.v > SLDDIR/fieimgv0.v
	rm junka.v junkb.v Junk

SLDDIR/fiermg.v : DATDIR/fiermg.H  DATDIR/fieimgzos.H
	Taplot tpow=0.4 < DATDIR/fiermg.H | Ta2vplot title="Residual depth migration " label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junka.v > Junk
	Taplot tpow=0. < DATDIR/fieimgzos.H | Ta2vplot title="Full depth migration " label1="Depth (km)" label2="Midpoint (km)" wheretitle=t TLSZFAT color=i d1num=0.2 d2num=0.5 o1num=0.6 o2num=1. grid=n polarity=1 xinch=10.24 yinch=7.5 out=junkb.v > Junk
        vp_Movie junka.v junkb.v > SLDDIR/fiermg.v
	rm junka.v junkb.v Junk

clean&:
	-RM_CMD *.o *.x *.dvi *.aux *.log
	-RM_CMD a.out

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
