R/EstimateNPP.r

In drGeorgeXenakis/fr3PGD: The Physiological Principles for Predicting Growth Model (3PG)

EstimateNPP <-
function (state, parms) 
{
    alpha <- parms[["alpha"]]
    gDM_mol <- parms[["gDM_mol"]]
    Y <- parms[["Y"]]
    phi.pa <- state[["phi.pa"]]
    fT <- state[["fT"]]
    fF <- state[["fF"]]
    fCalpha <- state[["fCalpha"]]
    fN <- state[["fN"]]
    PhysMod <- state[["PhysMod"]]
    alphaC <- alpha * fT * fF * fCalpha * fN * PhysMod
    GPP <- alphaC * gDM_mol * phi.pa/100
    NPP <- Y * GPP
    state[c("alphaC", "GPP", "NPP")] <- c(alphaC, GPP, NPP)
    return(state)
}