Benoît Pasquier

Gentle bump! :)

Hi @jonathansharp, No worries! I don't think merging would have reverted anything on your side (unless there is a conflict) but I just updated my branch (I pulled your changes...

Sorry for digging this up, but wanted to link this issue to that one on Unitful.jl: https://github.com/PainterQubits/Unitful.jl/issues/240, in case someone wants to add the functionality here into Unitful.jl.

FWIW, I find that in this case something like ```julia sum(x -> !isnan(x) * x, arr, dims=1) ``` (i.e., without using NaNMath) works just as well. (And is maybe a...

I can think of a workaround but FWIW `flatten` is actually not an AIBECS function but one from the aptly named [Flatten.jl](https://github.com/rafaqz/Flatten.jl) package! (AIBECS uses it to go from parameter...

Done (https://github.com/briochemc/AIBECS.jl/commit/d5e55585f7892c5170b47240904ed01696e307fc)

Maybe rename it as the Broecker model? http://weatherclasses.com/uploads/3/6/2/3/36231461/matlab_example_ode45.pdf

Thank you @zhenwu0728 for your comments. Sorry for the delayed response! I have a question regarding the "better `show`" comment: > * A better `show()` is needed for `OCIM2.load()`, `F_and_∇ₓF`,...

I gave a bit of thought to the issue with `F_and_∇ₓF` and think I have a decent solution. There is the `AIBECSFunction` that one can use instead, which returns a...

@zhenwu0728 FYI, there is now [an ideal-age tutorial deployed from the linked PR](https://juliaocean.github.io/AIBECS.jl/previews/PR85/tutorials/1_ideal_age/) which shows how this looks when using `AIBECSFunction` instead of `F_and_∇ₓF`. Does that seem like a good...