How to have vertical drift in only one speices of a biogeochemistry model? #4952
Unanswered
francispoulin
asked this question in
Computational science!
Replies: 1 comment 4 replies
-
|
I think the easiest way to do this would be with If you wanted to use the built in biogeochemistry infrastructure you can define your model like: using Oceananigans.Biogeochemistry: AbstractContinuousFormBiogeochemistry
struct NPZModel <: AbstractContinuousFormBiogeochemistry
parameters...
end
import Oceananigans.Biogeochemistry: biogeochemical_drift_velocity
biogeochemical_drift_velocity(bgc::NPZModel, ::Val{:D}) =
(u = ZeroField(), v = ZeroField(), w = ConstantField(-0.005)) |
Beta Was this translation helpful? Give feedback.
4 replies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
Below is a simple minimum working example that shows how you can add vertical advection to nutrients using a built in advection scheme. This uses first order upwinding to make sure that it preserves positivity and does not go beyond the maximum.
What I actually want to do is have an NPZ model say, and only have vertical drift of one species on top of the regular advection.
I think to do this i need to specify forcing to one of the species. Can someone show me how to do this using a built in advection scheme, rather than having to rebuild my own?
Beta Was this translation helpful? Give feedback.
All reactions