Merge pull request #245 from derpycode/master.bio_POCflux

fixed accidental double-DOM corection in prescribed POC flux scheme 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' &
  & )