DATA1 = ${CGG_DATA}/cgg.H 
#
# H files
#

cgg001data.H: ${DATA1}
	cp $< junk1.H
	echo in="/mnt/sepdata/1b/cgg_green_canyon/cgg.H@;/mnt/sepdata/1b/cgg_green_canyon/cgg.H@1;/mnt/sepdata/1b/cgg_green_canyon/cgg.H@2;/mnt/sepdata/1b/cgg_green_canyon/cgg.H@3;/mnt/sepdata/1b/cgg_green_canyon/cgg.H@4" >> junk1.H
	cp $<@@ junk1.H@@
	echo hff="junk1.H@@" >> junk1.H
	cp $<@@@ /scr1/claudio/junk1.H@@@
	echo in="/scr1/claudio/junk1.H@@@" >> junk1.H@@
	Window3d n2=1497000 < junk1.H > $@
#	Window3d n2=1000 < junk1.H > $@
	rm -f junk1.H*
	rm -f /scr1/claudio/junk1.H*

cgg002math1.H: cgg001data.H 
	<$< Headermath par=math1.P > $@

cggsou.H: cgg.H
	Window3d < $<@@ n1=1 f1=0 > sx.H
	Window3d < $<@@ n1=1 f1=1 > sy.H
	Cmplx sx.H sy.H maxsize=1000 > $@
	Rm sx.H sy.H

cggrec.H: cgg.H
	Window3d < $<@@ n1=1 f1=2 > rx.H
	Window3d < $<@@ n1=1 f1=3 > ry.H
	Cmplx rx.H ry.H  maxsize=1000 > $@
	Rm rx.H ry.H

cgg001data_test.H: math1.H
	Window3d <$< n2=100 > $@

cggsoutst.H: cgg001data_test.H
	Window3d < $<@@ n1=1 f1=37 > sx.H
	Window3d < $<@@ n1=1 f1=38 > sy.H
	Cmplx sx.H sy.H > $@
	Rm sx.H sy.H

cggrectst.H: cgg001data_test.H
	Window3d < $<@@ n1=1 f1=39 > rx.H
	Window3d < $<@@ n1=1 f1=40 > ry.H
	Cmplx rx.H ry.H > $@
	Rm rx.H ry.H

cgg003sort_shot.H: cgg002math1.H
	<$< Headermath par=math1.P > junk1.H
	Headermath <junk1.H key1=aoffset 'eqn1=@SQRT((s_x-g_x)^2+(s_y-g_y)^2)' key2=cmp_x 'eqn2=(s_x+g_x)/2)' \
	key3=cmp_y 'eqn3=(s_y+g_y)/2)' key4=srazim 'eqn4=@ATAN((g_y-s_y)/(g_x-s_x))*180/3.14159' \
	key5=offset_x 'eqn5=g_x-s_x' key6=offset_y 'eqn6=g_y-s_y' > junk2.H
	Sort3d par=./Par/sortshot_CGG.P max_memory=10000 <junk2.H > junk3.H
	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 > $@

cgg003sort_cmp.H: cgg002math1.H
	<$< Headermath par=math1.P > junk1.H
	Headermath <junk1.H key1=aoffset 'eqn1=@SQRT((s_x-g_x)^2+(s_y-g_y)^2)' key2=cmp_x 'eqn2=(s_x+g_x)/2)' \
	key3=cmp_y 'eqn3=(s_y+g_y)/2)' key4=srazim 'eqn4=@ATAN((g_y-s_y)/(g_x-s_x))*180/3.14159' \
	key5=offset_x 'eqn5=g_x-s_x' key6=offset_y 'eqn6=g_y-s_y' > junk2.H
	Sort3d par=./Par/sortcmp_CGG.P max_memory=10000 <junk2.H > junk4.H
#	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 > $@

cgg003sort_cmp2.H: cgg002math1.H
	<$< Headermath par=math1.P > junk1.H
	Headermath <junk1.H key1=aoffset 'eqn1=@SQRT((s_x-g_x)^2+(s_y-g_y)^2)' key2=cmp_x 'eqn2=(s_x+g_x)/2)' \
	key3=cmp_y 'eqn3=(s_y+g_y)/2)' key4=srazim 'eqn4=@ATAN((g_y-s_y)/(g_x-s_x))*180/3.14159' \
	key5=offset_x 'eqn5=g_x-s_x' key6=offset_y 'eqn6=g_y-s_y' > junk2.H
	Sort3d par=./Par/sortcmp_CGG1.P max_memory=10000 <junk2.H > junk4.H
#	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 > $@


cgg004sort_recv.H:
	Sort3d par=./Par/sortrecv_CGG.P max_memory=10000 <junk2.H > junk3.H
	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 maxsize=10000 > $@

cgg003sort_shota.H: 
	Sort3d par=./Par/sortshot_CGG.P max_memory=10000 <junk2.H > junk3a.H
	<junk3a.H Infill3d > $@

cgg004sort_recva.H:
	Sort3d par=./Par/sortrecv_CGG.P max_memory=10000 <junk2.H > junk3.H
	<junk3.H Infill3d > $@

cgg005sort_suprecv.H:
	Sort3d par=./Par/sortsuprecv_CGG.P max_memory=10000 <junk2.H > junk3.H
	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 maxsize=10000 > $@

cgg006sort_suprecv.H:
	Sort3d par=./Par/sortsuprecv_CGG.P max_memory=10000 <junk2.H > junk3.H
	<junk3.H Infill3d >junk4.H
	<junk4.H Stack3d normalize=1 maxsize=10000 > $@

