ProcessRoutingHBV.h

#ifndef HYDROBRICKS_PROCESS_ROUTING_HBV_H
#define HYDROBRICKS_PROCESS_ROUTING_HBV_H
 
#include "Includes.h"
#include "ProcessOutflow.h"
 
class ProcessRoutingHBV : public ProcessOutflow {
  public:
    explicit ProcessRoutingHBV(WaterContainer* container);
 
    ~ProcessRoutingHBV() override = default;
 
    static void RegisterProcessSettings(SettingsModel* modelSettings);
 
    [[nodiscard]] bool IsValid() const override;
 
    void SetParameters(const ProcessSettings& processSettings) override;
 
    void Reset() override;
 
    void Finalize() override;
 
    [[nodiscard]] const vecDouble& GetChangeRates() override {
        return GetRates();
    }
 
    double* GetValuePointer(std::string_view name) override;
 
  protected:
    const float* _maxbas;  // base of the triangular weighting function [d]
 
    // Delivery schedule (slot j: water due j timesteps from now) and pre-computed ordinates
    vecDouble _stuh;
    vecDouble _uhOrd;
    double _lastMaxbas;
    // Time step the ordinates were built for. The structure is built before the timer is
    // initialized, so the first computation falls back on a daily step; the ordinates are
    // rebuilt as soon as the real step is known.
    double _lastTimeStep;
 
    // End-of-step container content of the previous timestep (to measure the committed inflow)
    double _previousContent;
 
    // In-transit water (scheduled, not yet delivered) for water balance logging
    double _processStorage;
 
    const vecDouble& GetRates() override;
 
    void _recomputeUH();
 
    static double _cumulativeWeight(double t, double maxbas);
};
 
#endif  // HYDROBRICKS_PROCESS_ROUTING_HBV_H