#!/usr/bin/env python

import sys
import platform
import numpy as np
if platform.system() == 'Windows':
    print '...importing windows dummy'
    import windummy as dum
elif platform.system() == 'Linux':
    print '...importing linux dummy'
    import lindummy as dum
else:
    print 'ERROR:   Cannot import suitable binary modules.'
    sys.exit()

# load
dum.rst.load(True)

# dimensionalise
nstep, nx, ny, nz = dum.rst.dummy.shape

# transpose to local copy of dummy array
locdum = np.transpose(dum.rst.dummy.copy(), axes=(3,2,1,0)).flatten()
locdum.shape = (nstep, nx, ny, nz)

# output in the order I would expect post-transpose
for t in range(nstep):
    for x in range(nx):
        for y in range(ny):
            print  '\t%2d\t%2d\t%2d\t' % (t,x,y),
            check = []
            for z in range(nz):
                check.append(locdum[t,x,y,z])

            print  ''.join(check)