Skip to content

Commit a63a6e9

Browse files
zhenwu0728zhenwu0728
authored
add IronEnergyMode (#75)
* add `IronEnergyMode` - IronEnergyMode resolves intracellular iron quota and energy budget. - Photosynthesis is separated into light reaction (light harvesting) and dark reaction (carbon fixation). - Intracellular iron is allocated to photosynthesis and NO3 reduction. - Iron allocation is dynamic and regulated by light. * remove unused cell division strategies * bug fix * update dependencies
1 parent 831f17c commit a63a6e9

38 files changed

+783
-116
lines changed

‎Manifest.toml

+49-41
Original file line numberDiff line numberDiff line change
@@ -1,8 +1,8 @@
11
# This file is machine-generated - editing it directly is not advised
22

3-
julia_version = "1.10.3"
3+
julia_version = "1.10.4"
44
manifest_format = "2.0"
5-
project_hash = "1078913c9829ff096c73958654eaca165e82cef9"
5+
project_hash = "5e0f2ff9bdf8afe520c669148f1610daa9a8ddf0"
66

77
[[deps.AbstractFFTs]]
88
deps = ["LinearAlgebra"]
@@ -57,35 +57,37 @@ version = "0.5.0"
5757

5858
[[deps.CUDA]]
5959
deps = ["AbstractFFTs", "Adapt", "BFloat16s", "CEnum", "CUDA_Driver_jll", "CUDA_Runtime_Discovery", "CUDA_Runtime_jll", "Crayons", "DataFrames", "ExprTools", "GPUArrays", "GPUCompiler", "KernelAbstractions", "LLVM", "LLVMLoopInfo", "LazyArtifacts", "Libdl", "LinearAlgebra", "Logging", "NVTX", "Preferences", "PrettyTables", "Printf", "Random", "Random123", "RandomNumbers", "Reexport", "Requires", "SparseArrays", "StaticArrays", "Statistics"]
60-
git-tree-sha1 = "4e33522a036b39fc6f5cb7447ae3b28eb8fbe99b"
60+
git-tree-sha1 = "6e945e876652f2003e6ca74e19a3c45017d3e9f6"
6161
uuid = "052768ef-5323-5732-b1bb-66c8b64840ba"
62-
version = "5.3.3"
62+
version = "5.4.2"
6363

6464
[deps.CUDA.extensions]
6565
ChainRulesCoreExt = "ChainRulesCore"
66+
EnzymeCoreExt = "EnzymeCore"
6667
SpecialFunctionsExt = "SpecialFunctions"
6768

6869
[deps.CUDA.weakdeps]
6970
ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
71+
EnzymeCore = "f151be2c-9106-41f4-ab19-57ee4f262869"
7072
SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"
7173

7274
[[deps.CUDA_Driver_jll]]
7375
deps = ["Artifacts", "JLLWrappers", "LazyArtifacts", "Libdl", "Pkg"]
74-
git-tree-sha1 = "dc172b558adbf17952001e15cf0d6364e6d78c2f"
76+
git-tree-sha1 = "c48f9da18efd43b6b7adb7ee1f93fe5f2926c339"
7577
uuid = "4ee394cb-3365-5eb0-8335-949819d2adfc"
76-
version = "0.8.1+0"
78+
version = "0.9.0+0"
7779

7880
[[deps.CUDA_Runtime_Discovery]]
7981
deps = ["Libdl"]
80-
git-tree-sha1 = "38f830504358e9972d2a0c3e5d51cb865e0733df"
82+
git-tree-sha1 = "f3b237289a5a77c759b2dd5d4c2ff641d67c4030"
8183
uuid = "1af6417a-86b4-443c-805f-a4643ffb695f"
82-
version = "0.2.4"
84+
version = "0.3.4"
8385

8486
[[deps.CUDA_Runtime_jll]]
8587
deps = ["Artifacts", "CUDA_Driver_jll", "JLLWrappers", "LazyArtifacts", "Libdl", "TOML"]
86-
git-tree-sha1 = "4ca7d6d92075906c2ce871ea8bba971fff20d00c"
88+
git-tree-sha1 = "bcba305388e16aa5c879e896726db9e71b4942c6"
8789
uuid = "76a88914-d11a-5bdc-97e0-2f5a05c973a2"
88-
version = "0.12.1+0"
90+
version = "0.14.0+1"
8991

9092
[[deps.ColorTypes]]
9193
deps = ["FixedPointNumbers", "Random"]
@@ -95,9 +97,9 @@ version = "0.11.5"
9597

9698
[[deps.Colors]]
9799
deps = ["ColorTypes", "FixedPointNumbers", "Reexport"]
98-
git-tree-sha1 = "fc08e5930ee9a4e03f84bfb5211cb54e7769758a"
100+
git-tree-sha1 = "362a287c3aa50601b0bc359053d5c2468f0e7ce0"
99101
uuid = "5ae59095-9a9b-59fe-a467-6f913c188581"
100-
version = "0.12.10"
102+
version = "0.12.11"
101103

102104
[[deps.Compat]]
103105
deps = ["TOML", "UUIDs"]
@@ -180,19 +182,19 @@ uuid = "7b1f6079-737a-58dc-b8bc-7a2ca5c1b5ee"
180182

181183
[[deps.FixedPointNumbers]]
182184
deps = ["Statistics"]
183-
git-tree-sha1 = "335bfdceacc84c5cdf16aadc768aa5ddfc5383cc"
185+
git-tree-sha1 = "05882d6995ae5c12bb5f36dd2ed3f61c98cbb172"
184186
uuid = "53c48c17-4a7d-5ca2-90c5-79b7896eea93"
185-
version = "0.8.4"
187+
version = "0.8.5"
186188

187189
[[deps.Future]]
188190
deps = ["Random"]
189191
uuid = "9fa8497b-333b-5362-9e8d-4d0656e87820"
190192

191193
[[deps.GPUArrays]]
192194
deps = ["Adapt", "GPUArraysCore", "LLVM", "LinearAlgebra", "Printf", "Random", "Reexport", "Serialization", "Statistics"]
193-
git-tree-sha1 = "68e8ff56a4a355a85d2784b94614491f8c900cde"
195+
git-tree-sha1 = "c154546e322a9c73364e8a60430b0f79b812d320"
194196
uuid = "0c68f7d7-f131-5f86-a1c3-88cf8149b2d7"
195-
version = "10.1.0"
197+
version = "10.2.0"
196198

197199
[[deps.GPUArraysCore]]
198200
deps = ["Adapt"]
@@ -202,15 +204,21 @@ version = "0.1.6"
202204

203205
[[deps.GPUCompiler]]
204206
deps = ["ExprTools", "InteractiveUtils", "LLVM", "Libdl", "Logging", "Scratch", "TimerOutputs", "UUIDs"]
205-
git-tree-sha1 = "1600477fba37c9fc067b9be21f5e8101f24a8865"
207+
git-tree-sha1 = "518ebd058c9895de468a8c255797b0c53fdb44dd"
206208
uuid = "61eb1bfa-7361-4325-ad38-22787b887f55"
207-
version = "0.26.4"
209+
version = "0.26.5"
208210

209211
[[deps.InlineStrings]]
210212
deps = ["Parsers"]
211-
git-tree-sha1 = "9cc2baf75c6d09f9da536ddf58eb2f29dedaf461"
213+
git-tree-sha1 = "86356004f30f8e737eff143d57d41bd580e437aa"
212214
uuid = "842dd82b-1e85-43dc-bf29-5d0ee9dffc48"
213-
version = "1.4.0"
215+
version = "1.4.1"
216+
217+
[deps.InlineStrings.extensions]
218+
ArrowTypesExt = "ArrowTypes"
219+
220+
[deps.InlineStrings.weakdeps]
221+
ArrowTypes = "31f734f8-188a-4ce0-8406-c8a06bd891cd"
214222

215223
[[deps.InteractiveUtils]]
216224
deps = ["Markdown"]
@@ -227,10 +235,10 @@ uuid = "82899510-4779-5014-852e-03e436cf321d"
227235
version = "1.0.0"
228236

229237
[[deps.JLD2]]
230-
deps = ["FileIO", "MacroTools", "Mmap", "OrderedCollections", "Pkg", "PrecompileTools", "Printf", "Reexport", "Requires", "TranscodingStreams", "UUIDs"]
231-
git-tree-sha1 = "5ea6acdd53a51d897672edb694e3cc2912f3f8a7"
238+
deps = ["FileIO", "MacroTools", "Mmap", "OrderedCollections", "Pkg", "PrecompileTools", "Reexport", "Requires", "TranscodingStreams", "UUIDs", "Unicode"]
239+
git-tree-sha1 = "bdbe8222d2f5703ad6a7019277d149ec6d78c301"
232240
uuid = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
233-
version = "0.4.46"
241+
version = "0.4.48"
234242

235243
[[deps.JLLWrappers]]
236244
deps = ["Artifacts", "Preferences"]
@@ -246,9 +254,9 @@ version = "0.2.1+0"
246254

247255
[[deps.KernelAbstractions]]
248256
deps = ["Adapt", "Atomix", "InteractiveUtils", "LinearAlgebra", "MacroTools", "PrecompileTools", "Requires", "SparseArrays", "StaticArrays", "UUIDs", "UnsafeAtomics", "UnsafeAtomicsLLVM"]
249-
git-tree-sha1 = "ed7167240f40e62d97c1f5f7735dea6de3cc5c49"
257+
git-tree-sha1 = "8e5a339882cc401688d79b811d923a38ba77d50a"
250258
uuid = "63c18a36-062a-441e-b654-da1e3ab1ce7c"
251-
version = "0.9.18"
259+
version = "0.9.20"
252260

253261
[deps.KernelAbstractions.extensions]
254262
EnzymeExt = "EnzymeCore"
@@ -258,9 +266,9 @@ version = "0.9.18"
258266

259267
[[deps.LLVM]]
260268
deps = ["CEnum", "LLVMExtra_jll", "Libdl", "Preferences", "Printf", "Requires", "Unicode"]
261-
git-tree-sha1 = "839c82932db86740ae729779e610f07a1640be9a"
269+
git-tree-sha1 = "389aea28d882a40b5e1747069af71bdbd47a1cae"
262270
uuid = "929cbde3-209d-540e-8aea-75f648917ca0"
263-
version = "6.6.3"
271+
version = "7.2.1"
264272
weakdeps = ["BFloat16s"]
265273

266274
[deps.LLVM.extensions]
@@ -405,9 +413,9 @@ version = "1.4.3"
405413

406414
[[deps.PrettyTables]]
407415
deps = ["Crayons", "LaTeXStrings", "Markdown", "PrecompileTools", "Printf", "Reexport", "StringManipulation", "Tables"]
408-
git-tree-sha1 = "88b895d13d53b5577fd53379d913b9ab9ac82660"
416+
git-tree-sha1 = "66b20dd35966a748321d3b2537c4584cf40387c7"
409417
uuid = "08abe8d2-0d0c-5749-adfa-8a2ac140af0d"
410-
version = "2.3.1"
418+
version = "2.3.2"
411419

412420
[[deps.Printf]]
413421
deps = ["Unicode"]
@@ -456,9 +464,9 @@ version = "1.2.1"
456464

457465
[[deps.SentinelArrays]]
458466
deps = ["Dates", "Random"]
459-
git-tree-sha1 = "0e7508ff27ba32f26cd459474ca2ede1bc10991f"
467+
git-tree-sha1 = "90b4f68892337554d31cdcdbe19e48989f26c7e6"
460468
uuid = "91c51154-3ec4-41a3-a24f-3f23e20d615c"
461-
version = "1.4.1"
469+
version = "1.4.3"
462470

463471
[[deps.Serialization]]
464472
uuid = "9e88b42a-f829-5b0c-bbe9-9e923198166b"
@@ -479,9 +487,9 @@ version = "1.10.0"
479487

480488
[[deps.StaticArrays]]
481489
deps = ["LinearAlgebra", "PrecompileTools", "Random", "StaticArraysCore"]
482-
git-tree-sha1 = "bf074c045d3d5ffd956fa0a461da38a44685d6b2"
490+
git-tree-sha1 = "6e00379a24597be4ae1ee6b2d882e15392040132"
483491
uuid = "90137ffa-7385-5640-81b9-e52037218182"
484-
version = "1.9.3"
492+
version = "1.9.5"
485493

486494
[deps.StaticArrays.extensions]
487495
StaticArraysChainRulesCoreExt = "ChainRulesCore"
@@ -492,9 +500,9 @@ version = "1.9.3"
492500
Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
493501

494502
[[deps.StaticArraysCore]]
495-
git-tree-sha1 = "36b3d696ce6366023a0ea192b4cd442268995a0d"
503+
git-tree-sha1 = "192954ef1208c7019899fbf8049e717f92959682"
496504
uuid = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
497-
version = "1.4.2"
505+
version = "1.4.3"
498506

499507
[[deps.Statistics]]
500508
deps = ["LinearAlgebra", "SparseArrays"]
@@ -553,14 +561,14 @@ uuid = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
553561

554562
[[deps.TimerOutputs]]
555563
deps = ["ExprTools", "Printf"]
556-
git-tree-sha1 = "f548a9e9c490030e545f72074a41edfd0e5bcdd7"
564+
git-tree-sha1 = "5a13ae8a41237cff5ecf34f73eb1b8f42fff6531"
557565
uuid = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
558-
version = "0.5.23"
566+
version = "0.5.24"
559567

560568
[[deps.TranscodingStreams]]
561-
git-tree-sha1 = "5d54d076465da49d6746c647022f3b3674e64156"
569+
git-tree-sha1 = "a947ea21087caba0a798c5e494d0bb78e3a1a3a0"
562570
uuid = "3bb67fe8-82b1-5028-8e26-92a6c54297fa"
563-
version = "0.10.8"
571+
version = "0.10.9"
564572
weakdeps = ["Random", "Test"]
565573

566574
[deps.TranscodingStreams.extensions]
@@ -580,9 +588,9 @@ version = "0.2.1"
580588

581589
[[deps.UnsafeAtomicsLLVM]]
582590
deps = ["LLVM", "UnsafeAtomics"]
583-
git-tree-sha1 = "323e3d0acf5e78a56dfae7bd8928c989b4f3083e"
591+
git-tree-sha1 = "d9f5962fecd5ccece07db1ff006fb0b5271bdfdd"
584592
uuid = "d80eeb9a-aca5-4d75-85e5-170c8b632249"
585-
version = "0.1.3"
593+
version = "0.1.4"
586594

587595
[[deps.Zlib_jll]]
588596
deps = ["Libdl"]

‎src/Biogeochemistry/nutrient_fields.jl

+5-3
Original file line numberDiff line numberDiff line change
@@ -3,7 +3,9 @@ function nutrients_init(arch, g, FT = Float32)
33
Field(arch, g, FT), Field(arch, g, FT),
44
Field(arch, g, FT), Field(arch, g, FT),
55
Field(arch, g, FT), Field(arch, g, FT),
6-
Field(arch, g, FT), Field(arch, g, FT))
6+
Field(arch, g, FT), Field(arch, g, FT),
7+
Field(arch, g, FT), Field(arch, g, FT),
8+
Field(arch, g, FT))
79

