bug fix in cell division again

JuliaOcean · Oct 19, 2024 · 6446c36 · 6446c36
1 parent 44814bb
commit 6446c36
Show file tree

Hide file tree

Showing 8 changed files with 6 additions and 6 deletions.
diff --git a/src/Plankton/CarbonMode/division_death.jl b/src/Plankton/CarbonMode/division_death.jl
@@ -83,5 +83,6 @@ function divide!(plank, nuts, deactive_ind, arch::Architecture)
     get_tind!(plank.idx, plank.dvid, nuts.idc_int, deactive_ind, arch)
     copy_daughter_individuals!(plank, plank.dvid, plank.idx, arch)
     divide_to_half!(plank, arch)
+    plank.idx .= 0
     return nothing
 end
diff --git a/src/Plankton/CarbonMode/plankton_update.jl b/src/Plankton/CarbonMode/plankton_update.jl
@@ -30,7 +30,6 @@ function plankton_update!(plank, nuts, rnd, p, plk, diags_spcs, ΔT, t, arch::Ar
             throw(ArgumentError("number of individual exceeds the capacity at timestep $(t/86400.0) days"))
         end
         divide!(plank, nuts, deactive_ind, arch)
-        nuts.idc .= 0
     end
 
     ##### diagnostic for individual distribution

diff --git a/src/Plankton/IronEnergyMode/division_death.jl b/src/Plankton/IronEnergyMode/division_death.jl
@@ -96,6 +96,7 @@ function divide!(plank, nuts, deactive_ind, arch::Architecture)
     nuts.idc_int .= unsafe_trunc.(Int, nuts.idc)
     get_tind!(plank.idx, plank.dvid, nuts.idc_int, deactive_ind, arch)
     copy_daughter_individuals!(plank, plank.dvid, plank.idx, arch)
-    divide_to_half!(plank, arch)    
+    divide_to_half!(plank, arch)
+    plank.idx .= 0
     return nothing
 end
diff --git a/src/Plankton/IronEnergyMode/plankton_update.jl b/src/Plankton/IronEnergyMode/plankton_update.jl
@@ -31,7 +31,6 @@ function plankton_update!(plank, nuts, rnd, p, plk, diags_spcs, ΔT, t, arch::Ar
             throw(ArgumentError("number of individual exceeds the capacity at timestep $(t/86400.0) days"))
         end
         divide!(plank, nuts, deactive_ind, arch)
-        nuts.idc .= 0
     end
 
     ##### diagnostic for individual distribution

diff --git a/src/Plankton/MacroMolecularMode/division_death.jl b/src/Plankton/MacroMolecularMode/division_death.jl
@@ -92,5 +92,6 @@ function divide!(plank, nuts, deactive_ind, arch::Architecture)
     get_tind!(plank.idx, plank.dvid, nuts.idc_int, deactive_ind, arch)
     copy_daughter_individuals!(plank, plank.dvid, plank.idx, arch)
     divide_to_half!(plank, arch)
+    plank.idx .= 0
     return nothing
 end
diff --git a/src/Plankton/MacroMolecularMode/plankton_update.jl b/src/Plankton/MacroMolecularMode/plankton_update.jl
@@ -31,7 +31,6 @@ function plankton_update!(plank, nuts, rnd, p, plk, diags_spcs, ΔT, t, arch::Ar
             throw(ArgumentError("number of individual exceeds the capacity at timestep $(t/86400.0) days"))
         end
         divide!(plank, nuts, deactive_ind, arch)
-        nuts.idc .= 0
     end
 
     ##### diagnostic for individual distribution

diff --git a/src/Plankton/QuotaMode/division_death.jl b/src/Plankton/QuotaMode/division_death.jl
@@ -90,6 +90,7 @@ function divide!(plank, nuts, deactive_ind, arch::Architecture)
     nuts.idc_int .= unsafe_trunc.(Int, nuts.idc)
     get_tind!(plank.idx, plank.dvid, nuts.idc_int, deactive_ind, arch)
     copy_daughter_individuals!(plank, plank.dvid, plank.idx, arch)
-    divide_to_half!(plank, arch)    
+    divide_to_half!(plank, arch)
+    plank.idx .= 0
     return nothing
 end
diff --git a/src/Plankton/QuotaMode/plankton_update.jl b/src/Plankton/QuotaMode/plankton_update.jl
@@ -31,7 +31,6 @@ function plankton_update!(plank, nuts, rnd, p, plk, diags_spcs, ΔT, t, arch::Ar
             throw(ArgumentError("number of individual exceeds the capacity at timestep $(t/86400.0) days"))
         end
         divide!(plank, nuts, deactive_ind, arch)
-        nuts.idc .= 0
     end
 
     ##### diagnostic for individual distribution