# Makefile for figures

FIG      = >/dev/null out=
FAT      = axisfat=2 legendfat=2 titlefat=2 labelfat=2 symbolsz=2 plotfat=5
EIKGRLOC =     xll=5     yll=1.4  xinch=6.1364 yinch=7.5 ${FAT}



###########################
##  Laura, October 2008  ##
###########################

LREIKGRPARSV  = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=250 eclip=400   ${EIKGRLOC} ${FAT}
LREIKGRPARSQa = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=10    ${EIKGRLOC} ${FAT}
LREIKGRPARSQb = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=1200  ${EIKGRLOC} ${FAT}
LREIKGRPARSQc = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=600   ${EIKGRLOC} ${FAT}

# Station Maps
${R}/Laura_V_xy.v: Laura_V_xy.H
	Graph < Laura_V_xy.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

${R}/Laura_V_xy_alive.v: Laura_V_xy_alive.H
	Graph < Laura_V_xy_alive.H symbol='o' symbolsz=3 plotcol=7 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

# Final Tomo Maps
${R}/Laura_V_Etom_v.v: Laura_V_Etom.H ${R}/Laura_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_V_Etom.H n3=1 f3=0 | sdrMaskValhall.x value=1000 header=Laura_V_xy_alive.H unit=1000 rmax=125 bclip=255 eclip=395 angmin=180 angrmax=1500 | Grey ${LREIKGRPARSV} ${FIG}${R}/Laurav.v
	Window < Laura_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=250 eclip=400.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Laurav.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Laurav2.v
	Window < Laura_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='P Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Laurav2.v ${R}/Laura_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Laura_V_Etom_q.v: Laura_V_Etom.H ${R}/Laura_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_V_Etom.H n3=1 f3=1 | sdrMaskValhall.x value=1000 header=Laura_V_xy_alive.H unit=1000 rmax=125 bclip=.1 eclip=9.9  angmin=180 angrmax=1500 | Grey ${LREIKGRPARSQa} ${FIG}${R}/Lauraq.v
	Window < Laura_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=10.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Lauraq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Lauraq2.v
	Window < Laura_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='P Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Lauraq2.v ${R}/Laura_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Laura_V_Etom_pd.v: Laura_V_Etom.H ${R}/Laura_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_V_Etom.H n3=1 f3=2 | sdrMaskValhall.x value=1500 header=Laura_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=1190  angmin=180 angrmax=1500 | Grey ${LREIKGRPARSQb} ${FIG}${R}/Lauraq.v
	Window < Laura_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=1200.01  \
		title=' ' color=j plotcol=0 barlabel='accepted picks [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Lauraq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Lauraq2.v
	Window < Laura_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='P Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Lauraq2.v ${R}/Laura_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Laura_V_Etom_ne.v: Laura_V_Etom.H ${R}/Laura_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_V_Etom.H n3=1 f3=3 | sdrMaskValhall.x value=1000 header=Laura_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=595  angmin=180 angrmax=1500 | Grey ${LREIKGRPARSQc} ${FIG}${R}/Lauraq.v
	Window < Laura_V_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=600.01  \
		title=' ' color=j plotcol=0 barlabel='estimate density [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Lauraq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Lauraq2.v
	Window < Laura_V_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='P Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Lauraq2.v ${R}/Laura_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Laura-V-Etom.v: ${R}/Laura_V_Etom_v.v ${R}/Laura_V_Etom_q.v ${R}/Laura_V_Etom_pd.v ${R}/Laura_V_Etom_ne.v
	vppen < ${R}/Laura_V_Etom_v.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/1.v
	vppen < ${R}/Laura_V_Etom_q.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/2.v
	vppen < ${R}/Laura_V_Etom_pd.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/3.v
	vppen < ${R}/Laura_V_Etom_ne.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/4.v
	vp_SideBySideIso ${R}/1.v ${R}/2.v > ${R}/top.v
	vp_SideBySideIso ${R}/3.v ${R}/4.v > ${R}/bot.v
	vp_OverUnderIso  ${R}/top.v ${R}/bot.v > $@



