thermal_time_bilinear.h

Go to the documentation of this file.

#ifndef THERMAL_TIME_BILINEAR_H
#define THERMAL_TIME_BILINEAR_H
 
#include "../framework/module.h"
#include "../framework/state_map.h"
 
namespace standardBML
{
class thermal_time_bilinear : public differential_module
{
   public:
    thermal_time_bilinear(
        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")},
          topt{get_input(input_quantities, "topt")},
          tmax{get_input(input_quantities, "tmax")},
 
          // 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_bilinear"; }
 
   private:
    // References to input quantities
    double const& fractional_doy;
    double const& sowing_fractional_doy;
    double const& temp;
    double const& tbase;
    double const& topt;
    double const& tmax;
 
    // Pointers to output quantities
    double* TTc_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector thermal_time_bilinear::get_inputs()
{
    return {
        "fractional_doy",         // days
        "sowing_fractional_doy",  // days
        "temp",                   // degrees C
        "tbase",                  // degrees C
        "topt",                   // degrees C
        "tmax"                    // degrees C
    };
}
 
string_vector thermal_time_bilinear::get_outputs()
{
    return {
        "TTc"  // degrees C * day / hr
    };
}
 
void thermal_time_bilinear::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 <= topt                         ? temp - tbase
        : temp <= tmax                         ? (tmax - temp) * (topt - tbase) / (tmax - topt)
                                               : 0.0;  // 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