
#define MARMDATA ../../Marmdata

#define CMBINDIR ../../Bin/CM

#define BINDIR ../../Bin/MTYPE

#include <SEP.defs>

#ifdef ISsun4


%.TT1: CMBINDIR/TTtowave3 %.TT ttpar_%
	time CMBINDIR/TTtowave3 <%.TT par=ttpar_% nt=160 > %.TT1

%.GR2: CMBINDIR/Omgreens greenpar2_% %.vel
	time CMBINDIR/Omgreens <%.vel par=greenpar2_% > %.GR2 outform=1 

CMBINDIR/% : if not exist CMBINDIR/% and % in TTtowave3 Omgreens
	pushd ../../TTparam/Src ; cake %
#else


%.TT1: CMBINDIR/TTtowave2 %.TT ttpar_%
	echo "You need a CM5 to create %.TT1"

%.GR2: CMBINDIR/Omgreens greenpar2_% %.vel
	echo "You need a CM5 to create %.GR2"
#endif


$======================================
$ finite difference eikonal traveltimes
$======================================

marm%.ETT :  BINDIR/TTmaps ettpar_marm% marm%.slow
	BINDIR/TTmaps <marm%.slow par=ettpar_marm% >marm%.ETT

BINDIR/TTmaps : if not exist BINDIR/TTmaps
	cd ../Src; cake TTmaps

$======================================
$ band-limited traveltimes
$======================================

marm%.TT :  BINDIR/MTTmaps ttpar_marm% marm%.slow
	BINDIR/MTTmaps <marm%.slow par=ttpar_marm% >marm%.TT

BINDIR/MTTmaps : if not exist BINDIR/MTTmaps
	cd ../Src; cake MTTmaps


$========================================
$ subsampled velocity models
$========================================

#define INV_ARGS log=yes scale=-1. exp=yes

truemarm.slow : MARMDATA/marmvel.H
	Add MARMDATA/marmvel.H INV_ARGS >truemarm.slow out=stdout

marm%.slow : truemarm.slow smoothpar_marm% velpar_marm%
	Smooth < truemarm.slow par=smoothpar_marm% | Window par=velpar_marm% >marm%.slow

marm%.vel: marm%.slow 
	Add marm%.slow  INV_ARGS  >marm%.vel

$=================================

clean&:
	-Rm marm*.slow marm*.vel
	-Rm *.ETT *.TT *.TT1 *.GR2