#############################
##  Jianhua, January 2005  ##
#############################

JHEIKGRPARSV  = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=250 eclip=400   ${EIKGRLOC} ${FAT}
JHEIKGRPARSQa = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=10    ${EIKGRLOC} ${FAT}
JHEIKGRPARSQb = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=1200   ${EIKGRLOC} ${FAT}
JHEIKGRPARSQc = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=600   ${EIKGRLOC} ${FAT}
${R}/Jianhua_V_xy.v: Jianhua_V_xy.H
	Graph < Jianhua_V_xy.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

${R}/Jianhua_V_xy_alive.v: Jianhua_V_xy_alive.H
	Graph < Jianhua_V_xy_alive.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2


# Final Tomo Maps
${R}/Jianhua_V_Etom_v.v: Jianhua_V_Etom.H ${R}/Jianhua_V_xy_alive.v Jianhua_V_xy_alive.H
	Window < Jianhua_V_Etom.H n3=1 f3=0 | sdrMaskValhall.x value=1000 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 bclip=255 eclip=395 angmin=180 angrmax=1500 | Grey ${JHEIKGRPARSV} ${FIG}${R}/Jianhuav.v
	Window < Jianhua_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=250 eclip=400.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Jianhuav.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Jianhuav2.v
	Window < Jianhua_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Jianhuav2.v ${R}/Jianhua_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Jianhua_V_Etom_q.v: Jianhua_V_Etom.H ${R}/Jianhua_V_xy_alive.v Jianhua_V_xy_alive.H
	Window < Jianhua_V_Etom.H n3=1 f3=1 | sdrMaskValhall.x value=1000 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 bclip=.1 eclip=9.9  angmin=180 angrmax=1500 | Grey ${JHEIKGRPARSQa} ${FIG}${R}/Jianhuaq.v
	Window < Jianhua_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=10.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Jianhuaq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Jianhuaq2.v
	Window < Jianhua_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Jianhuaq2.v ${R}/Jianhua_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Jianhua_V_Etom_pd.v: Jianhua_V_Etom.H ${R}/Jianhua_V_xy_alive.v Jianhua_V_xy_alive.H
	Window < Jianhua_V_Etom.H n3=1 f3=2 | sdrMaskValhall.x value=1500 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=1190  angmin=180 angrmax=1500 | Grey ${JHEIKGRPARSQb} ${FIG}${R}/Jianhuaq.v
	Window < Jianhua_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=1200.01  \
		title=' ' color=j plotcol=0 barlabel='accepted picks [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Jianhuaq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Jianhuaq2.v
	Window < Jianhua_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Jianhuaq2.v ${R}/Jianhua_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Jianhua_V_Etom_ne.v: Jianhua_V_Etom.H ${R}/Jianhua_V_xy_alive.v Jianhua_V_xy_alive.H
	Window < Jianhua_V_Etom.H n3=1 f3=3 | sdrMaskValhall.x value=1000 header=Jianhua_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=595  angmin=180 angrmax=1500 | Grey ${JHEIKGRPARSQc} ${FIG}${R}/Jianhuaq.v
	Window < Jianhua_V_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=600.01  \
		title=' ' color=j plotcol=0 barlabel='estimate density [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Jianhuaq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Jianhuaq2.v
	Window < Jianhua_V3_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Jianhuaq2.v ${R}/Jianhua_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Jianhua-V-Etom.v: ${R}/Jianhua_V_Etom_v.v ${R}/Jianhua_V_Etom_q.v ${R}/Jianhua_V_Etom_pd.v ${R}/Jianhua_V_Etom_ne.v
	vppen < ${R}/Jianhua_V_Etom_v.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/1.v
	vppen < ${R}/Jianhua_V_Etom_q.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/2.v
	vppen < ${R}/Jianhua_V_Etom_pd.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/3.v
	vppen < ${R}/Jianhua_V_Etom_ne.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/4.v
	vp_SideBySideIso ${R}/1.v ${R}/2.v > ${R}/top.v
	vp_SideBySideIso ${R}/3.v ${R}/4.v > ${R}/bot.v
	vp_OverUnderIso  ${R}/top.v ${R}/bot.v > $@




