#define LINDIR ./../..
#define FIGDIR LINDIR/Figs/Findif2
#define CHPDIR LINDIR/Findif2
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB
#define FIGLIST vppolar vpgrid1 vpgrid2 ttmgra raygra ttmrec waverec wavettmfie
#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/vppolar.v: vppolar.x
        vppolar.x outtype=0 title=" " out=junk.v > Junk
        vppen junk.v vpstyle=n scale=0.528 > FIGDIR/vppolar.v
        rm junk.v Junk

FIGDIR/vpgrid1.v: vpgrid.x
	vpgrid.x outtype=0 out=junk.v > Junk
	vppen junk.v vpstyle=n scale=0.493 > FIGDIR/vpgrid1.v 
	rm junk.v Junk

FIGDIR/vpgrid2.v: vpgrid.x
	vpgrid.x outtype=1 out=junk.v > Junk
	vppen junk.v vpstyle=n scale=0.493 > FIGDIR/vpgrid2.v
	rm junk.v Junk

FIGDIR/ttmgra.v : DATDIR/ttm_gra.H DATDIR/strvel_gra.H DATDIR/errttm_gra.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=3. clip=1.4 | 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=4.01 > tempbar.H
	Taplot < tempbar.H gpow=1 bias=3. clip=1.4 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.5 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/ttm_gra.H > ttm_gra.temp
        Contour nc=20 c0=0 dc=0.1 < ttm_gra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.1 < ttm_gra.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/errttm_gra.H > errttm_gra.temp
        Taplot < errttm_gra.temp bias=0.0005 clip=0.0007 | 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.001 > tempbar.H
	Taplot < tempbar.H bias=0.0005 clip=0.0007 | Ta2vplot title=" " label2="Error (s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.00025 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/ttmgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H strvel_gra.temp ttm_gra.temp errttm_gra.temp

FIGDIR/raygra.v : DATDIR/ray_gra.H DATDIR/strvel_gra.H DATDIR/errray_gra.H
	Window j1=2 j2=2 < DATDIR/strvel_gra.H > strvel_gra.temp
        Taplot < strvel_gra.temp gpow=1 bias=3. clip=1.4 | 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=4.01 > tempbar.H
	Taplot < tempbar.H gpow=1 bias=3. clip=1.4 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.5 o2num=2. grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3. > bar.v
	Window j1=2 j2=2 < DATDIR/ray_gra.H > ray_gra.temp
        Contour nc=39 c0=-1.4835299 dc=0.087266463 < ray_gra.temp THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=39 c0=-1.4835299 dc=0.087266463 < ray_gra.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/errray_gra.H > errray_gra.temp
        Taplot clip=0.00896 < errray_gra.temp | 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.4 maxval=0.4 > tempbar.H
	Taplot < tempbar.H clip=0.55 | Ta2vplot title=" " label2="Error (degrees)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.2 o2num=-0.4 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 > 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/raygra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H strvel_gra.temp ray_gra.temp errray_gra.temp

FIGDIR/ttmrec.v: DATDIR/ttm_rec.H DATDIR/strvel_rec.H
	Window j1=2 j2=2 < DATDIR/strvel_rec.H > strvel_rec.temp
        Taplot < strvel_rec.temp gpow=1 bias=4.0 pclip=100 | 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=junk.v > Junk
        vppen junk.v yshift=1. > back.v
	Bar n2=200 type=h minval=0.99 maxval=5.01 > tempbar.H
 	Taplot < tempbar.H gpow=1 bias=4.0 pclip=100 | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=1. o2num=1. grid=n polarity=1 xinch=7.5 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-3.2 > bar.v
	Window j1=2 j2=2 < DATDIR/ttm_rec.H > ttm_rec.temp
        Contour nc=41 c0=0 dc=0.05 < ttm_rec.temp THICK FAST pad=n xinch=7.5 yinch=7.5 out=junk.v > Junk
        vppen junk.v yshift=1. > junk1.v
        Contour nc=41 c0=0 dc=0.05 < ttm_rec.temp THIN FAST pad=n xinch=7.5 yinch=7.5 out=junk.v > Junk
        vppen junk.v yshift=1. > junk2.v
        vppen back.v bar.v junk1.v junk2.v erase=o vpstyle=n scale=0.698 > FIGDIR/ttmrec.v
	rm back.v bar.v junk.v junk1.v junk2.v Junk
	Rm tempbar.H strvel_rec.temp ttm_rec.temp

