#define BINDIR /r3/martin/bin/MTYPE
#include <SEP.defs>
#define ALIGN xll=0.25 yll=0.25 xur=5. yur=5.
#include <SEP.defs>
#define MING min2=0 max2=3.998 min1=-2.06 max1=7.2
#define MINT min1=0 max1=3.998 min2=-2.06 max2=7.2
#define VSLAN /r3/martin/bin/MTYPE/Vslan
/* note that py-quantities are not generated in 2D modeling */

/* subsample data since frequency content is <50 Hz ->   6 min per transform*/
$last&: sxseis.C.90.hudwet.cmp.SL szseis.C.90.hudwet.cmp.SL
last&: syseis.C.90.hudwet.cmp.SL sxseis.C.0.hudwet.cmp.SL szseis.C.0.hudwet.cmp.SL syseis.C.0.hudwet.cmp.SL sxseis.B.0.hudwet.cmp.SL szseis.B.0.hudwet.cmp.SL syseis.B.0.hudwet.cmp.SL sxseis.B.90.hudwet.cmp.SL szseis.B.90.hudwet.cmp.SL syseis.B.90.hudwet.cmp.SL sxseis.A.0.hudwet.cmp.SL szseis.A.0.hudwet.cmp.SL syseis.A.0.hudwet.cmp.SL sxseis.A.90.hudwet.cmp.SL szseis.A.90.hudwet.cmp.SL syseis.A.90.hudwet.cmp.SL

paper-sl&:   
	Taplot 	< pxseis.f10.0.cmp.SUB | Ta2vplot title="CMP p-src x-comp (inline)"  out=FIGDIR/pxseis.f10.0.cmp.SUB.v head=/dev/null
	Taplot 	< pzseis.f10.0.cmp.SUB | Ta2vplot title="CMP p-src z-comp (inline)"  out=FIGDIR/pzseis.f10.0.cmp.SUB.v head=/dev/null
	Taplot 	< pxseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src x-comp (inline)" label2="Horiz. Slowness [s/km]" out=FIGDIR/pxseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pzseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src z-comp (inline)" label2="Horiz. Slowness [s/km]" out=FIGDIR/pzseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pxseis.f10.SL.DIFF | Ta2vplot title="Slant Stack p-src z-comp (inline)"  label2="Horiz. Slowness [s/km]" out=FIGDIR/pxseis.f10.SL.DIFF.v head=/dev/null
	Taplot 	< pzseis.f10.SL.DIFF | Ta2vplot title="Slant Stack p-src z-comp (inline)"  label2="Horiz. Slowness [s/km]" out=FIGDIR/pzseis.f10.SL.DIFF.v head=/dev/null

all.vslan&: all.vslan.0 all.vslan.90

all.hudC&: pzseis.C.90.huddry.cmp.SUB.SL pxseis.C.90.huddry.cmp.SUB.SL                     pzseis.C.0.huddry.cmp.SUB.SL  pxseis.C.0.huddry.cmp.SUB.SL    

all.hudCC&: pzseis.C.90.hudwet.cmp.SUB.SL pxseis.C.90.hudwet.cmp.SUB.SL                     pzseis.C.0.hudwet.cmp.SUB.SL  pxseis.C.0.hudwet.cmp.SUB.SL    

all.hudB&: pzseis.B.0.hudwet.cmp.SUB.SL  pxseis.B.0.hudwet.cmp.SUB.SL
        /* pzseis.B.90.hudwet.cmp.SUB.SL pxseis.B.90.hudwet.cmp.SUB.SL       */ 

all.hudA&: pzseis.A.90.hudwet.cmp.SUB.SL pxseis.A.90.hudwet.cmp.SUB.SL                  /*   pzseis.A.0.hudwet.cmp.SUB.SL  pxseis.A.0.hudwet.cmp.SUB.SL    */


all.huddry&: pzseis.f10.0.huddry.cmp.SUB.SL pxseis.f10.0.huddry.cmp.SUB.SL                   pzseis.f10.90.huddry.cmp.SUB.SL pxseis.f10.90.huddry.cmp.SUB.SL

all.hudwet&: szseis.f10.0.hudwet.cmp.SUB.SL sxseis.f10.0.hudwet.cmp.SUB.SL                   szseis.f10.90.hudwet.cmp.SUB.SL sxseis.f10.90.hudwet.cmp.SUB.SL

