varying_Jmax25.h

Go to the documentation of this file.

#ifndef VARYING_JMAX25_H
#define VARYING_JMAX25_H
 
#include "../framework/module.h"
#include "../framework/state_map.h"
 
namespace standardBML
{
class varying_Jmax25 : public direct_module
{
   public:
    varying_Jmax25(
        state_map const& input_quantities,
        state_map* output_quantities)
        : direct_module{},
 
          // Get references to input quantities
          DVI{get_input(input_quantities, "DVI")},
          Jmax_at_25_mature{get_input(input_quantities, "Jmax_at_25_mature")},
          DVI0_jmax{get_input(input_quantities, "DVI0_jmax")},
          D_DVI{get_input(input_quantities, "Delta_DVI_jmax")},
          sf_jmax{get_input(input_quantities, "sf_jmax")},
 
          // Get pointers to output quantities
          Jmax_at_25_op{get_op(output_quantities, "Jmax_at_25")}
    {
    }
    static string_vector get_inputs();
    static string_vector get_outputs();
    static std::string get_name() { return "varying_Jmax25"; }
 
   private:
    // Pointers to input quantities
    double const& DVI;
    double const& Jmax_at_25_mature;
    double const& DVI0_jmax;
    double const& D_DVI;
    double const& sf_jmax;
 
    // Pointers to output quantities
    double* Jmax_at_25_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector varying_Jmax25::get_inputs()
{
    return {
        "DVI",               // dimensionless
        "sf_jmax",           // dimensionless
        "DVI0_jmax",         // dimensionless
        "Delta_DVI_jmax",    // dimensionless
        "Jmax_at_25_mature"  // micromol / m^2 / s
    };
}
 
string_vector varying_Jmax25::get_outputs()
{
    return {
        "Jmax_at_25"  // micromol / m^2 / s
    };
}
 
void varying_Jmax25::do_operation() const
{
    double const k = -log(1.0 / 0.95 - 1.0);  // dimensionless
 
    double const C = Jmax_at_25_mature * sf_jmax;        // micromol / m^2 / s
    double const L = Jmax_at_25_mature * (1 - sf_jmax);  // micromol / m^2 / s
 
    double const Jmax_at_25 = L / (1.0 + exp(-k * (DVI - DVI0_jmax) / D_DVI)) + C;  // micromol / m^2 / s
 
    update(Jmax_at_25_op, Jmax_at_25);
}
 
}  // namespace standardBML
#endif