
UPSTEXFLAGS = color=n fat=1 fatmult=1.5 invras=y force=y
#RESDIR = /net/server/biondo/Fig/2011/SEP143
RESDIR = ../../Fig_test


ifndef OVER
OVER=n
endif

#TEMPVP=/net/glad/scr1/biondo/SCRATCH/VP-
TEMPVP=${TEMPDIR}/VP_

ifeq ($(OVER),y)  
ERASE=y
else

STACK-TITL="a)"
STACK-XSHIFT=0
STACK-YSHIFT=0

MIG-X0-TITL="b)"
MIG-X550-TITL="c)"
MIG-X1100-TITL="d)"
MIG-X0-XSHIFT=6
MIG-X550-XSHIFT=8
MIG-X1100-XSHIFT=10
MIG-YSHIFT=0

CIG-0-TITL="a)"
CIG-0-WANTAXIS1=y
CIG-0-XSHIFT=1

CIG-0-VLowFreq-TITL="b)"
CIG-0-VLowFreq-WANTAXIS1=n
CIG-0-VLowFreq-XSHIFT=3

CIG-YSHIFT=0

RMO-X0-RHO=1.06
RMO-X0-O1NUM=2020
RMO-X550-RHO=.965
RMO-X550-O1NUM=2000

RMO-TITL="b)"
RMO-WANTAXIS1=n
RMO-XSHIFT=3
RMO-YSHIFT=0

POWER-XSHIFT=0
POWER-X0-XSHIFT=5.4
POWER-X550-XSHIFT=9
POWER-YSHIFT=0
POWER-X0-RHO0=.94
POWER-X0-DRHO=.06
POWER-X550-RHO0=.93
POWER-X550-DRHO=.035

endif

Der-Power-all-over$(OVER).Vp : \
	Der-Power-Ang-Mig-0-$(MYHOST).H \
	Der-Smooth-Power-Ang-Mig-0-$(MYHOST).H  \
	Der-Power-Ang-Mig-0-VLowFreq-$(MYHOST).H 
	Merge axis=2 space=n $^ | \
	Window3d min1=-1600 max1=1600 | \
	Scale axis=1 | \
	Graph grid=y \
	screenht=8 screenratio=.65 \
	o1num=0 d1num=550 grid2=y \
	label1="Horizontal location [m]"  \
	label2="Derivative of stack power"  \
	label2="d Power/d\F10 r" \
	title=" " wheretitle=b wherexlabel=t dash=0,4,2   \
	legend=y legendloc="bl" legendsz=5 legendtitle="Central Freq." \
	curvelabel="25 Hz:25 Hz - Smooth:5 Hz" \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=0 yshift=0 scale=1 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@


Smooth-Power-all-%-over$(OVER).Vp :  \
	Smooth-Power-Ang-Mig-X550-%.Vp \
	Smooth-Power-Ang-Mig-X0-%.Vp \
	Smooth-Power-Ang-Mig-%.Vp
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@

Avg-Power-all-%-over$(OVER).Vp :  \
	Avg-Power-Ang-Mig-X550-%.Vp \
	Avg-Power-Ang-Mig-X0-%.Vp \
	Avg-Power-Ang-Mig-%.Vp
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@


