#define LINDIR ./../..
#define SLDDIR LINDIR/Slds/Findif2
#define CHPDIR LINDIR/Findif2
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB
#define SLDLIST 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>

SLDDIR/vppolar.v: vppolar.x
	vppolar.x outtype=0 title="Polar coordinates" out=junk.v > Junk
	vppen junk.v erase=o vpstyle=n scale=1. > SLDDIR/vppolar.v 
	rm junk.v Junk

SLDDIR/vpgrid1.v: vpgrid.x 
	vpgrid.x outtype=0 title="Grid and Notation" out=junk.v > Junk
	vppen junk.v erase=o vpstyle=n scale=1. > SLDDIR/vpgrid1.v 
	rm junk.v Junk

SLDDIR/vpgrid2.v: vpgrid.x
	vpgrid.x outtype=1 title="Grid and Notation" out=junk.v > Junk
	vppen junk.v erase=o vpstyle=n scale=1. > SLDDIR/vpgrid2.v
	rm junk.v Junk

SLDDIR/ttmgra.v : DATDIR/ttm_gra.H DATDIR/strvel_gra.H DATDIR/anattm_gra.H DATDIR/errttm_gra.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=3. clip=1.4 | 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 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=G 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
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttm_gra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=20 c0=0 dc=0.1 < DATDIR/ttm_gra.H 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
        Taplot < DATDIR/errttm_gra.H bias=0.0005 clip=0.0007 | 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 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=G 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
	vp_Movie junk1.v junk2.v > SLDDIR/ttmgra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H

SLDDIR/raygra.v : DATDIR/ray_gra.H DATDIR/strvel_gra.H DATDIR/anaray_gra.H DATDIR/errray_gra.H
        Taplot < DATDIR/strvel_gra.H gpow=1 bias=3. clip=1.4 | 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 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=G 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
        Contour nc=39 c0=-1.4835299 dc=0.087266463 < DATDIR/ray_gra.H THICK FAST pad=n xinch=10.24 yinch=5.12 out=junka.v > Junk
        Contour nc=39 c0=-1.4835299 dc=0.087266463 < DATDIR/ray_gra.H 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
        Taplot clip=0.00896 < DATDIR/errray_gra.H | 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 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=G 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.0 > bar.v
	vp_Overlay back.v bar.v > junk2.v
	vp_Movie junk1.v junk2.v > SLDDIR/raygra.v
	rm back.v bar.v junka.v junkb.v junk1.v junk2.v junk.v Junk
	Rm tempbar.H

SLDDIR/ttmrec.v: DATDIR/ttm_rec.H DATDIR/strvel_rec.H
        Taplot < DATDIR/strvel_rec.H gpow=1 bias=4.0 pclip=100 | Ta2vplot title="Traveltime field" 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=junk.v > Junk
        vppen junk.v yshift=0.8 > 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=G 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.4 > bar.v
        Contour nc=41 c0=0 dc=0.05 < DATDIR/ttm_rec.H THICK FAST pad=n xinch=7.5 yinch=7.5 out=junk.v > Junk
        vppen junk.v yshift=0.8 > junk1.v
        Contour nc=41 c0=0 dc=0.05 < DATDIR/ttm_rec.H THIN FAST pad=n xinch=7.5 yinch=7.5 out=junk.v > Junk
        vppen junk.v yshift=0.8 > junk2.v
        vp_Overlay back.v bar.v junk1.v junk2.v > SLDDIR/ttmrec.v
	rm back.v bar.v junk.v junk1.v junk2.v Junk
	Rm tempbar.H

SLDDIR/waverec.v : DATDIR/ttm_rec.H DATDIR/snap_rec.H
	Window 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=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.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 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=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.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 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=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.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 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=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.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
	vp_Movie junk1.v junk2.v junk3.v junk4.v > SLDDIR/waverec.v
	rm back.v junka.v junkb.v junk1.v junk2.v junk3.v junk4.v Junk
	Rm temprec.H

SLDDIR/wavettmfie.v : DATDIR/ttm_fie.H DATDIR/snap_fie.H DATDIR/strvel_fie.H 
        Taplot < DATDIR/strvel_fie.H bias=2 clip=1. | Ta2vplot title="Traveltime field" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G 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=G 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
        Contour nc=16 c0=0. dc=0.1 < DATDIR/ttm_fie.H THICK FAST pad=n xinch=10.24 yinch=4.28 out=junka.v > Junk
        Contour nc=16 c0=0. dc=0.1 < DATDIR/ttm_fie.H 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 f3=9 n3=1 < DATDIR/snap_fie.H > tempsnap.H
	Taplot < tempsnap.H | Ta2vplot title="Wave field" label1="Depth (km)" label2="Surface (km)" wheretitle=t TLSZFAT color=G 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
	vp_Movie junk1.v junk2.v > SLDDIR/wavettmfie.v
	rm back.v bar.v junk.v junk1.v junk2.v junka.v junkb.v Junk
	Rm tempbar.H
	Rm tempsnap.H

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

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
