#define DEBUG
#define CMLIB -v -lcmssl -lcmtiming 
#define XCMLIB -lcmfortgraph_slice -lcmsr_slice -lX11
#define XCMLIBOLD -lcmfortgraph -lcmsr -lX11
#define BINDIR /r3/martin/bin/MTYPE
#define ULTIMOD Ultimod


#define SEQ -S0-1
#include <SEP.defs>

#define CMFCOPTS  -extend_source -cmsim -cmdebug
test.x: test.o
	CMFC CMFCOPTS -c test.fcm
	CMFC CMFCOPTS -o test.x test.o SEPLIB SYSLIB

Test.x: Test.o
	CMFC CMFCOPTS -o Test.x Test.o SEPLIB SYSLIB

#define RebuildableFortran  ( cando %.f and exist %.f )
#define RebuildableFortran  (not {{cake  -q junk.f }}) and exist junk.f

%.fclean&: if     RebuildableFortran
	@-/bin/rm %.f
	
%.fclean&: if not RebuildableFortran
	@echo %.f not removed since it cannot be rebuilt

#define RebuildableCMFortran  cando %.fcm and exist %.fcm 
#define RebuildableCMFortran (not {{ cake -q %.fcm  }}) and exist %.fcm 

%.fcmclean&: if     RebuildableCMFortran
	@-/bin/rm %.fcm
	
%.fcmclean&: if not RebuildableCMFortran
	@echo %.fcm not removed since it cannot be rebuilt
	
dave&:
	Testsec.x <homoiso.h  >junk 

/* reciprocal traces and source positions are by one, due to staggered grid */
testreci1&:  homoiso.h wavelet.vert   /* looks ok */
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoiso.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=52 src_depth=104 geo_depth=104 xseis=xseis1 zseis=zseis1
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoiso.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=156 src_depth=104 geo_depth=104 xseis=xseis2 zseis=zseis2
	Window f1=155 n1=1 < zseis1 > jtrace1
	Window f1=51  n1=1 < zseis2 > jtrace2
	Merge jtrace1  jtrace2 axis=2 space=n > junk
	Graph <junk |Tube

testreci2&:  homoisorho.h wavelet.vert   /* looks ok */
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoisorho.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=52 src_depth=104 geo_depth=104 xseis=xseis1 zseis=zseis1
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoisorho.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=156 src_depth=104 geo_depth=104 xseis=xseis2 zseis=zseis2
	Window f1=155 n1=1 < zseis1 > jtrace1
	Window f1=51  n1=1 < zseis2 > jtrace2
	Merge jtrace1  jtrace2 axis=2 space=n > junk
	Graph <junk |Tube

testreci3&:  homoisoc.h wavelet.vert   /* looks ok */
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoisoc.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=52 src_depth=104 geo_depth=104 xseis=xseis1 zseis=zseis1
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoisoc.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 src_1=156 src_depth=104 geo_depth=104 xseis=xseis2 zseis=zseis2
	Window f1=155 n1=1 < zseis1 > jtrace1
	Window f1=51  n1=1 < zseis2 > jtrace2
	Merge jtrace1  jtrace2 axis=2 space=n > junk
	Graph <junk |Tube





test&:  homoiso.h wavelet.vert 
	time  BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoiso.h par=parfile.test  <wavelet.vert >junk  ncoef=0 coeff=ceoff.H dtseis=0.008 nsrc=2 src_inc=10 


split&: homoiso3c.h wavelet.vert
	time  BINDIR/Ultimod.split type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  xsnap=xsnap.split ysnap=ysnap.split zsnap=zsnap.split xseis=xseis.split yseis=yseis.split zseis=zseis.split <wavelet.vert >junk


bound&: /* homoac.h wavelet.vert */
	Window n1=200 n2=200 <homoac.h >jjj
	time BINDIR/Ultimod  oneway=0 damping=0.002 type=1 modeling=1 moduli=jjj par=parfile.bound <wavelet.vert >junk nt=400



testsuite higher&: ISO TI ORTHO MONO TRIC

iso ISO&: homojustiso.h wavelet.vert
        time BINDIR/Ultimod  type=2 src_type=0 ncomp=2 modeling=1 moduli=homojustiso.h par=parfile  xsnap=xsnap.ISO ysnap=ysnap.ISO zsnap=zsnap.ISO xseis=xseis.ISO yseis=yseis.ISO zseis=zseis.ISO<wavelet.vert >junk

