respiration.cpp

Go to the documentation of this file.

#include "respiration.h"
#include "temperature_response_functions.h"  // for Q10_temperature_response
 
double growth_resp(
    double const base_rate,  // any flux units
    double const grc         // dimensionless
)
{
    // Check for error conditions
    if (grc < 0.0 || grc > 1.0) {
        throw std::range_error("Thrown in growth_resp: grc is outside [0, 1].");
    }
 
    return base_rate < 0 ? 0.0 : base_rate * grc;
}
 
double growth_resp_Q10(
    double const base_rate,  // any flux units
    double const grc0,       // dimensionless
    double const Tleaf,      // degrees C
    double const Tref        // degrees C
)
{
    // Get the growth respiration coefficient at the current temperature
    double const grc{grc0 * Q10_temperature_response(Tleaf, Tref)};  // dimensionless
 
    return growth_resp(base_rate, grc);
}
 
double maintenance_resp_Q10(
    double const tissue_mass,  // Mg / ha
    double const mrc0,         // kg / kg / hr
    double const Tleaf,        // degrees C
    double const Tref          // degrees C
)
{
    // Get the maintenance respiration coefficient at the current temperature
    double const mrc{mrc0 * Q10_temperature_response(Tleaf, Tref)};  // kg / kg / hr
 
    return tissue_mass * mrc;
}