Fabian Bernhard

Note that if the ODE doesn't return `0.`, the callback at `t=0.` seems to work. ```Julia # With callback at 0.0, but f not returning 0.0, it solves without error...

Thanks, @ChrisRackauckas . Glad to hear that you are aware of the cause of this behaviour. Indeed, `dt=1e-3` solves the issue, and should be perfectly fine in my use case....

Thanks for this helpful answer, Jonnie! The first part solves my problem perfectly well. For the second part regarding multidimensional `ComponentMatrix`: in your example the second axis is simply `FlatAxis()`....

The underlying code in the package `LWFBrook90.jl` simplifies to the following MWE: ```julia using ComponentArrays # version: [b0b7db55] ComponentArrays v0.15.2 u0 = ComponentArray( GWAT = (mm = 0.0, d18O =...

Actually there might be two ways to fix: - either set `mins` and `maxs` based on `plot_vecs` as mentioned at bottom of first post - or refrain from setting `xlims`...