How to plot the parallelized example in "Parallelized Morris and Sobol Sensitivity Analysis of an ODE"

Hello!
I am trying GolbalSensitivity.jl for the first time, so please forgive me if this is a rookie mistake.

I have trouble with the example in the doc :

I reproduce the example until there :
sobol_result = gsa(f1, Sobol(), A, B, batch = true)
but I do not understand how I can plot the result obtained.

In this case, the sobol_result.ST returns a 4004 matrix instead of a 24 matrix for the non-parallized case.
Am I doing something wrong? Is something deprecated? How can I extract the result I want?

Thank you for your help!

The doc looks like it has a bug because it doesn’t pass the output_func so it doesn’t reduce the solution, it gives back the full 200 times timeseries for 2 species, hence the 400. I’ll fix that.