#############################
##  Artman, February 2004  ##
#############################

AREIKGRPARSV  = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=250 eclip=400   ${EIKGRLOC} ${FAT}
AREIKGRPARSQa = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=10    ${EIKGRLOC} ${FAT}
AREIKGRPARSQb = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=1200  ${EIKGRLOC} ${FAT}
AREIKGRPARSQc = color=j newclip=1 min1=521000 max1=530000 min2=6231000 max2=6242000 yreverse=0 xreverse=0 transp=0 title=' '  wantaxis=0   bclip=0   eclip=600   ${EIKGRLOC} ${FAT}


${R}/Artman_V_xy.v: Artman_V_xy.H
	Graph < Artman_xy.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

${R}/Artman_V_xy_alive.v: Artman_V_xy_alive.H
	Graph < Artman_V_xy_alive.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

# Final Tomo Maps
${R}/Artman_V_Etom_v.v: Artman_V_Etom.H ${R}/Artman_V_xy_alive.v Artman_V_xy_alive.H
	Window < Artman_V_Etom.H n3=1 f3=0 | sdrMaskValhall.x value=1000 header=Artman_V_xy_alive.H unit=1000 rmax=125 bclip=255 eclip=395 angmin=180 angrmax=1500 | Grey ${AREIKGRPARSV} ${FIG}${R}/Artmanv.v
	Window < Artman_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=250 eclip=400.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Artmanv.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Artmanv2.v
	Window < Artman_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Artmanv2.v ${R}/Artman_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Artman_V_Etom_q.v: Artman_V_Etom.H ${R}/Artman_V_xy_alive.v Artman_V_xy_alive.H
	Window < Artman_V_Etom.H n3=1 f3=1 | sdrMaskValhall.x value=1000 header=Artman_V_xy_alive.H unit=1000 rmax=125 bclip=.1 eclip=9.9  angmin=180 angrmax=1500 | Grey ${AREIKGRPARSQa} ${FIG}${R}/Artmanq.v
	Window < Artman_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=10.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Artmanq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Artmanq2.v
	Window < Artman_V_Etom.H n3=1 f3=1 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Artmanq2.v ${R}/Artman_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Artman_V_Etom_pd.v: Artman_V_Etom.H ${R}/Artman_V_xy_alive.v Artman_V_xy_alive.H
	Window < Artman_V_Etom.H n3=1 f3=2 | sdrMaskValhall.x value=1500 header=Artman_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=1190  angmin=180 angrmax=1500 | Grey ${AREIKGRPARSQb} ${FIG}${R}/Artmanq.v
	Window < Artman_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=1200.01  \
		title=' ' color=j plotcol=0 barlabel='accepted picks [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Artmanq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Artmanq2.v
	Window < Artman_V_Etom.H n3=1 f3=2 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Artmanq2.v ${R}/Artman_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Artman_V_Etom_ne.v: Artman_V_Etom.H ${R}/Artman_V_xy_alive.v Artman_V_xy_alive.H
	Window < Artman_V_Etom.H n3=1 f3=3 | sdrMaskValhall.x value=1000 header=Artman_V_xy_alive.H unit=1000 rmax=125 bclip=5 eclip=595  angmin=180 angrmax=1500 | Grey ${AREIKGRPARSQc} ${FIG}${R}/Artmanq.v
	Window < Artman_V_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=600.01  \
		title=' ' color=j plotcol=0 barlabel='estimate density [-]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Artmanq.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Artmanq2.v
	Window < Artman_V_Etom.H n3=1 f3=3 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Artmanq2.v ${R}/Artman_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

