simple example code: save_and_load_particles.pyΒΆ

(Source code, png, hires.png, pdf)

Saves and loads a set of particles.
The 'write_set_to_file' and 'read_set_from_file' functionso can hadle
files in a variety of formats.

from amuse.units import nbody_system
from amuse.ic.plummer import new_plummer_model
from amuse.plot import plot, native_plot

# the imports needed to read and write a set
from import write_set_to_file
from import read_set_from_file

import os.path

if __name__ == "__main__":
    if os.path.exists('plummer128.hdf'):

    # generate a particle set

    # write the set to file 'testfile'
    # the third argument is the file format, 'amuse' is an hdf5 based format
    # that saves all information. Other formats
    # available are e.g. csv, txt, gadget, starlab

    # reading back the file
    # we close the file, causing all data to be copied to memory
    # instead of being read from an open file
    particles=read_set_from_file('plummer128.hdf','amuse', close_file = True)

    # plotting
    plot(particles.x, particles.y,'r.')

Keywords: python, amuse, astrophysics, matplotlib, pylab, example, codex (see how-to-search-examples)