Status

Ongoing
Python front end complete, together with routines for cooling, star formation (including star formation histories) and feedback. Have begun to convert the physics routines back into C, using ctypes to enable them to be called from the python front-end. This seems to work fine, so simply needs to be extended to cover all the required physics modules.

Collaborators

Peter Thomas, Will Roper

Description

The purpose of this project is create a python front-end for L-galaxies to make the reading in of parameters much easier (and definable at run time) and to simply the I/O. At the same time, the code has been restructured so as to process (sub)halos from past to present, one snapshot at a time, so that a maximum of 2 snapshots need to be stored in memory at any given time.
The halos and subhalos are stored as python lists of instances of halo and subhalo classes. The galaxies are stored in numpy arrays within each snapshot. All data is output in HDF files.
The majority of the physics is carried out by the usual L-galaxies C-routines that will be given headers to convert them into python modules that can be imported at run-time (under development).