#define LINDIR ./../..
#define FIGDIR LINDIR/Figs/Wftprg3
#define CHPDIR LINDIR/Wftprg3
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB LINDIR/../MTYPE/lib/libWftprg.a
#define FIGLIST vptrigrid vpinit vpgus vpupdate vpraytrc vpbound \
	        vpwftprg1 vpwftprg2 vpwftprg3 vpwptype1 vpwptype2 \
	        vpwptype3 ttmgra thgra rgra jjgra ttmlay wavelay \
		twpn twlr twmar1 twmar2
#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>

FIGDIR/vpupdate.v: vpupdate.P vpupdate.x
	vpupdate.x par=vpupdate.P neconst=4 step=1 title=" " out=junk1.v > Junk
	vpupdate.x par=vpupdate.P neconst=4 step=2 title=" " out=junk2.v > Junk
	vpupdate.x par=vpupdate.P neconst=4 step=3 title=" " out=junk3.v > Junk
	vppen junk1.v yshift=9.2 > junka.v
	vppen junk2.v yshift=3.4 > junkb.v
	vppen junk3.v yshift=-2.4 > junkc.v
	vppen junka.v junkb.v junkc.v erase=o vpstyle=n scale=0.44 > FIGDIR/vpupdate.v
	rm junk1.v junk2.v junk3.v junka.v junkb.v junkc.v Junk

FIGDIR/vpinit.v: vpinit.P vpinit.x
	vpinit.x par=vpinit.P neconst=3 title=" " out=junk.v > Junk
	vppen junk.v vpstyle=n scale=0.5758 > FIGDIR/vpinit.v
	rm junk.v Junk

FIGDIR/vpgus.v: vpgus.P vpgus.x
	vpgus.x par=vpgus.P neconst=4 title=" " out=junk.v > Junk
	vppen junk.v vpstyle=n scale=0.5758 > FIGDIR/vpgus.v
	rm junk.v Junk

FIGDIR/vptrigrid.v: vptrigrid.P vptrigrid.x
	vptrigrid.x par=vptrigrid.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n scale=0.5758 > FIGDIR/vptrigrid.v
	rm junk1.v Junk

FIGDIR/vpraytrc.v: vpraytrc.P vpraytrc.x
	vpraytrc.x par=vpraytrc.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=-0.8 scale=0.36 > FIGDIR/vpraytrc.v
	rm junk1.v Junk

FIGDIR/vpbound.v: vpbound.P vpbound.x
	vpbound.x par=vpbound.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=-0.5 scale=0.285 > FIGDIR/vpbound.v
	rm junk1.v Junk

FIGDIR/vpwftprg.v: vpwftprg.P vpwftprg.x
	vpwftprg.x par=vpwftprg.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=0. scale=1. > FIGDIR/vpwftprg.v
	rm junk1.v Junk

FIGDIR/vpwftprg1.v: vpwftprg1.P vpwftprg1.x
	vpwftprg1.x par=vpwftprg1.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=0. scale=0.325 > FIGDIR/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

FIGDIR/vpwftprg2.v: vpwftprg2.P vpwftprg2.x
	vpwftprg2.x par=vpwftprg2.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=0. scale=0.325 > FIGDIR/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

FIGDIR/vpwftprg3.v: vpwftprg3.P vpwftprg3.x
	vpwftprg3.x par=vpwftprg3.P title=" " out=junk1.v > Junk 
	vppen junk1.v vpstyle=n yshift=0. scale=0.325 > FIGDIR/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

FIGDIR/vpwptype%.v: vpwptype%.P vpwptype%.x
	vpwptype%.x par=vpwptype%.P outtype=0 title=" " out=junk1.v > Junk 
	vpwptype%.x par=vpwptype%.P outtype=1 title=" " out=junk.v > Junk 
	vppen junk.v xshift=9.8 > junk2.v
	vppen junk1.v junk2.v vpstyle=n erase=n scale=0.315 > FIGDIR/vpwptype%.v
	rm junk1.v Junk

