#define DATAPATH /SDA/martin/AVO
###define DATAPATH /scr5/martin/AVO
#define DEBUG
#define BINDIR /r3/martin/bin/MTYPE
#define SEQ -S0-1
#include <SEP.defs>
#ifdef LARGE
#define SIZE  n1=488 n2=976 d1=0.020 d2=0.040
#else
#define SIZE  n1=232 n2=464 d1=0.020 d2=0.040
#endif

default&: /* model.calibrate  CMP/pzseis.iso.cmp CMP/pxseis.iso.cm CMP/pzseis.homoiso.cmp CMP/pxseis.homoiso.cmp */ CMP/pzseis.iso.cmp.SUB.SL CMP/pxseis.iso.cmp.SUB.SL CMP/pzseis.homoiso.cmp.SUB.SL CMP/pxseis.homoiso.cmp.SUB.SL


model.%&:
	cake -f Cake.model %


test&: frac10.90.layer wavelet
	time cmattach SEQ BINDIR/Ultimod  type=7 src_type=1 ncomp=3 modeling=1 moduli=frac10.90.layer par=parfile  <wavelet >junk




clean&:
	-RM_CMD *.x *.o  a.out core j*


/* cmp are in half offset ( m,h) */

last&: CMP/sxseis.C.90.hudwet.cmp CMP/szseis.C.90.hudwet.cmp                           CMP/syseis.C.90.hudwet.cmp CMP/sxseis.C.0.hudwet.cmp                            CMP/szseis.C.0.hudwet.cmp CMP/syseis.C.0.hudwet.cmp                             CMP/sxseis.B.0.hudwet.cmp CMP/szseis.B.0.hudwet.cmp                             CMP/syseis.B.0.hudwet.cmp CMP/sxseis.B.90.hudwet.cmp                            CMP/szseis.B.90.hudwet.cmp CMP/syseis.B.90.hudwet.cmp                           CMP/sxseis.A.0.hudwet.cmp CMP/szseis.A.0.hudwet.cmp                             CMP/syseis.A.0.hudwet.cmp CMP/sxseis.A.90.hudwet.cmp                            CMP/szseis.A.90.hudwet.cmp CMP/syseis.A.90.hudwet.cmp     


all.cmp&: all.cmp.10.y all.cmp.10.x all.cmp.10.z

all.cmp.10.y&: CMP/pyseis.f10.0.cmp CMP/pyseis.f10.90.cmp                                      CMP/syseis.f10.0.cmp CMP/syseis.f10.90.cmp

all.cmp.10.x&: CMP/pxseis.f10.0.cmp CMP/pxseis.f10.90.cmp		                       CMP/sxseis.f10.0.cmp CMP/sxseis.f10.90.cmp

all.cmp.10.z&: CMP/pzseis.f10.0.cmp CMP/pzseis.f10.90.cmp                                      CMP/szseis.f10.0.cmp CMP/szseis.f10.90.cmp

/* this is for calibration a isotropic model */
all.cmp.iso&: CMP/pzseis.iso.cmp CMP/pxseis.iso.cmp

/* this is for calibration a isotropic model */
all.cmp.homoiso&: CMP/pzseis.homoiso.cmp CMP/pxseis.homoiso.cmp
          /*    CMP/szseis.f10.0.cmp CMP/szseis.f10.90.cmp */

/* this is the reference top layer  */
all.cmp.homoti&: CMP/pzseis.homoti.cmp CMP/pxseis.homoti.cmp                                     /* CMP/szseis.homoti.cmp CMP/sxseis.homoti.cmp */

/* Hudson Cracks  */
all.cmp.C&:  CMP/pxseis.C.90.huddry.cmp  CMP/pzseis.C.90.huddry.cmp                          CMP/pxseis.C.0.huddry.cmp   CMP/pzseis.C.0.huddry.cmp 

$all.cmp.C&:  CMP/pxseis.C.90.hudwet.cmp  CMP/pzseis.C.90.hudwet.cmp                          CMP/pxseis.C.0.hudwet.cmp   CMP/pzseis.C.0.hudwet.cmp 

$all.cmp.A&:  CMP/pxseis.A.90.hudwet.cmp  CMP/pzseis.A.90.hudwet.cmp                      /*    CMP/pxseis.A.0.hudwet.cmp  CMP/pzseis.A.0.hudwet.cmp */
all.cmp.B&:  CMP/pxseis.B.90.hudwet.cmp  CMP/pzseis.B.90.hudwet.cmp                          CMP/pxseis.B.0.hudwet.cmp   CMP/pzseis.B.0.hudwet.cmp 

$all.cmp.A&:  CMP/pxseis.A.90.hudwet.cmp  CMP/pzseis.A.90.hudwet.cmp                      /*    CMP/pxseis.A.0.hudwet.cmp  CMP/pzseis.A.0.hudwet.cmp */

all.cmp.A&:   CMP/pxseis.A.0.hudwet.cmp  CMP/pzseis.A.0.hudwet.cmp

$all.cmp.hud&:  CMP/pzseis.f10.0.huddry.cmp  CMP/pxseis.f10.0.huddry.cmp                        CMP/pzseis.f10.90.huddry.cmp  CMP/pxseis.f10.90.huddry.cmp
all.cmp.hud&:  CMP/syseis.f10.0.hudwet.cmp  CMP/syseis.f10.90.hudwet.cmp                     /*   CMP/szseis.f10.0.hudwet.cmp  CMP/sxseis.f10.0.hudwet.cmp                        CMP/szseis.f10.90.hudwet.cmp  CMP/sxseis.f10.90.hudwet.cmp */

$all.cmp.hud&:  CMP/pzseis.f10.0.huddry.cmp  CMP/pxseis.f10.0.huddry.cmp                        CMP/pzseis.f10.90.huddry.cmp  CMP/pxseis.f10.90.huddry.cmp

CMP/%1.%2.cmp.%3.%4:
	cd CMP ; cake %1.%2.cmp.%3.%4

CMP/%1.cmp.%2:
	cd CMP ; cake %1.cmp.%2

/* CMP/pzseis.f10.90.hudwet.cmp  CMP/pxseis.f10.90.hudwet.cmp */
/* CMP/pzseis.f10.0.hudwet.cmp  CMP/pxseis.f10.0.hudwet.cmp */

CMP/%.cmp:
	Transp plane=12 <% >jtemp 
	Reverse which=2 <jtemp >CMP/%.cmp out=DATAPATH/%.cmp@
	echo "o2=-4.12 d2=0.04 " >> CMP/%.cmp


all.snap&: all.snap.10.y all.snap.10.x all.snap.10.z

all.snap.10.y&: SNAP/pysnap.f10.0 SNAP/pysnap.f10.90                                            SNAP/sysnap.f10.0 SNAP/sysnap.f10.90

all.snap.10.x&: SNAP/pxsnap.f10.0 SNAP/pxsnap.f10.90		                                SNAP/sxsnap.f10.0 SNAP/sxsnap.f10.90

all.snap.10.z&: SNAP/pzsnap.f10.0 SNAP/pzsnap.f10.90                                            SNAP/szsnap.f10.0 SNAP/szsnap.f10.90



SNAP/%:
	Cp <% >SNAP/% out=DATAPATH/%@
	echo "o1=0. d2=0.02 o2=-4.12 d2=0.04 o3=0 d3=0.2" >>SNAP/% 
	echo 'label1="depth [km]" label2="distance [km]"' >>SNAP/%


showC&: ../Curves/ptigoup.h
	TTticurve <zseis.unfrac.t velocity=0 vel=ptigroup.h >pcurve 


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