${R}/Artman-V-Etom.v: ${R}/Artman_V_Etom_v.v ${R}/Artman_V_Etom_q.v ${R}/Artman_V_Etom_pd.v ${R}/Artman_V_Etom_ne.v
	vppen < ${R}/Artman_V_Etom_v.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/1.v
	vppen < ${R}/Artman_V_Etom_q.v  xscale=.675 yscale=.675 xshift=-2.1 yshift=0   > ${R}/2.v
	vppen < ${R}/Artman_V_Etom_pd.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/3.v
	vppen < ${R}/Artman_V_Etom_ne.v xscale=.675 yscale=.675 xshift=-2.1 yshift=3.9 > ${R}/4.v
	vp_SideBySideIso ${R}/1.v ${R}/2.v > ${R}/top.v
	vp_SideBySideIso ${R}/3.v ${R}/4.v > ${R}/bot.v
	vp_OverUnderIso  ${R}/top.v ${R}/bot.v > $@



##############################
##  Example shot for Laura  ##
##############################
${R}/Laura_SHOT_V_xy_alive.v: Laura_V_xy_alive.H
	Graph < Laura_V_xy_alive.H symbol='o' symbolsz=3 plotcol=8 min1=521 max1=530 min2=6231 max2=6242 \
		title=' ' label1='UTM easting [km]' label2='UTM northing [km]' d2num=2 d1num=2 grid=1 ${EIKGRLOC} ${FIG}$@ plotfat=2

# Interpolated surface
${R}/Laura_SHOT_V_tsurf.v: Laura_SHOT_V_tsurf.H ${R}/Laura_SHOT_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_SHOT_V_tsurf.H | Grey ${LREIKGRPARSQa} ${FIG}${R}/Laurav.v
	Window < Laura_SHOT_V_tsurf.H | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=0   eclip=10.01  \
		title=' ' color=j plotcol=0 barlabel='Travel time [s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Laurav.v | sed -e 's/^C 256 0 0 128$$/C 256 255 255 255/' | plas | vppen >${R}/Laurav2.v
	#pldb  < ${R}/Laurav.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Laurav2.v
	Window < Artman_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Laurav2.v ${R}/Laura_SHOT_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@

# Velocity estimations
${R}/Laura_SHOT_V_Etom_v.v: Laura_SHOT_V_Etom.H ${R}/Laura_SHOT_V_xy_alive.v Laura_V_xy_alive.H
	Window < Laura_SHOT_V_Etom.H n3=1 f3=0  | Grey ${LREIKGRPARSV} ${FIG}${R}/Laurav.v
	Window < Laura_SHOT_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=1 wantaxis=0 wantframe=n newclip=1 bclip=250 eclip=400.01  \
		title=' ' color=j plotcol=0 barlabel='Velocity [m/s]' >/dev/null out=${R}/bar_.v
	vppen < ${R}/bar_.v xwmin=11.25 vpstyle=n yscale=.98 vshift=.13 > ${R}/bar.v
	pldb  < ${R}/Laurav.v | sed -e 's/^C 256 0 0 128$$/C 256 255 255 255/' | plas | vppen >${R}/Laurav2.v
	#pldb  < ${R}/Laurav.v | sed -e 's/^C 511 131 0 0$$/C 511 255 255 255/' | plas | vppen >${R}/Laurav2.v
	Window < Artman_V_Etom.H n3=1 f3=0 | Grey ${FAT} wantscalebar=0 wantaxis=0 wantframe=n newclip=1 bclip=1.76 eclip=2.0   color=j plotcol=0   title=' ' wheretitle=t \
		barlabel='Velocity [km/s]' >/dev/null out=${R}/title_.v; vppen < ${R}/title_.v ywmin=9 vpstyle=n xshift=.8 vshift=-.5 > ${R}/title.v
	vp_Overlay ${R}/Laurav2.v ${R}/Laura_SHOT_V_xy_alive.v ${R}/title.v ${R}/bar.v > $@



