Hi there-- The covering_grid returns arrays that are in Enzo (thus Fortran) ordering, correct? So if I do open(filename,'wb').write( h.covering_grid(0,left,resolution)[ 'Density'] ) I should wind up with a binary file in the same order as the original enzo data? d. -- Sent from my Stone Tablet and carried by my Pterodactyl.