TerrainRadiationSimple.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 TERRAINRADIATIONSIMPLE_H
#define TERRAINRADIATIONSIMPLE_H
 
#include <alpine3d/MPIControl.h>
#include <meteoio/MeteoIO.h>
#include <alpine3d/ebalance/TerrainRadiationAlgorithm.h>
#include <alpine3d/ebalance/ViewFactorsAlgorithm.h>
 
class TerrainRadiationSimple : public TerrainRadiationAlgorithm {
 
        public:
                TerrainRadiationSimple(const mio::Config& i_cfg, const mio::DEMObject &dem_in, const std::string& method);
 
                virtual void getRadiation(mio::Array2D<double>& direct, mio::Array2D<double>& diffuse,
                                    mio::Array2D<double>& terrain, const mio::Array2D<double>& direct_unshaded_horizontal,
                                    const mio::Array2D<double>& total_ilwr, mio::Array2D<double>& sky_ilwr,
                                    mio::Array2D<double>& terrain_ilwr, double solarAzimuth, double solarElevation) override;
                virtual void setMeteo(const mio::Array2D<double>& albedo, const mio::Array2D<double>& ta) override;
                void getSkyViewFactor(mio::Array2D<double> &o_sky_vf) override;
 
        private:
                double getAlbedo(const size_t& ii, const size_t& jj);
                void initSkyViewFactor(const mio::DEMObject &dem);
 
                mio::Array2D<double> albedo_grid, sky_vf;
 
                const size_t dimx, dimy;
                size_t startx, endx;
};
 
#endif