finally decided on the function of thermal history

JuliaOcean · Sep 6, 2023 · 82c8865 · 82c8865
1 parent f7d8c2d
commit 82c8865
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/Plankton/CarbonMode/growth_kernels.jl b/src/Plankton/CarbonMode/growth_kernels.jl
@@ -75,9 +75,8 @@ end
 ##### track temperature history
 @kernel function calc_thermal_history_kernel!(plank, nuts, p, ΔT)
     i = @index(Global)
-    # @inbounds plank.Th[i] += ΔT/3600 * exp(0.22*(nuts.T[i] - p.Topt)) * isless(p.Topt, nuts.T[i])
     @inbounds plank.Th[i] = plank.Th[i] * (1.0 - 1e-1 * ΔT/3600) +
-                            ΔT/3600 * exp(0.22*(nuts.T[i] - p.Topt)) * isless(p.Topt, nuts.T[i])
+                            ΔT/3600 * (nuts.T[i] - p.Topt) * isless(p.Topt, nuts.T[i])
 end
 function calc_thermal_history!(plank, nuts, p, ΔT, arch)
     kernel! = calc_thermal_history_kernel!(device(arch), 256, (size(plank.ac,1)))