include ${SEPINC}/SEP.top
RESDIR    = ../Fig
RESULTSER = pef3 pef2 pef1

tplt = Taplot > junky.A; < junky.A Ta2vplot titlesz=18 titlefat=3 >/dev/null 
grph = Dots titlesize=24 > junk.V

jshft = geom.shift

pef1=3
${RESDIR}/pef1.v : Sinus.java Pef.java 
	         ${JX}.Sinus n0=20                               > in.H
	< in.H   ${JX}.Pef     pef=${pef1}                       > ft.H
	< in.H   ${JX}.IcaiPef pef=${pef1}                       > ot.H
	< in.H   ${grph} title="Input"                           out=in.v
	< ft.H   ${grph} title="PEF"                             out=ft.v
	< ot.H   ${grph} title="Output (10^-7)"                  out=ot.v
	vp_OverUnderAniso in.v ot.v ft.v                         > $@ 

pef2=5,3
${RESDIR}/pef2.v : Pef.java
	${JAVAX} ${jshft}.SynSatTrac                             > in.H
	< in.H  ${JX}.Pef     pef=${pef2}                        > ft.H
	< in.H  ${JX}.IcaiPef pef=${pef2}                        > ot.H
	< in.H  ${tplt} title="Input"                            out=in.v
	< ft.H  ${tplt} title="PEF"                              out=ft.v
	< ot.H  ${tplt} title="Output (10^-6)"                   out=ot.v
	vp_SideBySideAniso in.v ft.v ot.v                        > $@ 

cplt = Taplot pclip=100 allpos=n > junky.A; < junky.A Cubeplot wheretitle="b" titlesz=18 titlefat=3 flat=n wantaxes=n wantcoordlabel=n point1=.66 point2=.66 >/dev/null frame1=0 frame2=5 frame3=0

aplt = Taplot gainpanel=all pclip=100  
win0 = Window n3=1 f3=0 
win1 = Window n3=1 f3=1 
vplt = Ta2vplot titlesz=18 titlefat=3 > junk.V 

pef3=5,3,2
${RESDIR}/pef3.v : Pef.java 
	${JAVAX} ${jshft}.PlaneWaveVol n0=6 n1=6 n2=6            > in.H
	< in.H   ${JX}.Pef     pef=${pef3}                       > ft.H 
	< in.H   ${JX}.IcaiPef pef=${pef3}                       > ot.H 
	< in.H   ${cplt} title="Input"                           out=in.v
	< ot.H   ${cplt} title="Output (10^-6)"                  out=ot.v
	< ft.H   ${aplt}                                         > ft.A
	< ft.A   ${win0} > j1.A; < j1.A ${vplt} title="PEF f3=0" out=j1.v
	< ft.A   ${win1} > j2.A; < j2.A ${vplt} title="PEF f3=1" out=j2.v
	vp_OverUnderAniso  j1.v j2.v                             > ft.v 
	vp_SideBySideAniso in.v ot.v ft.v                        > $@ 

clean : jclean
include ${SEPINC}/SEP.bottom
