GitHub

Box model

In this example we setup a LOBSTER biogeochemical model in a single box configuration. This example demonstrates:

  • How to setup OceanBioME's biogeochemical models as a stand-alone box model

Install dependencies

First we check we have the dependencies installed

using Pkg
pkg"add OceanBioME"

Model setup

Load the packages and setup the initial and forcing conditions

using OceanBioME, Oceananigans, Oceananigans.Units
using Oceananigans.Fields: FunctionField

const year = years = 365day

grid = BoxModelGrid()
clock = Clock(time = zero(grid))

This is forced by a prescribed time-dependent photosynthetically available radiation (PAR)

PAR⁰(t) = 60 * (1 - cos((t + 15days) * 2π / year)) * (1 / (1 + 0.2 * exp(-((mod(t, year) - 200days) / 50days)^2))) + 2

const z = -10 # specify the nominal depth of the box for the PAR profile
PAR_func(t) = PAR⁰(t) * exp(0.2z) # Modify the PAR based on the nominal depth and exponential decay

PAR = FunctionField{Center, Center, Center}(PAR_func, grid; clock)

Set up the model. Here, first specify the biogeochemical model, followed by initial conditions and the start and end times

model = BoxModel(; biogeochemistry = LOBSTER(; grid, light_attenuation = PrescribedPhotosyntheticallyActiveRadiation(PAR)),
                   clock)

set!(model, NO₃ = 10.0, NH₄ = 0.1, P = 0.1, Z = 0.01)

simulation = Simulation(model; Δt = 5minutes, stop_time = 5years)

simulation.output_writers[:fields] = JLD2Writer(model, model.fields; filename = "box.jld2", schedule = TimeInterval(10days), overwrite_existing = true)

prog(sim) = @info "$(prettytime(time(sim))) in $(prettytime(simulation.run_wall_time))"

simulation.callbacks[:progress] = Callback(prog, IterationInterval(1000000))
Callback of prog on IterationInterval(1000000)

Run the model (should only take a few seconds)

