$#define DEBUG
#define BINDIR /r3/martin/bin/MTYPE
#define FIGDIR ../Fig
#define REALLOG ../RealLog
#include <SEP.defs>
#define GRAPHMAX min2=1400 max2=5200
#define GRAPHMAXSMO min2=1400 max2=3400

default&:

reflectz.H:
	Add REALLOG/velp.h REALLOG/rho.h mode=product >jimp.H
	echo "+1. -1." >jfilt
	atoF <jfilt >jfilt.H
	Conv filtin=jfilt.H <jimp.H >jnum.H edge=1
	echo "+1. +1." >jfilt
	atoF <jfilt >jfilt.H
	Conv filtin=jfilt.H <jimp.H >jdenom.H edge=1
	Add  jnum.H jdenom.H mode=divide >reflectz.H

reflectt.H: reflectz.H
	Map2time.x velmap=REALLOG/velzp.pi.h <reflectz.H >jreft.H
	Real <jreft.H >j1
	Imag <jreft.H >j2
	Merge j1 j2 axis=2 space=n >jref.H
	Interpwell nzo=900 dzo=0.004 zo1=0 <jref.H >j.H
	Window n2=1 f2=1 <j.H >reflectt.H

band.H: reflectt.H
	Bandpass fhigh=50. <reflectt.H >band.H
	Graph <band.H |Tube

syntrace.H: reflectt.H
	Spectrum <../Stack/nmo.H >jspec.H
	<../Velan/cdp641.h Wavelet wavelet=data >jwavelet.H n1=1024 phase=none
	Window <jwavelet.H n1=100 >wavelet.H
	Conv filtin=wavelet.H <reflectt.H | Bandpass flo=1. fhi=50 >syntrace.H

comp&: syntrace.H
	 <../Stack/stack641.H Scale >j1 dscale=1
	Scale <syntrace.H >j2 dscale=1
	Merge j1 j2 axis=3 space=n >jcomp

autocor.H:
	Pad n1=0 <reflectz.H |Rtoc | Ft3d  sign1=1 sign2=0 sign3=0 >jft.H
	Cabs <jft.H |Rtoc >jftcor.H
	Ft3d sign1=-1 sign2=0 sign3=0 <jftcor.H >autocor.H


pef:
	Wdec < ../Stack/stack641.H nlag=1 npef=10  >pef.H
$ Iner < ../Stack/stack641.H padin=padin.H lambda=0. >pef.H


veryclean&: jclean
	RM_CMD *.h  *log j* 

#include <SEP.obj.rules>
#include <SEP.prog.rules>
#include <SEP.idoc.rules>
