TechSnowA3D.h

Go to the documentation of this file.

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

    Alpine3D 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 Lesser General Public License for more details.

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

class SnowpackInterfaceWorker;

#include <meteoio/MeteoIO.h>
#include <snowpack/libsnowpack.h>
#include <alpine3d/MeteoObj.h>

class TechSnowA3D
{
        public:
                TechSnowA3D(const mio::Config& cfg, const mio::DEMObject& dem);

                void setMeteo(const mio::Grid2DObject& ta,
                              const mio::Grid2DObject& rh,
                              const mio::Grid2DObject& hs,
                              const mio::Date& timestamp);
                
                mio::Grid2DObject getGrid(const SnGrids::Parameters& param) const;
                
                static std::string getGridsRequirements() { return "TA RH HS"; }

        private:
                typedef struct CONDITION {
                        CONDITION() 
                           : slope_area(0.), number_snowguns(0), priority(mio::IOUtils::unodata), min_height(0.)  {}
                        
                        void setUserProperties(const unsigned int& nr_snowguns, const unsigned int& i_priority, const double& i_min_height, const double& wb_thresh) {
                                number_snowguns = nr_snowguns;
                                priority = i_priority;
                                min_height = i_min_height;
                                wet_bulb_thresh = wb_thresh;
                        }
                        
                        const std::string toString() {std::ostringstream os; os << "[" << slope_area << "m², " << number_snowguns << " snowguns, " << priority << " priority, " << min_height << "m min, min wbt: " << wet_bulb_thresh << "°C]"; return os.str();}
                        
                        double slope_area;              //area of the slope section
                        unsigned int number_snowguns;   //number of snow guns per slope section
                        unsigned int priority;  //snow production priority for this slope
                        double min_height;              //minimum snow height for technical snow production
                        double wet_bulb_thresh; //wet bulb temperature threshold (°C) below which snow can be produced
                } condition;
                
                typedef struct SNOWSTRATEGY {
                        SNOWSTRATEGY() : startProd(), endProd(), startAim(0.), endAim(0.), gunOperation(0.) {}
                        const std::string toString() {std::ostringstream os; os << "[" << startProd.toString(mio::Date::ISO) << "-" << endProd.toString(mio::Date::ISO) << " aim: " << startAim << " -> " << endAim << " @ " << gunOperation << "]"; return os.str();}
                        
                        mio::Date startProd;    //start date of snow production
                        mio::Date endProd;              //end date of snow production
                        double startAim;                //snow production aim in snow height [m] until ski resort opening (factor x snow_prod_min)
                        double endAim;                  //snow production aim in snow height [m] between ski resort opening and end date of snow production (factor x snow_prod_min)
                        double gunOperation;    //percentage of snow gun operation
                } snowStrategy;
                
                static TechSnowA3D::snowStrategy setSnowStrategy(const mio::Config& cfg, const double& TZ, const unsigned int& nr);
                void setSlopeConditions(const std::string& filename);
                static size_t getSlopeNumber(const double& dbl_code);
                double setPriority(const mio::Date& date, const TechSnowA3D::snowStrategy &ppt, const double& snow_height, const TechSnowA3D::condition& slope, const int date_hour) const;
                
                mio::Grid2DObject skiRunsMap; 
                mio::Grid2DObject psum_tech;
                
                std::map<size_t, condition> slope_conditions;
                mio::Date start_season, end_season;             
                mio::Date earliest_production;                  
                std::vector<snowStrategy> priorities;
                
                double max_snowgun_water;       
                const double mean_elevation;
                int slope_open, slope_closed; 
};

#endif