thermal_time_beta.h

Go to the documentation of this file.

#ifndef THERMAL_TIME_BETA_H
#define THERMAL_TIME_BETA_H
 
#include <cmath>  // for pow
#include "../framework/module.h"
#include "../framework/state_map.h"
 
namespace standardBML
{
class thermal_time_beta : public differential_module
{
   public:
    thermal_time_beta(
        state_map const& input_quantities,
        state_map* output_quantities)
        : differential_module{},
 
          // Get references to input quantities
          fractional_doy{get_input(input_quantities, "fractional_doy")},
          sowing_fractional_doy{get_input(input_quantities, "sowing_fractional_doy")},
          temp{get_input(input_quantities, "temp")},
          tbase{get_input(input_quantities, "tbase")},
          tmax{get_input(input_quantities, "tmax")},
          talpha{get_input(input_quantities, "talpha")},
          tbeta{get_input(input_quantities, "tbeta")},
          ttc_scale{get_input(input_quantities, "ttc_scale")},
 
          // Get pointers to output quantities
          TTc_op{get_op(output_quantities, "TTc")}
    {
    }
    static string_vector get_inputs();
    static string_vector get_outputs();
    static std::string get_name() { return "thermal_time_beta"; }
 
   private:
    // References to input quantities
    double const& fractional_doy;
    double const& sowing_fractional_doy;
    double const& temp;
    double const& tbase;
    double const& tmax;
    double const& talpha;
    double const& tbeta;
    double const& ttc_scale;
 
    // Pointers to output quantities
    double* TTc_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector thermal_time_beta::get_inputs()
{
    return {
        "fractional_doy",         // days
        "sowing_fractional_doy",  // days
        "temp",                   // degrees C
        "tbase",                  // degrees C
        "tmax",                   // degrees C
        "talpha",                 // dimensionless
        "tbeta",                  // dimensionless
        "ttc_scale"               // degrees C * day / hr
    };
}
 
string_vector thermal_time_beta::get_outputs()
{
    return {
        "TTc"  // degrees C * day / hr
    };
}
 
void thermal_time_beta::do_operation() const
{
    // Find the rate of change
    double const rate =
        fractional_doy < sowing_fractional_doy ? 0.0
        : temp <= tbase                        ? 0.0
        : temp <= tmax                         ? ttc_scale * pow(temp - tbase, talpha) * pow(tmax - temp, tbeta)
                                               : 0.0;  // degrees C * day / hr
 
    // Update the output quantity list
    update(TTc_op, rate);
}
 
}  // namespace standardBML
#endif