FIGDIR/waverec.v : DATDIR/ttm_rec.H DATDIR/snap_rec.H
	Window j1=2 j2=2 n3=1 f3=3 < DATDIR/snap_rec.H | Scale dscale=0 axis=123 > temprec.H
        Taplot < temprec.H | Ta2vplot title="t=0.395 second" label1="Depth (km)" label2="Surface (km)" wheretitle=b 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
        Contour nc=1 c0=0.395 < DATDIR/ttm_rec.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.395 < DATDIR/ttm_rec.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
	Window j1=2 j2=2 n3=1 f3=4 < DATDIR/snap_rec.H | Scale dscale=0 axis=123 > temprec.H
        Taplot < temprec.H | Ta2vplot title="t=0.494 second" label1="Depth (km)" label2="Surface (km)" wheretitle=b 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
        Contour nc=1 c0=0.494 < DATDIR/ttm_rec.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.494 < DATDIR/ttm_rec.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 n3=1 f3=5 < DATDIR/snap_rec.H | Scale dscale=0 axis=123 > temprec.H
        Taplot < temprec.H | Ta2vplot title="t=0.593 second" label1="Depth (km)" label2="Surface (km)" wheretitle=b 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
        Contour nc=1 c0=0.593 < DATDIR/ttm_rec.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.593 < DATDIR/ttm_rec.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 > junk3.v
	Window j1=2 j2=2 n3=1 f3=6 < DATDIR/snap_rec.H | Scale dscale=0 axis=123 > temprec.H
        Taplot < temprec.H | Ta2vplot title="t=0.692 second" label1="Depth (km)" label2="Surface (km)" wheretitle=b 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
        Contour nc=1 c0=0.692 < DATDIR/ttm_rec.H THICK d1num=0 d2num=0 FAST pad=n xinch=7.5 yinch=7.5 out=junka.v > Junk
        Contour nc=1 c0=0.692 < DATDIR/ttm_rec.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=9.6 > junka.v
        vppen junk2.v xshift=7 yshift=9.6 > junkb.v
	vppen junk3.v xshift=-2 yshift=-0.6 > junkc.v
        vppen junk4.v xshift=7 yshift=-0.6 > junkd.v
	vppen junka.v junkb.v junkc.v junkd.v erase=o vpstyle=n scale=0.34 > FIGDIR/waverec.v
	rm back.v junka.v junkb.v junkc.v junkd.v junk1.v junk2.v junk3.v junk4.v Junk
	Rm temprec.H

FIGDIR/wavettmfie.v : DATDIR/ttm_fie.H DATDIR/snap_fie.H DATDIR/strvel_fie.H 
	Window j1=2 j2=2 < DATDIR/strvel_fie.H > strvel_fie.temp
        Taplot < strvel_fie.temp bias=2 clip=1. | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i d1num=1. d2num=1. o1num=0. o2num=0.5 grid=n polarity=1 xinch=10.24 yinch=4.28 out=back.v > Junk
        Bar n2=200 type=h minval=1.459 maxval=2.61 > tempbar.H
	Taplot < tempbar.H bias=2. clip=1. | Ta2vplot title=" " label2="Velocity (km/s)" wherexlabel=b TLSZFAT color=i wantaxis1=n d2num=0.3 o2num=1.46 grid=n polarity=1 xinch=10.24 yinch=0.5 out=junk.v > Junk
        vppen junk.v yshift=-2.8 > bar.v
	Window j1=2 j2=2 < DATDIR/ttm_fie.H > ttm_fie.temp
        Contour nc=16 c0=0. dc=0.1 < ttm_fie.temp THICK FAST pad=n xinch=10.24 yinch=4.28 out=junka.v > Junk
        Contour nc=16 c0=0. dc=0.1 < ttm_fie.temp THIN FAST pad=n xinch=10.24 yinch=4.28 out=junkb.v > Junk
        vp_Overlay back.v bar.v junka.v junkb.v > junk1.v
	Window j1=2 j2=2 f3=9 n3=1 < DATDIR/snap_fie.H > tempsnap.H
	Taplot < tempsnap.H | Ta2vplot title=" " label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=i d1num=1. d2num=1. o1num=0. o2num=0.5 grid=n polarity=1 xinch=10.24 yinch=4.28 out=back.v > Junk
        Contour nc=1 c0=0.999 < DATDIR/ttm_fie.H THICK FAST pad=n xinch=10.24 yinch=4.28 out=junka.v > Junk
        Contour nc=1 c0=0.999 < DATDIR/ttm_fie.H THIN FAST pad=n xinch=10.24 yinch=4.28 out=junkb.v > Junk
        vp_Overlay back.v junka.v junkb.v > junk2.v
	vppen junk1.v yshift=5. > junka.v
	vppen junk2.v yshift=-2.6 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.53 > FIGDIR/wavettmfie.v
	rm back.v bar.v junk.v junk1.v junk2.v junka.v junkb.v Junk
	Rm tempbar.H strvel_fie.temp ttm_fie.temp
	Rm tempsnap.H

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

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
