diff --git a/src/Individuals/Plankton/IronEnergyMode/division_death.jl b/src/Individuals/Plankton/IronEnergyMode/division_death.jl
index f8d100a9..287561e6 100644
--- a/src/Individuals/Plankton/IronEnergyMode/division_death.jl
+++ b/src/Individuals/Plankton/IronEnergyMode/division_death.jl
@@ -41,6 +41,9 @@ end
         @inbounds plank.qNH4[idx[i]] = plank.qNH4[i]
         @inbounds plank.qP[idx[i]]   = plank.qP[i]
         @inbounds plank.qFe[idx[i]]  = plank.qFe[i]
+        @inbounds plank.qFePS[idx[i]] = plank.qFePS[i]
+        @inbounds plank.qFeNR[idx[i]] = plank.qFeNR[i]
+        @inbounds plank.qFeNF[idx[i]] = plank.qFeNF[i]
         @inbounds plank.exEn[idx[i]] = plank.exEn[i]
         @inbounds plank.Chl[idx[i]]  = plank.Chl[i]
         @inbounds plank.gen[idx[i]]  = plank.gen[i]
@@ -65,6 +68,9 @@ end
     @inbounds plank.qNH4[i] *= (2.0f0 - plank.dvid[i]) / 2.0f0 
     @inbounds plank.qP[i]   *= (2.0f0 - plank.dvid[i]) / 2.0f0 
     @inbounds plank.qFe[i]  *= (2.0f0 - plank.dvid[i]) / 2.0f0 
+    @inbounds plank.qFePS[i]*= (2.0f0 - plank.dvid[i]) / 2.0f0 
+    @inbounds plank.qFeNR[i]*= (2.0f0 - plank.dvid[i]) / 2.0f0 
+    @inbounds plank.qFeNF[i]*= (2.0f0 - plank.dvid[i]) / 2.0f0 
     @inbounds plank.exEn[i] *= (2.0f0 - plank.dvid[i]) / 2.0f0 
     @inbounds plank.Chl[i]  *= (2.0f0 - plank.dvid[i]) / 2.0f0 
     @inbounds plank.gen[i]  += plank.dvid[i]