#define LINDIR ./../..
#define SLDDIR LINDIR/Slds/Wftprg3
#define CHPDIR LINDIR/Wftprg3
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB LINDIR/../MTYPE/lib/libWftprg.a
#define SLDLIST vpupdate vptrigrid vpraytrc vpinit vpgus vpbound vpwftprg1 vpwftprg2 vpwftprg3 vpwptype1 vpwptype2 vpwptype3 ttmgra thgra rgra jjgra ttmlay wavelay twpn twlr twmar1 twmar2 tmmovie
#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/vpupdate.v: vpupdate.P vpupdate.x
	vpupdate.x par=vpupdate.P neconst=1 step=1 title="Phase 1" out=junk1.v > Junk
	vpupdate.x par=vpupdate.P neconst=1 step=2 title="Phase 2" out=junk2.v > Junk
	vpupdate.x par=vpupdate.P neconst=1 step=3 title="Phase 3" out=junk3.v > Junk
	vpupdate.x par=vpupdate.P neconst=2 step=1 title="Phase 1" out=junk4.v > Junk
	vpupdate.x par=vpupdate.P neconst=2 step=2 title="Phase 2" out=junk5.v > Junk
	vpupdate.x par=vpupdate.P neconst=2 step=3 title="Phase 3" out=junk6.v > Junk
	vpupdate.x par=vpupdate.P neconst=3 step=1 title="Phase 1" out=junk7.v > Junk
	vpupdate.x par=vpupdate.P neconst=3 step=2 title="Phase 2" out=junk8.v > Junk
	vpupdate.x par=vpupdate.P neconst=3 step=3 title="Phase 3" out=junk9.v > Junk
	vpupdate.x par=vpupdate.P neconst=4 step=1 title="Phase 1" out=junk10.v > Junk
	vpupdate.x par=vpupdate.P neconst=4 step=2 title="Phase 2" out=junk11.v > Junk
	vpupdate.x par=vpupdate.P neconst=4 step=3 title="Phase 3" out=junk12.v > Junk
	vpupdate.x par=vpupdate.P neconst=5 step=1 title="Phase 1" out=junk13.v > Junk
	vpupdate.x par=vpupdate.P neconst=5 step=2 title="Phase 2" out=junk14.v > Junk
	vpupdate.x par=vpupdate.P neconst=5 step=3 title="Phase 3" out=junk15.v > Junk
	vp_Overlay junk1.v > junk.v
	vp_Movie junk.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v junk9.v junk10.v junk11.v junk12.v junk13.v junk14.v junk15.v > SLDDIR/vpupdate.v
	rm junk.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v junk9.v junk10.v junk11.v junk12.v junk13.v junk14.v junk15.v

SLDDIR/vptrigrid.v: vptrigrid.P vptrigrid.x
	vptrigrid.x par=vptrigrid.P detail=0 out=junk1.v > Junk 
	vptrigrid.x par=vptrigrid.P detail=1 out=junk2.v > Junk 
	vp_Overlay junk1.v > junk.v
	vp_Movie junk.v junk2.v > SLDDIR/vptrigrid.v
	rm junk.v junk1.v junk2.v Junk

SLDDIR/vpinit.v: vpinit.P vpinit.x
	vpinit.x par=vpinit.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpinit.v
	rm junk1.v Junk

SLDDIR/vpgus.v: vpgus.P vpgus.x
	vpgus.x par=vpgus.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpgus.v
	rm junk1.v Junk

SLDDIR/vpraytrc.v: vpraytrc.P vpraytrc.x
	vpraytrc.x par=vpraytrc.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n yshift=-1.25 > SLDDIR/vpraytrc.v
	rm junk1.v Junk

SLDDIR/vpbound.v: vpbound.P vpbound.x
	vpbound.x par=vpbound.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpbound.v
	rm junk1.v Junk

SLDDIR/vpwftprg1.v: vpwftprg1.P vpwftprg1.x
	vpwftprg1.x par=vpwftprg1.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpwftprg1.v
	rm junk1.v Junk

vpwftprg1.x : vpwftprg1.o vpwppnt1.o vpwpq.o vpwpwur.o
        LD LDOPTS vpwftprg1.o vpwppnt1.o vpwpq.o vpwpwur.o USRLIBS SEPLIB VLIB FLIB2 SYSLIB -o vpwftprg1.x