$all.hudwet&: pzseis.f10.0.hudwet.cmp.SUB.SL pxseis.f10.0.hudwet.cmp.SUB.SL                   pzseis.f10.90.hudwet.cmp.SUB.SL pxseis.f10.90.hudwet.cmp.SUB.SL

all.homoti&: pzseis.homoti.cmp.SUB.SL pxseis.homoti.cmp.SUB.SL

all.vslan.0&: pxseis.f10.0.cmp.SUB.SL                                                         pzseis.f10.0.cmp.SUB.SL  /* sxseis.f10.0.cmp.SUB.SL                                syseis.f10.0.cmp.SUB.SL  szseis.f10.0.cmp.SUB.SL */

/* for SL for the calibration  */
all.iso&:  pxseis.iso.cmp.SUB.SL  pzseis.iso.cmp.SUB.SL                                    pxseis.homoiso.cmp.SUB.SL  pzseis.homoiso.cmp.SUB.SL

all.vslan.90&: pxseis.f10.90.cmp.SUB.SL                                                        pzseis.f10.90.cmp.SUB.SL  /* sxseis.f10.90.cmp.SUB.SL                              syseis.f10.90.cmp.SUB.SL  szseis.f10.90.cmp.SUB.SL */


pcurve&: ../Curves/ptigroup.h
	$echo "d2=0.02" >> zseis.unfrac.cmp
	TTticurve <zseis.unfrac.cmp  d2=0.02 tshift=0.08 z0=1.0 velocity=0 vel=../Curves/ptigroup.h >pcurve 


/* forward least squares slant stack */
%.SL:  %
        time VSLAN <% option=3 p0=-0.675 dp=.00375 np=360 noalias=n trim=0 >%.SL


/* Inverse Slant stack  (transpose)   */
%.In:
	time VSLAN <%.SL transp=1 option=2 x0=-4.12 dx=0.04 nx=464 noalias=y trim=0  >junk
        Window n1=2000 <junk >%.In

/* Subsample the seismograms  */
%.SUB: % 
	Window <% >%.SUB  j1=4
	
/* Bandpass           */
%.B: %
	Bandpass <% >%.B  flo=1. fhi=50.
	
/* Subtract Slant Stacks*/
%seis.f10.SL.DIFF: %seis.f10.0.cmp.SUB.SL %seis.f10.90.cmp.SUB.SL
        Add scale=1,-1 %seis.f10.0.cmp.SUB.SL %seis.f10.90.cmp.SUB.SL >%seis.f10.SL.DIFF

/* Subtract Seismograms   */
%seis.f10.DIFF: %seis.f10.0.cmp %seis.f10.90.cmp
        Add scale=1,-1 %seis.f10.0.cmp %seis.f10.90.cmp >%seis.f10.DIFF


/* plot everything  */

#define DIFFLIST pxseis.f10 pzseis.f10 sxseis.f10 syseis.f10 szseis.f10

plot&: seismo slants sldiffs diffs


/* how to plot subtracted quantitites */

sldiffs&: [[sub -i X X.SL.DIFF.v DIFFLIST]]

%.SL.DIFF.v: %.SL.DIFF
	Taplot 	< %.SL.DIFF | Ta2vplot title="% Slant Stack Difference"  out=FIGDIR/%.SL.DIFF.v head=/dev/null

diffs&: [[sub -i X X.DIFF.v DIFFLIST]]

%.DIFF.v:    %.DIFF
	Taplot 	< %.DIFF | Ta2vplot title="% Seismogram Difference"  out=FIGDIR/%.DIFF.v head=/dev/null
	


