FESDIAbudget {FESDIA}R Documentation

Budgets for the CNPDIA or MPBDIA model.

Description

FESDIAbudgetO2, FESDIAbudgetC, FESDIAbudgetN, FESDIAbudgetP estimate mass budgets from FESDIA model solutions.

Usage

  FESDIAbudgetO2(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 
  FESDIAbudgetC(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 
  FESDIAbudgetN(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 
  FESDIAbudgetP(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 
  FESDIAbudgetS(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 
  FESDIAbudgetFe(out, ..., which = c("All", "Rates", "Fluxes", "Losses", "Fluxmat")) 

Arguments

out

an output object returned by FESDIAsolve or FESDIAdyna.

which

if not NULL, a vector with names of the items to return.

...

unused.

Value

FESDIAbudgetx returns the element budget (C, N, P, S, Fe, O2) of the solution as a list, with the following items.

If more than one FESDIA object is passed to these functions, a matrix is returned, one column for each object.

For dynamic runs, the budget is taken over the mean of the run; Delta is then the integrated mean rate of change; for steady-state runs, Delta should be very small.

Author(s)

Karline Soetaert

References

Soetaert K, PMJ Herman and JJ Middelburg, 1996a. A model of early diagenetic processes from the shelf to abyssal depths. Geochimica Cosmochimica Acta, 60(6):1019-1040.

Soetaert K, PMJ Herman and JJ Middelburg, 1996b. Dynamic response of deep-sea sediments to seasonal variation: a model. Limnol. Oceanogr. 41(8): 1651-1668.

Examples


# some runs to work with  
  defsteady <- FESDIAsolve()
  defdyn    <- FESDIAdyna()

# altered steady-state run
  out <- FESDIAsolve(parms = list(Cflux = 1000))
  cbind(default = FESDIAparms(), altered = FESDIAparms(out))

# budgets
  B1 <- FESDIAbudgetO2(out)  
  B1$Fluxmat
  colSums(B1$Fluxmat) - rowSums(B1$Fluxmat)  # Small only for O2
  FESDIAbudgetO2(out, defsteady)$Rates  
  
  B2 <- FESDIAbudgetC(out)
  colSums(B2$Fluxmat) - rowSums(B2$Fluxmat)  # Small for all states
  
  FESDIAbudgetC(out,defsteady)
  FESDIAbudgetC(out,defsteady,defdyn, which = "Rates")

[Package FESDIA version 1.0 Index]