thermal_time_linear.h

Go to the documentation of this file.

#ifndef THERMAL_TIME_LINEAR_H
#define THERMAL_TIME_LINEAR_H
 
#include "../framework/module.h"
#include "../framework/state_map.h"
 
namespace standardBML
{
class thermal_time_linear : public differential_module
{
   public:
    thermal_time_linear(
        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")},
 
          // 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_linear"; }
 
   private:
    // References to input quantities
    double const& fractional_doy;
    double const& sowing_fractional_doy;
    double const& temp;
    double const& tbase;
 
    // Pointers to output quantities
    double* TTc_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector thermal_time_linear::get_inputs()
{
    return {
        "fractional_doy",         // days
        "sowing_fractional_doy",  // days
        "temp",                   // degrees C
        "tbase"                   // degrees C
    };
}
 
string_vector thermal_time_linear::get_outputs()
{
    return {
        "TTc"  // degrees C * day / hr
    };
}
 
void thermal_time_linear::do_operation() const
{
    // Find the rate of change on a daily basis
    double const rate_per_day =
        fractional_doy < sowing_fractional_doy ? 0.0
        : temp <= tbase                        ? 0.0
                                               : temp - tbase;  // degrees C
 
    // Convert to an hourly rate
    double const rate_per_hour = rate_per_day / 24.0;  // degrees C * day / hr
 
    // Update the output quantity list
    update(TTc_op, rate_per_hour);
}
 
}  // namespace standardBML
#endif