Metamorphism.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 METAMORPHISM_H
#define METAMORPHISM_H
 
#include <snowpack/DataClasses.h>
 
#include <meteoio/MeteoIO.h>
#include <map>
#include <string>
 
#define CALL_MEMBER_FN(object,ptrToMember)  ((object).*(ptrToMember))
 
class SnowStation;
class Metamorphism;
typedef void (Metamorphism::*MetaModelFn)(const CurrentMeteo&, SnowStation&) const;
typedef double (Metamorphism::*MetaSpRateFn)(const ElementData&) const;
 
class Metamorphism {
        public:
                Metamorphism(const SnowpackConfig& i_cfg);
 
                void runMetamorphismModel(const CurrentMeteo& Mdata, SnowStation& Xdata) const noexcept;
 
                static double csPoreArea(const ElementData& Edata);
 
                static double getCoordinationNumberN3(const double& Rho);
 
                static double ddRate(const ElementData& Edata);
 
                static const double mm_tg_dpdz, ba_g_fudge, sa_g_fudge, max_grain_growth, bond_size_stop;
                static const double max_grain_bond_ratio, wind_slab_enhance, wind_slab_vw, wind_slab_depth;
                static const double melt_crust_relaxation;
 
        private:
                static double TGBondRate(const ElementData& Edata);
 
                double LatticeConstant0(const double& th_ice) const;
 
                double TGGrainRate(const ElementData& Edata, const double& Tbot, const double& Ttop,
                                   const double& gradTSub, const double& gradTSup) const;
 
                static double ETBondRate(ElementData& Edata);
                static double ETGrainRate(const ElementData& Edata);
 
                static double PressureSintering(ElementData& Edata);
 
                void metamorphismDEFAULT(const CurrentMeteo& Mdata, SnowStation& Xdata) const;
                void metamorphismNIED(const CurrentMeteo& Mdata, SnowStation& Xdata) const;
 
 
                double spRateDEFAULT(const ElementData& Edata) const;
                double spRateNIED(const ElementData& Edata) const;
 
                static const bool __init;    
                static bool initStaticData();
                static std::map<std::string, MetaModelFn> mapMetamorphismModel;
                static std::map<std::string, MetaSpRateFn> mapSpRate;
 
                const std::string metamorphism_model;
                const double sn_dt, new_snow_grain_size;
                const double melt_crust_relaxation_sp_threshold;
};
 
#endif //End of Metamorphism.h