From 5d45a18ab1ea9ee97de7db9080109b0ed30b555a Mon Sep 17 00:00:00 2001 From: Zhen Wu Date: Wed, 25 Oct 2023 16:52:09 -0400 Subject: [PATCH] use biomass normalized f_T2B instead --- src/Parameters/param_default.jl | 2 +- src/Plankton/CarbonMode/division_death.jl | 1 - src/Plankton/CarbonMode/growth_kernels.jl | 2 +- 3 files changed, 2 insertions(+), 3 deletions(-) diff --git a/src/Parameters/param_default.jl b/src/Parameters/param_default.jl index 79b7e1f0..ab33da80 100644 --- a/src/Parameters/param_default.jl +++ b/src/Parameters/param_default.jl @@ -217,7 +217,7 @@ function phyt_params_default(N::Int64, mode::CarbonMode) "Ea" => [5.3e4], # Free energy "PCmax" => [4.2e-5], # Maximum primary production rate (per second) "respir" => [1.2e-6], # Respiration rate(per second) - "f_T2B" => [1.0e-16], # Thermal damage rate (mmolC/K/s) + "f_T2B" => [2.7e-6], # Thermal damage rate (1.0/K/s) "grz_P" => [0.0], # Grazing probability per second "dvid_P" => [5e-5], # Probability of cell division per second. "dvid_type" => [1], # The type of cell division, 1:sizer, 2:adder. diff --git a/src/Plankton/CarbonMode/division_death.jl b/src/Plankton/CarbonMode/division_death.jl index 1c256747..3f5a1e6f 100644 --- a/src/Plankton/CarbonMode/division_death.jl +++ b/src/Plankton/CarbonMode/division_death.jl @@ -83,6 +83,5 @@ function divide!(plank, deactive_ind, arch::Architecture) get_tind!(plank.idx, plank.dvid, Int.(con_ind), deactive_ind, arch) copy_daughter_individuals!(plank, plank.dvid, Int.(plank.idx), arch) divide_to_half!(plank, arch) - # println("time step ends") return nothing end diff --git a/src/Plankton/CarbonMode/growth_kernels.jl b/src/Plankton/CarbonMode/growth_kernels.jl index d8a7fe98..05c511c7 100644 --- a/src/Plankton/CarbonMode/growth_kernels.jl +++ b/src/Plankton/CarbonMode/growth_kernels.jl @@ -71,7 +71,7 @@ end ##### keep thermal damage between 0.0 and healthy biomass (Bm-Bd) @kernel function calc_thermal_damage_kernel!(plank, T, p, ΔT) i = @index(Global) - @inbounds plank.TD[i] = (T[i] - p.Topt) * p.f_T2B * + @inbounds plank.TD[i] = (T[i] - p.Topt) * p.f_T2B * (plank.Bm[i] - plank.Bd[i]) isless(p.Topt, T[i]) * isless(plank.Bd[i], plank.Bm[i]) * isless(0.0, p.thermal)