Update aed_phytoplankton.F90

Note that deposited phytoplankton mass in zones where mpb is not set as active was missing a -ve.
AquaticEcoDynamics · Mar 12, 2023 · 58ea02f · 58ea02f
1 parent 36dcc7f
commit 58ea02f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/aed_phytoplankton.F90 b/src/aed_phytoplankton.F90
@@ -1341,7 +1341,7 @@ SUBROUTINE aed_calculate_benthic_phytoplankton(data,column,layer_idx)
  fI = zero_
  mpb_prod = zero_
  mpb_resp = zero_
- mpb_burial = Psed_phy_c/secs_per_day ! all deposited material is buried
+ mpb_burial = -Psed_phy_c/secs_per_day ! set all deposited material to be buried (noting that Psed is -ve downawards)
  ENDIF
  ENDIF
  mpb_flux = (mpb_prod-mpb_resp)*mpb