ProcessRoutingGR6J.h

#ifndef HYDROBRICKS_PROCESS_ROUTING_GR6J_H
#define HYDROBRICKS_PROCESS_ROUTING_GR6J_H
 
#include "Includes.h"
#include "ProcessOutflow.h"
 
class ProcessRoutingGR6J : public ProcessOutflow {
  public:
    explicit ProcessRoutingGR6J(WaterContainer* container);
 
    ~ProcessRoutingGR6J() override = default;
 
    static void RegisterProcessSettings(SettingsModel* modelSettings);
 
    void SetParameters(const ProcessSettings& processSettings) override;
 
    void Reset() override;
 
    void Finalize() override;
 
    double* GetValuePointer(std::string_view name) override;
 
    [[nodiscard]] vecDouble GetChangeRates() override;
 
  protected:
    const float* _exchangeFactor;     // X2 [mm/d]
    const float* _routingCapacity;    // X3 [mm]
    const float* _uhBaseTime;         // X4 [d]
    const float* _exchangeThreshold;  // X5 [-]
    const float* _expStoreCoeff;      // X6 [mm]
 
    // Unit-hydrograph convolution buffers (length changes with X4)
    vecDouble _stuh1;
    vecDouble _stuh2;
 
    // Pre-computed UH ordinates
    vecDouble _uh1Ord;
    vecDouble _uh2Ord;
 
    // Store levels [mm] — updated in Finalize()
    double _r;     // power routing store
    double _rexp;  // exponential store (may be negative)
 
    // Logged discharge components exposed via GetValuePointer
    double _qr;
    double _qrexp;
    double _qd;
 
    // Cached total process-internal water storage for water balance logging
    double _processStorage;
 
    vecDouble GetRates() override;
 
    void _recomputeUH();
};
 
#endif  // HYDROBRICKS_PROCESS_ROUTING_GR6J_H