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