diff --git a/generic_tracers/generic_COBALT.F90 b/generic_tracers/generic_COBALT.F90
index b6bba8f..c3eb4c7 100644
--- a/generic_tracers/generic_COBALT.F90
+++ b/generic_tracers/generic_COBALT.F90
@@ -7702,10 +7702,10 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
     !
     ! Sinking phytoplankton: Iron
     !
+    !> Iron flux to the sediment is removed, and flux from the sediment is
+    !! handled separately later using a relationship based on Dale et al., 2015. 
     call g_tracer_get_values(tracer_list,'fedi','btm_reservoir',phyto(DIAZO)%ffe_btm,isd,jsd)
     phyto(DIAZO)%ffe_btm = phyto(DIAZO)%ffe_btm/dt
-    !call g_tracer_get_pointer(tracer_list,'fedi_btf','field',temp_field)
-    !temp_field(:,:,1) = phyto(DIAZO)%ffe_btm(:,:)
     call g_tracer_set_values(tracer_list,'fedi','btm_reservoir',0.0)
     if (phyto(DIAZO)%id_ffe_btm .gt. 0)           &
          used = g_send_data(phyto(DIAZO)%id_ffe_btm,phyto(DIAZO)%ffe_btm, &
@@ -7714,8 +7714,6 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'felg','btm_reservoir',phyto(LARGE)%ffe_btm,isd,jsd)
     phyto(LARGE)%ffe_btm = phyto(LARGE)%ffe_btm/dt
-    !call g_tracer_get_pointer(tracer_list,'felg_btf','field',temp_field)
-    !temp_field(:,:,1) = phyto(LARGE)%ffe_btm(:,:)
     call g_tracer_set_values(tracer_list,'felg','btm_reservoir',0.0)
     if (phyto(LARGE)%id_ffe_btm .gt. 0)           &
          used = g_send_data(phyto(LARGE)%id_ffe_btm,phyto(LARGE)%ffe_btm, &
@@ -7724,8 +7722,6 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'femd','btm_reservoir',phyto(MEDIUM)%ffe_btm,isd,jsd)
     phyto(MEDIUM)%ffe_btm = phyto(MEDIUM)%ffe_btm/dt
-    !call g_tracer_get_pointer(tracer_list,'femd_btf','field',temp_field)
-    !temp_field(:,:,1) = phyto(MEDIUM)%ffe_btm(:,:)
     call g_tracer_set_values(tracer_list,'femd','btm_reservoir',0.0)
     if (phyto(MEDIUM)%id_ffe_btm .gt. 0)           &
          used = g_send_data(phyto(MEDIUM)%id_ffe_btm,phyto(MEDIUM)%ffe_btm, &
@@ -7734,8 +7730,6 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'fesm','btm_reservoir',phyto(SMALL)%ffe_btm,isd,jsd)
     phyto(SMALL)%ffe_btm = phyto(SMALL)%ffe_btm/dt
-    !call g_tracer_get_pointer(tracer_list,'fesm_btf','field',temp_field)
-    !temp_field(:,:,1) = phyto(SMALL)%ffe_btm(:,:)
     call g_tracer_set_values(tracer_list,'fesm','btm_reservoir',0.0)
     if (phyto(SMALL)%id_ffe_btm .gt. 0)           &
          used = g_send_data(phyto(SMALL)%id_ffe_btm,phyto(SMALL)%ffe_btm, &
@@ -7746,6 +7740,8 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
     !
     call g_tracer_get_values(tracer_list,'pdi','btm_reservoir',phyto(DIAZO)%fp_btm,isd,jsd)
     phyto(DIAZO)%fp_btm = phyto(DIAZO)%fp_btm/dt
+    call g_tracer_get_pointer(tracer_list,'pdi_btf','field',temp_field)
+    temp_field(:,:,1) = phyto(DIAZO)%fp_btm(:,:)
     call g_tracer_set_values(tracer_list,'pdi','btm_reservoir',0.0)
     if (phyto(DIAZO)%id_fp_btm .gt. 0)           &
          used = g_send_data(phyto(DIAZO)%id_fp_btm,phyto(DIAZO)%fp_btm, &
@@ -7754,6 +7750,8 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'plg','btm_reservoir',phyto(LARGE)%fp_btm,isd,jsd)
     phyto(LARGE)%fp_btm = phyto(LARGE)%fp_btm/dt
+    call g_tracer_get_pointer(tracer_list,'plg_btf','field',temp_field)
+    temp_field(:,:,1) = phyto(LARGE)%fp_btm(:,:)
     call g_tracer_set_values(tracer_list,'plg','btm_reservoir',0.0)
     if (phyto(LARGE)%id_fp_btm .gt. 0)           &
          used = g_send_data(phyto(LARGE)%id_fp_btm,phyto(LARGE)%fp_btm, &
@@ -7762,6 +7760,8 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'pmd','btm_reservoir',phyto(MEDIUM)%fp_btm,isd,jsd)
     phyto(MEDIUM)%fp_btm = phyto(MEDIUM)%fp_btm/dt
+    call g_tracer_get_pointer(tracer_list,'pmd_btf','field',temp_field)
+    temp_field(:,:,1) = phyto(MEDIUM)%fp_btm(:,:)
     call g_tracer_set_values(tracer_list,'pmd','btm_reservoir',0.0)
     if (phyto(MEDIUM)%id_fp_btm .gt. 0)           &
          used = g_send_data(phyto(MEDIUM)%id_fp_btm,phyto(MEDIUM)%fp_btm, &
@@ -7770,6 +7770,8 @@ subroutine generic_COBALT_update_from_bottom(tracer_list, dt, tau, model_time)
 
     call g_tracer_get_values(tracer_list,'psm','btm_reservoir',phyto(SMALL)%fp_btm,isd,jsd)
     phyto(SMALL)%fp_btm = phyto(SMALL)%fp_btm/dt
+    call g_tracer_get_pointer(tracer_list,'psm_btf','field',temp_field)
+    temp_field(:,:,1) = phyto(SMALL)%fp_btm(:,:)
     call g_tracer_set_values(tracer_list,'psm','btm_reservoir',0.0)
     if (phyto(SMALL)%id_fp_btm .gt. 0)           &
          used = g_send_data(phyto(SMALL)%id_fp_btm,phyto(SMALL)%fp_btm, &
@@ -9900,7 +9902,7 @@ subroutine generic_COBALT_update_from_source(tracer_list,Temp,Salt,rho_dzt,dzt,h
 
           ! iron from sediment (Elrod)
           !cobalt%ffe_sed(i,j) = cobalt%fe_2_n_sed * cobalt%f_ndet_btf(i,j,1)
-          ! iron from sediment (Dale)
+          ! iron from sediment (Dale, 2015)
           cobalt%ffe_sed(i,j) = cobalt%ffe_sed_max * tanh( (cobalt%fntot_btm(i,j)*cobalt%c_2_n*sperd*1.0e3)/ &
                                 max(cobalt%btm_o2(i,j)*1.0e6,epsln) )
           cobalt%ffe_geotherm(i,j) = cobalt%ffe_geotherm_ratio*internal_heat(i,j)*4184.0/dt