Segfault with large values of mesh_side_length
I'm trying to run the 100Mpc and 200Mpc EAGLE DMONLY boxes on Irene. Initially I assumed that I should scale mesh_side_length with the simulation volume, but I find that if I set it too large the code segfaults despite having plenty of memory available. Looking at mesh_gravity.c I think the problem is that the mesh array (which has size mesh_side_length^3 ) is indexed using ints. This means we crash if mesh_side_length > ~1290.
Maybe the answer is that I shouldn't be trying to use such huge mesh sizes, but if so we should at least stop with an error message if mesh_side_length is too large.