# Noise
${R}/Laura_V_bpC.v: Laura_Z3d_bpC.H
	Window < $< f2=117 n2=96  > part0.H 
	Pad    < part0.H end2=7   > part1.H ; echo 'hff=-1' >> part1.H
	Window < $< f2=216 n2=118 > part2.H ; echo 'hff=-1' >> part2.H
	Cat part1.H part2.H axis=2 | sdrMute.x kill=19 > temp_snA3.H
	echo "o2=0 d2=50 label1='time [s]'" >> temp_snA3.H
	echo "o3=0 d3=1  label2='g_x [m]"   >> temp_snA3.H
	Window < temp_snA3.H n1=600 n3=1  | Grey o2=0 d2=50 max2=11000 max1=30 o1=0 plotfat=2 axisfat=2 legendfat=2 titlefat=2 labelfat=2 \
		${FAT} ${FIG}$@ title=' ' label1='Time [s]' label2='Approximate Distance [m]'
	Rm part1.H; Rm part2.H; Rm temp_snA3.H

# Virtual Shot (at 190)
${R}/Laura_SHOT_V_bpC.v: Laura_Z3d_bpC.H
	Window < $< f2=117 n2=96  > part0.H 
	Pad    < part0.H end2=7   > part1.H ; echo 'hff=-1' >> part1.H
	Window < $< f2=216 n2=118 > part2.H ; echo 'hff=-1' >> part2.H
	Cat part1.H part2.H axis=2 | sdrMute.x kill=19 > temp_snA3.H
	echo "o2=0 d2=50 label1='Time [s]'" >> temp_snA3.H
	echo "o2=-9500 d2=50  label2='Approximate Offset [m]'"   >> temp_snA3.H
	sdrSI2d_self1d_t_omp.x < temp_snA3.H  ns=1 fs=190 js=1 ntau=512 sym_t=1    > temp_snA4.H
	Window <  temp_snA4.H f1=511 n1=302 n3=1  | Grey min1=0 max1=30  plotfat=2 axisfat=2 legendfat=2 titlefat=2 labelfat=2 \
		${FAT} ${FIG}$@ title=' '

