Stability.h

Go to the documentation of this file.

/*
 *  SNOWPACK stand-alone
 *
 *  Copyright WSL Institute for Snow and Avalanche Research SLF, DAVOS, SWITZERLAND
*/
/*  This file is part of Snowpack.
    Snowpack is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    Snowpack is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with Snowpack.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef STABILITY_H
#define STABILITY_H

#include <snowpack/StabilityAlgorithms.h>
#include <snowpack/DataClasses.h>

#include <map>
#include <string>

class InstabilityData;

typedef double (*StabMemFn)(const ElementData&, const double&);
typedef bool (*StabFnShearStrength)(const double&, const double&, const mio::Date&,
                                               ElementData&, NodeData&, StabilityData&);

class Stability {
        public:
                Stability (const SnowpackConfig& i_cfg, const bool& i_classify_profile=false);
                
                void checkStability(const CurrentMeteo& Mdata, SnowStation& Xdata);

                static const double psi_ref, max_stability, minimum_slab, ground_rough;
                static const double min_depth_ssi, skier_depth, min_thick_crust;
                static const int sh_mod, prof_classi;
                static const size_t nmax_lemon;

        private:
                void initStability(SnowStation& Xdata);
                static double initStructuralStabilityIndex(const ElementData& Edata_lower, const ElementData& Edata_upper,
                                              const double& Sk, unsigned short &n_lemon);
                
                static void findWeakLayer(const double& Pk, std::vector<unsigned short>& n_lemon, SnowStation& Xdata, double &Swl_ssi, double &Swl_Sk38, size_t &Swl_lemon);
                
                static const bool __init;    
                static bool initStaticData();
                static std::map<std::string, StabMemFn> mapHandHardness;
                static std::map<std::string, StabFnShearStrength> mapShearStrength;

                std::string strength_model, hardness_parameterization;
                double hoar_density_buried;
                bool plastic;
                bool classify_profile, multi_layer_sk38, RTA_ssi;
};

#endif //End of Stability.h