#!/usr/bin/perl

$file = $ARGV[0];
($name = $file)=~s/\.dat//g;

open (IN, "< $file") || die("Cannot find $file\n");
open (OUT, "> $name.dx");

print OUT "# Water density on a regular grid\n";
print OUT "object 1 class gridpositions counts 17 51 51\n";
print OUT "origin  18.907 16.989 15.881\n";
print OUT "delta  0.5 0 0\n";
print OUT "delta  0 0.5 0\n";
print OUT "delta  0 0 0.5\n";
print OUT "object 2 class gridconnections counts 17 51 51\n";
print OUT "object 3 class array type double rank 0 items 44217 data follows\n\n";

$a=0; $t=0;
while (<IN>) {
	chomp $_;
	$line[$a] = $_;
	$a++; $t++;
	if ($a==3) {
		printf OUT $line[$a-3] .' '. $line[$a-2] .' '. $line[$a-1] ."\n";
		$a=0;
	}
}
close(IN);
$n = $t%3;
for ($a=0; $a<$n; $a++) {
	print OUT "$line[$a]\n";
}
close(OUT);