From 476ad87224a5ae8d26ab61ea76445fbe216c4936 Mon Sep 17 00:00:00 2001 From: Jago Stong-Wright Date: Wed, 8 Nov 2023 12:22:25 +0000 Subject: [PATCH] removed `Relaxation` since it temporarily doesn't work --- docs/src/model_implementation.md | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/docs/src/model_implementation.md b/docs/src/model_implementation.md index 0dd37bc50..11c5428c9 100644 --- a/docs/src/model_implementation.md +++ b/docs/src/model_implementation.md @@ -235,7 +235,7 @@ using OceanBioME.Sediments: InstantRemineralisation @inline surface_PAR(t) = 200 * (1 - cos((t + 15days) * 2π / year)) * (1 / (1 + 0.2 * exp(-((mod(t, year) - 200days) / 50days)^2))) + 2 -@inline temp(z, t) = cos(t * 2π / year + 50days) + 28 +@inline ∂ₜT(z, t) = - 2π / year * sin(t * 2π / year + 50days) @inline κₜ(z, t) = 1e-2 * (1 + tanh((z - 50) / 10)) / 2 + 1e-4 @@ -267,7 +267,7 @@ biogeochemistry = Biogeochemistry(NutrientPhytoplankton(; sinking_velocity); model = NonhydrostaticModel(; grid, biogeochemistry, closure = ScalarDiffusivity(ν = κₜ, κ = κₜ), - forcing = (T = Relaxation(rate = 1/day, target = temp), )) + forcing = (; T = ∂ₜT)) set!(model, P = 0.01, N = 15, T = 28)