ti TI&: homoiso3c.h wavelet.vert
	time BINDIR/ULTIMOD type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  xsnap=xsnap.TI ysnap=ysnap.TI zsnap=zsnap.TI xseis=xseis.TI yseis=yseis.TI zseis=zseis.TI <wavelet.vert >junk

ortho ORTHO&: ortho.h wavelet.vert
	time BINDIR/Ultimod  type=5 src_type=0 ncomp=3 modeling=1 moduli=ortho.h par=parfile  xsnap=xsnap.ORTHO ysnap=ysnap.ORTHO zsnap=zsnap.ORTHO xseis=xseis.ORTHO yseis=yseis.ORTHO zseis=zseis.ORTHO <wavelet.vert >junk

mono MONO&: mono.h wavelet.vert
	time BINDIR/Ultimod  type=6 src_type=0 ncomp=3 modeling=1 moduli=mono.h par=parfile  xsnap=xsnap.MONO ysnap=ysnap.MONO zsnap=zsnap.MONO xseis=xseis.MONO yseis=yseis.MONO zseis=zseis.MONO <wavelet.vert >junk

tric TRIC&: triclinic.h wavelet.vert
	time BINDIR/Ultimod  type=7 src_type=0 ncomp=3 modeling=1 moduli=triclinic.h par=parfile  xsnap=xsnap.TRIC ysnap=ysnap.TRIC zsnap=zsnap.TRIC xseis=xseis.TRIC yseis=yseis.TRIC zseis=zseis.TRIC <wavelet.vert >junk


/* derivative coefficients*/
coeff.H:
$        very poor
        echo "-1. 1." |atoF >coeff.H
$
$        Francis' central limit
$        echo "0.0006975446 -0.0095703125 0.079752604 -1.1962891 1.1962891 -0.07
$ 9752604 0.0095703125 0.0006975446" |atoF >coeff.H
$
$       Holberg's coefficients
$       echo "0.1326224603581672E-04 -0.1260422755482438E-03 0.6766549997685445E
$-03 -0.2713077307504363E-02 0.9277622829117149E-02 -0.3045387209576234E-01 0.117
$8572129899180E+00 -0.1247830531633298E+01 0.1247830531633298E+01 -0.117857212989
$9180E+00  0.3045387209576234E-01 -0.9277622829117149E-02  0.2713077307504363E-02
$ -0.6766549997685445E-03  0.1260422755482438E-03 -0.1326224603581672E-04" | atoF
$ >coeff.H


ovacoust:
	Window f3=3 n3=1 <zsnap >jj
$   4th snapshot:     0.8 * 4 kms = 3.2
	Circle <jj radius=3.2 x0=4.12 z0=0.8 >junk
	Graph <junk yreverse=1 symbol="*" out=jj2.v head=/dev/null min1=0 min2=0. max2=8.32 max1=8.32
	Byte <jj |Ta2vplot out=jj.v head=/dev/null
	tube jj.v jj2.v erase=once

ovacoust2:
	Window f3=7 n3=1 <zsnap2 >jj
$   4th snapshot:     (0.8-0.04) * 4 kms = 3.04
$   4th snapshot:     (2*0.8-0.04) * 4 kms = 6.16
	Circle <jj radius=6.4 x0=4.12 z0=0.00 >junk
	Graph <junk yreverse=1 symbol="*" out=jj2.v head=/dev/null min1=0 min2=0. max2=8.32 max1=8.32
	Byte <jj |Ta2vplot out=jj.v head=/dev/null
	tube jj.v jj2.v erase=once

acoust&: homoac.h wavelet.vert
	time BINDIR/Ultimod  type=1 modeling=1 moduli=homoac.h par=parfile  <wavelet.vert >junk

acoust2&: homoac.h wavelet.vert60
	time BINDIR/Ultimod  type=1 modeling=1 moduli=homoac.h par=parfile  src_depth=0 surf_type=0 <wavelet.vert60 >junk zsnap=zsnap2 nt=400

