#define PRINT | Tube
#include <SEP.defs>
/* note that py-quantities are not generated in 2D modeling */


#define SNAPLIST pzsnap.f10.0   sxsnap.f10.90  szsnap.f10.0  pxsnap.f10.0                        pzsnap.f10.90  sysnap.f10.0   szsnap.f10.90 pxsnap.f10.90                       sxsnap.f10.0   sysnap.f10.90 

#define DIFFLIST pxsnap.diff pzsnap.diff sxsnap.diff sysnap.diff szsnap.diff

#define DIFFPOST px pz sx sy sz
poster&: [[sub X X.post DIFFPOST ]]

%.post&: %snap.f10.0 %snap.f10.90
	Window f3=5 j3=5 n3=3 <%snap.f10.0 datapath=/scr6/martin/ >j%.1
	Window f3=5 j3=5 n3=3 <%snap.f10.90 datapath=/scr6/martin/ >j%.2
	Add scale=1,-1 j%.1 j%.2 datapath=/scr6/martin/ >%snap.diff
	echo 'title="Wave Field Difference %"' >>%snap.diff
	Taplot <%snap.diff |Ta2vplot | Pspen  label=n
$	-RM_CMD j%.1 j%.2 %snap.diff


plot&: [[ sub -i X X.v SNAPLIST]]  [[ sub -i X X.v DIFFLIST]]

snapplot&: [[ sub -i X X.v SNAPLIST]]

/* make snapshot difference plots */
%.diff.v: %.diff
	Taplot <%.diff gainpanel=every |Ta2vplot out=FIGDIR/%.diff.v head=/dev/null

/* make snapshot plots */
%.v: % 
	Window f3=5 j3=5 n3=3  <% >j%$$
	Taplot <j%$$ gainpanel=every |Ta2vplot out=FIGDIR/%.v head=/dev/null
	-RM_CMD j%$$

/* calculate snapshot difference  */
%snap.diff: %snap.f10.0 %snap.f10.90
	Window f3=5 j3=5 n3=3 <%snap.f10.0 >j%$$.1
	Window f3=5 j3=5 n3=3 <%snap.f10.90 >j%$$.2
	Add scale=1,-1 jsub1jj jsub2jj >%snap.diff
	echo 'title="Wave Field Difference %"' >>%snap.diff
	-RM_CMD j%$$.1 j%$$.2

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