SLDDIR/vpwftprg2.v: vpwftprg2.P vpwftprg2.x
	vpwftprg2.x par=vpwftprg2.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpwftprg2.v
	rm junk1.v Junk

vpwftprg2.x : vpwftprg2.o vpwppnt2.o vpwpq.o vpwpwur.o
        LD LDOPTS vpwftprg2.o vpwppnt2.o vpwpq.o vpwpwur.o USRLIBS SEPLIB VLIB FLIB2 SYSLIB -o vpwftprg2.x

SLDDIR/vpwftprg3.v: vpwftprg3.P vpwftprg3.x
	vpwftprg3.x par=vpwftprg3.P out=junk1.v > Junk 
	vppen junk1.v erase=o vpstyle=n > SLDDIR/vpwftprg3.v
	rm junk1.v Junk

vpwftprg3.x : vpwftprg3.o vpwppnt3.o vpwpq.o vpwpwur.o
        LD LDOPTS vpwftprg3.o vpwppnt3.o vpwpq.o vpwpwur.o USRLIBS SEPLIB VLIB FLIB2 SYSLIB -o vpwftprg3.x

vpwptype%.x : vpwptype%.o 
        LD LDOPTS vpwptype%.o USRLIBS SEPLIB VLIB FLIB2 SYSLIB -o vpwptype%.x

SLDDIR/vpwptype%.v: vpwptype%.P vpwptype%.x
	vpwptype%.x par=vpwptype%.P outtype=0 out=junk1.v > Junk 
	vpwptype%.x par=vpwptype%.P outtype=1 title=" " out=junk2.v > Junk 
	vp_Overlay junk1.v > junk.v
	vp_Movie junk.v junk2.v > SLDDIR/vpwptype%.v
	rm junk.v junk1.v junk2.v Junk

SLDDIR/ttmgra.v : DATDIR/ttmgra.H DATDIR/strvel_gra.H DATDIR/ttmgra_err.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=5. clip=4.2 | Ta2vplot title="Traveltime field" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=1.99 maxval=8.01 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1. o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
        Contour nc=21 dc=0.04 < DATDIR/ttmgra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junk1.v > Junk
        Contour nc=21 dc=0.04 < DATDIR/ttmgra.H THIN FAST pad=n xinch=10.24 yinch=5.12 out=junk2.v > Junk
	vp_Overlay back.v bar.v junk1.v junk2.v > junka.v
	Window j1=2 < DATDIR/ttmgra_err.H > ttmgra.temp
        Taplot < ttmgra.temp bias=0.0 clip=0.000014 | Ta2vplot title="Traveltime error" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=-0.01 maxval=0.01> tempbar.H
        Taplot < tempbar.H gpow=1 bias=0. clip=0.014 | Ta2vplot title=" " label2="Error (ms)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.005 o2num=-0.01 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
	vp_Overlay back.v bar.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/ttmgra.v
        rm back.v bar.v junk1.v junk2.v junka.v junkb.v junk.v Junk
        Rm tempbar.H ttmgra.temp

SLDDIR/thgra.v : DATDIR/thgra.H DATDIR/strvel_gra.H DATDIR/thgra_ana.H DATDIR/thgra_err.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=5. clip=4.2 | Ta2vplot title="Take-off-angle field" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=1.99 maxval=8.01 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1. o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
        Contour nc=37 dc=0.087266462 c0=-1.570796327 < DATDIR/thgra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junk1.v > Junk
        Contour nc=37 dc=0.087266462 c0=-1.570796327 < DATDIR/thgra.H THIN FAST pad=n xinch=10.24 yinch=5.12 out=junk2.v > Junk
	vp_Overlay  back.v bar.v junk1.v junk2.v > junka.v
	Window j1=2 < DATDIR/thgra_err.H > thgra.temp
        Taplot < thgra.temp bias=0.0 clip=0.004547 | Ta2vplot title="Take-off-angle error" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=-0.2 maxval=0.2 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=0. clip=0.28 | Ta2vplot title=" " label2="Error (degrees)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.1 o2num=-0.2 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
	vp_Overlay back.v bar.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/thgra.v
        rm back.v bar.v junk1.v junk2.v junka.v junkb.v junk.v Junk
        Rm tempbar.H thgra.temp

