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 (3.916 seconds)
[ Info: Executing initial time step...
[ Info:     ... initial time step complete (185.762 ms).
[ Info: Simulation is stopping after running for 24.512 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.040019989013672, 10.065675735473633, 10.083023071289062, 10.095608711242676, 10.105438232421875, 10.113603591918945, 10.120672225952148, 10.126932144165039, 10.132525444030762, 10.137530326843262, 10.14199447631836, 10.14596176147461, 10.1494722366333, 10.152569770812988, 10.155299186706543, 10.15770435333252, 10.159826278686523, 10.161704063415527, 10.163373947143555, 10.16486644744873, 10.166208267211914, 10.1674222946167, 10.168529510498047, 10.169546127319336, 10.170485496520996, 10.171359062194824, 10.172176361083984, 10.172945022583008, 10.17367172241211, 10.174360275268555, 10.175017356872559, 10.175643920898438, 10.176244735717773, 10.176820755004883, 10.177374839782715, 10.177909851074219, 10.178424835205078, 10.178922653198242, 10.179404258728027, 10.17987060546875, 10.180322647094727, 10.180760383605957, 10.181185722351074, 10.181597709655762, 10.181998252868652, 10.182388305664062, 10.18276596069336, 10.183135032653809, 10.183492660522461, 10.183841705322266, 10.184181213378906, 10.1845121383667, 10.184834480285645, 10.185149192810059, 10.185455322265625, 10.185754776000977, 10.186046600341797, 10.186331748962402, 10.186609268188477, 10.186881065368652, 10.18714714050293, 10.187406539916992, 10.187660217285156, 10.187908172607422, 10.188150405883789, 10.188386917114258, 10.188619613647461, 10.188846588134766, 10.189068794250488, 10.189286231994629, 10.189499855041504, 10.189708709716797, 10.189913749694824, 10.19011402130127, 10.19031047821045, 10.190503120422363, 10.190691947937012, 10.190877914428711, 10.191059112548828, 10.191237449645996, 10.191412925720215, 10.191584587097168, 10.191753387451172, 10.19191837310791, 10.192081451416016, 10.192240715026855, 10.192398071289062, 10.19255256652832, 10.192703247070312, 10.192852020263672, 10.192998886108398, 10.193142890930176, 10.193284034729004, 10.1934232711792, 10.193559646606445, 10.193694114685059, 10.193826675415039, 10.19395637512207, 10.194085121154785, 10.19421100616455, 10.194334983825684, 10.194457054138184, 10.19457721710205, 10.194695472717285, 10.194812774658203, 10.194927215576172, 10.195040702819824, 10.195151329040527, 10.19526195526123, 10.195369720458984, 10.195476531982422, 10.195581436157227, 10.195684432983398, 10.195786476135254, 10.195887565612793, 10.1959867477417, 10.196084976196289, 10.196181297302246, 10.19627571105957, 10.196370124816895, 10.196462631225586, 10.196554183959961, 10.196643829345703, 10.196732521057129, 10.196820259094238, 10.196907043457031, 10.196991920471191, 10.197076797485352, 10.197159767150879, 10.19724178314209, 10.197322845458984, 10.197402954101562, 10.197482109069824, 10.19756031036377, 10.197637557983398, 10.197713851928711, 10.197789192199707, 10.197863578796387, 10.19793701171875, 10.198009490966797, 10.198081016540527, 10.198151588439941, 10.198221206665039, 10.198290824890137, 10.198359489440918, 10.198426246643066, 10.198493003845215, 10.198559761047363, 10.198624610900879, 10.198689460754395, 10.198753356933594, 10.198816299438477, 10.198878288269043, 10.19894027709961, 10.19900131225586, 10.199061393737793, 10.199121475219727, 10.199179649353027, 10.199238777160645, 10.199295997619629, 10.199353218078613, 10.199409484863281, 10.19946575164795, 10.1995210647583, 10.199576377868652, 10.199629783630371, 10.199684143066406, 10.199736595153809, 10.199789047241211, 10.199841499328613, 10.1998929977417, 10.199943542480469, 10.199994087219238, 10.200044631958008, 10.200093269348145, 10.200142860412598, 10.200191497802734, 10.200239181518555, 10.200286865234375, 10.200333595275879, 10.200380325317383, 10.20042610168457], NH₄ = [0.10000000149011612, 0.0629749745130539, 0.04138703644275665, 0.028993042185902596, 0.021883627399802208, 0.017698973417282104, 0.015069189481437206, 0.013233933597803116, 0.011794044636189938, 0.010555262677371502, 0.009433391503989697, 0.00839900877326727, 0.00744673004373908, 0.006579192355275154, 0.005799497477710247, 0.005108272191137075, 0.004503076430410147, 0.003978857304900885, 0.003528757719323039, 0.003144931746646762, 0.002819227520376444, 0.0025436843279749155, 0.002310865791514516, 0.0021140470635145903, 0.0019472995772957802, 0.0018055003602057695, 0.0016842937329784036, 0.0015800247201696038, 0.001489659771323204, 0.0014107036404311657, 0.0013411182444542646, 0.0012792490888386965, 0.001223759027197957, 0.0011735704028978944, 0.0011278170859441161, 0.0010858031455427408, 0.0010469689732417464, 0.0010108641581609845, 0.000977124902419746, 0.0009454560349695385, 0.0009156167507171631, 0.0008874092600308359, 0.0008606698829680681, 0.0008352617733180523, 0.000811069505289197, 0.0007879948243498802, 0.000765953038353473, 0.0007448705146089196, 0.0007246825261972845, 0.0007053317385725677, 0.0006867667543701828, 0.0006689414149150252, 0.0006518136360682547, 0.0006353450589813292, 0.0006195004098117352, 0.0006042472086846828, 0.00058955536223948, 0.0005753968143835664, 0.0005617454298771918, 0.0005485768197104335, 0.0005358681082725525, 0.0005235979333519936, 0.0005117460386827588, 0.0005002935067750514, 0.0004892224678769708, 0.0004785160126630217, 0.0004681582504417747, 0.0004581341054290533, 0.00044842937495559454, 0.00043903052574023604, 0.0004299248394090682, 0.0004211002087686211, 0.00041254510870203376, 0.0004042486834805459, 0.00039620057214051485, 0.0003883909375872463, 0.0003808104374911636, 0.00037345013697631657, 0.0003663016250357032, 0.000359356781700626, 0.00035260795266367495, 0.00034604783286340535, 0.00033966940827667713, 0.00033346604323014617, 0.00032743136398494244, 0.0003215592878405005, 0.00031584405223838985, 0.00031028006924316287, 0.0003048620419576764, 0.00029958493541926146, 0.00029444388928823173, 0.0002894342178478837, 0.0002845515264198184, 0.00027979150763712823, 0.00027515014517121017, 0.00027062345179729164, 0.00026620770222507417, 0.0002618993166834116, 0.0002576948027126491, 0.0002535908715799451, 0.0002495843218639493, 0.00024567206855863333, 0.00024185123038478196, 0.0002381189406150952, 0.00023447247804142535, 0.00023090923787094653, 0.00022742671717423946, 0.00022402245667763054, 0.0002206941571785137, 0.00021743956313002855, 0.00021425650629680604, 0.00021114292030688375, 0.00020809676789212972, 0.00020511614275164902, 0.0002021991676883772, 0.0001993440237129107, 0.0001965489936992526, 0.00019381240417715162, 0.00019113261078018695, 0.0001885080710053444, 0.0001859372714534402, 0.0001834187569329515, 0.00018095110135618597, 0.00017853296594694257, 0.00017616301192902029, 0.0001738399878377095, 0.00017156264220830053, 0.00016932979633565992, 0.00016714030061848462, 0.0001649930200073868, 0.0001628868922125548, 0.0001608208694960922, 0.00015879393322393298, 0.00015680510841775686, 0.00015485343465115875, 0.0001529380097053945, 0.00015105793136171997, 0.00014921232650522143, 0.00014740036567673087, 0.0001456212339689955, 0.00014387414557859302, 0.00014215832925401628, 0.00014047304284758866, 0.00013881757331546396, 0.00013719122216571122, 0.0001355932909063995, 0.0001340231392532587, 0.00013248012692201883, 0.0001309636136284098, 0.0001294730172958225, 0.00012800772674381733, 0.00012656720355153084, 0.0001251508656423539, 0.00012375818914733827, 0.00012238863564562052, 0.00012104171037208289, 0.00011971691128564999, 0.00011841375089716166, 0.00011713176354533061, 0.00011587049812078476, 0.00011462951079010963, 0.00011340835771989077, 0.00011220663145650178, 0.00011102390999440104, 0.00010985980043187737, 0.00010871391714317724, 0.00010758587450254709, 0.0001064753087121062, 0.00010538185597397387, 0.00010430517431814224, 0.00010324491449864581, 0.00010220074909739196, 0.00010117235069628805, 0.00010015940642915666, 9.916159615386277e-5, 9.817864338401705e-5, 9.721022797748446e-5, 9.625608799979091e-5, 9.53159251366742e-5, 9.438947745366022e-5, 9.347648301627487e-5, 9.257667261408642e-5, 9.168979886453599e-5], P = [0.10000000149011612, 0.08886265754699707, 0.07981985807418823, 0.07253268361091614, 0.06666655838489532, 0.06192740052938461, 0.05807190760970116, 0.05490536615252495, 0.052274834364652634, 0.05006144195795059, 0.048173341900110245, 0.04653981700539589, 0.04510650411248207, 0.04383166879415512, 0.042683277279138565, 0.04163672775030136, 0.04067312553524971, 0.03977792337536812, 0.03893988952040672, 0.038150325417518616, 0.037402454763650894, 0.036690954118967056, 0.03601161763072014, 0.03536108508706093, 0.03473663330078125, 0.034136030822992325, 0.03355742245912552, 0.03299923986196518, 0.03246014565229416, 0.03193897008895874, 0.03143468499183655, 0.030946381390094757, 0.03047322854399681, 0.030014479532837868, 0.029569443315267563, 0.029137490317225456, 0.028718028217554092, 0.028310509398579597, 0.027914419770240784, 0.0275292806327343, 0.02715463750064373, 0.026790061965584755, 0.02643515355885029, 0.026089530438184738, 0.025752831250429153, 0.025424713268876076, 0.02510485239326954, 0.02479294128715992, 0.02448868751525879, 0.024191809818148613, 0.0239020437002182, 0.023619139567017555, 0.023342851549386978, 0.023072954267263412, 0.0228092260658741, 0.022551460191607475, 0.022299453616142273, 0.02205301821231842, 0.021811969578266144, 0.021576132625341415, 0.021345341578125954, 0.021119436249136925, 0.020898262038826942, 0.020681671798229218, 0.020469525828957558, 0.02026168815791607, 0.020058028399944305, 0.01985842175781727, 0.019662749022245407, 0.01947089470922947, 0.0192827470600605, 0.019098201766610146, 0.018917156383395195, 0.01873951032757759, 0.018565168604254723, 0.018394041806459427, 0.01822604238986969, 0.018061082810163498, 0.017899081110954285, 0.01773996092379093, 0.01758364588022232, 0.017430059611797333, 0.017279135063290596, 0.01713079959154129, 0.01698499172925949, 0.016841642558574677, 0.01670069433748722, 0.0165620855987072, 0.016425758600234985, 0.016291657462716103, 0.016159728169441223, 0.01602991856634617, 0.01590217836201191, 0.015776457265019417, 0.01565270870923996, 0.015530885197222233, 0.015410944819450378, 0.01529284194111824, 0.01517653651535511, 0.015061985701322556, 0.014949151314795017, 0.014837995171546936, 0.014728480018675327, 0.014620569534599781, 0.014514228329062462, 0.014409423805773258, 0.01430612150579691, 0.014204289764165878, 0.014103896915912628, 0.014004914090037346, 0.013907310552895069, 0.01381105836480856, 0.013716128654778004, 0.013622495345771313, 0.013530132360756397, 0.013439012691378593, 0.013349112123250961, 0.01326040644198656, 0.013172872364521027, 0.013086485676467419, 0.01300122495740652, 0.012917067855596542, 0.012833993881940842, 0.012751980684697628, 0.012671009637415409, 0.012591060250997543, 0.012512113898992538, 0.012434150092303753, 0.012357152998447418, 0.012281103059649467, 0.012205984443426132, 0.012131777592003345, 0.012058468535542488, 0.011986040510237217, 0.011914476752281189, 0.01184376236051321, 0.011773882433772087, 0.0117048230022192, 0.011636568233370781, 0.011569105088710785, 0.011502419598400593, 0.01143649872392416, 0.011371329426765442, 0.01130689773708582, 0.011243192479014397, 0.011180201545357704, 0.011117911897599697, 0.01105631235986948, 0.010995391756296158, 0.010935138911008835, 0.01087554357945919, 0.010816593654453754, 0.010758278891444206, 0.0107005899772048, 0.010643516667187214, 0.010587048716843128, 0.010531176812946796, 0.010475891642272472, 0.01042118389159441, 0.010367044247686863, 0.01031346432864666, 0.010260435752570629, 0.010207949206233025, 0.01015599723905325, 0.010104571469128132, 0.0100536635145545, 0.010003265924751759, 0.009953371249139309, 0.009903972037136555, 0.009855060838162899, 0.009806630201637745, 0.009758672676980495, 0.009711182676255703, 0.00966415274888277, 0.009617575444281101, 0.009571445174515247, 0.009525755420327187, 0.009480499662458897, 0.00943567231297493, 0.009391266852617264, 0.009347276762127876, 0.009303697384893894, 0.009260522201657295], Z = [0.009999999776482582, 0.009559744037687778, 0.008661212399601936, 0.007543097250163555, 0.006377417128533125, 0.005271041300147772, 0.00428025471046567, 0.0034270575270056725, 0.0027125615160912275, 0.002126515144482255, 0.0016534717287868261, 0.0012764973798766732, 0.0009792278287932277, 0.0007468870608136058, 0.0005666868528351188, 0.0004278743581380695, 0.00032159656984731555, 0.00024068157654255629, 0.00017939390090759844, 0.00013319564459379762, 9.852897346718237e-5, 7.262658618856221e-5, 5.335130845196545e-5, 3.906314304913394e-5, 2.851102726708632e-5, 2.074576150334906e-5, 1.5050852198328357e-5, 1.0888043107115664e-5, 7.854798241169192e-6, 5.651374976878287e-6, 4.05548462367733e-6, 2.902918367908569e-6, 2.072824827337172e-6, 1.4765802234251169e-6, 1.0494159141671844e-6, 7.441537377417262e-7, 5.265373488327896e-7, 3.717690333360224e-7, 2.619507029066881e-7, 1.8420109881844837e-7, 1.2927486636726826e-7, 9.055393235257725e-8, 6.331301705131409e-8, 4.4186712955252005e-8, 3.07838519120196e-8, 2.1409459094456906e-8, 1.4864771635814122e-8, 1.0303815578538433e-8, 7.130865231630423e-9, 4.927276808786019e-9, 3.3994460579833685e-9, 2.341860261267925e-9, 1.610942934604509e-9, 1.1065727223424915e-9, 7.590564288229018e-10, 5.199670694011616e-10, 3.5571168233161643e-10, 2.4302637680051475e-10, 1.6582639428275314e-10, 1.130083804312676e-10, 7.691935283160589e-11, 5.2292517538354844e-11, 3.550856136902425e-11, 2.4083880723058115e-11, 1.6316604317667682e-11, 1.1042125547255921e-11, 7.464566291470742e-12, 5.040743430301253e-12, 3.400412763865379e-12, 2.2915235247528143e-12, 1.5426983692237717e-12, 1.0375520971897534e-12, 6.97138147990467e-13, 4.679673532362316e-13, 3.1383947632530884e-13, 2.1028201424057313e-13, 1.4076878497495254e-13, 9.415147968986487e-14, 6.291746502051429e-14, 4.200926648103946e-14, 2.8025618189543057e-14, 1.8681326898286645e-14, 1.244253198330646e-14, 8.280644914334875e-15, 5.506551872737787e-15, 3.658988361593656e-15, 2.4294814986549392e-15, 1.611917357031208e-15, 1.0686947991662911e-15, 7.080298113517453e-16, 4.687499388346611e-16, 3.101184819756026e-16, 2.050290894432317e-16, 1.3545948092348804e-16, 8.943631631134005e-17, 5.90110582820881e-17, 3.89109896048007e-17, 2.5640981023865704e-17, 1.6885919440869698e-17, 1.1113387423789803e-17, 7.309772739254567e-18, 4.805081838562881e-18, 3.1567557232060124e-18, 2.0726604638146804e-18, 1.3600845877575325e-18, 8.919854572580352e-19, 5.846652156253168e-19, 3.8301667686078824e-19, 2.5077973257833404e-19, 1.6410988057004454e-19, 1.0733657187201476e-19, 7.016726920808129e-20, 4.584567119483608e-20, 2.9939324057565204e-20, 1.9541975566090488e-20, 1.2749133640681646e-20, 8.313451922874194e-21, 5.4184292349669184e-21, 3.529875730106055e-21, 2.298487131608203e-21, 1.4959738082864275e-21, 9.732126248178199e-22, 6.328427469332375e-22, 4.113302765810408e-22, 2.672358581607652e-22, 1.7354431257350573e-22, 1.1265226551724067e-22, 7.309465872009026e-23, 4.740779832327626e-23, 3.073509687368534e-23, 1.9917835399954975e-23, 1.290252140223439e-23, 8.354761160601508e-24, 5.407823732159142e-24, 3.498985050929198e-24, 2.26305261046183e-24, 1.4631275366146407e-24, 9.455981007575642e-25, 6.108995056371012e-25, 3.945240805391943e-25, 2.5469456543614716e-25, 1.6436524629695834e-25, 1.0603427018911983e-25, 6.838017657183583e-26, 4.4082222229602653e-26, 2.84084650817949e-26, 1.8301415845244012e-26, 1.178625676985443e-26, 7.58792362482865e-27, 4.883457300466273e-27, 3.141888540749461e-27, 2.0207592681749634e-27, 1.2992727850124614e-27, 8.351211647062099e-28, 5.366157937252907e-28, 3.4470188218495806e-28, 2.213561081567907e-28, 1.4210466872423928e-28, 9.12001495487729e-29, 5.851327856969488e-29, 3.753067205791252e-29, 2.4065368846757332e-29, 1.5426742488892062e-29, 9.88627687776024e-30, 6.333872764264525e-30, 4.056814946043731e-30, 2.5976550004584122e-30, 1.6628740569783017e-30, 1.0641921509971649e-30, 6.80870684239218e-31, 4.355061748892554e-31, 2.784903398242139e-31, 1.7803821146593815e-31, 1.1379014537358153e-31, 7.270855150054969e-32, 4.644690567722589e-32, 2.966330606941272e-32, 1.8939775100063906e-32, 1.2089924868060672e-32, 7.715549658822031e-33, 4.922725691481975e-33, 3.140080794278299e-33, 2.0025040341226998e-33], sPOM = [0.0, 0.002583166118711233, 0.0036341671366244555, 0.003868874628096819, 0.0037018738221377134, 0.0033643203787505627, 0.002979036420583725, 0.002606767462566495, 0.002273865742608905, 0.001988436095416546, 0.001749521354213357, 0.0015521722380071878, 0.001390137942507863, 0.0012572066625580192, 0.0011478046653792262, 0.0010571989696472883, 0.000981504563242197, 0.0009176013991236687, 0.0008630204829387367, 0.0008158276905305684, 0.0007745178882032633, 0.0007379245944321156, 0.0007051456486806273, 0.0006754834321327507, 0.0006483977194875479, 0.0006234689499251544, 0.0006003694725222886, 0.0005788420676253736, 0.0005586827755905688, 0.0005397287313826382, 0.0005218481528572738, 0.0005049330648034811, 0.000488894060254097, 0.00047365593491122127, 0.0004591546021401882, 0.00044533476466313004, 0.0004321481683291495, 0.00041955214692279696, 0.0004075086326338351, 0.0003959833993576467, 0.0003849453933071345, 0.0003743662964552641, 0.0003642201190814376, 0.0003544829087331891, 0.0003451324882917106, 0.0003361483686603606, 0.000327511370414868, 0.00031920368201099336, 0.0003112086560577154, 0.00030351069290190935, 0.00029609521152451634, 0.00028894850402139127, 0.0002820576773956418, 0.0002754106535576284, 0.00026899605290964246, 0.00026280313613824546, 0.00025682177511043847, 0.0002510424528736621, 0.0002454561472404748, 0.00024005438899621367, 0.0002348290872760117, 0.0002297726459801197, 0.00022487787646241486, 0.00022013795387465507, 0.00021554638806264848, 0.00021109706722199917, 0.0002067841705866158, 0.00020260218298062682, 0.00019854583661071956, 0.00019461015472188592, 0.00019079039338976145, 0.00018708205607254058, 0.00018348083540331572, 0.00017998265684582293, 0.00017658360593486577, 0.00017327998648397624, 0.00017006824782583863, 0.00016694501391611993, 0.0001639070687815547, 0.0001609513274161145, 0.00015807485033292323, 0.0001552748290123418, 0.00015254857135005295, 0.0001498935162089765, 0.00014730720431543887, 0.00014478727825917304, 0.00014233148249331862, 0.0001399376633344218, 0.0001376037544105202, 0.00013532776210922748, 0.00013310778012964875, 0.0001309420185862109, 0.00012882868759334087, 0.0001267661282327026, 0.00012475272524170578, 0.00012278692156542093, 0.00012086722563253716, 0.0001189922186313197, 0.00011716051085386425, 0.00011537079262780026, 0.00011362178338458762, 0.00011191225348738953, 0.00011024102423107252, 0.00010860696784220636, 0.00010700897837523371, 0.00010544600809225813, 0.00010391704563517123, 0.00010242110147373751, 0.00010095723700942472, 9.952453547157347e-5, 9.812213102122769e-5, 9.674916509538889e-5, 9.540481551084667e-5, 9.408830374013633e-5, 9.279885853175074e-5, 9.15357522899285e-5, 9.029825741890818e-5, 8.908569725463167e-5, 8.789740240899846e-5, 8.673272532178089e-5, 8.559105481253937e-5, 8.44717797008343e-5, 8.337432518601418e-5, 8.229812374338508e-5, 8.124262240016833e-5, 8.020730456337333e-5, 7.919166091596708e-5, 7.819518214091659e-5, 7.721740257693455e-5, 7.625784928677604e-5, 7.531607116106898e-5, 7.439163891831413e-5, 7.348412327701226e-5, 7.259312405949458e-5, 7.17182265361771e-5, 7.085905235726386e-5, 7.001523044891655e-5, 6.918639701325446e-5, 6.83721955283545e-5, 6.757227674825117e-5, 6.678632780676708e-5, 6.601400673389435e-5, 6.525500793941319e-5, 6.45090258331038e-5, 6.377576937666163e-5, 6.305494025582448e-5, 6.234626926016062e-5, 6.164947262732312e-5, 6.096430297475308e-5, 6.0290487454039976e-5, 5.9627785958582535e-5, 5.8975954743800685e-5, 5.833475006511435e-5, 5.7703950005816296e-5, 5.7083329011220485e-5, 5.647266880259849e-5, 5.587175110122189e-5, 5.5280379456235096e-5, 5.4698350140824914e-5, 5.4125463066156954e-5, 5.3561529057333246e-5, 5.3006366215413436e-5, 5.245978900347836e-5, 5.192162643652409e-5, 5.139170025358908e-5, 5.086985038360581e-5, 5.035590584157035e-5, 4.98497138323728e-5, 4.935111792292446e-5, 4.885996168013662e-5, 4.8376106860814616e-5, 4.7899404307827353e-5, 4.7429712140001357e-5, 4.696689575212076e-5, 4.651082053896971e-5, 4.606135917128995e-5, 4.561838068184443e-5, 4.51817613793537e-5, 4.4751384848495945e-5, 4.432712739799172e-5, 4.390887261251919e-5, 4.349651135271415e-5, 4.308993447921239e-5], bPOM = [0.0, 0.002584952162578702, 0.0036416808143258095, 0.0038819925393909216, 0.0037178059574216604, 0.003380199894309044, 0.002993041882291436, 0.002618159167468548, 0.0022826273925602436, 0.001994913909584284, 0.0017541770357638597, 0.0015554510755464435, 0.001392413629218936, 0.001258770003914833, 0.001148870913311839, 0.0010579228401184082, 0.0009819945553317666, 0.0009179325425066054, 0.0008632441167719662, 0.0008159787976182997, 0.0007746200426481664, 0.0007379937451332808, 0.0007051924476400018, 0.0006755152135156095, 0.0006484193145297468, 0.0006234836182557046, 0.0006003794842399657, 0.0005788488197140396, 0.0005586874322034419, 0.0005397319328039885, 0.0005218503065407276, 0.0005049345199950039, 0.0004888950497843325, 0.00047365660429932177, 0.00045915503869764507, 0.00044533508480526507, 0.0004321483720559627, 0.00041955229244194925, 0.0004075087490491569, 0.0003959834575653076, 0.0003849454515147954, 0.00037436632555909455, 0.00036422014818526804, 0.0003544829087331891, 0.00034513251739554107, 0.0003361483686603606, 0.000327511370414868, 0.00031920368201099336, 0.0003112086560577154, 0.00030351069290190935, 0.00029609521152451634, 0.00028894850402139127, 0.0002820576773956418, 0.0002754106535576284, 0.00026899605290964246, 0.00026280313613824546, 0.00025682177511043847, 0.0002510424528736621, 0.0002454561472404748, 0.00024005438899621367, 0.0002348290872760117, 0.0002297726459801197, 0.00022487787646241486, 0.00022013795387465507, 0.00021554638806264848, 0.00021109706722199917, 0.0002067841705866158, 0.00020260218298062682, 0.00019854583661071956, 0.00019461015472188592, 0.00019079039338976145, 0.00018708205607254058, 0.00018348083540331572, 0.00017998265684582293, 0.00017658360593486577, 0.00017327998648397624, 0.00017006824782583863, 0.00016694501391611993, 0.0001639070687815547, 0.0001609513274161145, 0.00015807485033292323, 0.0001552748290123418, 0.00015254857135005295, 0.0001498935162089765, 0.00014730720431543887, 0.00014478727825917304, 0.00014233148249331862, 0.0001399376633344218, 0.0001376037544105202, 0.00013532776210922748, 0.00013310778012964875, 0.0001309420185862109, 0.00012882868759334087, 0.0001267661282327026, 0.00012475272524170578, 0.00012278692156542093, 0.00012086722563253716, 0.0001189922186313197, 0.00011716051085386425, 0.00011537079262780026, 0.00011362178338458762, 0.00011191225348738953, 0.00011024102423107252, 0.00010860696784220636, 0.00010700897837523371, 0.00010544600809225813, 0.00010391704563517123, 0.00010242110147373751, 0.00010095723700942472, 9.952453547157347e-5, 9.812213102122769e-5, 9.674916509538889e-5, 9.540481551084667e-5, 9.408830374013633e-5, 9.279885853175074e-5, 9.15357522899285e-5, 9.029825741890818e-5, 8.908569725463167e-5, 8.789740240899846e-5, 8.673272532178089e-5, 8.559105481253937e-5, 8.44717797008343e-5, 8.337432518601418e-5, 8.229812374338508e-5, 8.124262240016833e-5, 8.020730456337333e-5, 7.919166091596708e-5, 7.819518214091659e-5, 7.721740257693455e-5, 7.625784928677604e-5, 7.531607116106898e-5, 7.439163891831413e-5, 7.348412327701226e-5, 7.259312405949458e-5, 7.17182265361771e-5, 7.085905235726386e-5, 7.001523044891655e-5, 6.918639701325446e-5, 6.83721955283545e-5, 6.757227674825117e-5, 6.678632780676708e-5, 6.601400673389435e-5, 6.525500793941319e-5, 6.45090258331038e-5, 6.377576937666163e-5, 6.305494025582448e-5, 6.234626926016062e-5, 6.164947262732312e-5, 6.096430297475308e-5, 6.0290487454039976e-5, 5.9627785958582535e-5, 5.8975954743800685e-5, 5.833475006511435e-5, 5.7703950005816296e-5, 5.7083329011220485e-5, 5.647266880259849e-5, 5.587175110122189e-5, 5.5280379456235096e-5, 5.4698350140824914e-5, 5.4125463066156954e-5, 5.3561529057333246e-5, 5.3006366215413436e-5, 5.245978900347836e-5, 5.192162643652409e-5, 5.139170025358908e-5, 5.086985038360581e-5, 5.035590584157035e-5, 4.98497138323728e-5, 4.935111792292446e-5, 4.885996168013662e-5, 4.8376106860814616e-5, 4.7899404307827353e-5, 4.7429712140001357e-5, 4.696689575212076e-5, 4.651082053896971e-5, 4.606135917128995e-5, 4.561838068184443e-5, 4.51817613793537e-5, 4.4751384848495945e-5, 4.432712739799172e-5, 4.390887261251919e-5, 4.349651135271415e-5, 4.308993447921239e-5], DOM = [0.0, 0.0034145123790949583, 0.007179975975304842, 0.010157641023397446, 0.012044407427310944, 0.012919594533741474, 0.013003219850361347, 0.012536533176898956, 0.011730199679732323, 0.010747780092060566, 0.009706081822514534, 0.008682388812303543, 0.007723473943769932, 0.006854043807834387, 0.006083735730499029, 0.005412491504102945, 0.004834474064409733, 0.004340791143476963, 0.003921308554708958, 0.0035657924599945545, 0.0032645834144204855, 0.003008951898664236, 0.002791242441162467, 0.002604888752102852, 0.0024443536531180143, 0.002305028960108757, 0.002183118835091591, 0.0020755219738930464, 0.001979721011593938, 0.0018936814740300179, 0.0018157653976231813, 0.0017446561250835657, 0.0016792957903817296, 0.001618834096007049, 0.0015625861706212163, 0.0015099983429536223, 0.001460620784200728, 0.001414086320437491, 0.0013700926210731268, 0.001328388461843133, 0.0012887639459222555, 0.0012510411906987429, 0.0012150682741776109, 0.0011807138798758388, 0.0011478636879473925, 0.0011164168827235699, 0.0010862837079912424, 0.0010573839535936713, 0.0010296450927853584, 0.0010030008852481842, 0.000977391260676086, 0.0009527604561299086, 0.0009290574234910309, 0.0009062346071004868, 0.0008842477691359818, 0.0008630557567812502, 0.0008426200947724283, 0.0008229046943597496, 0.0008038756786845624, 0.0007855013245716691, 0.0007677516550756991, 0.0007505984976887703, 0.0007340153097175062, 0.0007179769454523921, 0.0007024597725830972, 0.000687441264744848, 0.0006729002925567329, 0.0006588167743757367, 0.0006451716180890799, 0.0006319468375295401, 0.0006191252032294869, 0.0006066905916668475, 0.0005946274031884968, 0.0005829210858792067, 0.0005715576116926968, 0.0005605237674899399, 0.000549806805793196, 0.0005393947358243167, 0.0005292760906741023, 0.0005194399273023009, 0.0005098758265376091, 0.0005005737766623497, 0.0004915242898277938, 0.0004827183729503304, 0.0004741472948808223, 0.000465802731923759, 0.0004576768260449171, 0.000449761952040717, 0.0004420508921612054, 0.0004345366614870727, 0.00042721262434497476, 0.00042007240699604154, 0.00041310989763587713, 0.0004063192754983902, 0.00039969486533664167, 0.0003932313120458275, 0.00038692346424795687, 0.0003807663160841912, 0.0003747551527339965, 0.0003688853175844997, 0.00036315247416496277, 0.0003575523442123085, 0.0003520808822941035, 0.00034673415939323604, 0.0003415084502194077, 0.0003364000585861504, 0.00033140555024147034, 0.0003265215491410345, 0.0003217447956558317, 0.00031707220477983356, 0.0003125007206108421, 0.0003080274909734726, 0.0003036496927961707, 0.000299364619422704, 0.0002951697097159922, 0.0002910624025389552, 0.00028704031137749553, 0.00028310107882134616, 0.0002792424347717315, 0.00027546219644136727, 0.0002717582683544606, 0.00026812858413904905, 0.00026457119383849204, 0.00026108414749614894, 0.0002576656697783619, 0.0002543139271438122, 0.0002510272024665028, 0.00024780380772426724, 0.0002446421713102609, 0.00024154066340997815, 0.00023849781427998096, 0.00023551213962491602, 0.00023258219880517572, 0.00022970663849264383, 0.0002268841053592041, 0.00022411331883631647, 0.0002213929983554408, 0.00021872192155569792, 0.0002160989388357848, 0.00021352287149056792, 0.00021099262812640518, 0.00020850711734965444, 0.00020606527687050402, 0.0002036661171587184, 0.00020130863413214684, 0.000198991852812469, 0.00019671487098094076, 0.00019447675731498748, 0.0001922766532516107, 0.0001901136856758967, 0.0001879870251286775, 0.00018589585670270026, 0.00018383942369837314, 0.00018181692576035857, 0.00017982764984481037, 0.00017787083925213665, 0.0001759458245942369, 0.0001740519073791802, 0.00017218840366695076, 0.0001703546877251938, 0.00016855011926963925, 0.00016677408711984754, 0.00016502599464729428, 0.00016330524522345513, 0.0001616112858755514, 0.00015994356363080442, 0.0001583015255164355, 0.00015668464766349643, 0.00015509243530686945, 0.00015352437912952155, 0.00015197998436633497, 0.00015045879990793765, 0.0001489603309892118, 0.0001474841556046158, 0.0001460298226447776, 0.0001445968955522403, 0.00014318496687337756, 0.00014179362915456295, 0.00014042247494217008, 0.00013907112588640302, 0.0001377391890855506, 0.00013642631529364735, 0.00013513211160898209])

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.