vert2c.iso&: homojustiso.h wavelet.vert
	time BINDIR/Ultimod  type=2 src_type=0 ncomp=2 modeling=1 moduli=homojustiso.h par=parfile  <wavelet.vert >junk

vert2c&: homoiso.h wavelet.vert
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoiso.h par=parfile  <wavelet.vert >junk

inline2c&: homoiso.h wavelet.inline
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=2 modeling=1 moduli=homoiso.h par=parfile  <wavelet.inline >junk

pressure2c&: homoiso.h wavelet.vert
	time BINDIR/Ultimod  type=4 src_type=1 ncomp=2 modeling=1 moduli=homoiso.h par=parfile  <wavelet.vert >junk

vert3c&: homoiso3c.h wavelet.vert
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  xsnap=xsnap.vert3c ysnap=ysnap.vert3c zsnap=zsnap.vert3c xseis=xseis.vert3c yseis=yseis.vert3c zseis=zseis.vert3c <wavelet.vert >junk

inline3c&: homoiso3c.h wavelet.inline
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  xsnap=xsnap.inline ysnap=ysnap.inline zsnap=zsnap.inline xseis=xseis.inline yseis=yseis.inline zseis=zseis.inline <wavelet.inline >junk

xline3c&: homoiso3c.h wavelet.xline
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  <wavelet.xline >junk


pressure3c&: homoiso3c.h wavelet.vert
	time BINDIR/Ultimod  type=4 src_type=1 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile  <wavelet.vert >junk


inline+45&: wavelet+i45.h homoiso3c.h
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile xsnap=xsnap.+i45 ysnap=ysnap.+i45 zsnap=zsnap.+i45 xseis.+i45 yseis.+i45 zseis.+i45 <wavelet+i45.h >junk

inline-45&: wavelet-i45.h homoiso3c.h
	time BINDIR/Ultimod  type=4 src_type=0 ncomp=3 modeling=1 moduli=homoiso3c.h par=parfile xsnap=xsnap.-i45 ysnap=ysnap.-i45 zsnap=zsnap.-i45 xseis.-i45 yseis.-i45 zseis.-i45 <wavelet-i45.h >junk



subtraction&: /* inline+45 inline-45 vert3c inline3c */
	Add xsnap.+i45 xsnap.-i45 mode=sum >xsnap.add.+i45
	Add ysnap.+i45 ysnap.-i45 mode=sum >ysnap.add.+i45
	Add zsnap.+i45 zsnap.-i45 mode=sum >zsnap.add.+i45
	Add xsnap.+i45 xsnap.-i45 mode=sum scale=1.,-1. >xsnap.sub.+i45
	Add ysnap.+i45 ysnap.-i45 mode=sum scale=1.,-1. >ysnap.sub.+i45
	Add zsnap.+i45 zsnap.-i45 mode=sum scale=1.,-1. >zsnap.sub.+i45

viscoacoust&: /* homoac.h wavelet.vert */
$	Scale <homoac.h >homoviscoac.h dscale=0.
	time BINDIR/Ultivisco  type=1 modeling=1 moduli=homoac.h viscomoduli=homoviscoac.h par=parfile  <wavelet.vert >junk

alliso&: homoiso.h wavelet.vert 
	Scale <homoiso.h >homoviscoiso.h dscale=0.0
	time BINDIR/Ultiall visco=0  type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.iso xseis=xseis.iso nt=400
$	time BINDIR/Ultiall  type=2 modeling=1 moduli=homoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.iso.all xseis=xseis.iso.all nt=400
$	time BINDIR/Ultimod  type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.iso xseis=xseis.iso nt=400
$	Scale <homoiso.h >homoviscoiso.h dscale=0.005
	Scale <homoiso.h >homoviscoiso.h dscale=0.05
	time BINDIR/Ultiall visco=1 type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.viscoiso xseis=xseis.viscoiso nt=400

viscoiso&: homoiso.h wavelet.vert 
	Scale <homoiso.h >homoviscoiso.h dscale=0.0
$	time BINDIR/Ultivisco visco=1  type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.iso xseis=xseis.iso nt=400
	time BINDIR/Ultimod  type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.iso xseis=xseis.iso nt=400