FIGDIR/ttmgra.v : DATDIR/ttmgra.H DATDIR/strvel_gra.H DATDIR/ttmgra_err.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t 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=i 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.0 > bar.v
	Window j1=2 j2=2 < DATDIR/ttmgra.H > ttmgra.temp
        Contour nc=21 c0=0 dc=0.04 < ttmgra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=21 c0=0 dc=0.04 < ttmgra.temp THIN 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 j1=2 j2=2 < DATDIR/ttmgra_err.H > ttmgra.temp
        Taplot < ttmgra.temp bias=0.0 clip=0.000014 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t polarity=1 out=back.v > Junk
	Bar n2=200 type=h minval=-0.01 maxval=0.01 > tempbar.H
	Taplot < tempbar.H bias=0.0 clip=0.014 | Ta2vplot title=" " label2="Error (s)" wherexlabel=b TLSZFAT color=i 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.0 > bar.v
	vp_Overlay back.v bar.v > junk2.v
	vppen junk1.v yshift=8.2 > junka.v
	vppen junk2.v yshift=0 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.45 > FIGDIR/ttmgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H ttmgra.temp strvel_gra.temp ttmgra.temp

FIGDIR/thgra.v : DATDIR/thgra.H DATDIR/strvel_gra.H DATDIR/thgra_err.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t 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=i 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.0 > bar.v
	Window j1=2 j2=2 < DATDIR/thgra.H > thgra.temp
        Contour nc=37 dc=0.087266462 c0=-1.570796327 < thgra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=37 dc=0.087266462 c0=-1.570796327 < thgra.temp THIN 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 j1=2 j2=2 < DATDIR/thgra_err.H > thgra.temp
        Taplot < thgra.temp bias=0.0 clip=0.004547 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t polarity=1 out=back.v > Junk
	Bar n2=200 type=h minval=-0.2 maxval=0.2 > tempbar.H
	Taplot < tempbar.H bias=0.0 clip=0.28 | Ta2vplot title=" " label2="Error (s)" wherexlabel=b TLSZFAT color=i 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.0 > bar.v
	vp_Overlay back.v bar.v > junk2.v
	vppen junk1.v yshift=8.2 > junka.v
	vppen junk2.v yshift=0 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.45 > FIGDIR/thgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H thgra.temp strvel_gra.temp

FIGDIR/rgra.v : DATDIR/rgra.H DATDIR/strvel_gra.H DATDIR/rgra_err.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t 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=i 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.0 > bar.v
	Window j1=2 j2=2 < DATDIR/rgra.H > rgra.temp
        Contour nc=19 dc=0.2 c0=0. < rgra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=19 dc=0.2 c0=0. < rgra.temp THIN 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 j1=2 j2=2 < DATDIR/rgra_err.H > rgra.temp
        Taplot < rgra.temp bias=0.0025 clip=0.0063 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t polarity=1 out=back.v > Junk
	Bar n2=200 type=h minval=-0.002 maxval=0.007 > tempbar.H
	Taplot < tempbar.H bias=0.0025 clip=0.0063 | Ta2vplot title=" " label2="Error (s)" wherexlabel=b TLSZFAT color=i 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.0 > bar.v
	vp_Overlay back.v bar.v > junk2.v
	vppen junk1.v yshift=8.2 > junka.v
	vppen junk2.v yshift=0 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.45 > FIGDIR/rgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H rgra.temp strvel_gra.temp

FIGDIR/jjgra.v : DATDIR/jjgra.H DATDIR/strvel_gra.H DATDIR/jjgra_err.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=5. clip=4.2 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t 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=i 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.0 > bar.v
	Window j1=2 j2=2 < DATDIR/jjgra.H > jjgra.temp
        Contour nc=21 c0=0 dc=0.4 < jjgra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=21 c0=0 dc=0.4 < jjgra.temp THIN 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 j1=2 < DATDIR/jjgra_err.H > jjgra.temp
        Taplot < jjgra.temp bias=0.012 clip=0.0168 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n wheretitle=t polarity=1 out=back.v > Junk
	Bar n2=200 type=h minval=0.0 maxval=0.024 > tempbar.H
	Taplot < tempbar.H bias=0.012 clip=0.0168 | Ta2vplot title=" " label2="Error (s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.006 o2num=0.0 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3.0 > bar.v
	vp_Overlay back.v bar.v > junk2.v
	vppen junk1.v yshift=8.2 > junka.v
	vppen junk2.v yshift=0 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.45 > FIGDIR/jjgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H jjgra.temp strvel_gra.temp

