How to preserve reproducibility when applying perturbations in the event of a numerical instability crash

Hi Tammas,

I was wondering if you would be happy to run through the steps you take to get your version of the script working on Gadi? I’m working through some model crashes and so it would be great to use your version which lets you set a random seed.

I’ve copied the script over, but think I’m getting problems with python and numpy versions. When I load only the pythonlib/umfile_utils module, I receive an error about the ‘default_rng’ attribute not existing, while if I additionally load the conda/analysis3 module, numpy fails to import.

Many thanks,
Spencer