$	Scale <homoiso.h >homoviscoiso.h dscale=0.005
	Scale <homoiso.h >homoviscoiso.h dscale=0.05
	time BINDIR/Ultivisco visco=1 type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h par=parfile  <wavelet.vert >junk zseis=zseis.viscoiso xseis=xseis.viscoiso nt=400

rho&:
	Scale <homoiso.h >homoviscoiso.h dscale=0.0
	Window f3=2 n3=1 <homoiso.h >jJrho 
	Cat jJrho jJrho jJrho >rho.H
	time BINDIR/Ultiall viscorho=1 type=2 modeling=1 moduli=homoiso.h viscomoduli=homoviscoiso.h density=rho.H par=parfile  <wavelet.vert >junk zseis=zseis.rhoiso xseis=xseis.rhoiso nt=400

/* make a two component source */
wavelet.vert:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Spike n1=128 n2=1 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk junk2 space=n axis=2 >wavelet.vert

wavelet.vert60:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.040 fund=30 fhigh=60. | Scale dscale=1000000. >junk
	Spike n1=128 n2=1 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk junk2 space=n axis=2 >wavelet.vert60

wavelet+i45.h:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Scale <junk dscale=0.7071 >junk2
	Scale <junk dscale=0.7071 >junk3
	Spike n1=128 n2=1 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk4
	Merge junk2 junk3 junk4 space=n axis=2 >wavelet+i45.h


wavelet-i45.h:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Scale <junk dscale=0.7071 >junk2
	Scale <junk dscale=-0.7071 >junk3
	Spike n1=128 n2=1 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk4
	Merge junk2 junk3 junk4 space=n axis=2 >wavelet-i45.h

wavelet.inline:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Spike n1=128 n2=1 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk2 junk space=n axis=2 >wavelet.inline

wavelet.xline:
	Wavelet n1=128 d1=0.004 domain=time wavelet=ricker1 tdelay=0.080 fund=15 fhigh=30. | Scale dscale=1000000. >junk
	Spike n1=128 n2=2 n3=1 d1=0.004 k1=-1 k2=-1 k3=-1 >junk2
	Merge junk2 junk space=n axis=2 >wavelet.xline


homoiso.h:
	Geninput nx=208 nz=208 dx=40. dz=40.  symmetry=2 vp=4.0 vs=2.3  rho=2.5   >homoiso.h

homoisorho.h:
	Window f3=4 n3=1 <homoiso.h >rho.h
	Window f3=0 n3=4 <homoiso.h >rest.h
	Window f2=0   n2=104 <rho.h >junk1
	Window f2=104 n2=104 <rho.h >junk2
	Scale <junk2 dscale=0.9 >junk3
	Merge junk1 junk3 axis=2 space=n >rhonew.h
	Merge rest.h rhonew.h axis=3 space=n > homoisorho.h

homoisoc.h:
	Window f3=4 n3=1 <homoiso.h >rho.h
	Window f3=0 n3=4 <homoiso.h >rest.h
	Window f2=0   n2=104 <rest.h >junk1
	Window f2=104 n2=104 <rest.h >junk2
	Scale <junk2 dscale=0.8 >junk3
	Merge junk1 junk3 axis=2 space=n >new.h
	Merge new.h rho.h axis=3 space=n > homoisoc.h

homojustiso.h:
	Window n3=1 <homoiso.h >junk.c11
	Window f3=3 n3=1  <homoiso.h >junk.c55
	Window f3=4 <homoiso.h >junk.rho
	Merge junk.c11 junk.c55 junk.rho axis=3 space=n >homojustiso.h

homoac.h:
	Window n3=1 <homoiso.h >junk.c11
	Window f3=4 <homoiso.h >junk.rho
	Merge junk.c11 junk.rho axis=3 space=n >homoac.h

homoiso3c.h: homoiso.h
	Window n3=4 <homoiso.h >junk      /* get sitffnesses only */
	Window f3=4 <homoiso.h >junk2     /* get rho */
	Window f3=3 n3=1 <junk >junk3    /* get c55*/
	Merge junk junk3 junk2 axis=3 space=n >homoiso3c.h
	RM_CMD junk*
	
