diff --git a/genie-biogem/src/fortran/biogem_box.f90 b/genie-biogem/src/fortran/biogem_box.f90
index 94f776a07..d26fef917 100644
--- a/genie-biogem/src/fortran/biogem_box.f90
+++ b/genie-biogem/src/fortran/biogem_box.f90
@@ -883,10 +883,10 @@ SUBROUTINE sub_calc_bio_uptake(dum_i,dum_j,dum_k1,dum_dt)
        ! prescribed POC flux
        ! NOTE: force_restore_docn_nuts has been assigned a negative sign (in the forcing update of the main biogem subroutine) ...
        ! NOTE: allow depletion of PO4 < 0.0 so as to force export production to the prescribed value
-       ! NOTE: correct for DOM (given that force_restore_docn_nuts has been devied from a prescribed particulate flux)
-       ! NOTE: assume par_bio_red_DOMfrac for the DOM fraction 
-       !       (loc_bio_red_DOMtotal has not been set to a non-zero value yet because the code that did it was moved back ...)
-       loc_dPO4 = -force_restore_docn_nuts(io_PO4)/(1.0 - par_bio_red_DOMfrac)
+       ! NOTE: there is no correction for DOM production (this occurs later)
+       !       hence at this point, loc_dPO4 is the PO4 uptake associated with POP, and not total production (POP+DOP)
+       !       (and hence idffers from the other biological scheme calculations of loc_dPO4) 
+       loc_dPO4 = -force_restore_docn_nuts(io_PO4)
     CASE ( &
          & '2N2T_PN_Tdep' &
          & )