LIMCoralRockall: Linear inverse model specification for the Deep-water Coral...

In LIM: Linear Inverse Model Examples and Solution Methods

Linear inverse model specification for the Deep-water Coral food web at Rockall Bank

Description

Linear inverse model specification for the deep-water coral ecosystem at Rockall Bank, North-East Atlantic. See van Oevelen et al. (2009)

Units of the flows are mmol C/m2/day

The linear inverse model LIMCoralRockall is generated from the file ‘CWCRockall.input’ which can be found in subdirectory /examples/FoodWeb of the package directory

Usage

data(LIMCoralRockall)

Format

a list of matrices, vectors, names and values that specify the linear inverse model problem.

see the return value of Setup for more information about this list

A more complete description of this structure is in vignette("LIM")

Author(s)

Dick van Oevelen <[e-mail]>

Karline Soetaert <[e-mail]>

References

van Oevelen, Dick, Gerard Duineveld, Marc Lavaleye, Furu Mienis, Karline Soetaert, and Carlo H. R. Heip, 2009.

The cold-water coral community as hotspot of carbon cycling on continental margins: A food web analysis from Rockall Bank (northeast Atlantic). Limnology and Oceangraphy 54 : 1829 – 1844.

http://www.aslo.org/lo/toc/vol_54/issue_6/1829.pdf

See Also

browseURL(paste(system.file(package="LIM"), "/examples/Foodweb/", sep=""))

contains "CWCRockall.input", the input file; read this with Setup

LIMTakapoto, LIMRigaSummer and many others

Examples

Coral <- Flowmatrix(LIMCoralRockall)

plotweb(Coral, main = "Deep Water Coral Foodweb, Rockall Bank",
        sub = "mmolC/m2/day", lab.size = 0.8)


## Not run: 
xr <- LIMCoralRockall$NUnknowns
i1 <- 1:(xr/2)
i2 <- (xr/2+1):xr

pm <- par(mfrow = c(1, 1))
Simplest <- Ldei(LIMCoralRockall)$X
Ranges <- Xranges(LIMCoralRockall)
Plotranges(Ranges[i1, 1], Ranges[i1, 2], Simplest[i1], lab.cex = 0.7,
        main = "Deep Water Coral - ranges part 1")

Plotranges(Ranges[i2, 1], Ranges[i2, 2], Simplest[i2], lab.cex = 0.7,
        main = "Deep Water Coral - ranges part 2")
par(mfrow = pm)

## End(Not run)

LIM documentation built on May 11, 2022, 5:05 p.m.