Smooth-Power-Ang-Mig-X550-%.Vp: Smooth-Power-Ang-Mig-%-$(MYHOST).H 
	<Smooth-Power-Ang-Mig-$*-$(MYHOST).H Window3d n2=1 min2=550 | \
	Scale axis=1 |\
	Graph grid=y \
	screenht=8 screenratio=1.05 \
	o1num=$(POWER-X550-RHO0) d1num=$(POWER-X550-DRHO) d2num=.25 min2=0 \
	label1="Residual moveout parameter [s\_mig\^/s]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Stack power"  \
	title="c)" wheretitle=b wherexlabel=t  wantaxis2=n \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-X550-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Smooth-Power-Ang-Mig-X0-%.Vp: Smooth-Power-Ang-Mig-%-$(MYHOST).H 
	<Smooth-Power-Ang-Mig-$*-$(MYHOST).H Window3d n2=1 min2=0 | \
	Scale axis=1 |\
	Graph grid=y \
	screenht=8 screenratio=1.05 \
	o1num=$(POWER-X0-RHO0) d1num=$(POWER-X0-DRHO) d2num=.25 min2=0 \
	label1="Residual moveout parameter [s/s\_mig\^]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Stack power"  \
	title="b)" wheretitle=b wherexlabel=t  wantaxis2=y \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-X0-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Smooth-Power-Ang-Mig-%.Vp: Smooth-Power-Ang-Mig-%-$(MYHOST).H 
	<Smooth-Power-Ang-Mig-$*-$(MYHOST).H  \
	Grey pclip=100 grid=y color=i allpos=y \
	screenht=8 screenratio=.8 \
	o2num=0 d2num=550 min2=-1600 max2=1600 \
	label1="Residual moveout parameter [s/s\_mig\^]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Horizontal location [m]"  \
	title="a)" \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Avg-Power-Ang-Mig-X550-%.Vp: Avg-Power-Ang-Mig-%-$(MYHOST).H 
	<Avg-Power-Ang-Mig-$*-$(MYHOST).H Window3d n2=1 min2=550 | \
	Scale axis=1 |\
	Graph grid=y \
	screenht=8 screenratio=1.05 \
	o1num=$(POWER-X550-RHO0) d1num=$(POWER-X550-DRHO) d2num=.25 min2=0 \
	label1="Residual moveout parameter [s/s\_mig\^]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Stack power"  \
	title="c)" wheretitle=b wherexlabel=t  wantaxis2=n \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-X550-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Avg-Power-Ang-Mig-X0-%.Vp: Avg-Power-Ang-Mig-%-$(MYHOST).H 
	<Avg-Power-Ang-Mig-$*-$(MYHOST).H Window3d n2=1 min2=0 | \
	Scale axis=1 |\
	Graph grid=y \
	screenht=8 screenratio=1.05 \
	o1num=$(POWER-X0-RHO0) d1num=$(POWER-X0-DRHO) d2num=.25 min2=0 \
	label1="Residual moveout parameter [s/s\_mig\^]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Stack power"  \
	title="b)" wheretitle=b wherexlabel=t  wantaxis2=y \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-X0-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@


Avg-Power-Ang-Mig-%.Vp: Avg-Power-Ang-Mig-%-$(MYHOST).H 
	<Avg-Power-Ang-Mig-$*-$(MYHOST).H  \
	Grey pclip=100 grid=y color=i allpos=y \
	screenht=8 screenratio=.8 \
	o2num=0 d2num=550 min2=-1600 max2=1600 \
	label1="Residual moveout parameter [s/s\_mig\^]" \
	label1="Res. mov. par. [s/s\_mig\^]" \
	label1="\F10 r \F0 =[s\_mig\^/s]" \
	label2="Horizontal location [m]"  \
	title="a)" \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(POWER-XSHIFT) yshift=$(POWER-YSHIFT) scale=.55 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@


Rmo-all-X550-over$(OVER).Vp :  \
	Cig-0-X550-over$(OVER).Vp \
	Rmo-0-X550-over$(OVER).Vp  
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@

Rmo-all-X0-over$(OVER).Vp :  \
	Cig-0-X0-over$(OVER).Vp \
	Rmo-0-X0-over$(OVER).Vp  
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@