SLDDIR/rgra.v : DATDIR/rgra.H DATDIR/strvel_gra.H DATDIR/rgra_ana.H DATDIR/rgra_err.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=5. clip=4.2 | Ta2vplot title="Curvature-radius field" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=1.99 maxval=8.01 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1. o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
        Contour nc=19 dc=0.2 c0=0. < DATDIR/rgra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junk1.v > Junk
        Contour nc=19 dc=0.2 c0=0. < DATDIR/rgra.H THIN FAST pad=n xinch=10.24 yinch=5.12 out=junk2.v > Junk
	vp_Overlay  back.v bar.v junk1.v junk2.v > junka.v
	Window j1=2 < DATDIR/rgra_err.H > rgra.temp
        Taplot < rgra.temp bias=0.0025 clip=0.0063 | Ta2vplot title="Curvature-radius error" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=-0.002 maxval=0.007 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=0.0025 clip=0.0063 | Ta2vplot title=" " label2="Error (km)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.003 o2num=-0.002 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
	vp_Overlay back.v bar.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/rgra.v
        rm back.v bar.v junk1.v junk2.v junka.v junkb.v junk.v Junk
        Rm tempbar.H rgra.temp

SLDDIR/jjgra.v : DATDIR/jjgra.H DATDIR/strvel_gra.H DATDIR/jjgra_ana.H DATDIR/jjgra_err.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=5. clip=4.2 | Ta2vplot title="Spreading-factor field" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=1.99 maxval=8.01 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1. o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
        Contour nc=21 dc=0.4 c0=0. < DATDIR/jjgra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junk1.v > Junk
        Contour nc=21 dc=0.4 c0=0. < DATDIR/jjgra.H THIN FAST pad=n xinch=10.24 yinch=5.12 out=junk2.v > Junk
	vp_Overlay  back.v bar.v junk1.v junk2.v > junka.v
	Window j1=2 < DATDIR/jjgra_err.H > jjgra.temp
        Taplot < jjgra.temp bias=0.012 clip=0.0168 | Ta2vplot title="Spreading-factor error" label1="Depth (km)" label2="Surface (km)" TLSZFAT color=G xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 wheretitle=t grid=n polarity=1 out=back.v > Junk
        Bar n2=200 type=h minval=0.0 maxval=0.024 > tempbar.H
        Taplot < tempbar.H gpow=1 bias=0.012 clip=0.0168 | Ta2vplot title=" " label2="Error (km)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.012 o2num=-0.024 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
	vp_Overlay back.v bar.v > junkb.v
	vp_Movie junka.v junkb.v > SLDDIR/jjgra.v
        rm back.v bar.v junk1.v junk2.v junka.v junkb.v junk.v Junk
        Rm tempbar.H jjgra.temp

SLDDIR/ttmlay.v : DATDIR/strvel_lay.H DATDIR/ttmlay0.H DATDIR/ttmlay1.H DATDIR/ttmlay2.H
        Taplot < DATDIR/strvel_lay.H bias=3. clip=3. | Ta2vplot title="First arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Bar n2=100 type=h minval=1.499 maxval=4.501 > tempbar.H
	Taplot < tempbar.H bias=3. clip=3. | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1. o2num=1.5 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
	vppen junk.v yshift=-3 > bar.v
        Contour nc=24 c0=0 dc=0.1 < DATDIR/ttmlay0.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=n wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=24 c0=0 dc=0.1 < DATDIR/ttmlay0.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
        Taplot < DATDIR/strvel_lay.H bias=3. clip=3. | Ta2vplot title="Most energetic arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
        Contour nc=30 c0=0 dc=0.1 < DATDIR/ttmlay1.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=n wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=30 c0=0 dc=0.1 < DATDIR/ttmlay1.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk2.v
        Taplot < DATDIR/strvel_lay.H bias=3. clip=3. | Ta2vplot title="Upgoing arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
        Contour nc=30 c0=0 dc=0.1 < DATDIR/ttmlay2.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=n wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=30 c0=0 dc=0.1 < DATDIR/ttmlay2.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk3.v
	vp_Movie junk1.v junk2.v junk3.v > SLDDIR/ttmlay.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk3.v junk.v Junk
	Rm tempbar.H

