#!/usr/bin/perl
#a simple program to make sample data for use with conplot or gnuplot
#December 29, 1998
$outfile="squash.dat";
$im=16;
$jm=16;
$in=$im+1;
$jn=$jm+1;
open OUT,">$outfile";
$pi=3.14159265;
print OUT "# sample data made by makedata.pl\n";
print OUT "#-      x       |        y       |  intensity   |      u      |      v\n";
for ($j=0; $j<=$jm; $j++){
	for ($i=0; $i<=$im; $i++){
		$x=2*$i*$pi/$im;
		$y=2*$j*$pi/$jm;
		$z=sin($x)*sin($y);
		$y=$y*(1+.2*sin($x));
 		$x=$x*(1+.2*sin($y));
		$u=$y-3;
		$v=3-$x;
		print OUT sprintf("%14.6e  %14.6e  %14.6e  %14.6e  %14.6e\n",$x,$y,$z,$u,$v);
#		print OUT "$x\t$y\t$z\t$u\t$v\n";
	}
	print OUT "\n";
}
print "made data file $outfile\n";