MeteoObj.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 METEOOBJ_H
#define METEOOBJ_H
 
#include <meteoio/MeteoIO.h>
#include <snowpack/libsnowpack.h>
 
#include <alpine3d/MPIControl.h>
#include <alpine3d/Glaciers.h>
 
#include <iostream>
 
class SnGrids {
        public:
                //Make sure you also put the parameters in the same order in SnGrids::initStaticData()!!
                enum Parameters {firstparam=0,
                                          TA=firstparam, 
                                          RH, 
                                          VW, 
                                          DW, 
                                          ISWR, 
                                          ISWR_DIFF, 
                                          ISWR_DIR, 
                                          ILWR, 
                                          HS, 
                                          PSUM, 
                                          PSUM_PH, 
                                          PSUM_TECH, 
                                          GROOMING, 
                                          TSG, 
                                          TSS, 
                                          TS0, 
                                          TSNOW, 
                                          TSNOW_AVG, 
                                          RHOSNOW_AVG, 
                                          SWE, 
                                          RSNO, 
                                          TOP_ALB, 
                                          SURF_ALB, 
                                          SP, 
                                          RB, 
                                          RG, 
                                          N3, 
                                          MS_SNOWPACK_RUNOFF, 
                                          MS_SURFACE_MASS_FLUX, 
                                          MS_SOIL_RUNOFF, 
                                          MS_WATER, 
                                          MS_WATER_SOIL, 
                                          MS_ICE_SOIL, 
                                          SFC_SUBL, 
                                          MNS, 
                                          STORE, 
                                          GLACIER, 
                                          GLACIER_EXPOSED, 
                                          ET, 
                                          ISWR_TERRAIN, 
                                          ILWR_TERRAIN, 
                                          ISWR_BELOW_CAN,
                                          TSOIL1, TSOIL2, TSOIL3, TSOIL4, TSOIL_MAX, 
                                          SOIL_RUNOFF1, SOIL_RUNOFF2, SOIL_RUNOFF3, SOIL_RUNOFF4, SOIL_RUNOFF_MAX,
                                          lastparam=SOIL_RUNOFF_MAX};
 
 
                static const size_t nrOfParameters; 
                static const std::string& getParameterName(const size_t& parindex);
                static size_t getParameterIndex(const std::string& parname);
 
        private:
                static std::vector<std::string> paramname;
                static const bool __init;    
                static bool initStaticData();
};
 
class MeteoObj
{
        public:
                MeteoObj(const mio::Config& config, const mio::DEMObject& in_dem);
                ~MeteoObj();
 
                void setSkipWind(const bool& i_skipWind);
                void prepare(const mio::Date& in_date);
                void get(const mio::Date& in_date,
                         mio::Grid2DObject& ta,
                         mio::Grid2DObject& rh,
                         mio::Grid2DObject& psum,
                         mio::Grid2DObject& psum_ph,
                         mio::Grid2DObject& vw,
                         mio::Grid2DObject& dw,
                         mio::Grid2DObject& p,
                         mio::Grid2DObject& ilwr,
                         mio::Grid2DObject& iswr_dir,
                         mio::Grid2DObject& iswr_diff);
                void get(const mio::Date& in_date, std::vector<mio::MeteoData>& o_vecMeteo);
                bool fillPrecSplitting();
                void checkMeteoForcing(const mio::Date& calcDate);
                void setGlacierMask(const mio::Grid2DObject& glacierMask);
                void setDEM(const mio::DEMObject& in_dem);
                double getTiming() const;
 
        private:
                static void checkLapseRate(const std::vector<mio::MeteoData>& i_vecMeteo, const mio::MeteoData::Parameters& param);
                static void checkGridRange(const mio::Date& calcDate, const mio::Grid2DObject& grid, const mio::MeteoData::Parameters& param);
                static void checkInputsRequirements(std::vector<mio::MeteoData>& vecData);
                void fillMeteoGrids(const mio::Date& calcDate);
                void getMeteo(const mio::Date& calcDate);
 
                mio::Timer timer;
                const mio::Config &config;
                mio::IOManager io;
                mio::DEMObject dem;
                mio::Grid2DObject ta, rh, psum, psum_ph, vw, dw, p, ilwr, iswr_dir, iswr_diff;
                mio::Grid2DObject sum_ta, sum_rh, sum_rh_psum, sum_psum, sum_psum_ph, sum_vw, sum_ilwr;
                std::vector<mio::MeteoData> vecMeteo;
                mio::Date date;
                Glaciers *glaciers;
                unsigned int count_sums, count_precip;
                bool skipWind; 
                bool dataFromGrids; 
};
 
#endif