seismo&:
	Taplot 	< pxseis.f10.0.cmp | Ta2vplot title="Seismogram p-src x-comp (inline)"  out=FIGDIR/pxseis.f10.0.cmp.v head=/dev/null
	Taplot 	< pzseis.f10.0.cmp | Ta2vplot title="Seismogram p-src z-comp (inline)"  out=FIGDIR/pzseis.f10.0.cmp.v head=/dev/null
	Taplot 	< pxseis.f10.90.cmp  | Ta2vplot title="Seismogram p-src x-comp (xline)"  out=FIGDIR/pxseis.f10.90.cmp.v head=/dev/null
	Taplot 	< pzseis.f10.90.cmp | Ta2vplot title="Seismogram p-src z-comp (xline)"  out=FIGDIR/pzseis.f10.90.cmp.v head=/dev/null
	Taplot 	< sxseis.f10.0.cmp | Ta2vplot title="Seismogram s-src x-comp (inline)"  out=FIGDIR/sxseis.f10.0.cmp.v head=/dev/null
	Taplot 	< syseis.f10.0.cmp | Ta2vplot title="Seismogram s-src z-comp (inline)"  out=FIGDIR/syseis.f10.0.cmp.v head=/dev/null
	Taplot 	< szseis.f10.0.cmp | Ta2vplot title="Seismogram s-src z-comp (inline)"  out=FIGDIR/szseis.f10.0.cmp.v head=/dev/null
	Taplot 	< sxseis.f10.90.cmp | Ta2vplot title="Seismogram s-src x-comp (xline)"  out=FIGDIR/sxseis.f10.90.cmp.v head=/dev/null
	Taplot 	< syseis.f10.90.cmp | Ta2vplot title="Seismogram s-src z-comp (xline)"  out=FIGDIR/syseis.f10.90.cmp.v head=/dev/null
	Taplot 	< szseis.f10.90.cmp | Ta2vplot title="Seismogram s-src z-comp (xline)"  out=FIGDIR/szseis.f10.90.cmp.v head=/dev/null


slants&:
	Taplot 	< pxseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src x-comp (inline)"  out=FIGDIR/pxseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pzseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src z-comp (inline)"  out=FIGDIR/pzseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pxseis.f10.90.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src x-comp (xline)"  out=FIGDIR/pxseis.f10.90.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pzseis.f10.90.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src z-comp (xline)"  out=FIGDIR/pzseis.f10.90.cmp.SUB.SL.v head=/dev/null
	Taplot 	< sxseis.f10.0.cmp.SUB.SL | Ta2vplot title="Seismogram s-src x-comp (inline)"  out=FIGDIR/sxseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< syseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack s-src z-comp (inline)"  out=FIGDIR/syseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< szseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack s-src z-comp (inline)"  out=FIGDIR/szseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< sxseis.f10.90.cmp.SUB.SL | Ta2vplot title="Slant Stack s-src x-comp (xline)"  out=FIGDIR/sxseis.f10.90.cmp.SUB.SL.v head=/dev/null
	Taplot 	< syseis.f10.90.cmp.SUB.SL | Ta2vplot title="Slant Stack s-src z-comp (xline)"  out=FIGDIR/syseis.f10.90.cmp.SUB.SL.v head=/dev/null
	Taplot 	< szseis.f10.90.cmp.SUB.SL | Ta2vplot title="Slant Stack s-src z-comp (xline)"  out=FIGDIR/szseis.f10.90.cmp.SUB.SL.v head=/dev/null


poster&:   
	Taplot 	< pxseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src x-comp (inline)"  out=FIGDIR/pxseis.f10.0.cmp.SUB.SL.v head=/dev/null
	Taplot 	< pzseis.f10.0.cmp.SUB.SL | Ta2vplot title="Slant Stack p-src z-comp (inline)"  out=FIGDIR/pzseis.f10.0.cmp.SUB.SL.v head=/dev/null
        Taplot  < pxseis.f10.0.cmp.SUB.SL.DIFF | Ta2vplot title="% Seismogram Difference"  out=FIGDIR/%.DIFF.v head=/dev/null




extract&:
	TPisocurve velocity=3.6 z0=1.0 n1=512 n2=360 d1=0.008 o2=-0.675 d2=0.00375 > jcurvep
	TPisocurve velocity=2.06 z0=1.0 n1=512 n2=360 d1=0.008 o2=-0.675 d2=0.00375 > jcurves
	Graph < jcurvep min1=-0.675 max1=0.675 max2=4.096 yreverse=1 wantaxis=n wantframe=n wanttitle=n wantlabel=n >j2.V
	Graph < jcurves min1=-0.675 max1=0.675 max2=4.096 yreverse=1 wantaxis=n wantframe=n wanttitle=n wantlabel=n >j3.V
	Taplot <pzseis.iso.cmp.SUB.SL |Ta2vplot >j1.V
	Vppen <j1.V j2.V j3.V vpstyle=n erase=once |Tube




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