permafrost.h

Go to the documentation of this file.

#ifndef PERMAFROST_H
#define PERMAFROST_H
 
#include "layeredgrid.h"
 
 
class WaterCycle; // forward
struct ClimateDay; // forward
class ResourceUnit; // forward
class WaterOut; // forward
 
namespace Water {
 
class Permafrost
{
public:
    Permafrost();
    ~Permafrost();
    void setup(WaterCycle *wc);
 
    //const SStats &stats() const { return stats; }
 
    void newYear();
 
    void run(const ClimateDay *clim_day);
 
    void debugData(DebugList &out);
 
    double mossLayerThickness() const { return mMossBiomass / mosspar.bulk_density; } // kg/m2  / rho [kg/m3] = m
    double SOLLayerThickness() const {return mSOLDepth; }
    double mossBiomass() const { return mMossBiomass; }
    void burnMoss(const double biomass_kg) { mMossBiomass = mMossBiomass - biomass_kg/cRUArea;
                                           mMossBiomass = std::max(mMossBiomass, cMinMossBiomass); }
 
private:
    const double cMinMossBiomass = 0.0001; // kg/m2
    void setupMossLayer();
    void calculateMoss();
 
 
    void setupThermalConductivity();
 
    double thermalConductivity(bool from_below) const;
    double thermalConductivityFrozen() const;
 
    struct FTResult {
        FTResult(): delta_mm(0.), delta_soil(0.), new_depth(0.), orig_depth(0.) {}
        double delta_mm; // change of water (mm within iLand water bucket), freezing is negative
        double delta_soil; // change of ice layer (m) (within iLand water bucket), freezing is negative
        double new_depth; // final depth (m)
        double orig_depth;  // starting depth (m)
    };
    FTResult calcFreezeThaw(double at, double temp, bool lowerIceEdge, bool fromAbove);
    FTResult mResult; 
 
    WaterCycle *mWC;
    double mCurrentSoilFrozen;
    double mCurrentWaterFrozen;
    double mSoilDepth;
    double mPWP; 
    double mFC; 
    double mTop;
    bool mTopFrozen; 
    double mBottom;
    double mFreezeBack; 
 
    double mSOLDepth; 
 
    double mGroundBaseTemperature; 
 
    double mKdry;
    double mKsat;
    double mKice;
    bool mSoilIsCoarse; 
 
    // moss related variables
    double mMossBiomass; 
    // stats
    struct SStats {
        void reset() { maxSnowDepth=0.; daysSnowCover=0; maxFreezeDepth=0.; maxThawDepth=0.; }
        double maxSnowDepth; 
        int daysSnowCover; 
        double maxFreezeDepth; 
        double maxThawDepth; 
    };
    SStats stats;
 
    // parameters of the permafrost model are stored in a static struct
    // (and exist only once)
    struct SParam {
        void init() { EFusion = 0.333; // MJ/litre H20
                      maxPermafrostDepth = 2.; // 2m
                 }
        // physical constants
        double EFusion; 
        double maxPermafrostDepth; 
 
        double lambdaSnow; 
        double lambdaOrganicLayer; 
 
        double organicLayerDensity; 
        double SOLDefaultDepth = 0.1; 
 
        // parameters
        double groundBaseDepth; 
        double maxFreezeThawPerDay; 
 
        bool onlySimulate; 
    };
 
    struct SMossParam {
        const double SLA=1.; 
        const double AMax=0.3; 
        double light_k; 
        double light_comp; 
        double light_sat; 
        double respiration_q, respiration_b; 
        double bulk_density; 
        // carbon cycle
        double CNRatio; 
        double r_decomp; 
    };
 
    static SParam par;
    static SMossParam mosspar;
 
    friend class PermafrostLayers;
    friend class ::WaterOut;
};
 
class PermafrostLayers: public LayeredGrid<ResourceUnit*> {
  public:
    ~PermafrostLayers() {}
    void setGrid(const Grid<ResourceUnit*> &grid) { mGrid = &grid; }
    void clearGrid() { mGrid = nullptr; }
    double value(ResourceUnit* const &data, const int index) const;
    const QVector<LayeredGridBase::LayerElement> &names();
    //bool onClick(const QPointF &world_coord) const;
private:
    QVector<LayeredGridBase::LayerElement> mNames;
};
 
 
} // end namespace
 
#endif // PERMAFROST_H