EnergyBalance.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 ENERGYBALANCE_H
#define ENERGYBALANCE_H

#include <string>
#include <meteoio/MeteoIO.h>

class SnowpackInterface;

#include <alpine3d/ebalance/RadiationField.h>
#include <alpine3d/ebalance/TerrainRadiationAlgorithm.h>
#include <alpine3d/SnowpackInterface.h>

class EnergyBalance
{
        public:
                EnergyBalance( const unsigned int& i_nbworkers, const mio::Config& cfg, const mio::DEMObject &dem_in);
                EnergyBalance(const EnergyBalance&);
                ~EnergyBalance( );

                EnergyBalance& operator=(const EnergyBalance&); 

                void setSnowPack(SnowpackInterface &mysnowpack );
                void setAlbedo( const mio::Grid2DObject &in_albedo );

                void setMeteo(const mio::Grid2DObject& in_ilwr,
                              const mio::Grid2DObject& in_ta, const mio::Grid2DObject& in_rh, const mio::Grid2DObject& in_p, const mio::Date timestamp);

                void setStations(const std::vector<mio::MeteoData>& in_vecMeteo);
                double getTiming() const;
                void Destroy();
                std::string getGridsRequirements() const;

        private:
                SnowpackInterface *snowpack;
                TerrainRadiationAlgorithm *terrain_radiation;
                std::vector<RadiationField*> radfields;
                mio::DEMObject dem;
                std::vector<mio::MeteoData> vecMeteo;
                mio::Grid2DObject albedo;
                mio::Array2D<double> direct, diffuse, reflected;
                mio::Timer timer;

                size_t dimx, dimy;
                unsigned int nbworkers;
};

#endif