ProcessLateralSnowRedistributionFrey.h

Go to the documentation of this file.

 
#ifndef HYDROBRICKS_PROCESS_LATERAL_SNOW_REDISTRIBUTION_FREY_H
#define HYDROBRICKS_PROCESS_LATERAL_SNOW_REDISTRIBUTION_FREY_H
 
#include "ProcessLateral.h"
 
class FluxToBrick;
class Forcing;
 
class ProcessLateralSnowRedistributionFrey : public ProcessLateral {
  public:
    explicit ProcessLateralSnowRedistributionFrey(WaterContainer* container, bool dynamicDensity = false);
 
    ~ProcessLateralSnowRedistributionFrey() override = default;
 
    [[nodiscard]] bool IsValid() const override;
 
    static void RegisterProcessSettings(SettingsModel* modelSettings);
 
    static void RegisterProcessSettingsDynamic(SettingsModel* modelSettings);
 
    void SetHydroUnitProperties(HydroUnit* unit, Brick* brick) override;
 
    void SetParameters(const ProcessSettings& processSettings) override;
 
    void AttachForcing(Forcing* forcing) override;
 
    void Reset() override;
 
    void Finalize() override;
 
    [[nodiscard]] bool IsLateralProcess() const noexcept override {
        return true;
    }
 
  protected:
    bool _dynamicDensity;           // Whether the snow density evolves dynamically
    float _slopeDeg;                // Slope of the hydro unit [°]
    Forcing* _temperature;          // Air temperature forcing (dynamic variant only)
    const float* _correction;       // Correction coefficient C []
    const float* _holdingCapacity;  // Land-cover snow holding capacity H_v (as a snow depth) [mm]
    const float* _rhoMax;           // Maximum snow density for redistribution [kg/m3]
    const float* _maxSnowDepth;     // Maximum snow depth (over-accumulation guard) [mm]
    const float* _snowDensity;      // Constant snow density (constant variant) [kg/m3]
    const float* _rhoMin;           // Minimum snow density (dynamic variant) [kg/m3]
    const float* _rhoFreshMax;      // Maximum fresh-snow density (dynamic variant) [kg/m3]
    const float* _rhoSettling;      // Settling rate toward rhoMax (dynamic variant) [1/day]
    const float* _rhoScale;         // Density sigmoid slope coefficient (dynamic variant) []
    const float* _tScale;           // Density sigmoid shift coefficient (dynamic variant) [°C]
    double _density;                // Tracked bulk snow density (dynamic variant) [kg/m3]
    double _prevSwe;                // Committed SWE at the previous time step (dynamic variant) [mm]
 
    const vecDouble& GetRates() override;
 
  private:
    [[nodiscard]] double GetSnowDensity() const;
 
    [[nodiscard]] double ComputeFreshSnowDensity() const;
 
    double AvoidUnrealisticAccumulation(double rate, Flux* flux);
};
 
#endif  // HYDROBRICKS_PROCESS_LATERAL_SNOW_REDISTRIBUTION_FREY_H