Merge pull request #53 from MichaelBarryBMT/Phyto-N-P-sed-fluxes

Modifies phyto N and P sedimentation
AquaticEcoDynamics · Mar 7, 2022 · 284d8ad · 284d8ad
2 parents 2b71366 + ee56b57
commit 284d8ad
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/aed_phytoplankton.F90 b/src/aed_phytoplankton.F90
@@ -1415,11 +1415,15 @@ SUBROUTINE aed_mobility_phytoplankton(data,column,layer_idx,mobility)
  IF(data%phytos(phy_i)%simINDynamics>0) THEN
  mobility(data%id_in(phy_i)) = vvel
  _DIAG_VAR_(data%id_PhySEDn(phy_i)) = vvel*_STATE_VAR_(data%id_in(phy_i))* secs_per_day
+ ELSE
+ _DIAG_VAR_(data%id_PhySEDn(phy_i)) = vvel*_STATE_VAR_(data%id_p(phy_i))* data%phytos(phy_i)%X_ncon* secs_per_day
  ENDIF
  ! set global mobility array for phytoplankton IP pool, and record flux
  IF(data%phytos(phy_i)%simIPDynamics>0) THEN
  mobility(data%id_ip(phy_i)) = vvel
  _DIAG_VAR_(data%id_PhySEDp(phy_i)) = vvel*_STATE_VAR_(data%id_ip(phy_i))* secs_per_day
+ ELSE
+ _DIAG_VAR_(data%id_PhySEDp(phy_i)) = vvel*_STATE_VAR_(data%id_p(phy_i))* data%phytos(phy_i)%X_pcon* secs_per_day
  ENDIF
 
  ! cumulate the community sedimentation flux (mmmol/m2/s) for later use