#! /usr/bin/env perl
#
# simple program reading a 2dim array (written previously by simple_xy_wr.pl)
#
# Heiko Klein <Heiko.Klein@met.no> 2007-02-06
#

use strict;
use warnings;

use PDL::Lite;
use PDL::NetCDF;
use Fcntl;


my $ncfile = new PDL::NetCDF("simple_xy.nc", {REVERSE_DIMS => 1,
					      MODE => O_RDONLY});
my $dimNames = $ncfile->getdimensionnames('data');
my @dimSizes = map {$ncfile->dimsize($_);} @$dimNames;
my @start = map {0} @dimSizes;
my $data = $ncfile->get('data', \@start, \@dimSizes);
# for reading complete arrays, the following can be used, too:
#my $data = $ncfile->get('data');
$ncfile->close;

print "data(", join(',', @$dimNames), "):\n";
print $data;