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/ebalance/SolarPanel.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 compute(const mio::Grid2DObject& in_ilwr,
                             const mio::Grid2DObject& in_ta, const mio::Grid2DObject& in_rh, const mio::Grid2DObject& in_p,
                             const mio::Grid2DObject& in_iswr_dir, const mio::Grid2DObject& in_iswr_diff,
                             const mio::Date timestamp);
 
                void writeSP(const unsigned int max_steps); //write Solar Panels information
 
                bool hasSP(){return terrain_radiation->hasSP();} //are Solar Panels defined?
                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_unshaded_horizontal, direct, diffuse, reflected; //FELIX: direct_unshaded_horizontal
                mio::Timer timer;
                const mio::Config& cfg;
                size_t dimx, dimy;
                unsigned int nbworkers;
                bool dataFromGrids;
};
 
#endif