# 199: Implementation of a 6-pool Michaelis-Menten model

Description Usage Arguments Value Author(s) Examples

### Description

This function implements a 6-pool Michaelis-Meneten model with pairs of microbial biomass and substrate pools.

### Usage

 `1` ```ThreepairMMmodel(t, ks, kb, Km, r, Af = 1, ADD, ival) ```

### Arguments

 `t` vector of times to calculate a solution. `ks` a vector of length 3 representing SOM decomposition rate (m3 d-1 (gCB)-1) `kb` a vector of length 3 representing microbial decay rate (d-1) `Km` a vector of length 3 representing the Michaelis constant (g m-3) `r` a vector of length 3 representing the respired carbon fraction (unitless) `Af` a scalar representing the Activity factor; i.e. a temperature and moisture modifier (unitless) `ADD` a vector of length 3 representing the annual C input to the soil (g m-3 d-1) `ival`

### Value

An object of class NlModel that can be further queried.

### Author(s)

Carlos A. Sierra, Markus Mueller

### Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34``` ``` days=seq(0,1000) #Run the model with default parameter values MMmodel=ThreepairMMmodel(t=days,ival=rep(c(100,10),3),ks=c(0.1,0.05,0.01), kb=c(0.005,0.001,0.0005),Km=c(100,150,200),r=c(0.9,0.9,0.9), ADD=c(3,1,0.5)) # fixme mm: # the next line causes trouble on Rforge Windows patched build # Cpools=getC(MMmodel) # # #Time solution # matplot(days,Cpools,type="l",ylab="Concentrations",xlab="Days",lty=rep(1:2,3), # ylim=c(0,max(Cpools)*1.2),col=rep(1:3,each=2), # main="Multi-substrate microbial model") # legend("topright",c("Substrate 1", "Microbial biomass 1", # "Substrate 2", "Microbial biomass 2", # "Substrate 3", "Microbial biomass 3"), # lty=rep(1:2,3),col=rep(1:3,each=2), # bty="n") # # # #State-space diagram # plot(Cpools[,2],Cpools[,1],type="l",ylab="Substrate",xlab="Microbial biomass") # lines(Cpools[,4],Cpools[,3],col=2) # lines(Cpools[,6],Cpools[,5],col=3) # legend("topright",c("Substrate-Enzyme pair 1","Substrate-Enzyme pair 2", # "Substrate-Enzyme pair 3"),col=1:3,lty=1,bty="n") # # #Microbial biomass over time # plot(days,Cpools[,2],type="l",col=2,xlab="Days",ylab="Microbial biomass") # # ```

Search within the SoilR package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.