cggazimof.H: cgg003sort_shot.H
	Headermath <$< key1=srazim 'eqn1=@ATAN((g_y-s_y)/(g_x-s_x))*180/3.14159' > junk1.H 
	Window3d n1=1 f1=9 <junk1.H@@ > aof.H
	Window3d n1=1 f1=10 <junk1.H@@ > azi.H
	Cmplx aof.H azi.H > $@
	rm aof.H azi.H

cggsousort.H: cgg003sort_shot.H
	Window3d < $<@@ n1=1 f1=0 > sx.H
	Window3d < $<@@ n1=1 f1=1 > sy.H
	Cmplx sx.H sy.H > $@
	Rm sx.H sy.H

cggrecsort.H: cgg003sort_shot.H
	Window3d < $<@@ n1=1 f1=2 > rx.H
	Window3d < $<@@ n1=1 f1=3 > ry.H
	Cmplx rx.H ry.H > $@
	Rm rx.H ry.H
#
# Figs
#
cgg003sort_shot.v:
	< cgg003sort_shot.H Grey gainpanel=all eout=1 pclip=90  label1="time(s)"> cgg003sort_shot.T
	< cgg003sort_shot.T Cubeplot movie=0 frame1=1361 frame2=5 frame3=8 point1=.9 point2=.9 title=" "\
	wantcoordlabel=1 wantaxes=y popup=n screenratio=1.5 > /dev/null out=$R/$@

cgg004sort_recv.v:
	< cgg004sort_recv.H Grey gainpanel=all eout=1 pclip=90 label1="time(s)"> cgg004sort_recv.T
	< cgg004sort_recv.T Cubeplot movie=0 frame1=1361 frame2=1 frame3=7 point1=.9 point2=.9 title=" "\
	wantcoordlabel=1 wantaxes=y popup=n screenratio=1.5 > /dev/null out=$R/$@

cgg005sort_suprecv.v:
	< cgg005sort_suprecv.H Grey gainpanel=all eout=1 pclip=90 label1="time(s)"> cgg005sort_suprecv.T
	< cgg005sort_suprecv.T Cubeplot movie=0 frame1=1361 frame2=5 frame3=4 point1=.9 point2=.9 title=" "\
	wantcoordlabel=1 wantaxes=y popup=n screenratio=1.5 > /dev/null out=temp.v
	<temp.v vp_annotate batch=y text=label.txt | vppen vpstyle=n > $R/$@

cgg006sort_suprecv.v:
	< cgg006sort_suprecv.H Grey gainpanel=all eout=1 pclip=90 label1="time(s)"> cgg006sort_suprecv.T
	< cgg006sort_suprecv.T Cubeplot movie=0 frame1=1361 frame2=10 frame3=14 point1=.9 point2=.9 title=" "\
	wantcoordlabel=1 wantaxes=y popup=n screenratio=2. > /dev/null out=$R/$@

cgg005sup.v: cgg005sort_suprecv.H
	< cgg005sort_suprecv.H Grey gainpanel=all eout=1 pclip=90 label1="time(s)"> cgg005sort_suprecv.T
	< cgg005sort_suprecv.T Cubeplot movie=0 frame1=1361 frame2=5 frame3=4 point1=.9 point2=.9 title=" "\
	wantcoordlabel=1 wantaxes=y popup=n screenratio=2. > /dev/null out=temp.v
	vppen < temp.v xcenter=-0.5 ycenter=-0.65 vpstyle=n > $R/$@

cgg006sup.v: cgg006sort_suprecv.v
	< $R/cgg006sort_suprecv.v vppen xcenter=-6.0 ycenter=-0.65 vpstyle=n > ${R}/$@

suprec.v: cgg005sup.v cgg006sup.v
	vp_Overlay $R/cgg006sup.v $R/cgg005sup.v > $R/$@

cggsou.v: cggsou.H
	Graph <$< symbol="+" label2="Y(m)" label1="X(m)" title=" " \
	 > /dev/null out=$R/$@
	Rm cggsou.H
#	min2=7700. max2=9400. min1=21500. max1=36000.

cggrec.v: cggrec.H
	Graph <$< symbol="." label2="Y(m)" label1="X(m)" title=" " \
	> /dev/null out=$R/$@
	Rm cggrec.H
#	min2=7700. max2=9400. min1=21500. max1=36000. 

cggazimof.v: cggazimof.H
	Graph <$< symbol="." d2num=20. d1num=1500. label2="azimuth(deg)" label1="aoffset(m)" title=" " \
	min2=-75. max2=75. min1=0. max1=9000. > /dev/null out=$R/$@
	Rm cggazimof.H

cggsousort.v: cggsousort.H
	Graph <$< symbol="+" label2="Y(m)" label1="X(m)" title=" " \
	min2=7900. max2=9600. min1=28000. max1=52000. > /dev/null out=$R/$@
	Rm cggsousort.H

cggrecsort.v: cggrecsort.H
	Graph <$< symbol="." label2="Y(m)" label1="X(m)" title=" " \
	min2=7900. max2=9600. min1=28000. max1=52000. > /dev/null out=$R/$@
	Rm cggrecsort.H

cggsoutst.v: cggsoutst.H
	Graph <$< symbol="+" label2="Y(m)" label1="X(m)" title=" " \
	> /dev/null out=$R/$@
	Rm cggsoutst.H

cggrectst.v: cggrectst.H
	Graph <$< symbol="." label2="Y(m)" label1="X(m)" title=" " \
	> /dev/null out=$R/$@
	Rm cggrectst.H