Rmo-0-X%-over$(OVER).Vp: Rmo-Ang-Mig-0-$(MYHOST).H 
	Window3d <Rmo-Ang-Mig-0-$(MYHOST).H  \
	min1=1800 max1=2200 n2=1 min2=$* n4=1 min4=$(RMO-X$*-RHO) |  \
	Grey pclip=100 color=i \
	screenwd=5 \
	label1="Depth [m]" label2="Angl [deg]" wantaxis1=$(RMO-WANTAXIS1) \
	grid1=y o1num=$(RMO-X$*-O1NUM) d1num=10000 o2num=-30 d2num=30 \
	title=$(RMO-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(RMO-XSHIFT) yshift=$(RMO-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Cig-%-X0-over$(OVER).Vp: Ang-Mig-%-$(MYHOST).H 
	Window3d <Ang-Mig-$*-$(MYHOST).H  \
	min1=1800 max1=2200 n4=1 min4=0 |  \
	Grey pclip=100 color=i \
	screenwd=5 \
	label1="Depth [m]" label2="Angl [deg]" wantaxis1=$(CIG-$*-WANTAXIS1) \
	grid1=y o1num=2020 d1num=10000 o2num=-30 d2num=30 \
	title=$(CIG-$*-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(CIG-$*-XSHIFT) yshift=$(CIG-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Cig-%-X550-over$(OVER).Vp: Ang-Mig-%-$(MYHOST).H 
	Window3d <Ang-Mig-$*-$(MYHOST).H  \
	min1=1800 max1=2200 n4=1 min4=550 |  \
	Grey pclip=100 color=i \
	screenwd=5 \
	label1="Depth [m]" label2="Angl [deg]" wantaxis1=$(CIG-$*-WANTAXIS1) \
	grid1=y o1num=2000 d1num=10000 o2num=-30 d2num=30 \
	title=$(CIG-$*-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(CIG-$*-XSHIFT) yshift=$(CIG-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Mig-0-VLowFreq-stack-X-over$(OVER).Vp :  \
	Mig-0-VLowFreq-stack-over$(OVER).Vp \
	Mig-0-VLowFreq-X0-over$(OVER).Vp \
	Mig-0-VLowFreq-X550-over$(OVER).Vp \
	Mig-0-VLowFreq-X1100-over$(OVER).Vp 
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@

Mig-0-stack-X-over$(OVER).Vp :  \
	Mig-0-stack-over$(OVER).Vp \
	Mig-0-X0-over$(OVER).Vp \
	Mig-0-X550-over$(OVER).Vp \
	Mig-0-X1100-over$(OVER).Vp 
	Vppen big=n vpstyle=n erase=n \
	$^ \
	>$@ datapath=$(TEMPVP)  
#	Tube <$@

Mig-%-stack-over$(OVER).Vp: Mig-%-$(MYHOST).H 
	Window3d <Mig-$*-$(MYHOST).H n3=1 min3=0 \
	min1=1800 max1=2200 min2=-1600 max2=1600 |  \
	Grey pclip=100 color=i \
	label1="Depth [m]" label2="Horizontal location [m]" wantaxis1=y \
	o2num=0 d2num=550 grid2=y \
	title=$(STACK-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(STACK-XSHIFT) yshift=$(STACK-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@

Mig-0-X%-over$(OVER).Vp: Ang-Mig-0-$(MYHOST).H 
	Window3d <Ang-Mig-0-$(MYHOST).H  \
	min1=1800 max1=2200 n4=1 min4=$* |  \
	Grey pclip=100 color=i \
	screenwd=5 \
	label1="Depth [m]" label2="Angl [deg]" wantaxis1=n \
	grid1=y o1num=2000 d1num=10000 o2num=-30 d2num=30 \
	title=$(MIG-X$*-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(MIG-X$*-XSHIFT) yshift=$(MIG-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@


Mig-0-VLowFreq-X%-over$(OVER).Vp: Ang-Mig-0-VLowFreq-$(MYHOST).H 
	Window3d <Ang-Mig-0-VLowFreq-$(MYHOST).H  \
	min1=1800 max1=2200 n4=1 min4=$* |  \
	Grey pclip=100 color=i \
	screenwd=5 \
	label1="Depth [m]" label2="Angl [deg]" wantaxis1=n \
	grid1=y o1num=2000 d1num=10000 o2num=-30 d2num=30 \
	title=$(MIG-X$*-TITL) \
	datapath=$(TEMPVP) \
	>tmp.Vp 
	Vppen <tmp.Vp big=n vpstyle=n  \
	xshift=$(MIG-X$*-XSHIFT) yshift=$(MIG-YSHIFT) scale=.5 \
	datapath=$(TEMPVP)  \
	>$@
#	Tube <$@


${RESDIR}/%.v   :  %.Vp
	Vppen vpstyle=n big=n  <$*.Vp >junk.H out=$(RESDIR)/$*.v
