#define LINDIR ./../..
#define FIGDIR LINDIR/Figs/Resmig4
#define CHPDIR LINDIR/Resmig4
#define SRCDIR CHPDIR/Src
#define DATDIR CHPDIR/Dat
#define BINDIR .
#define USRLIBS VPLOTLIB
#define FIGLIST vprmig velsyn datspk imgspk datsyn imgsyn
#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/vprmig.v: vprmig.x
	vprmig.x out=junk.v > Junk
	vppen junk.v vpstyle=n scale=0.79 > FIGDIR/vprmig.v
	rm junk.v Junk

FIGDIR/velsyn.v : DATDIR/strvel_syn.H DATDIR/str_syn.H
        Taplot < DATDIR/strvel_syn.H gpow=1 bias=2.7 clip=1. | 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 polarity=1 out=back.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 out=junk.v > Junk
	vp_Overlay back.v junk.v > junka.v
	Bar n2=200 type=h minval=1.89 maxval=3.01 > tempbar.H
	Taplot < tempbar.H gpow=1 bias=2.7 clip=1. | 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=junkb.v > Junk
	vppen junka.v yshift=-0.3 > junk1.v
        vppen junkb.v yshift=-3.5 > junk2.v
	vppen junk1.v junk2.v erase=o vpstyle=n scale=0.52 > FIGDIR/velsyn.v 
	rm back.v junk.v junka.v junkb.v junk1.v junk2.v Junk
	Rm tempbar.H

FIGDIR/datsyn.v : DATDIR/syncos.H DATDIR/synimg0.H
	Window j1=2 j2=2 < DATDIR/syncos.H > syncos.temp
        Taplot < syncos.temp | Ta2vplot title=" " label1="Time (s)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=junk1.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THICK min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junka.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junkb.v > Junk
	vp_Overlay junka.v junkb.v > junk.v
	Window j1=2 j2=2 < DATDIR/synimg0.H > synimg0.temp
        Taplot < synimg0.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 polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junk2.v
	vppen junk1.v yshift=5 > junka.v
        vppen junk2.v yshift=-2.2 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.52 > FIGDIR/datsyn.v 
	rm back.v junk.v junka.v junkb.v junk1.v junk2.v Junk
	Rm syncos.temp synimg0.temp

FIGDIR/imgsyn.v : DATDIR/synimg.H DATDIR/str_syn.H DATDIR/synrmg.H 
	Graph < DATDIR/str_syn.H FAST yreverse=y THICK min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junka.v > Junk
	Graph < DATDIR/str_syn.H FAST yreverse=y THIN min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.9,1.9,1.9,1.9,1.9,1.9,1.9 out=junkb.v > Junk
	vp_Overlay junka.v junkb.v > junk.v
	Window j1=2 j2=2 < DATDIR/synrmg.H > synrmg.temp
        Taplot < synrmg.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 polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junka.v
	Window j1=2 j2=2 < DATDIR/synimg.H > synimg.temp
        Taplot < synimg.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 polarity=1 out=back.v > Junk
	vp_Overlay back.v junk.v > junkb.v
	vppen junka.v yshift=5 > junk1.v
        vppen junkb.v yshift=-2.2 > junk2.v
	vppen junk1.v junk2.v erase=o vpstyle=n scale=0.52 > FIGDIR/imgsyn.v 
	rm back.v junk.v junka.v junkb.v junk1.v junk2.v Junk
	Rm synrmg.temp synimg.temp

FIGDIR/datspk.v : DATDIR/spkcos.H DATDIR/spkimg0.H DATDIR/synoper.H
	Window j1=2 j2=2 < DATDIR/spkcos.H > spkcos.temp
        Taplot < spkcos.temp | Ta2vplot title=" " label1="Time (s)" label2="Surface (km)" TLSZFAT color=i xinch=10.24 yinch=5.12 d1num=1 d2num=1 o1num=0 o2num=0 grid=n polarity=1 out=junk1.v > Junk
	Window j1=2 j2=2 < DATDIR/spkimg0.H > spkimg0.temp
        Taplot pclip=99.8 < spkimg0.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 polarity=1 out=back.v > Junk
	Graph < DATDIR/synoper.H FAST yreverse=y THICK min1=0. max1=4. min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.5,1.5,1.5,1.5,1.5,1.5,1.5 out=junka.v > Junk
	Graph < DATDIR/synoper.H FAST yreverse=y THIN min1=0. max1=4. min2=0. max2=2. pad=n xinch=10.24 yinch=5.12 dash=1.5,1.5,1.5,1.5,1.5,1.5,1.5 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.2 > junkb.v
	vppen junka.v junkb.v erase=o vpstyle=n scale=0.52 > FIGDIR/datspk.v 
	rm back.v junka.v junkb.v junk1.v junk2.v Junk
	Rm spkcos.temp spkimg0.temp

FIGDIR/imgspk.v : DATDIR/spkimg0.H DATDIR/str_syn.H DATDIR/spkrmg.H 
	Window j1=2 j2=2 < DATDIR/spkrmg.H > spkrmg.temp
        Taplot pclip=99.8 < spkrmg.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 polarity=1 out=back.v > Junk
	vp_Overlay back.v > junka.v
	Window j1=2 j2=2 < DATDIR/spkimg.H > spkimg.temp
        Taplot pclip=99.8 < spkimg.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 polarity=1 out=back.v > Junk
	vp_Overlay back.v > junkb.v
	vppen junka.v yshift=5 > junk1.v
        vppen junkb.v yshift=-2.2 > junk2.v
	vppen junk1.v junk2.v erase=o vpstyle=n scale=0.52 > FIGDIR/imgspk.v 
	rm back.v junka.v junkb.v junk1.v junk2.v Junk
	Rm spkrmg.temp spkimg.temp

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

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

#include <SEP.obj.rules>

#include <SEP.prog.rules>