FIGDIR/ttmlay.v : DATDIR/strvel_lay.H DATDIR/ttmlay0.H DATDIR/ttmlay1.H DATDIR/ttmlay2.H
	Window j1=2 j2=2 < DATDIR/strvel_lay.H > strvel_lay.temp
        Taplot < strvel_lay.temp bias=3. clip=3. | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Window j1=2 j2=2 < DATDIR/ttmlay0.H > ttmlay0.temp
        Contour nc=24 c0=0 dc=0.1 < ttmlay0.temp 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 < ttmlay0.temp 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
	Window j1=2 j2=2 < DATDIR/strvel_lay.H > strvel_lay.temp
        Taplot < strvel_lay.temp bias=3. clip=3. | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i wantaxis2=n d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Window j1=2 j2=2 < DATDIR/ttmlay1.H > ttmlay1.temp
        Contour nc=30 c0=0 dc=0.1 < ttmlay1.temp 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 < ttmlay1.temp 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
	Window j1=2 j2=2 < DATDIR/strvel_lay.H > strvel_lay.temp
        Taplot < strvel_lay.temp bias=3. clip=3. | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i wantaxis2=n d1num=1. d2num=1. o1num=0 o2num=0 grid=n polarity=1 xinch=10.24 yinch=5.12 out=back.v > Junk
	Window j1=2 j2=2 < DATDIR/ttmlay2.H > ttmlay2.temp
        Contour nc=30 c0=0 dc=0.1 < ttmlay2.temp 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 < ttmlay2.temp 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
	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=i wantaxis1=n d2num=1. o2num=1.5 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
	vppen junk1.v yshift=10.4 > junka.v
	vppen junk2.v yshift=4.9 > junkb.v
	vppen junk3.v yshift=-0.6 > junkc.v
	vppen junk.v yshift=-3.8 > bar.v
	vppen junka.v junkb.v junkc.v bar.v erase=o vpstyle=n scale=0.40 > FIGDIR/ttmlay.v
	rm back.v junka.v junkb.v junkc.v junk1.v junk2.v junk3.v junk.v bar.v Junk
	Rm tempbar.H strvel_lay.temp ttmlay1.temp ttmlay2.temp ttmlay0.temp

FIGDIR/wavelay.v : DATDIR/ttmlay0.H DATDIR/ttmlay1.H DATDIR/ttmlay2.H DATDIR/snaplay.H
	Window j1=2 j2=2 n3=1 f3=6 < DATDIR/snaplay.H | Scale dscale=0 axis=123 > templay.H
	Taplot < templay.H | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i 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=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=n color=i 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=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT wantaxis2=n color=i 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
	vppen junk1.v yshift=8.6 > junka.v
	vppen junk2.v yshift=3.1 > junkb.v
	vppen junk3.v yshift=-2.4 > junkc.v
	vppen junka.v junkb.v junkc.v erase=n vpstyle=n scale=0.43 > FIGDIR/wavelay.v
	rm back.v junka.v junkb.v junkc.v junk1.v junk2.v junk3.v Junk
	Rm templay.H

