BOD: Benthic oxygen demand

In DanielReedOcean/COOLBEANS: Coastal Ocean Oxygen Linked to Benthic Exchange and Nutrient Supply model

Description

Calculate the benthic oxygen demand of a coastal sediment for given conditions.

Usage

BOD(z, NPP, C.flux.extra = 0)

Arguments

z	Water depth (m)
NPP	Net primary production (mol C/[m2 y])
C.flux.extra	Additional organic carbon fluxes to sediments (mol C/[m2 y])

Value

Returns benthic oxygen demand of sediment (mol O2/[m2 y])

Author(s)

Daniel C. Reed

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (z, NPP, T, C.flux.extra = 0) 
{
    C.flux.sed <- C.flux(NPP, z) + C.flux.extra
    Oxygen.demand <- (((77/60) * p.remin(z) - (17/60)) * C.flux.sed)
    if (any(NPP <= 0)) 
        warning(paste("Net primary production is ", NPP[NPP <= 
            0], " - must be positive.\n", sep = ""))
    if (any(z <= 0)) 
        warning(paste("Water depth is ", z[z <= 0], " - must be positive.\n", 
            sep = ""))
    return(Oxygen.demand)
  }

DanielReedOcean/COOLBEANS documentation built on May 6, 2019, 1:35 p.m.