DATA1 = $(ELF_DATA)/elf.H ####### This is the data header you are trying to find #########
#
# H files
#

data.H: ${DATA1}
	cp $< junk.H
	cp $<@ /scr1/claudio/junk.H@
	cp $<@@ junk.H@@
	cp $<@@@ /scr1/claudio/junk.H@@@
	echo in="/scr1/claudio/junk.H@" >>junk.H
	echo hff="junk.H@@">>junk.H
	echo in="/scr1/claudio/junk.H@@@" >>junk.H@@
	Window3d n2=100000 < junk.H > $@
	mv /scr1/claudio/junk.H@ /scr1/claudio/$@@
	echo in="/scr1/claudio/$@@" >>$@
	rm -f junk.H*
	rm -f /scr1/claudio/junk*

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

sou.H: math1.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

rec.H: math1.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

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

soutst.H: data_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

rectst.H: data_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

sort_shot.H: math1.H
	Headermath <$< 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' > junk1.H
	Sort3d par=./Par/sortshot_ELF.P <junk1.H >junk2.H 
	<junk2.H Infill3d >junk3.H
	<junk3.H Stack3d normalize=1 > $@
	
sort_recv.H: math1.H
	Headermath <$< 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' > rjunk1.H
	Sort3d par=./Par/sortrecv_ELF.P <rjunk1.H >rjunk2.H 
	<rjunk2.H Infill3d >rjunk3.H
	<rjunk3.H Stack3d normalize=1 > $@

azimof.H: sort_shot.H
	Window3d n1=1 f1=42 <sort_shot.H@@ > aof.H
	Window3d n1=1 f1=49 <sort_shot.H@@ > azi.H
	Cmplx aof.H azi.H > $@

sousort.H: sort_shot.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

recsort.H: sort_shot.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
#
# Figs
#
sort_shot.V:
	Grey <sort_shot.H eout=1 gainpanel=all pclip=95 > sort_shot.T
	Cubeplot movie=0 frame1=37 frame2=5 frame3=5 point1=.75 point2=.75 wantcoordlabel=1 wantaxes=y popup=n \
	screenratio=1.5 < sort_shot.T > /dev/null out=$R/$@

sou.v: sou.H
	Graph <sou.H symbol="+" label2="Y(m)" label1="X(m)" title=" " \
	min2=3000. max2=4000. min1=-1000. max1=13000. > /dev/null out=$R/$@
	Rm sou.H

rec.v: rec.H
	Graph <rec.H symbol="." label2="Y(m)" label1="X(m)" title=" " \
	min2=3000. max2=4000. min1=-1000. max1=13000. > /dev/null out=$R/$@
	Rm rec.H

sousort.v: sousort.H
	Graph <sousort.H symbol="+" label2="Y(m)" label1="X(m)" title=" " \
	min2=3000. max2=4000. min1=-1000. max1=13000. > /dev/null out=$R/$@
	Rm sousort.H

recsort.v: recsort.H
	Graph <recsort.H symbol="." label2="Y(m)" label1="X(m)" title=" " \
	min2=3000. max2=4000. min1=-1000. max1=13000. > /dev/null out=$R/$@
	Rm recsort.H

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

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

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