#define USRLIBS TDFLIB
#define TUBEARGS see_progress=y
#include <SEP.defs>
#define FIGLIST subsine3 subsine5 il.1.10 il.2.15 il.2.10 il.3.20


#define DOTS Dots dots=0 connect=0 gaineach=0 constsep=1 strings=0 clip=1.2>/dev/null radius=.1

#define GRAPH Graph clip=1 max2=1 min2=-1 fastplot=20 > junk.H

default& : tube

$FIGDIR/il.%1.%2.v                    : Model.x model.p Levint.x  
FIGDIR/il.%1.%2.v FIGDIR/il.%1.%2.v3 : Model.x model.p Levint.x  
	Zero movie.H
	Model.x par=model.p job=%1          ncont=1000  > cont.H
	Model.x par=model.p job=%1 npick=%2 ncont=1000  > neighbor.H
	Model.x par=model.p job=%1 npick=%2             > pairs.H
	<pairs.H Levint.x na=5              n1=100 > ilint.H
	<movie.H Transp plane=23 > transp.H
	<transp.H GRAPH                      out=FIGDIR/il.%1.%2.v3
	<cont.H     DOTS par=model.p  connect=1 out=cont.v
	<ilint.H    DOTS par=model.p  strings=1 out=ilint.v
	<neighbor.H DOTS par=model.p     dots=2 out=neighbor.v
	vp_Overlay cont.v ilint.v neighbor.v > FIGDIR/il.%1.%2.v

#define O1D1 `<alias.P Get o1 d1`
FIGDIR/subsine%.v!            : Alias.x Levint.x
	Zero movie.H
	Alias.x tmax=80. pairs=1   n1=200 > alias.P
	Alias.x tmax=80. pairs=0   n1=200 > alias.H
	<alias.P Levint.x n1=200 O1D1 na=%> levint.H
	< movie.H Transp plane=23     > transp.H
	<levint.H GRAPH out=levint.v
	< alias.H GRAPH out=alias.v
	<transp.H GRAPH out=subsine%.v3
        < alias.H DOTS  strings=1 out=alias.v
        <levint.H DOTS  connect=1 out=levint.v
	vp_Overlay alias.v levint.v > FIGDIR/subsine%.v

clean&: jclean
	@-/bin/rm \*.P

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