removed Relaxation since it temporarily doesn't work

OceanBioME · Nov 8, 2023 · 476ad87 · 476ad87
1 parent b62b819
commit 476ad87
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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)