Glaciers.h

Go to the documentation of this file.

/***********************************************************************************/
/*  Copyright 2009-2015 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 GLACIERS_H
#define GLACIERS_H
 
#include <meteoio/MeteoIO.h>
 
class Glaciers
{
        public:
                Glaciers(const mio::Config& cfg);
                Glaciers(const mio::Config& cfg, const mio::DEMObject& in_dem);
                ~Glaciers() {}
 
                void setDEM(const mio::DEMObject& in_dem);
                void setGlacierMap(const mio::Grid2DObject& glacierMask);
                const mio::Grid2DObject correctTemperatures(const mio::Grid2DObject& hs, const mio::Grid2DObject& tss, const mio::Grid2DObject& ta) const;
                void correctTemperatures(const mio::Grid2DObject& hs, const mio::Grid2DObject& tss, mio::Grid2DObject& ta) const;
                void correctTemperatures(mio::Grid2DObject& ta) const;
 
                void getGrids(mio::Grid2DObject &alt, mio::Grid2DObject &dist) const;
 
        private:
                void init(const mio::Config& cfg);
                static bool enableKatabatikFlows(const mio::Grid2DObject& hs, const mio::Grid2DObject& tss, const mio::Grid2DObject& ta, const mio::Grid2DObject& isGlacier);
                static bool hillslope_distribute_cell(const mio::Grid2DObject& dem, const mio::Grid2DObject& mask, const double& A, const size_t ii, const size_t jj, mio::Grid2DObject &flow, mio::Grid2DObject &src_altitude, mio::Grid2DObject &src_distance);
                void hillslope_flow(mio::Grid2DObject glacier_mask);
 
                mio::DEMObject dem;
                mio::Grid2DObject isGlacier; 
                mio::Grid2DObject flowpath, src_altitude, src_distance;
                double KBL, K, scale;
};
 
#endif