Synopsis
		sfpwd < in.rsf dip=dip.rsf > out.rsf both=n n4=2 order=1 nj1=1 nj2=1 drift=n

February 2013 program of the month:
http://ahay.org/blog/2013/02/09/program-of-the-month-sfpwd/

Parameters
		bool both=n [y/n] if y, compute both left and right predictions   file dip= auxiliary input file name   bool drift=n [y/n] if shift filter   int n4=2 what to compute in 3-D. 0: in-line, 1: cross-line, 2: both   int nj1=1 in-line aliasing   int nj2=1 cross-line aliasing   int order=1 accuracy

Used In
		CWP          geo2006TimeShiftImagingCondition/zicig geo2008InterferometricImagingCondition/circle   DATA          alaska/line31-81-m8r cant12/stk   GEE          lsq/puck   GEO384W          hw4/sigsbee   MILANO          taupvel/synth   SEP          pwd/blast pwd/lomo vc2/beivc   TCCS          anisodiff/barrolka anisodiff/pcable anisodiff/simple3d diffr/chan3d diffr/fault diffr/gom dlchan/testing enhemd/sigmoid flat/flat fpwd/sigmoid fpwd/teapot ovc/nankai pi/bei pmig/mig pwshape/qdome seislet/lena shape/lomo sobel/pari sobel/penobscot sparse/sparse-experim-shaping sparse/vg uncert/beivc   ZJU          diffrlrra/fault diffrlrra/nankai_new diffrlrra/synth