SLDDIR/wavelay.v : DATDIR/ttmlay0.H DATDIR/ttmlay1.H DATDIR/ttmlay2.H DATDIR/snaplay.H
	Window n3=1 f3=6 < DATDIR/snaplay.H | Scale dscale=0 axis=123 > templay.H
	Taplot < templay.H | Ta2vplot title="First arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay0.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay0.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk1.v
	Taplot < templay.H | Ta2vplot title="Most Energetic arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay1.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay1.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
	Taplot < templay.H | Ta2vplot title="Upgoing arrival" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay2.H THICK d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
	Contour nc=1 c0=1.399 < DATDIR/ttmlay2.H THIN d1num=0 d2num=0 title=" " label1=" " label2=" " wantframe=y wantaxis=n pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk3.v
	vp_Movie junk1.v junk2.v junk3.v > SLDDIR/wavelay.v
	rm back.v junka.v junkb.v junk1.v junk2.v junk3.v Junk
	Rm templay.H

SLDDIR/twpn.v : DATDIR/gauneg.H DATDIR/gaupos.H DATDIR/ttmneg.H DATDIR/ttmpos.H DATDIR/snapneg.H DATDIR/snappos.H
        Taplot < DATDIR/gauneg.H bias=3. clip=1.3 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
	Bar n2=100 type=h minval=1.49 maxval=3.01 > tempbar.H
	Taplot < tempbar.H bias=3. clip=1.3 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.5 o2num=1.5 grid=n polarity=1 xinch=7.5 yinch=0.4 out=junk.v > Junk
	vppen junk.v yshift=-4.2 > bar.v
        Contour nc=20 c0=0 dc=0.08 < DATDIR/ttmneg.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.08 < DATDIR/ttmneg.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
	Window n3=1 f3=3 < DATDIR/snapneg.H | Scale dscale=0 axis=123 > tempneg.H
	Taplot < tempneg.H | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
        Contour nc=1 c0=0.319 dc=0.08 < DATDIR/ttmneg.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.319 dc=0.08 < DATDIR/ttmneg.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
        Taplot < DATDIR/gaupos.H bias=3. clip=1.3 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
	Bar n2=100 type=h minval=2.99 maxval=4.51 > tempbar.H
	Taplot < tempbar.H bias=3. clip=1.3 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=0.5 o2num=3. grid=n polarity=1 xinch=7.5 yinch=0.4 out=junk.v > Junk
	vppen junk.v yshift=-4.2 > bar.v
        Contour nc=20 c0=0 dc=0.08 < DATDIR/ttmpos.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.08 < DATDIR/ttmpos.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk3.v
	Window n3=1 f3=3 < DATDIR/snappos.H | Scale dscale=0 axis=123 > temppos.H
	Taplot < temppos.H | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
        Contour nc=1 c0=0.319 dc=0.14 < DATDIR/ttmpos.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.319 dc=0.14 < DATDIR/ttmpos.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk4.v
	vp_Movie junk1.v junk2.v junk3.v junk4.v > junk.v
	vppen junk.v yshift=0.4 vpstyle=n > SLDDIR/twpn.v
	rm back.v bar.v junka.v junkb.v junk.v junk1.v junk2.v junk3.v junk4.v Junk
	Rm temppos.H tempneg.H tempbar.H

SLDDIR/twlr.v : DATDIR/ttmlft.H DATDIR/ttmrgt.H DATDIR/snapneg.H
	Window n3=1 f3=5 < DATDIR/snapneg.H | Scale dscale=0 axis=123 > tempneg.H
	Taplot < tempneg.H | Ta2vplot title="Left-going" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
        Contour nc=1 c0=0.479 < DATDIR/ttmlft.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.479 < DATDIR/ttmlft.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk1.v
	Taplot < tempneg.H | Ta2vplot title="Right-going" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=7.5 yinch=7.5 out=back.v > Junk
        Contour nc=1 c0=0.479 < DATDIR/ttmrgt.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.479 < DATDIR/ttmrgt.H THIN d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/twlr.v
	rm back.v junka.v junkb.v junk1.v junk2.v Junk
	Rm tempneg.H

SLDDIR/twmar1.v : DATDIR/marmvel1.H DATDIR/ttmmar0.H DATDIR/snapmar.H
        Taplot bias=3000 clip=2500 < DATDIR/marmvel1.H | Ta2vplot title="Traveltime field" label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT color=G d1num=1000 d2num=1000. o1num=0 o2num=2000 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Bar n2=100 type=h minval=1500 maxval=5500 | Taplot bias=3000 clip=2500 | Ta2vplot title=" " label2="Velocity (m/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1000 o2num=1500 grid=n polarity=1 xinch=10.24 yinch=0.3 out=junk.v > Junk
	vppen junk.v yshift=-3 > bar.v
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
	Window n3=1 f3=9 < DATDIR/snapmar.H | Scale dscale=0 axis=123 > tempsnap.H
	Taplot pclip=96 < tempsnap.H | Ta2vplot title="Wave field" label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1000. d2num=1000. o1num=0 o2num=2000 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Contour nc=1 c0=0.999 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
	Contour nc=1 c0=0.999 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/twmar1.v
	rm back.v junka.v junkb.v junk1.v junk2.v bar.v junk.v Junk
	Rm tempsnap.H

