aba_decay.h

Go to the documentation of this file.

#ifndef ABA_DECAY_H
#define ABA_DECAY_H
 
#include "../framework/module.h"
#include "../framework/state_map.h"
 
namespace standardBML
{
class aba_decay : public differential_module
{
   public:
    aba_decay(state_map const& input_quantities, state_map* output_quantities)
        : differential_module{},
 
          // Get pointers to input quantities
          soil_aba_concentration_ip{get_ip(input_quantities, "soil_aba_concentration")},
          aba_decay_constant_ip{get_ip(input_quantities, "aba_decay_constant")},
 
          // Get pointers to output quantities
          soil_aba_concentration_op{get_op(output_quantities, "soil_aba_concentration")}
    {
    }
    static string_vector get_inputs();
    static string_vector get_outputs();
    static std::string get_name() { return "aba_decay"; }
 
   private:
    // Pointers to input quantities
    const double* soil_aba_concentration_ip;
    const double* aba_decay_constant_ip;
 
    // Pointers to output quantities
    double* soil_aba_concentration_op;
 
    // Main operation
    void do_operation() const;
};
 
string_vector aba_decay::get_inputs()
{
    return {
        "aba_decay_constant",     //
        "soil_aba_concentration"  //
    };
}
 
string_vector aba_decay::get_outputs()
{
    return {
        "soil_aba_concentration"  //
    };
}
 
void aba_decay::do_operation() const
{
    // Collect input quantities and make calculations
    double aba_decay_constant = *aba_decay_constant_ip;
    double soil_aba_concentration = *soil_aba_concentration_ip;
 
    double decay_rate = -aba_decay_constant * soil_aba_concentration;
 
    // Update the output quantity list
    update(soil_aba_concentration_op, decay_rate);
}
 
}  // namespace standardBML
#endif