@info "Running the model..."
run!(simulation)
[ Info: Running the model...
[ Info: Initializing simulation...
[ Info: 0 seconds in 0 seconds
[ Info:     ... simulation initialization complete (1.255 seconds)
[ Info: Executing initial time step...
[ Info:     ... initial time step complete (501.908 ms).
[ Info: Simulation is stopping after running for 15.215 seconds.
[ Info: Simulation time 1825 days equals or exceeds stop time 1825 days.

Load the output

times = FieldTimeSeries("box.jld2", "P").times

timeseries = NamedTuple{keys(model.fields)}(FieldTimeSeries("box.jld2", "$field")[1, 1, 1, :] for field in keys(model.fields))
(NO₃ = [10.0, 10.015174865722656, 9.961457252502441, 9.81340503692627, 9.481049537658691, 8.78665542602539, 7.772372245788574, 6.71930456161499, 5.456292152404785, 3.7530651092529297, 1.6817322969436646, 0.22698144614696503, 0.014157782308757305, 0.0024202847853302956, 0.0012515581911429763, 0.0010644716676324606, 0.001099957968108356, 0.0012859976850450039, 0.0016993563622236252, 0.002630256349220872, 0.005181585904210806, 0.01364092342555523, 0.04050086438655853, 0.11302825808525085, 0.28080010414123535, 0.6053680181503296, 1.1175591945648193, 1.8015879392623901, 2.6145544052124023, 3.5032308101654053, 4.415060997009277, 5.30452299118042, 6.1361494064331055, 6.885622024536133, 7.53949499130249, 8.093769073486328, 8.551620483398438, 8.92063045501709, 9.209507942199707, 9.422353744506836, 9.535173416137695, 9.2625150680542, 8.173480033874512, 6.555201053619385, 5.160356044769287, 3.6473162174224854, 1.7884249687194824, 0.3061993420124054, 0.02045947127044201, 0.003275787690654397, 0.0014364334056153893, 0.0011037595104426146, 0.0010669095208868384, 0.0011754705337807536, 0.0014561830321326852, 0.0020661738235503435, 0.0035602941643446684, 0.008155648596584797, 0.023509694263339043, 0.06854546070098877, 0.18103069067001343, 0.4203503429889679, 0.8381097316741943, 1.440375804901123, 2.1951937675476074, 3.0528271198272705, 3.95931077003479, 4.865062236785889, 5.729351043701172, 6.5222392082214355, 7.22495698928833, 7.829020977020264, 8.334342002868652, 8.746667861938477, 9.074613571166992, 9.32568645477295, 9.495933532714844, 9.501469612121582, 8.822708129882812, 7.363541603088379, 5.838147163391113, 4.446770191192627, 2.7479753494262695, 0.8989087343215942, 0.07879916578531265, 0.006828206591308117, 0.001966425683349371, 0.0012069109361618757, 0.001064614742062986, 0.0011036444921046495, 0.001288687577471137, 0.0017012226162478328, 0.002631847746670246, 0.005183500237762928, 0.013643880374729633, 0.04050518199801445, 0.11303344368934631, 0.2808052599430084, 0.6053726077079773, 1.117563009262085, 1.8015910387039185, 2.6145570278167725, 3.503232717514038, 4.41506290435791, 5.3045244216918945, 6.136150360107422, 6.885622978210449, 7.539495468139648, 8.093769073486328, 8.551620483398438, 8.92063045501709, 9.209508895874023, 9.422353744506836, 9.535173416137695, 9.2625150680542, 8.173480033874512, 6.555201053619385, 5.160356521606445, 3.6473164558410645, 1.7884252071380615, 0.3061993718147278, 0.020459473133087158, 0.0032757879234850407, 0.0014364334056153893, 0.0011037595104426146, 0.0010669095208868384, 0.0011754705337807536, 0.0014561830321326852, 0.0020661738235503435, 0.0035602941643446684, 0.008155648596584797, 0.023509694263339043, 0.06854546070098877, 0.18103069067001343, 0.4203503429889679, 0.8381097316741943, 1.440375804901123, 2.1951937675476074, 3.0528271198272705, 3.95931077003479, 4.865062236785889, 5.729351043701172, 6.5222392082214355, 7.22495698928833, 7.829020977020264, 8.334342002868652, 8.746667861938477, 9.074613571166992, 9.32568645477295, 9.495933532714844, 9.501469612121582, 8.822708129882812, 7.363541603088379, 5.838147163391113, 4.446770191192627, 2.7479753494262695, 0.8989087343215942, 0.07879916578531265, 0.006828206591308117, 0.001966425683349371, 0.0012069109361618757, 0.001064614742062986, 0.0011036444921046495, 0.001288687577471137, 0.0017012226162478328, 0.002631847746670246, 0.005183500237762928, 0.013643880374729633, 0.04050518199801445, 0.11303344368934631, 0.2808052599430084, 0.6053726077079773, 1.117563009262085, 1.8015910387039185, 2.6145570278167725, 3.503232717514038, 4.41506290435791, 5.3045244216918945, 6.136150360107422, 6.885622978210449, 7.539495468139648, 8.093769073486328, 8.551620483398438], NH₄ = [0.10000000149011612, 0.043240949511528015, 0.0001223123981617391, 0.0001139868181780912, 0.00011856096534756944, 0.0001727785711409524, 0.0004039824998471886, 0.0007761779124848545, 0.0008759932243265212, 0.0008587001939304173, 0.0009447065531276166, 0.0019731158390641212, 0.009732905775308609, 0.009632930159568787, 0.007402525749057531, 0.007039217744022608, 0.00754015427082777, 0.008929948322474957, 0.011863507330417633, 0.018496019765734673, 0.03655117750167847, 0.08803946524858475, 0.1986335664987564, 0.3714835047721863, 0.6054473519325256, 0.9028093814849854, 1.230083703994751, 1.519329309463501, 1.719581961631775, 1.8148984909057617, 1.8119916915893555, 1.7280194759368896, 1.5844063758850098, 1.4034514427185059, 1.20576012134552, 1.008233904838562, 0.8228419423103333, 0.656194806098938, 0.5092864036560059, 0.37491533160209656, 0.22074398398399353, 0.00015944176993798465, 0.00014487933367490768, 0.00040241319220513105, 0.0009838792029768229, 0.0011213222751393914, 0.0011272826232016087, 0.0019841876346617937, 0.010275823064148426, 0.0116951335221529, 0.007991655729711056, 0.007115870248526335, 0.007228013128042221, 0.008122209459543228, 0.010139813646674156, 0.014464893378317356, 0.025152310729026794, 0.05617465823888779, 0.13523367047309875, 0.2776438891887665, 0.4804040789604187, 0.7470194101333618, 1.066762924194336, 1.383507251739502, 1.6322953701019287, 1.7802892923355103, 1.8248002529144287, 1.778852939605713, 1.6622618436813354, 1.4972418546676636, 1.3055355548858643, 1.1061041355133057, 0.9134812951087952, 0.7369732856750488, 0.5804613828659058, 0.441533625125885, 0.3047516345977783, 0.09762074798345566, 0.00013055410818196833, 0.00021895721147302538, 0.0007007925305515528, 0.0011153913801535964, 0.0011003970867022872, 0.0013022986240684986, 0.0044281017035245895, 0.014429138042032719, 0.009204195812344551, 0.007384237367659807, 0.007076871581375599, 0.007567816413938999, 0.008946670219302177, 0.011874630115926266, 0.01850549504160881, 0.03656217083334923, 0.08805251121520996, 0.19864404201507568, 0.37148839235305786, 0.6054481267929077, 0.902808427810669, 1.230082392692566, 1.5193281173706055, 1.719580888748169, 1.8148976564407349, 1.8119908571243286, 1.7280187606811523, 1.584405779838562, 1.4034510850906372, 1.2057597637176514, 1.008233666419983, 0.8228417634963989, 0.6561946868896484, 0.5092862844467163, 0.3749152719974518, 0.22074392437934875, 0.00015944175538606942, 0.0001448793482268229, 0.00040241325041279197, 0.0009838792029768229, 0.0011213222751393914, 0.0011272825067862868, 0.00198418740183115, 0.010275822132825851, 0.0116951335221529, 0.007991655729711056, 0.007115870248526335, 0.007228013128042221, 0.008122209459543228, 0.010139813646674156, 0.014464893378317356, 0.025152310729026794, 0.05617465823888779, 0.13523367047309875, 0.2776438891887665, 0.4804040789604187, 0.7470194101333618, 1.066762924194336, 1.383507251739502, 1.6322953701019287, 1.7802892923355103, 1.8248002529144287, 1.778852939605713, 1.6622618436813354, 1.4972418546676636, 1.3055355548858643, 1.1061041355133057, 0.9134812951087952, 0.7369732856750488, 0.5804613828659058, 0.441533625125885, 0.3047516345977783, 0.09762074798345566, 0.00013055410818196833, 0.00021895721147302538, 0.0007007925305515528, 0.0011153913801535964, 0.0011003970867022872, 0.0013022986240684986, 0.0044281017035245895, 0.014429138042032719, 0.009204195812344551, 0.007384237367659807, 0.007076871581375599, 0.007567816413938999, 0.008946670219302177, 0.011874630115926266, 0.01850549504160881, 0.03656217083334923, 0.08805251121520996, 0.19864404201507568, 0.37148839235305786, 0.6054481267929077, 0.902808427810669, 1.230082392692566, 1.5193281173706055, 1.719580888748169, 1.8148976564407349, 1.8119908571243286, 1.7280187606811523, 1.584405779838562, 1.4034510850906372, 1.2057597637176514, 1.008233666419983, 0.8228417634963989], P = [0.10000000149011612, 0.12984344363212585, 0.2040421962738037, 0.3034106194972992, 0.49908363819122314, 0.7335052490234375, 0.6938731670379639, 0.6320943832397461, 0.7436588406562805, 0.9307419061660767, 1.0742682218551636, 0.8847528696060181, 0.681115984916687, 0.7011634111404419, 0.7250727415084839, 0.7335047125816345, 0.7359057068824768, 0.7364087104797363, 0.7360563278198242, 0.7345408797264099, 0.7294687628746033, 0.7148680090904236, 0.6881376504898071, 0.6497069597244263, 0.589862585067749, 0.5041258335113525, 0.41236838698387146, 0.333420991897583, 0.2688766121864319, 0.21477864682674408, 0.1688195914030075, 0.13038446009159088, 0.09934178739786148, 0.07543353736400604, 0.05814423784613609, 0.04681691154837608, 0.040997464209795, 0.04116232320666313, 0.050592243671417236, 0.08218029886484146, 0.19273090362548828, 0.6263865232467651, 1.092787504196167, 0.9423794746398926, 0.7490721344947815, 0.8639734387397766, 1.0260825157165527, 0.9116142392158508, 0.6912965774536133, 0.7003745436668396, 0.7243351340293884, 0.7332423329353333, 0.7358993291854858, 0.7365544438362122, 0.7364272475242615, 0.7355548739433289, 0.732766330242157, 0.7237504720687866, 0.7028493285179138, 0.6707512736320496, 0.6232634782791138, 0.5494474768638611, 0.45737573504447937, 0.37090009450912476, 0.2996038794517517, 0.24072124063968658, 0.19084355235099792, 0.14867070317268372, 0.11394920200109482, 0.0865224152803421, 0.06600203365087509, 0.05177439749240875, 0.04322350397706032, 0.04023028537631035, 0.04429246857762337, 0.06195245683193207, 0.11941591650247574, 0.3545417785644531, 0.8746058940887451, 1.1194695234298706, 0.7914396524429321, 0.7849479913711548, 0.9554653167724609, 1.0271373987197876, 0.7615671157836914, 0.6843312382698059, 0.7148023247718811, 0.7300383448600769, 0.7349802851676941, 0.7363613843917847, 0.7365670204162598, 0.7361159324645996, 0.7345635294914246, 0.7294764518737793, 0.7148701548576355, 0.688139021396637, 0.6497082114219666, 0.5898632407188416, 0.5041259527206421, 0.4123683571815491, 0.33342093229293823, 0.2688765823841095, 0.21477863192558289, 0.1688195914030075, 0.13038447499275208, 0.09934179484844208, 0.07543353736400604, 0.05814424157142639, 0.04681691154837608, 0.040997464209795, 0.04116232693195343, 0.050592247396707535, 0.08218030631542206, 0.19273091852664948, 0.6263865232467651, 1.092787504196167, 0.942379355430603, 0.7490721344947815, 0.8639734387397766, 1.0260825157165527, 0.9116142392158508, 0.6912965774536133, 0.7003745436668396, 0.7243351340293884, 0.7332423329353333, 0.7358993291854858, 0.7365544438362122, 0.7364272475242615, 0.7355548739433289, 0.732766330242157, 0.7237504720687866, 0.7028493285179138, 0.6707512736320496, 0.6232634782791138, 0.5494474768638611, 0.45737573504447937, 0.37090009450912476, 0.2996038794517517, 0.24072124063968658, 0.19084355235099792, 0.14867070317268372, 0.11394920200109482, 0.0865224152803421, 0.06600203365087509, 0.05177439749240875, 0.04322350397706032, 0.04023028537631035, 0.04429246857762337, 0.06195245683193207, 0.11941591650247574, 0.3545417785644531, 0.8746058940887451, 1.1194695234298706, 0.7914396524429321, 0.7849479913711548, 0.9554653167724609, 1.0271373987197876, 0.7615671157836914, 0.6843312382698059, 0.7148023247718811, 0.7300383448600769, 0.7349802851676941, 0.7363613843917847, 0.7365670204162598, 0.7361159324645996, 0.7345635294914246, 0.7294764518737793, 0.7148701548576355, 0.688139021396637, 0.6497082114219666, 0.5898632407188416, 0.5041259527206421, 0.4123683571815491, 0.33342093229293823, 0.2688765823841095, 0.21477863192558289, 0.1688195914030075, 0.13038447499275208, 0.09934179484844208, 0.07543353736400604, 0.05814424157142639, 0.04681691154837608, 0.040997464209795], Z = [0.009999999776482582, 0.010352029465138912, 0.013263306580483913, 0.022448604926466942, 0.05595647543668747, 0.19885462522506714, 0.5087844729423523, 0.6716187000274658, 0.7686651945114136, 0.9158943891525269, 1.0677695274353027, 1.0844931602478027, 0.918433666229248, 0.8443203568458557, 0.8449364900588989, 0.8554612994194031, 0.8615809679031372, 0.8640700578689575, 0.8647333383560181, 0.8642415404319763, 0.8618202805519104, 0.854006290435791, 0.8362618088722229, 0.8071979284286499, 0.7619893550872803, 0.6911334991455078, 0.5971585512161255, 0.4966892600059509, 0.4025420844554901, 0.31879058480262756, 0.2460085153579712, 0.18444396555423737, 0.13428464531898499, 0.0951780304312706, 0.06603315472602844, 0.045207373797893524, 0.03087635524570942, 0.021376250311732292, 0.015431412495672703, 0.012384152971208096, 0.013257419690489769, 0.033119335770606995, 0.1937560737133026, 0.7000173330307007, 0.8611879348754883, 0.910651683807373, 1.029584288597107, 1.0738544464111328, 0.9310296773910522, 0.848743200302124, 0.8455564975738525, 0.855301558971405, 0.861463189125061, 0.8640960454940796, 0.8649256229400635, 0.8647755980491638, 0.8634608387947083, 0.858941376209259, 0.8465222716331482, 0.8232482075691223, 0.787247359752655, 0.7300936579704285, 0.6461127996444702, 0.5466909408569336, 0.44842293858528137, 0.3593061864376068, 0.28101643919944763, 0.21380497515201569, 0.1579435020685196, 0.11340132355690002, 0.07945415377616882, 0.05469467118382454, 0.03734097629785538, 0.025613293051719666, 0.018019165843725204, 0.013549317605793476, 0.012099391780793667, 0.017797915264964104, 0.07649126648902893, 0.4360930323600769, 0.8291451930999756, 0.8763893842697144, 0.9664925336837769, 1.0759012699127197, 1.0114351511001587, 0.8743828535079956, 0.842917263507843, 0.8505656123161316, 0.8589423298835754, 0.863096296787262, 0.8646619319915771, 0.8649584054946899, 0.864327073097229, 0.8618524074554443, 0.8540177345275879, 0.836266040802002, 0.8072001338005066, 0.7619907259941101, 0.6911341547966003, 0.5971588492393494, 0.4966893196105957, 0.4025420844554901, 0.3187905550003052, 0.24600850045681, 0.18444395065307617, 0.13428464531898499, 0.0951780304312706, 0.06603315472602844, 0.045207373797893524, 0.03087635524570942, 0.021376250311732292, 0.015431413426995277, 0.01238415390253067, 0.013257420621812344, 0.03311933949589729, 0.193756103515625, 0.7000173330307007, 0.8611878752708435, 0.9106516242027283, 1.029584288597107, 1.0738544464111328, 0.9310296773910522, 0.848743200302124, 0.8455564975738525, 0.855301558971405, 0.861463189125061, 0.8640960454940796, 0.8649256229400635, 0.8647755980491638, 0.8634608387947083, 0.858941376209259, 0.8465222716331482, 0.8232482075691223, 0.787247359752655, 0.7300936579704285, 0.6461127996444702, 0.5466909408569336, 0.44842293858528137, 0.3593061864376068, 0.28101643919944763, 0.21380497515201569, 0.1579435020685196, 0.11340132355690002, 0.07945415377616882, 0.05469467118382454, 0.03734097629785538, 0.025613293051719666, 0.018019165843725204, 0.013549317605793476, 0.012099391780793667, 0.017797915264964104, 0.07649126648902893, 0.4360930323600769, 0.8291451930999756, 0.8763893842697144, 0.9664925336837769, 1.0759012699127197, 1.0114351511001587, 0.8743828535079956, 0.842917263507843, 0.8505656123161316, 0.8589423298835754, 0.863096296787262, 0.8646619319915771, 0.8649584054946899, 0.864327073097229, 0.8618524074554443, 0.8540177345275879, 0.836266040802002, 0.8072001338005066, 0.7619907259941101, 0.6911341547966003, 0.5971588492393494, 0.4966893196105957, 0.4025420844554901, 0.3187905550003052, 0.24600850045681, 0.18444395065307617, 0.13428464531898499, 0.0951780304312706, 0.06603315472602844, 0.045207373797893524, 0.03087635524570942], sPOM = [0.0, 0.003560015931725502, 0.009505986236035824, 0.021664997562766075, 0.0552557148039341, 0.1581156849861145, 0.3376777470111847, 0.4167248010635376, 0.4905034899711609, 0.6372584104537964, 0.8074907064437866, 0.7920577526092529, 0.5860166549682617, 0.5341794490814209, 0.5428787469863892, 0.5532804131507874, 0.5583259463310242, 0.560195803642273, 0.5605785250663757, 0.5599197149276733, 0.5572159290313721, 0.5487362146377563, 0.5305110812187195, 0.5023579001426697, 0.4603208899497986, 0.39875328540802, 0.3260424733161926, 0.2578658163547516, 0.20072641968727112, 0.15397393703460693, 0.11582940071821213, 0.08512669801712036, 0.061069779098033905, 0.042857278138399124, 0.029560208320617676, 0.020177893340587616, 0.01375873014330864, 0.009509485214948654, 0.006881908979266882, 0.005789491813629866, 0.008371444419026375, 0.045999035239219666, 0.22992269694805145, 0.5465846061706543, 0.5722912549972534, 0.6174810528755188, 0.757769763469696, 0.7865648865699768, 0.6000315546989441, 0.5370551347732544, 0.5430235266685486, 0.5530583262443542, 0.558229386806488, 0.560244083404541, 0.5607921481132507, 0.5605002641677856, 0.5590180158615112, 0.5540472269058228, 0.5408978462219238, 0.5177380442619324, 0.48358720541000366, 0.4318816661834717, 0.3627024292945862, 0.2907193899154663, 0.22788219153881073, 0.1761741042137146, 0.13391301035881042, 0.099600650370121, 0.07231486588716507, 0.05128813907504082, 0.035656511783599854, 0.024441786110401154, 0.016652239486575127, 0.011403979733586311, 0.00801451038569212, 0.006112963892519474, 0.006207654718309641, 0.01625591143965721, 0.11090543121099472, 0.40198373794555664, 0.5829227566719055, 0.5785907506942749, 0.6831624507904053, 0.8082001805305481, 0.6909435391426086, 0.5519042611122131, 0.5377330780029297, 0.5486077666282654, 0.5561856627464294, 0.5594980716705322, 0.5606406927108765, 0.5607466101646423, 0.559983491897583, 0.5572396516799927, 0.5487443804740906, 0.5305140614509583, 0.5023595690727234, 0.4603218734264374, 0.39875373244285583, 0.3260425925254822, 0.2578658163547516, 0.20072640478610992, 0.15397392213344574, 0.11582939326763153, 0.08512669801712036, 0.061069779098033905, 0.042857278138399124, 0.029560208320617676, 0.020177895203232765, 0.01375873014330864, 0.009509486146271229, 0.006881909444928169, 0.005789491813629866, 0.00837144535034895, 0.04599904268980026, 0.22992269694805145, 0.5465845465660095, 0.5722911953926086, 0.617480993270874, 0.757769763469696, 0.7865648865699768, 0.6000315546989441, 0.5370551347732544, 0.5430235266685486, 0.5530583262443542, 0.558229386806488, 0.560244083404541, 0.5607921481132507, 0.5605002641677856, 0.5590180158615112, 0.5540472269058228, 0.5408978462219238, 0.5177380442619324, 0.48358720541000366, 0.4318816661834717, 0.3627024292945862, 0.2907193899154663, 0.22788219153881073, 0.1761741042137146, 0.13391301035881042, 0.099600650370121, 0.07231486588716507, 0.05128813907504082, 0.035656511783599854, 0.024441786110401154, 0.016652239486575127, 0.011403979733586311, 0.00801451038569212, 0.006112963892519474, 0.006207654718309641, 0.01625591143965721, 0.11090543121099472, 0.40198373794555664, 0.5829227566719055, 0.5785907506942749, 0.6831624507904053, 0.8082001805305481, 0.6909435391426086, 0.5519042611122131, 0.5377330780029297, 0.5486077666282654, 0.5561856627464294, 0.5594980716705322, 0.5606406927108765, 0.5607466101646423, 0.559983491897583, 0.5572396516799927, 0.5487443804740906, 0.5305140614509583, 0.5023595690727234, 0.4603218734264374, 0.39875373244285583, 0.3260425925254822, 0.2578658163547516, 0.20072640478610992, 0.15397392213344574, 0.11582939326763153, 0.08512669801712036, 0.061069779098033905, 0.042857278138399124, 0.029560208320617676, 0.020177895203232765, 0.01375873014330864], bPOM = [0.0, 0.0035621989518404007, 0.00952228531241417, 0.021755993366241455, 0.055876024067401886, 0.16628307104110718, 0.44872337579727173, 0.850450336933136, 1.2632336616516113, 1.7961138486862183, 2.5194859504699707, 3.1563198566436768, 3.1747355461120605, 2.906909704208374, 2.7138304710388184, 2.6202754974365234, 2.581129550933838, 2.565303087234497, 2.558298110961914, 2.553694009780884, 2.546907901763916, 2.529226064682007, 2.485206127166748, 2.4017462730407715, 2.2668633460998535, 2.060858726501465, 1.7794359922409058, 1.45612370967865, 1.1372742652893066, 0.8540084958076477, 0.619322657585144, 0.43491697311401367, 0.29655328392982483, 0.19704748690128326, 0.128180131316185, 0.08207203447818756, 0.05203966423869133, 0.03293241932988167, 0.021108902990818024, 0.01439428050071001, 0.013571138493716717, 0.04924270138144493, 0.23887471854686737, 0.7620310187339783, 1.430733561515808, 1.9383189678192139, 2.508795976638794, 3.0824551582336426, 3.1554105281829834, 2.9113919734954834, 2.7195627689361572, 2.623603343963623, 2.5827910900115967, 2.5662684440612793, 2.559204339981079, 2.5551931858062744, 2.5506434440612793, 2.5401811599731445, 2.511448860168457, 2.448976516723633, 2.3418476581573486, 2.173860549926758, 1.9282578229904175, 1.619956374168396, 1.2937089204788208, 0.990032434463501, 0.7302997708320618, 0.5210034251213074, 0.3603852391242981, 0.24245107173919678, 0.1592971235513687, 0.10273317992687225, 0.06540511548519135, 0.04138174280524254, 0.0262824185192585, 0.017181172966957092, 0.012895346619188786, 0.020310241729021072, 0.11425092816352844, 0.45103153586387634, 1.1159526109695435, 1.6929022073745728, 2.2051730155944824, 2.825829267501831, 3.192734718322754, 3.040071725845337, 2.801870346069336, 2.6620209217071533, 2.5986878871917725, 2.572704792022705, 2.562072277069092, 2.557058095932007, 2.553215265274048, 2.5467209815979004, 2.5291500091552734, 2.485173225402832, 2.4017319679260254, 2.2668581008911133, 2.060857057571411, 1.779435634613037, 1.45612370967865, 1.1372742652893066, 0.8540084362030029, 0.6193225979804993, 0.4349169433116913, 0.29655328392982483, 0.19704747200012207, 0.128180131316185, 0.08207202702760696, 0.05203966423869133, 0.03293241560459137, 0.021108902990818024, 0.01439428050071001, 0.013571139425039291, 0.04924270883202553, 0.23887471854686737, 0.7620310187339783, 1.430733561515808, 1.9383188486099243, 2.508795976638794, 3.0824551582336426, 3.1554105281829834, 2.9113919734954834, 2.7195627689361572, 2.623603343963623, 2.5827910900115967, 2.5662684440612793, 2.559204339981079, 2.5551931858062744, 2.5506434440612793, 2.5401811599731445, 2.511448860168457, 2.448976516723633, 2.3418476581573486, 2.173860549926758, 1.9282578229904175, 1.619956374168396, 1.2937089204788208, 0.990032434463501, 0.7302997708320618, 0.5210034251213074, 0.3603852391242981, 0.24245107173919678, 0.1592971235513687, 0.10273317992687225, 0.06540511548519135, 0.04138174280524254, 0.0262824185192585, 0.017181172966957092, 0.012895346619188786, 0.020310241729021072, 0.11425092816352844, 0.45103153586387634, 1.1159526109695435, 1.6929022073745728, 2.2051730155944824, 2.825829267501831, 3.192734718322754, 3.040071725845337, 2.801870346069336, 2.6620209217071533, 2.5986878871917725, 2.572704792022705, 2.562072277069092, 2.557058095932007, 2.553215265274048, 2.5467209815979004, 2.5291500091552734, 2.485173225402832, 2.4017319679260254, 2.2668581008911133, 2.060857057571411, 1.779435634613037, 1.45612370967865, 1.1372742652893066, 0.8540084362030029, 0.6193225979804993, 0.4349169433116913, 0.29655328392982483, 0.19704747200012207, 0.128180131316185, 0.08207202702760696, 0.05203966423869133], DOM = [0.0, 0.004266202915459871, 0.01208703313022852, 0.02720102295279503, 0.06265975534915924, 0.16641317307949066, 0.4481649696826935, 0.9190309643745422, 1.486770510673523, 2.17606782913208, 3.0583088397979736, 4.063421726226807, 4.825807571411133, 5.211373805999756, 5.374627590179443, 5.4393744468688965, 5.464417457580566, 5.473806381225586, 5.476770877838135, 5.47647762298584, 5.472854137420654, 5.461483001708984, 5.43074893951416, 5.364479064941406, 5.244716167449951, 5.0469512939453125, 4.74735164642334, 4.344983100891113, 3.8664441108703613, 3.3503191471099854, 2.8329670429229736, 2.342585325241089, 1.8981947898864746, 1.510410189628601, 1.1828272342681885, 0.9137228727340698, 0.6978655457496643, 0.5281946659088135, 0.39719080924987793, 0.2979830205440521, 0.22615158557891846, 0.1925777941942215, 0.28103384375572205, 0.7033841013908386, 1.4353750944137573, 2.231137275695801, 3.098215103149414, 4.047327518463135, 4.801496505737305, 5.197464466094971, 5.368093967437744, 5.436574935913086, 5.463322162628174, 5.473539352416992, 5.477054595947266, 5.477445125579834, 5.475399017333984, 5.468749523162842, 5.449538230895996, 5.403096675872803, 5.312619686126709, 5.157346725463867, 4.910678863525391, 4.557850360870361, 4.112893104553223, 3.610649585723877, 3.0898163318634033, 2.58300518989563, 2.1137945652008057, 1.6968560218811035, 1.339097499847412, 1.0412306785583496, 0.7995544075965881, 0.6077295541763306, 0.4583166241645813, 0.3439839780330658, 0.2586964964866638, 0.2020040601491928, 0.21090812981128693, 0.43766137957572937, 1.0516916513442993, 1.8292840719223022, 2.6506311893463135, 3.572720766067505, 4.470092296600342, 5.038052558898926, 5.301506519317627, 5.4101762771606445, 5.453062057495117, 5.469667911529541, 5.475822925567627, 5.477545261383057, 5.476773262023926, 5.472964763641357, 5.461521148681641, 5.430758476257324, 5.36447811126709, 5.244712829589844, 5.046947956085205, 4.747349262237549, 4.3449811935424805, 3.8664426803588867, 3.350318193435669, 2.8329663276672363, 2.3425848484039307, 1.898194432258606, 1.5104098320007324, 1.1828269958496094, 0.9137226939201355, 0.6978654265403748, 0.5281946063041687, 0.39719074964523315, 0.29798296093940735, 0.22615155577659607, 0.1925777792930603, 0.28103384375572205, 0.7033841013908386, 1.4353750944137573, 2.231137275695801, 3.098215103149414, 4.047327518463135, 4.801496505737305, 5.197464466094971, 5.368093967437744, 5.436574935913086, 5.463322162628174, 5.473539352416992, 5.477054595947266, 5.477445125579834, 5.475399017333984, 5.468749523162842, 5.449538230895996, 5.403096675872803, 5.312619686126709, 5.157346725463867, 4.910678863525391, 4.557850360870361, 4.112893104553223, 3.610649585723877, 3.0898163318634033, 2.58300518989563, 2.1137945652008057, 1.6968560218811035, 1.339097499847412, 1.0412306785583496, 0.7995544075965881, 0.6077295541763306, 0.4583166241645813, 0.3439839780330658, 0.2586964964866638, 0.2020040601491928, 0.21090812981128693, 0.43766137957572937, 1.0516916513442993, 1.8292840719223022, 2.6506311893463135, 3.572720766067505, 4.470092296600342, 5.038052558898926, 5.301506519317627, 5.4101762771606445, 5.453062057495117, 5.469667911529541, 5.475822925567627, 5.477545261383057, 5.476773262023926, 5.472964763641357, 5.461521148681641, 5.430758476257324, 5.36447811126709, 5.244712829589844, 5.046947956085205, 4.747349262237549, 4.3449811935424805, 3.8664426803588867, 3.350318193435669, 2.8329663276672363, 2.3425848484039307, 1.898194432258606, 1.5104098320007324, 1.1828269958496094, 0.9137226939201355, 0.6978654265403748])

And plot

using CairoMakie

fig = Figure(size = (1200, 1200), fontsize = 24)

axs = []
for (name, tracer) in pairs(timeseries)
    idx = (length(axs))
    push!(axs, Axis(fig[floor(Int, idx/2), Int(idx%2)], ylabel = "$name", xlabel = "Year", xticks=(0:10)))
    lines!(axs[end], times / year, tracer, linewidth = 3)
end

fig

This page was generated using Literate.jl.