#define FIGDIR ../Fig
#define FIGSUP ../Fig2
#define FIGLIST marine fourcdp
#include <SEP.defs>
#define OSEBERG ../../Data/Oseberg

marine.action&:
	tube FIGDIR/marine.v3 ; exit(0)

FIGDIR/marine.v: menergycdp.H
	        Taplot allpos=y pclip=98. <menergycdp.H |Ta2vplot title="Marine Variability" out=FIGDIR/marine.v head=/dev/null label1="Offset [m]" label2="Midpoint Distance [m]"
	cake FIGDIR/marine.v3

FIGDIR/fourcdp.v: tpowcdp.H
	Window f3=101 n3=1 <tpowcdp.H | Taplot |Ta2vplot title="CMP Gather 5050" out=j1.v head=/dev/null label1="Time [m]" label2="Offset [m]"
	Window f3=111 n3=1 <tpowcdp.H | Taplot |Ta2vplot title="CMP Gather 5550" out=j2.v head=/dev/null label1="Time [m]" label2="Offset [m]"
	Window f3=121 n3=1 <tpowcdp.H | Taplot |Ta2vplot title="CMP Gather 6050" out=j3.v head=/dev/null label1="Time [m]" label2="Offset [m]"
	Window f3=131 n3=1 <tpowcdp.H | Taplot |Ta2vplot title="CMP Gather 6550" out=j4.v head=/dev/null label1="Time [m]" label2="Offset [m]"
	vp_SideBySideAniso j1.v j2.v >ja.v 
	vp_SideBySideAniso j3.v j4.v >jb.v 
	vp_OverUnderAniso  ja.v jb.v >FIGDIR/fourcdp.v
	vppen j1.v j2.v j3.v j4.v >FIGDIR/fourcdp.v3


FIGDIR/marine.v3: tpowcdpmov.H
	Taplot <tpowcdpmov.H |Ta2vplot title="CDP Gather" out=FIGDIR/marine.v3 >/dev/null


FIGDIR/marine1.v: menergy.H overlay.v
	Taplot allpos=y pclip=98. <menergy.H |Ta2vplot title="Marine Variability" out=marine1.v head=/dev/null label1="Offset [m]" label2="Source Distance [m]"
	vp_Overlay marine1.v overlay.v >FIGDIR/marine1.v
	cake FIGDIR/marine1.v3

FIGDIR/marine1.v3: tpowmov.H
	Taplot <tpowmov.H |Ta2vplot title="Raw Shot Gather" out=FIGDIR/marine1.v3 >/dev/null 

FIGDIR/marine2.v: menergy.H
	Thplot allpos=y <menergy.H sz=2 alpha=30 label3="Energy" title="Marine Variability" out=FIGDIR/marine2.v head=/dev/null


/* a rotated style plot !!! */
FIGDIR/marine2.ps: FIGDIR/marine2.v
	vppen big=n <FIGDIR/marine2.v |vppen align=lb >junk.v
		pstexpen junk.v FIGDIR/marine2.ps color=COLOR fatmult=FATMULT invras=INVRAS big=y

FIGDIR/marine3.v: menergy.H
	Graph <menergy.H title="Marine Energy Decay" label2="Energy" label1="Offset" out=FIGDIR/marine3.v head=/dev/null

FIGDIR/marine4.v: menergy.H
	Add log=y add=0.001 menergy.H >menergylog.H
	Graph <menergylog.H title="Marine Energy Decay" label2="Log(Energy)" label1="Offset" out=FIGDIR/marine4.v head=/dev/null


sgose.H:
	Window n      n1=4 <../../Data/Oseberg/Line_124.HDR >beforeshotno.H
	Window f1=5  n1=1 <../../Data/Oseberg/Line_124.HDR >shotno.H
	Window f1=6 n1=4 <../../Data/Oseberg/Line_124.HDR >beforeoffset.H
	Window f1=10 n1=1 <../../Data/Oseberg/Line_124.HDR >offset.H
	Window f1=12     <../../Data/Oseberg/Line_124.HDR >afteroffset.H
	Add cdpno

/* shot spacing is 50 m */

/* apply Tpow to raw data */

tpow.H: ../../Data/Oseberg/Line_124.H
	Tpow tpow=2 <OSEBERG/Line_124.H > tpow.H
	echo "d2=50 o2=0" >>tpow.H

menergy.H: tpow.H
	Energy <tpow.H lwind=1499 j1=1 >je.H o1=0
	Window n1=1 <je.H >menergy.H

midp.H: midp.ascii
	atoF <midp.ascii >midp.H
	echo esize=8 n1=2 n2=1 >>midp.H

overlay.v: midp.H
	Graph <midp.H plotfat=1 plotcol=6 yreverse=1 min2=185 max2=2585 wantaxis1=n wanttitle=n >/dev/null out=midp.v
	Graph <midp.H plotfat=2 plotcol=7 yreverse=1 min2=185 max2=2585 wantaxis1=n wanttitle=n >/dev/null out=midp2.v
	vp_Overlay midp2.v midp.v >overlay.v


/* show every 10th shot */
tpowmov.H: tpow.H
	Window j3=10 <tpow.H >tpowmov.H

/* show every 10th shot */
tpowcdpmov.H: tpowcdp.H
	Window j1=10 f3=97 j3=10 n3=20 <tpowcdp.H >tpowcdpmov.H

tpowcdp.H: cdpcube.H
	Tpow tpow=2 <cdpcube.H > tpowcdp.H

menergycdp.H: tpowcdp.H
	Energy <tpowcdp.H lwind=1499 j1=1 >je.H o1=0
	Window n1=1 <je.H >menergycdp.H

/* this sorts correctly, but is not a "Cube"  */
cdp.H: if not exist cdp.H
	Sort <OSEBERG/Line_124.H attrfile=OSEBERG/Line_124.HDR hdrkey1=cdp_no hdrkey2=sg_dist >cdp.H attrout=cdp.HDR 

/* with this we get an even odd effect, cannot be right  */
/*
cdpcube.H: cdp.H
	Cubify <cdp.H axis2=sg_dist axis3=cdp_no >cdpcube.H attrout=cdpcube.HDR d2=25 n2=48
*/

/* even , odd pattern is in near offset trace */
cdpcube.H: cdp.H
	Cubify <cdp.H axis2=sg_dist axis3=cdp_no >cdpcube.H attrout=cdpcube.HDR d2=50 n2=48
	echo "o3=0 d3=50" >> cdpcube.H
	echo "o3=0 d3=50" >> cdpcube.HDR



srcube.H:
	Cubify <OSEBERG/Line_124.H attrfile=OSEBERG/Line_124.HDR axis2=g_x_coord axis3=s_x_coord >srcube.H attrout=srcube.HDR d2=2500 n2=251 o2=49601781 

shotcube.H:
	Cubify <OSEBERG/Line_124.H attrfile=OSEBERG/Line_124.HDR axis2=sg_dist axis3=s_x_coord >shotcube.H attrout=hshotcube.H o2=185 n2=96 d2=25 

showkeys&:
	Dumpkeys <OSEBERG/Line_124.H attrfile=OSEBERG/Line_124.HDR 

clean&:
	RM_CMD *.H *.HDR *.v j*

#include <SEP.idoc.rules>
