Snowpack.h

#ifndef HYDROBRICKS_SNOWPACK_H
#define HYDROBRICKS_SNOWPACK_H
 
#include <memory>
 
#include "Includes.h"
#include "SnowContainer.h"
#include "SurfaceComponent.h"
 
class Snowpack : public SurfaceComponent {
  public:
    Snowpack();
 
    void Reset() override;
 
    void SaveAsInitialState() override;
 
    void SetParameters(const BrickSettings& brickSettings) override;
 
    void AttachFluxIn(Flux* flux) override;
 
    [[nodiscard]] bool IsValid(bool checkProcesses = true) const override;
 
    [[nodiscard]] WaterContainer* GetSnowContainer() const;
 
    void Finalize() override;
 
    void SetInitialState(double value, ContentType type) override;
 
    [[nodiscard]] double GetContent(ContentType type) const override;
 
    void UpdateContent(double value, ContentType type) override;
 
    void UpdateContentFromInputs(double timeStepInDays = 1.0) override;
 
    void ResetInputBooking() override;
 
    void ApplyConstraints(double timeStep) override;
 
    vecDoublePt GetDynamicContentChanges() override;
 
    double* GetValuePointer(std::string_view name) override;
 
    [[nodiscard]] bool HasSnow() const;
 
    [[nodiscard]] double GetSnowAge() const;
 
  protected:
    // Snowfall rate [mm/d] from which the surface counts as fresh snow (age reset).
    static constexpr double kFreshSnowfallRate = 0.01;
    // Snow content [mm] below which the snowpack counts as empty (age reset).
    static constexpr double kEmptySnowContent = 0.01;
 
    std::unique_ptr<SnowContainer> _snow;  // owning
    double _snowAge = 0;                   // age of the snow surface [d] since the last snowfall
    double _initialSnowAge = 0;            // snow surface age saved with the initial state [d]
    double _snowfallInput = 0;             // snow inflow of the current time step [mm] (for the age reset)
};
 
#endif  // HYDROBRICKS_SNOWPACK_H