ProcessRoutingGR4J.h

#ifndef HYDROBRICKS_PROCESS_ROUTING_GR4J_H
#define HYDROBRICKS_PROCESS_ROUTING_GR4J_H
 
#include "Includes.h"
#include "ProcessOutflow.h"
 
class ProcessRoutingGR4J : public ProcessOutflow {
  public:
    explicit ProcessRoutingGR4J(WaterContainer* container);
 
    ~ProcessRoutingGR4J() override = default;
 
    static void RegisterProcessSettings(SettingsModel* modelSettings);
 
    [[nodiscard]] bool RequiresDailyTimeStep() const override {
        return true;
    }
 
    void SetParameters(const ProcessSettings& processSettings) override;
 
    void Reset() override;
 
    void Finalize() override;
 
    double* GetValuePointer(std::string_view name) override;
 
  protected:
    const float* _exchangeFactor;   // X2 [mm/d]
    const float* _routingCapacity;  // X3 [mm]
    const float* _uhBaseTime;       // X4 [d]
 
    // Unit-hydrograph convolution buffers (length changes with X4)
    vecDouble _stuh1;
    vecDouble _stuh2;
 
    // Pre-computed UH ordinates
    vecDouble _uh1Ord;
    vecDouble _uh2Ord;
 
    // Routing store level [mm] — updated in Finalize()
    double _r;
 
    // Logged discharge components exposed via GetValuePointer
    double _qr;
    double _qd;
 
    // Cached total process-internal water storage for water balance logging
    double _processStorage;
 
    const vecDouble& GetRates() override;
 
    void _recomputeUH();
 
    static double _sh1(double t, double x4);
    static double _sh2(double t, double x4);
};
 
#endif  // HYDROBRICKS_PROCESS_ROUTING_GR4J_H