SLDDIR/twmar2.v : DATDIR/marmvel1.H DATDIR/ttmmar1.H DATDIR/snapmar.H
        Taplot bias=3000 clip=2500 < DATDIR/marmvel1.H | Ta2vplot title="Traveltime field" label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT color=G d1num=1000 d2num=1000. o1num=0 o2num=2000 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Bar n2=100 type=h minval=1500 maxval=5500 | Taplot bias=3000 clip=2500 | Ta2vplot title=" " label2="Velocity (m/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1000 o2num=1500 grid=n polarity=1 xinch=10.24 yinch=0.3 out=junk.v > Junk
	vppen junk.v yshift=-3 > bar.v
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttmmar1.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttmmar1.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
	Window n3=1 f3=9 < DATDIR/snapmar.H | Scale dscale=0 axis=123 > tempsnap.H
	Taplot pclip=96 < tempsnap.H | Ta2vplot title="Wave field" label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT wantaxis2=y color=G d1num=1000. d2num=1000. o1num=0 o2num=2000 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Contour nc=1 c0=0.999 < DATDIR/ttmmar1.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
	Contour nc=1 c0=0.999 < DATDIR/ttmmar1.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v junka.v junkb.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/twmar2.v
	rm back.v junka.v junkb.v junk1.v junk2.v bar.v junk.v Junk
	Rm tempsnap.H

SLDDIR/tmmovie.v : DATDIR/marmvel1.H DATDIR/ttmmar0.H
        Taplot bias=3000 clip=2500 < DATDIR/marmvel1.H | Ta2vplot title="Marmousi model" label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT color=G d1num=1000 d2num=1000. o1num=0 o2num=2000 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Bar n2=100 type=h minval=1500 maxval=5500 | Taplot bias=3000 clip=2500 | Ta2vplot title=" " label2="Velocity (m/s)" wherexlabel=b TLSZFAT color=G wantaxis1=n d2num=1000 o2num=1500 grid=n polarity=1 xinch=10.24 yinch=0.3 out=junk.v > Junk
	vppen junk.v yshift=-3 > bar.v
        Contour nc=1 c0=0.01 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.01 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
        Contour nc=1 c0=0.1 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.1 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk2.v
        Contour nc=1 c0=0.2 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.2 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk2.v
        Contour nc=1 c0=0.3 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.3 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk3.v
        Contour nc=1 c0=0.4 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.4 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk4.v
        Contour nc=1 c0=0.5 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.5 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk5.v
        Contour nc=1 c0=0.6 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.6 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk6.v
        Contour nc=1 c0=0.7 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.7 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk7.v
        Contour nc=1 c0=0.8 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.8 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk8.v
        Contour nc=1 c0=0.9 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=0.9 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk9.v
        Contour nc=1 c0=1. dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1. dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk10.v
        Contour nc=1 c0=1.1 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.1 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk11.v
        Contour nc=1 c0=1.2 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.2 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk12.v
        Contour nc=1 c0=1.3 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.3 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk13.v
        Contour nc=1 c0=1.4 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.4 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk14.v
        Contour nc=1 c0=1.5 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.5 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk15.v
        Contour nc=1 c0=1.6 dc=0.1 < DATDIR/ttmmar0.H THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=1 c0=1.6 dc=0.1 < DATDIR/ttmmar0.H THIN d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junkb.v > Junk
	vp_Overlay back.v bar.v junka.v junkb.v > junk16.v
	vp_Movie junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v junk9.v junk10.v junk11.v junk12.v junk13.v junk14.v junk15.v junk16.v > SLDDIR/tmmovie.v
	rm back.v junka.v junkb.v junk1.v junk2.v junk3.v junk4.v junk5.v junk6.v junk7.v junk8.v junk9.v junk10.v junk11.v junk12.v junk13.v junk14.v junk15.v junk16.v bar.v junk.v Junk

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

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