810
nut = NamedTuple{nut_names}(fields)
911
return nut
@@ -14,8 +16,8 @@ end
1416
Generate defalut nutrient initial conditions.
1517
"""
1618
function default_nut_init()
17-
init = (DIC=20.0, NH4=0.5, NO3=0.8, PO4=0.10, DOC=10.0, DON=0.1, DOP=0.05, POC=0.0, PON=0.0, POP=0.0)
18-
rand_noise = (DIC=0.1, NH4=0.1, NO3=0.1, PO4=0.1, DOC=0.1, DON=0.1, DOP=0.1, POC=0.1, PON=0.1, POP=0.1)
19+
init = (DIC=20.0, NH4=0.5, NO3=0.8, PO4=0.10, FeT=1.0e-6, DOC=10.0, DON=0.1, DOP=0.05, DOFe=0.0, POC=0.0, PON=0.0, POP=0.0, POFe=0.0)
20+
rand_noise = (DIC=0.1, NH4=0.1, NO3=0.1, PO4=0.1, FeT=0.1, DOC=0.1, DON=0.1, DOP=0.1, DOFe=0.1, POC=0.1, PON=0.1, POP=0.1, POFe=0.1)
1921
return (initial_condition = init, rand_noise = rand_noise)
2022
end
2123

‎src/Biogeochemistry/nutrient_forcings.jl

+7
Original file line numberDiff line numberDiff line change
@@ -28,6 +28,13 @@ function nut_forcing!(F, nut_temp, nut, params, ΔT)
2828

2929
@inbounds F.POP.data .= F.POP.data .- nut_temp.POP.data .* params["kPOP"] .* ΔT
3030

31+
@inbounds F.FeT.data .= F.FeT.data .+ nut_temp.DOFe.data .* params["kDOFe"] .* ΔT
32+
33+
@inbounds F.DOFe.data .= F.DOFe.data .- nut_temp.DOFe.data .* params["kDOFe"] .* ΔT .+
34+
nut_temp.POFe.data .* params["kPOFe"] .* ΔT
35+
36+
@inbounds F.POFe.data .= F.POFe.data .- nut_temp.POFe.data .* params["kPOFe"] .* ΔT
37+
3138
return nothing
3239
end
3340

‎src/Diagnostics/Diagnostics.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -11,7 +11,7 @@ using PlanktonIndividuals.Architectures: device, Architecture, array_type
1111
using PlanktonIndividuals.Grids
1212
using PlanktonIndividuals.Fields
1313

14-
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode
14+
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode, IronEnergyMode
1515

1616
import Base: show
1717

‎src/Diagnostics/diagnostics_kernels.jl

+2
Original file line numberDiff line numberDiff line change
@@ -19,6 +19,8 @@ function diags_spcs!(diags_sp, plank, ac, x, y, z, mode::AbstractMode, arch::Arc
1919
diags = (:PS, :BS, :VDOC, :VNH4, :VNO3, :VPO4, :resp, :exu, :Bm, :Cq, :Nq, :Pq, :Chl)
2020
elseif isa(mode, MacroMolecularMode)
2121
diags = (:PS, :VDOC, :VHN4, :VNO3, :VPO4, :S_PRO, :S_DNA, :S_RNA, :resp, :ρChl, :CH, :NST, :PST, :PRO, :DNA, :RNA, :Chl)
22+
elseif isa(mode, IronEnergyMode)
23+
diags = (:PS, :CF, :ECF, :RS, :ERS, :NR, :ENR, :BS, :VNH4, :VNO3, :VPO4, :VFe, :Bm, :En, :CH, :qNO3, :qNH4, :qP, :qFe, :Chl)
2224
end
2325

2426
for diag in keys(diags_sp)

‎src/Diagnostics/diagnostics_struct.jl

+3-1
Original file line numberDiff line numberDiff line change
@@ -98,7 +98,7 @@ function diag_avail(tracer, plank, mode::AbstractMode)
9898
end
9999

100100
function tracer_avail_diags()
101-
return (:PAR, :DIC, :DOC, :POC, :NH4, :NO3, :DON, :PON, :PO4, :DOP, :POP)
101+
return (:PAR, :DIC, :DOC, :POC, :NH4, :NO3, :DON, :PON, :PO4, :DOP, :POP, :FeT, :DOFe, :POFe)
102102
end
103103

104104
function plank_avail_diags(mode::AbstractMode)
@@ -109,6 +109,8 @@ function plank_avail_diags(mode::AbstractMode)
109109
plank_avail = (:num, :graz, :mort, :dvid, :PS, :BS, :VDOC, :VNH4, :VNO3, :VPO4, :resp, :exu, :Bm, :Cq, :Nq, :Pq, :Chl)
110110
elseif isa(mode, MacroMolecularMode)
111111
plank_avail = (:num, :graz, :mort, :dvid, :PS, :VDOC, :VNH4, :VNO3, :VPO4, :resp, :ρChl, :S_PRO, :S_DNA, :S_RNA, :exu, :CH, :NST, :PST, :PRO, :DNA, :RNA, :Chl)
112+
elseif isa(mode, IronEnergyMode)
113+
plank_avail = (:num, :graz, :mort, :dvid, :PS, :CF, :ECF, :RS, :ERS, :NR, :ENR, :BS, :VNH4, :VNO3, :VPO4, :VFe, :Bm, :En, :CH, :qNO3, :qNH4, :qP, :qFe, :Chl)
112114
end
113115
return plank_avail
114116
end

‎src/Fields/Fields.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -43,6 +43,6 @@ function zero_fields!(a)
4343
end
4444
end
4545

46-
const nut_names=(:DIC,:NH4,:NO3,:PO4,:DOC,:DON,:DOP,:POC,:PON,:POP)
46+
const nut_names=(:DIC,:NH4,:NO3,:PO4,:FeT,:DOC,:DON,:DOP,:DOFe,:POC,:PON,:POP,:POFe)
4747

4848
end

‎src/Model/Model.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -15,7 +15,7 @@ using PlanktonIndividuals.Parameters
1515
using PlanktonIndividuals.Plankton
1616
using PlanktonIndividuals.Diagnostics
1717

18-
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode
18+
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode, IronEnergyMode
1919

2020
import Base: show
2121

‎src/Model/time_step.jl

+8-6
Original file line numberDiff line numberDiff line change
@@ -43,7 +43,8 @@ function TimeStep!(model::PlanktonModel, ΔT, diags::PlanktonDiagnostics)
4343
model.individuals.phytos[sp].data.yi, model.individuals.phytos[sp].data.zi,
4444
model.individuals.phytos[sp].data.ac, model.nutrients.NH4.data,
4545
model.nutrients.NO3.data, model.nutrients.PO4.data, model.nutrients.DOC.data,
46-
model.timestepper.par, model.timestepper.temp, model.timestepper.pop, model.arch)
46+
model.nutrients.FeT.data, model.timestepper.par, model.timestepper.temp,
47+
model.timestepper.pop, model.arch)
4748

4849
plankton_update!(model.individuals.phytos[sp].data, model.timestepper.nuts,
4950
model.timestepper.rnd, model.individuals.phytos[sp].p,
@@ -65,11 +66,12 @@ function TimeStep!(model::PlanktonModel, ΔT, diags::PlanktonDiagnostics)
6566
model.individuals.phytos[sp].data.xi, model.individuals.phytos[sp].data.yi,
6667
model.individuals.phytos[sp].data.zi, model.arch)
6768

68-
find_NPT!(model.timestepper.nuts, model.individuals.phytos[sp].data.xi,
69-
model.individuals.phytos[sp].data.yi, model.individuals.phytos[sp].data.zi,
70-
model.individuals.phytos[sp].data.ac, model.nutrients.NH4.data,
71-
model.nutrients.NO3.data, model.nutrients.PO4.data, model.nutrients.DOC.data,
72-
model.timestepper.par, model.timestepper.temp, model.timestepper.pop, model.arch)
69+
find_NPT!(model.timestepper.nuts, model.individuals.phytos[sp].data.xi,
70+
model.individuals.phytos[sp].data.yi, model.individuals.phytos[sp].data.zi,
71+
model.individuals.phytos[sp].data.ac, model.nutrients.NH4.data,
72+
model.nutrients.NO3.data, model.nutrients.PO4.data, model.nutrients.DOC.data,
73+
model.nutrients.FeT.data, model.timestepper.par, model.timestepper.temp,
74+
model.timestepper.pop, model.arch)
7375

7476
plankton_update!(model.individuals.phytos[sp].data, model.timestepper.nuts,
7577
model.timestepper.rnd, model.individuals.phytos[sp].p,

‎src/Model/timestepper.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ function timestepper(arch::Architecture, FT::DataType, g::AbstractGrid, maxN)
4141
velos_d = replace_storage(array_type(arch), velos)
4242

4343
nuts = StructArray(NH4 = zeros(FT, maxN), NO3 = zeros(FT, maxN), PO4 = zeros(FT, maxN),
44-
DOC = zeros(FT, maxN), idc = zeros(Int, maxN),
44+
DOC = zeros(FT, maxN), FeT = zeros(FT, maxN), idc = zeros(Int,maxN),
4545
par = zeros(FT, maxN), T = zeros(FT, maxN), pop = zeros(FT, maxN))
4646
nuts_d = replace_storage(array_type(arch), nuts)
4747

‎src/Output/Output.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -12,7 +12,7 @@ using PlanktonIndividuals.Fields
1212
using PlanktonIndividuals.Diagnostics
1313
using PlanktonIndividuals.Model
1414

15-
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode
15+
using PlanktonIndividuals: AbstractMode, CarbonMode, QuotaMode, MacroMolecularMode, IronEnergyMode
1616

1717
import Base: show
1818

0 commit comments

Comments
 (0)