ortho.h: homoiso.h
	Window n3=4 <homoiso.h >junk      /* get sitffnesses only */
	Window f3=4 <homoiso.h >junk2     /* get rho */
	Window f3=3 n3=1 <junk >junk3    /* get c55 */
	Cp <junk3 >junk4
	Merge junk junk3 junk4 junk2 axis=3 space=n >ortho.h
	RM_CMD junk*

mono.h: homoiso.h
	Window f3=0 n3=1 <homoiso.h >junk.c11
	Window f3=1 n3=1 <homoiso.h >junk.c13
	Window f3=2 n3=1 <homoiso.h >junk.c33
	Window f3=3 n3=1 <homoiso.h >junk.c55
	Window f3=4 n3=1 <homoiso.h >junk.rho
	Spike n1=208 n2=208 d1=0.040 d2=0.040 nsp=0 >junk.c14
	Cp <junk.c13 >junk.c12
	Cp <junk.c14 >junk.c15
	Cp <junk.c11 >junk.c22
	Cp <junk.c13 >junk.c23
	Cp <junk.c14 >junk.c25
	Cp <junk.c14 >junk.c26
	Cp <junk.c14 >junk.c35
	Cp <junk.c55 >junk.c44
	Cp <junk.c14 >junk.c46
	Cp <junk.c55 >junk.c66
	Merge junk.c11 junk.c12 junk.c13 junk.c15 junk.c22 junk.c23 junk.c25 junk.c33 junk.c35 junk.c44 junk.c46 junk.c55 junk.c66 junk.rho axis=3 space=n >mono.h
	RM_CMD junk*

triclinic.h: homoiso.h
	Window f3=0 n3=1 <homoiso.h >junk.c11
	Window f3=1 n3=1 <homoiso.h >junk.c13
	Window f3=2 n3=1 <homoiso.h >junk.c33
	Window f3=3 n3=1 <homoiso.h >junk.c55
	Window f3=4 n3=1 <homoiso.h >junk.rho
	Cp <junk.c13 >junk.c12
	Spike n1=208 n2=208 d1=0.040 d2=0.040 nsp=0 >junk.c14
	Cp <junk.c14 >junk.c15
	Cp <junk.c14 >junk.c16
	Cp <junk.c11 >junk.c22
	Cp <junk.c13 >junk.c23
	Cp <junk.c14 >junk.c24
	Cp <junk.c14 >junk.c25
	Cp <junk.c14 >junk.c26
	Cp <junk.c14 >junk.c34
	Cp <junk.c14 >junk.c35
	Cp <junk.c14 >junk.c36
	Cp <junk.c55 >junk.c44
	Cp <junk.c14 >junk.c45
	Cp <junk.c14 >junk.c46
	Cp <junk.c14 >junk.c56
	Cp <junk.c55 >junk.c66
	Merge junk.c11 junk.c12 junk.c13 junk.c14 junk.c15 junk.c16 junk.c22 junk.c23 junk.c24 junk.c25 junk.c26 junk.c33  axis=3 space=n >jj
	Merge jj junk.c34 junk.c35 junk.c36 junk.c44 junk.c45 junk.c46 junk.c55 junk.c56 junk.c66 junk.rho axis=3 space=n >triclinic.h
	RM_CMD junk* jj*
	

isolayer.h:
	Geninput nx=208 nz=70 dx=20. dz=20.  symmetry=2 vp=4.0 vs=2.3  rho=2.5   >junk0
	Geninput nx=208 nz=138 dx=20. dz=20.  symmetry=2 vp=4.5 vs=3.2  rho=2.0   >junk1
	Merge junk0 junk1  space=n axis=2 |Transp plane=12 >isolayer.h



layer.h: 
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 vp=4.0 vs=3.3  rho=2.5   >junk1
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 >junk0
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 vp=3.0 vs=2.0  rho=2.5     >junk3
	Geninput nx=208 nz=50 dx=40. dz=15. symmetry=2 c11=34.3 c13=10.7 c33=22.7 c44=5.4 c66=10.6 rho=2.5 >junk2
	Geninput nx=208 nz=8 dx=40. dz=15.  symmetry=2 vp=4.0 vs=2.3  rho=2.5     >junk4
	Merge junk0 junk1 junk2 junk3 junk4 space=n axis=2 |Transp plane=12 >layer.h


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

%.T:	
	Taplot <% gainpanel=every  >%.T

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