I believe I fixed the memory leak. It's looking good now. I'll have it run it for a while and check then however. If it's not fixed it is dramatically reduced. The memory leak was fixed by writing destroy_model before copy_model. Putting destroy_model in copy_model didn't work - I have no clue why. OpenMP is currently on by the way.

I also implemented the multislice algorithm, although I haven't actually tried to use it yet. Setting use_multislice to true in rmc.f90 enables multislice every 10k moves.