PhaseChange.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 PHASE_CHANGE_H
#define PHASE_CHANGE_H
 
#include <snowpack/DataClasses.h>
#include <meteoio/MeteoIO.h>
 
class PhaseChange {
        public:
                PhaseChange(const SnowpackConfig& i_cfg);
                
                void reset();
                void initialize(SnowStation& Xdata);                                                            //Call before first call to compPhaseChange in a time step
                void finalize(const SurfaceFluxes& Sdata, SnowStation& Xdata, const mio::Date& date_in);        //Call after last call to compPhaseChange in a time step
                double compPhaseChange(SnowStation& Xdata, const mio::Date& date_in, const bool& verbose=true); //Call to do a phase change in a time step, returning the temperature of the top node (K)
 
                static const double RE_theta_r;         
                static const double RE_theta_threshold; 
                static const double theta_r;            
 
        private:
                //To prevent string comparisons, we define an enumerated list:
                enum watertransportmodels{BUCKET, NIED, RICHARDSEQUATION};
                watertransportmodels iwatertransportmodel_snow, iwatertransportmodel_soil;
 
                std::string watertransportmodel_snow;
                std::string watertransportmodel_soil;
                void compSubSurfaceMelt(ElementData& Edata, const unsigned int nSolutes, const double& dt,
                                        double& ql_Rest, const mio::Date& date_in);
                void compSubSurfaceFrze(ElementData& Edata, const unsigned int nSolutes, const double& dt,
                                        const mio::Date& date_in);
 
                const double sn_dt; 
 
                double cold_content_in;         
                double cold_content_soil_in;    
                double cold_content_out;        
                double cold_content_soil_out;   
 
                const bool alpine3d;                    
                double t_crazy_min, t_crazy_max;
                double max_theta_ice;           
 
                static const double theta_s;    
 
                bool enable_ice_reservoir;          
};
 
#endif