Follow-up to this: I opened an issue on github page for GlobalSensitivity: Handling NaNs for divergent trajectories in Sobol Method · Issue #185 · SciML/GlobalSensitivity.jl. I made a proposal that can possibly remedy the problem but don’t know if my proposal is correct or not.