FIGDIR/twpn.v : DATDIR/gauneg.H DATDIR/gaupos.H DATDIR/ttmneg.H DATDIR/ttmpos.H DATDIR/snapneg.H DATDIR/snappos.H
	Window j1=2 j2=2 < DATDIR/gauneg.H > gauneg.temp
        Taplot < gauneg.temp bias=3.0 clip=1.3 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i 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.0 clip=1.3 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i 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
	Window j1=2 j2=2 < DATDIR/ttmneg.H > ttmneg.temp
        Contour nc=20 c0=0 dc=0.08 < ttmneg.temp 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 < ttmneg.temp 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 j1=2 j2=2 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=i 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 < 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 < 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
	Window j1=2 j2=2 < DATDIR/gaupos.H > gaupos.temp
        Taplot < gaupos.temp bias=3.0 clip=1.3 | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i 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.0 clip=1.3 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.5 o2num=3.0 grid=n polarity=1 xinch=7.5 yinch=0.4 out=junk.v > Junk
	vppen junk.v yshift=-4.2 > bar.v
	Window j1=2 j2=2 < DATDIR/ttmpos.H > ttmpos.temp
        Contour nc=20 c0=0 dc=0.08 < ttmpos.temp 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 < ttmpos.temp 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 j1=2 j2=2 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=i 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 < 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 < 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
        vppen junk1.v xshift=-2 yshift=10. > junka.v
	vppen junk2.v xshift=-2 yshift=-1.2 > junkb.v
        vppen junk3.v xshift=7 yshift=10. > junkc.v
        vppen junk4.v xshift=7 yshift=-1.2 > junkd.v
	vppen junka.v junkb.v junkc.v junkd.v erase=o vpstyle=n scale=0.34 > FIGDIR/twpn.v
	rm back.v bar.v junka.v junkb.v junkc.v junkd.v junk.v junk1.v junk2.v junk3.v junk4.v Junk
	Rm temppos.H tempneg.H tempbar.H gauneg.temp gaupos.temp ttmneg.temp ttmpos.temp

FIGDIR/twlr.v : DATDIR/ttmlft.H DATDIR/ttmrgt.H DATDIR/snapneg.H 
	Window j1=2 j2=2 n3=1 f3=5 < 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=i 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
        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
	vppen junk1.v xshift=-2 yshift=0 > junka.v
        vppen junk2.v xshift=7 yshift=0 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.34 > FIGDIR/twlr.v
	rm back.v junka.v junkb.v junk1.v junk2.v Junk
	Rm tempneg.H

FIGDIR/twmar1.v : DATDIR/marmvel1.H DATDIR/ttmmar0.H DATDIR/snapmar.H
	Window j1=2 j2=2 < DATDIR/marmvel1.H > marmvel1.temp
        Taplot bias=3000 clip=2500 < marmvel1.temp | Ta2vplot title=" " label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT color=i 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=i 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
	Window j1=2 j2=2 < DATDIR/ttmmar0.H > ttmmar0.temp
        Contour nc=20 c0=0 dc=0.1 < ttmmar0.temp 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 < ttmmar0.temp 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 j1=2 j2=2 n3=1 f3=9 < DATDIR/snapmar.H | Scale dscale=0 axis=123 > snapmar.temp
	Taplot pclip=96 < snapmar.temp | Ta2vplot title=" " label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT wantaxis2=y color=i 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
	vppen junk1.v yshift=7. > junka.v
	vppen junk2.v yshift=-1. > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.51 > FIGDIR/twmar1.v
	rm back.v junka.v junkb.v junk1.v junk2.v bar.v junk.v Junk
	Rm snapmar.temp ttmmar0.temp marmvel1.temp

FIGDIR/twmar2.v : DATDIR/marmvel1.H DATDIR/ttmmar1.H DATDIR/snapmar.H
	Window j1=2 j2=2 < DATDIR/marmvel1.H > marmvel1.temp
        Taplot bias=3000 clip=2500 < marmvel1.temp | Ta2vplot title=" " label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT color=i 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=i 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
	Window j1=2 j2=2 < DATDIR/ttmmar1.H > ttmmar1.temp
        Contour nc=24 c0=0 dc=0.1 < ttmmar1.temp THICK d1num=0 d2num=0 FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=24 c0=0 dc=0.1 < ttmmar1.temp 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 j1=2 j2=2 n3=1 f3=9 < DATDIR/snapmar.H | Scale dscale=0 axis=123 > tempsnap.H
	Taplot pclip=96 < tempsnap.H | Ta2vplot title=" " label1="Depth (m)" label2="Surface (m)" wheretitle=t TLSZFAT wantaxis2=y color=i 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
	vppen junk1.v yshift=7. > junka.v
	vppen junk2.v yshift=-1. > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.51 > FIGDIR/twmar2.v
	rm back.v junka.v junkb.v junk1.v junk2.v bar.v junk.v Junk
	Rm tempsnap.H marmvel1.temp ttmmar1.temp

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

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
