use biomass normalized f_T2B instead

Author: zhenwu0728

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.

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

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)