PICKFIGPARS2 = min2=0 max2=30 yreverse=1 label1=' ' label2=' ' title=' ' wantaxis=1 labelsz=11 ${FAT}
# Virtual, picked, shot
${R}/Laura_SHOT_picked_V_bpC.v: Laura_Z3d_bpC.H #Laura_V_SNR.H
	Window < Laura_Z3d_bpC.H f2=117 n2=96  > part0.H 
	Pad    < part0.H end2=7   > part1.H ; echo 'hff=-1' >> part1.H
	Window < Laura_Z3d_bpC.H f2=216 n2=118 > part2.H ; echo 'hff=-1' >> part2.H
	Cat part1.H part2.H axis=2 | sdrMute.x kill=19 > temp_snA3.H

	Window < Laura_Z3d_bpC.H@@ f2=117 n2=96  > part0.H 
	Pad    < part0.H end2=7   > part1.H ; echo 'hff=-1' >> part1.H
	Window < Laura_Z3d_bpC.H@@ f2=216 n2=118 > part2.H ; echo 'hff=-1' >> part2.H
	Cat part1.H part2.H axis=2 | sdrMute.x kill=19 > temp_snA3.H@@

	echo "o2=0 d2=50 label1='Time [s]'" >> temp_snA3.H
	echo "o2=-9500 d2=50  label2='Approximate Offset [m]'"   >> temp_snA3.H
	sdrSI2d_self1d_t_omp.x < temp_snA3.H  ns=1 fs=190 js=1 ntau=512 sym_t=1    | Envelope > temp_snA4.H
	sdrSIexpand_headers.x  < temp_snA3.H@@  ns=1 fs=190 js=1 rx=4 ry=5 sx=1 sy=2  > temp_snA4.H@@; echo "hff=temp_snA4.H@@" >>temp_snA4.H

	Window <  temp_snA4.H f1=511 n1=302 n3=1  | Envelope | Grey min1=0 max1=30 plotfat=2 axisfat=2 legendfat=2 titlefat=2 labelfat=2 \
		${FAT} ${FIG}${R}/Laura_shot_env.v title=' ' allpos=1

	SNRcalc.x < temp_snA4.H header=temp_snA4.H@@ v=350 dtau=2  > temp_snA3_snr.H

	Window3d <  temp_snA3_snr.H n1=1 f1=0  > pickst.H synch=1;  echo "hff= " >> pickst.H;
	Window3d <  temp_snA3_snr.H n1=2 f1=10 >  windt.H synch=1;  echo "hff= " >> windt.H;  Transp <  windt.H >  wind.H; echo "o2=1" >>  wind.H
	echo "o2=0 d1=50 label1='Time [s]'" >> pickst.H
	echo "o1=-9500 d1=50  label2='Approximate Offset [m]'"   >> pickst.H
	echo "o1=0 d1=50 label1='Time [s]'" >> wind.H
	echo "o1=-9500 d1=50  label2='Approximate Offset [m]'"   >> wind.H
	Graph  < pickst.H   ${FIG}${R}/Laura_shot_picks.v  min1=-9500 max1=1550 min2=0 max2=30 yreverse=1 label1=' ' label2=' ' title=' ' wantaxis=0 plotcol=5 plotfat=10 symbol='+' symbolsz=2
	Graph  < wind.H     ${FIG}${R}/Laura_shot_wind.v   min1=-9500 max1=1550 min2=0 max2=30 yreverse=1 label1=' ' label2=' ' title=' ' wantaxis=0 plotcol=6 plotfat=10
	vp_Overlay ${R}/Laura_shot_env.v ${R}/Laura_shot_picks.v ${R}/Laura_shot_wind.v > $@

${R}/Laura-V-picked.v: Fig/Laura_V_bpC.v Fig/Laura_SHOT_V_bpC.v Fig/Laura_SHOT_picked_V_bpC.v
	vp_OverUnderAniso Fig/Laura_V_bpC.v Fig/Laura_SHOT_V_bpC.v Fig/Laura_SHOT_picked_V_bpC.v > $@

${R}/Valhall-V-Etom.v: ${R}/Laura-V-Etom.v ${R}/Jianhua-V-Etom.v ${R}/Artman-V-Etom.v
	cat ${R}/Laura-V-Etom.v ${R}/Jianhua-V-Etom.v ${R}/Artman-V-Etom.v > $@

# Scatter plot of SNR
${R}/snr-scat.v: Artman_V_SNR.H
	Window < Artman_V_SNR.H n1=1 f1=1 > snr.H
	Window < Artman_V_SNR.H n1=1 f1=6 > rad.H
	echo "n1=5225791 n2=1" >> snr.H
	echo "n1=5225791 n2=1" >> rad.H
	sdrLog.x < snr.H > snrl.H
	Cmplx rad.H snrl.H > snr.H
	Window < snr.H j1=100 | Graph min2=-10 max2=40   title=' ' label1='radius [m]' label2='Log\_10\^(SNR)' symbol=. ${FIG}${R}/temp.v
	sdrTapers.x o1=-10 n1=20 scale=0.3011,1,1 type=2,2,2 t1=-10 t2=10 t3=10 t4=10 | \
	    Graph min2=-10 max2=40  wantaxis=0 title=' ' ${FIG}${R}/snrbnd.v plotcol=5 plotfat=5
	vp_Overlay ${R}/temp.v ${R}/snrbnd.v > $@





