Benoît Pasquier
Benoît Pasquier
See doi: 10.1029/GB001i001p00015
Maybe worth reworking the interface to be simpler and yielding more efficient code. I was thinking of something along the lines of 1. Define tracers, e.g., `DIP`, `DOP`, `POP`, `DFe`,...
This would allow me to spot the failing runs quicker.
I.e., the partition of the tracer equation into local sources and sinks and transport.
I find I could sometimes use a very simple 2-box model, so let's make one. Should probably go with the values in the Sarmiento and Gruber book or the reference...
Assuming one wants to simulate two independent tracers at the same time, there is no need to construct the full matrix. A better approach is to split the model into...
 [AIBECS_logo_Cartopy.pdf](https://github.com/briochemc/AIBECS.jl/files/3292836/AIBECS_logo_Cartopy.pdf)
Probably add something like ```julia print_type(io, f, val::Float64, ppu, s, n) = @printf io "%$(n)s = %8.2e [%s] %s\n" f val ppu s ``` where `n` is the maximum length...