include ${SEPINC}/SEP.top
RESDIR = ./Fig
R = ${RESDIR}
PARDIR = ./par
P = ${PARDIR}

RESULTSER = 
RESULTSCR = 
RESULTSNR = 1

n = >/dev/null
sizes  =  symbolsz=7 screenratio=1. 
noframe= wantframe=n wantaxis=n
graph  = Graph title=" " label2="Y (m)" label1="X (m)"           
grey   = Grey  title=" " label1="Y (m)" label2="X (m)" pclip=100

ROTATE= -43.655792 
include EXXONmakefile
include MIGmakefile
include BPmakefile

EXXON311=/data/exxon-311/


Chapter5-repro-test:
	make clean
	rm -f ./Bin/*.[x] ./Src/*.o ./Src/*.mod 
	@-cd ./Src;make default; cd ..
			    # NR 5.1
			    # NR 5.2
	make -n CrossIn.v   # CR 5.3
	make -n Depth1300.v # CR 5.4
	make VELCUT.v       # ER 5.5
			    # NR 5.6
			    # NR 5.7
			    # NR 5.8
			    # NR 5.9
			    # NR 5.10
	make XOM-VEL.v      # ER 5.11
	make -n XOM-SRC.v   # CR 5.12
	make -n XOM-RCV.v   # CR 5.13
	make -n RFIG1.v     # CR 5.14
	make -n RFIG3.v     # CR 5.15
	make -n RFIG5.v     # CR 5.16




######################################
#
# . . Make floating point headers
# . . All traces
EXXON311a.H:
	< $(EXXON311)/EXXON311.HH Headermath \
		key1=s_x_old eqn1="1.0*sx" type1=scalar_float \
		key2=s_y_old eqn2="1.0*sy" type2=scalar_float \
		key3=g_x_old eqn3="1.0*gx" type3=scalar_float \
		key4=g_y_old eqn4="1.0*gy" type4=scalar_float \
		key5=cmp_x_old eqn5="0.5*(sx+gx)" type5=scalar_float \
		key6=cmp_y_old eqn6="0.5*(sy+gy)" type6=scalar_float \
		key7=off_x_old eqn7="(sx-gx)" type7=scalar_float \
		key8=off_y_old eqn8="(gy-sy)" type8=scalar_float \
		key9=absoff_old eqn9="@SQRT((sx-gx)*(sx-gx)+(gy-sy)*(gy-sy))" type9=scalar_float\
		key10=azimuth_old  eqn10="@ATAN((gy-sy)/(gx-sx))" type10=scalar_float\
		key11=s_x eqn11="(sx-497634.3106)*0.7235-0.6903*(sy-3183172.481)" type11=scalar_float\
		key12=s_y eqn12="-(sx-497634.3106)*0.6903-0.7235*(sy-3183172.481)" type12=scalar_float\
		key13=g_x eqn13="(gx-497634.3106)*0.7235-0.6903*(gy-3183172.481)" type13=scalar_float\
		key14=g_y eqn14="-(gx-497634.3106)*0.6903-0.7235*(gy-3183172.481)" type14=scalar_float\
		key15=off_x eqn15="(gx-497634.3106)*0.7235-0.6903*(gy-3183172.481)\
			-(sx-497634.3106)*0.7235+0.6903*(sy-3183172.481)" label15=scalar_float\
		key16=off_y eqn16="-(gx-497634.3106)*0.6903-0.7235*(gy-3183172.481)\
			+(sx-497634.3106)*0.6903+0.7235*(sy-3183172.481)" label16=scalar_float > $@

######################################
#
# . . Rotations
#
# . . Look at all grids
# . . Source position
G0.H:     EXXON311a.H
	< EXXON311a.H Grid_fold par=sort.p > $@

# . . Source Positions
XOM-SRC.v: G0.H
	<  G0.H Transp plane=13 | Stack axes=3 maxsize=1000 normalize=y  > G1.H
	echo "label3='Source Fold'" >> G1.H
	Math file1=$< exp='file1*10'| Grey newclip=1 bclip=0 eclip=4 out=$@ $(n) wantscalebar=y\
		label1="Inline Source (m)" label2="Crossline Source (m)" title=" " transp=n yreverse=n
	rm -f G1.H

# . . Receiver Position
XOM-RCV.v: G0.H
	<  G0.H Transp plane=24 | Stack axes=3 maxsize=1000 normalize=y > G2.H	
	echo "label3='Receiver Fold'" >> G2.H
	Math file1=$< exp='file1*10000' > azaz.H
	< azaz.H Grey newclip=1 bclip=0 eclip=2000 out=$@ $(n) wantscalebar=y\
		label1="Inline Receiver (m)" label2="Crossline Receiver (m)" title=" " transp=n yreverse=n
	rm -f G2.H azaz.H

# . . ER Figure
XOM-VEL.v: 
	< $(EXXON311)/Vel/exxon-velocity.H Window3d j1=2 n1=225 > vcut.H
	echo "d1=30" >>  vcut.H
	Byte < vcut.H pclip=95 allpos=y gainpanel=a|\
		Window3d min2=$(min2) max2=$(max2) min3=$(min3) max3=$(max3) max2=14000|\
		Cubeplot frame1=83 frame2=183 frame3=188 wantcoordlabel=-1 out=$@ $(dn) \
		label1="Depth (m)" label2="Inline (m)" label3="Crossline (m)" title=" "\
		labelsz=8 point1=0.6 point2=0.66 labelsz=6
	rm -f vcut.H
#######################################
#
# . . Velocity section

XOFF=497634.3106
YOFF=3183172.481

exxon-velocity.H:
	Math    file1=$(EXXON311)/Vel//p_vertical.H \
		file2=${EXXON311}/Vel/delta.anistrpy.H \
		exp="file1*@SQRT(1+2*file2)" > $@
	echo "o1=0 o2=0 o3=0" >> $@

NoSaltVelocity.H:
	Math    file1=$(EXXON311)/Vel/NoSaltVelocity \
		file2=$(EXXON311)/Vel/delta.anistrpy.H \
		exp="file1*@SQRT(1+2*file2)" > $@
	echo "o1=0 o2=0 o3=0 " >> $@


########################################
#
# . . Make the entire infilled cube

exxon311-%.H: EXXON311-%.H
	< EXXON311-$*.H Infill3d synch=1 maxsize=100 axes=2 verb=y > $@

exxon-dec-%.H: exxon311-0%0.H exxon311-0%1.H exxon311-0%2.H exxon311-0%3.H exxon311-0%4.H\
               exxon311-0%5.H exxon311-0%6.H exxon311-0%7.H exxon311-0%8.H exxon311-0%9.H
	Cat axis=5 space=n exxon311-0$*[0123456789].H > $@ 
	cp /scr2/jeff/exxon-dec-$*.H@ /data/exxon-311/exxon-dec-$*.H@
	rm -f /scr2/jeff/exxon-dec-$*.H@
	rm -f EXXON311-0$*[0123456789].H
	rm -f EXXON311-0$*[0123456789].H@@
	rm -f EXXON311-0$*[0123456789].H@@@@
	rm -f exxon311-0$*[0123456789].H
	rm -f /scr2/jeff/exxon311-0$*[0123456789].H@
	rm -f /scr2/jeff/exxon311-0$*[0123456789].H@@@
	rm -f /scr2/jeff/exxon311-0$*[0123456789].H@@@@@
	echo 'sets next: in="/data/exxon-311/exxon-dec-$*.H@"' >> $@

# . . Infill data
DEC: exxon-dec-0.H exxon-dec-1.H exxon-dec-2.H exxon-dec-3.H exxon-dec-4.H \
	   exxon-dec-5.H exxon-dec-6.H exxon-dec-7.H exxon-dec-8.H exxon-dec-9.H
	echo @shit

clean: jclean klean.3d

include ${SEPINC}/SEP.bottom