UF77LIBS = ${VPLOTLIB} ${TDFLIB}
INVRAS = y
include ${SEPINC}/SEP.top
RESULTSER = vec16 ras16 vec18 ras18 vec20 ras20 ras21
RESULTSNR = tenthou
RESULTSCR = last18 lap18 lint18
RESULTSER = #
RESULTSCR = last18 lap18 lint18 compare

TA2VPLOT = Ta2vplot crowd=.85 >/dev/null transp=n yreverse=n
TA2VPLOT = Ta2vplot Xcrowd    >/dev/null
TA2VPLOT = Ta2vplot crowd=.82 >/dev/null
#PSARGS = invras=y color=n
DATA =  ../../Data

default: tube

wind%.H : ${DATA}/apr%.H	
	< ${DATA}/apr$*.H                        > wind$*.H   Window    Xj2=13

COMPARE =  ${RESDIR}/last18.v ${RESDIR}/lap18.v ${RESDIR}/lint18.v
${RESDIR}/compare.v ${RESDIR}/compare.v3: ${COMPARE}
	vp_Movie           ${COMPARE} >         ${RESDIR}/compare.v3
	vp_SideBySideAniso ${COMPARE} > junk.v
	<junk.v  vppen vpstyle=n txscale=2. >${RESDIR}/compare.v

${RESDIR}/last18.v ${RESDIR}/lap18.v : wind18.H interp2.x     interp.p
	Zero last.H
	< wind18.H interp2.x par=interp.p warmup=99999 > lap.H 
	cat lap.H last.H | Byte pclip=100. > last.A
	cat lap.H  lap.H | Byte pclip=100. > lap.A
	<last.A ${TA2VPLOT} title=April_18_LastInBin  out=${RESDIR}/last18.v
	< lap.A ${TA2VPLOT} title=April_18_Laplacian  out=${RESDIR}/lap18.v

${RESDIR}/lint18.v ${RESDIR}/lint18.v3 : wind18.H interp2.x     interp.p
	Zero last.H
	< wind18.H interp2.x par=interp.p warmup=20    > int.H 
	cat int.H last.H | Byte pclip=100. > last.A
	cat int.H  int.H | Byte pclip=100. > int.A
	< int.A ${TA2VPLOT} title=April_18_PEF        out=int.v
	<last.A ${TA2VPLOT} title=April_18_LastInBin  out=last.v
	vp_Movie           last.v int.v >         ${RESDIR}/lint18.v3
	cp int.v                             ${RESDIR}/lint18.v
#	vp_SideBySideAniso last.v int.v > junk.v
#	<junk.v  vppen vpstyle=n txscale=2. >${RESDIR}/lint18.v

#							raster plot
${RESDIR}/ras%.v : ${DATA}/apr%.H shift.x     shift.%.p
	<       ${DATA}/apr$*.H shift.x par=shift.$*.p  > junk.H
	<junk.H Byte pclip=99.99 > shift.$*.A
	<shift.$*.A ${TA2VPLOT} title=April_$*   out=${RESDIR}/ras$*.v

#blinker : blink16.v3 blink18.v3 blink20.v3 blink21.v3
blinker :
	tube blink16.v3
	tube blink18.v3
	tube blink20.v3
	tube blink21.v3

BLINK = ${RESDIR}/ras%.v ${RESDIR}/vec%.v
blink%.v3 :      ${BLINK}
	vp_Movie ${subst %,$*,${BLINK}} > blink$*.v3

peek :
	<mm.H Byte gainpanel=every pclip=100> mm.A
	< mm.A Ta2vplot crowd=1. fastplot=20 > mm.V out=junk.v
	tube junk.v

movie : movie.v3
	tube movie.v3

FRAMELIST = F/vec16.v F/ras16.v F/vec18.v F/ras18.v F/vec20.v F/ras20.v
movie.v3 :       ${FRAMELIST}
	vp_Movie ${FRAMELIST} > movie.v3

#							vector plot
vec% ${RESDIR}/vec%.v : ${DATA}/apr%.H track.x
	<       ${DATA}/apr$*.H track.x  date=$* > junk.V
	Cp junk.V junk.V1 out=${RESDIR}/vec$*.v


MERGE = ${DATA}/apr16.H ${DATA}/apr18.H ${DATA}/apr20.H
${DATA}/apr21.H :    ${MERGE}
	Merge <${MERGE} space=n axis=2 > ${DATA}/apr21.H out=${DATA}/apr21.F

#                                            	Original Alistair asci to seplib.
#${DATA}/apr%.H :  New.x apr%.p ${DATA}/apr%.xyz.Z
#		cp      ${DATA}/apr%.xyz.Z fort.8.Z
#		uncompress fort.8.Z
#		New.x par=apr%.p > ${DATA}/apr%.H out=${DATA}/apr%.F
#		/bin/rm fort.8

clean: jclean

include ${SEPINC}/SEP.bottom
