Time-plotting the derivative of a species (reaction velocity) in catalyst

Hello,
I can’t find an easy way to time-plot a derivative of a species, although I’m sure data are stored somewhere. Thanks for help.
Andre

sol(t, Val{1}) is how it’s stored in the solution. ``sol(t, Val{1}; idxs = A) for speciesA. You can interpolate with a dense t` to then make the plot.