include ${SEPINC}/SEP.top

RESULTSER = transmissibility heater impWaveSnap impFlowSnap impDiffSnap resFlowSnap resWave1Snap resWave2Snap resDiffSnap 

${RESDIR}/transmissibility.v: ReservoirFactory.java
	${MAKE} ReservoirFactory.jc 
	${JX}.ReservoirFactory > transmissibility.H
	Taplot allpos=y < transmissibility.H >junk.A; 
	<junk.A Ta2vplot title=Transmissibility out=$@ >j.V

${RESDIR}/heater.v: HeatOpTest.java
	${MAKE} HeatOpTest.jc 
	${JX}.HeatOpTest > heater.H
	Taplot < heater.H >junk.A; <junk.A Ta2vplot out=heater.v >j.V
	Taplot < source.H >junk.A; <junk.A Ta2vplot out=source.v >j.V
	vp_SideBySideAniso heater.v source.v > $@

plot = Ta2vplot titlesz=24 labelsz=20 > junk.V
${RESDIR}/%Snap.v : 
	${MAKE} $*.H 
	< $*.H Taplot gainpanel=all > snap.A 
	< snap.A ${plot} title="Snapshots"       out=$@3
	< snap.A Window f3=10  n3=1 > junk.A 
	< junk.A ${plot} title="10  time steps" out=j1.v 
	< snap.A Window f3=30 n3=1 > junk.A 
	< junk.A ${plot} title="30  time steps" out=j2.v 
	< snap.A Window f3=50 n3=1 > junk.A 
	< junk.A ${plot} title="50  time steps" out=j3.v 
	< snap.A Window f3=70 n3=1 > junk.A 
	< junk.A ${plot} title="70  time steps" out=j4.v 
	< snap.A Window f3=90 n3=1 > junk.A 
	< junk.A ${plot} title="90  time steps" out=j5.v 
	< snap.A Window f3=110 n3=1 > junk.A 
	< junk.A ${plot} title="110 time steps" out=j6.v 
	< snap.A Window f3=130 n3=1 > junk.A 
	< junk.A ${plot} title="130 time steps" out=j7.v 
	< snap.A  Window f3=149 n3=1 > junk.A 
	< junk.A ${plot} title="150 time steps" out=j8.v 
	vp_OverUnderAniso  j1.v j3.v j5.v j.7.v > jA.v
	vp_OverUnderAniso  j2.v j4.v j6.v j.8.v > jB.v
	vp_SideBySideAniso jA.v jB.v            > $@

impWave.H: 
	${MAKE} ImpReflectedWave.jc
	${JX}.ImpReflectedWave sigma=.166 nt=150      > $@
impDiff.H   : 
	${MAKE} ImpReflectedWave.jc 
	${MAKE} ImpAbsorbedWave.jc
	#${JX}.ImpReflectedWave sigma=.166 nt=150      > imp1.H
	${JX}.ImpAbsorbedWave  sigma=.166 nt=150      > imp2.H
	Add < imp1.H imp2.H scale=1,-1 > $@
impFlow.H: 
	${MAKE} ImpulseFlow.jc
	${JX}.ImpulseFlow sigma=.166 nt=150  > $@

resFlow.H: ReservoirFlow.java
	${MAKE} ReservoirFlow.jc
	${JX}.ReservoirFlow sigma=.166 nt=150  > $@

resFlow1.H: resFlow.H
	Window n3=1 f3=10    < resFlow.H  > junky.H
	Add junky.H  scale=0.3            > junkx.H
	Add junkx.H    add=0.3 out=stdout > $@

resFlow2.H: resFlow.H 
	Window n3=1 f3=145   < resFlow.H  > junky.H
	Add junky.H  scale=0.3            > junkx.H
	Add junkx.H    add=0.3 out=stdout > $@

resWave1.H: resFlow1.H ReservoirWave.java
	${MAKE} ReservoirWave.jc 
	${JX}.ReservoirWave vel=resFlow1.H sigma=.166 nt=150  > $@

resWave2.H: resFlow2.H ReservoirWave.java
	${MAKE} ReservoirWave.jc 
	${JX}.ReservoirWave vel=resFlow2.H sigma=.166 nt=150  > $@

resDiff.H: resWave1.H resWave2.H 
	Add resWave1.H resWave2.H scale=-1.,1. > resDiff.H


clean : texclean jclean

include ${SEPINC}/SEP.bottom
