Meteo.h

Go to the documentation of this file.

/*
 *  SNOWPACK stand-alone
 *
 *  Copyright WSL Institute for Snow and Avalanche Research SLF, DAVOS, SWITZERLAND
*/
/*  This file is part of Snowpack.
    Snowpack is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    Snowpack 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 General Public License for more details.
 
    You should have received a copy of the GNU General Public License
    along with Snowpack.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef METEO_H
#define METEO_H
 
#include <meteoio/MeteoIO.h>
 
#include <snowpack/SnowpackConfig.h>
#include <snowpack/snowpackCore/Canopy.h>
#include <snowpack/DataClasses.h>
 
class Meteo {
        public:
                //except Richardson and Neutral, all are standard MO iterations with the specified stability correction. In unstable conditions, they use Paulson and Stearns & Weidner, 1993, for scalars
                typedef enum {
                        RICHARDSON,  
                        NEUTRAL,  
                        MO_LOG_LINEAR, 
                        MO_HOLTSLAG, 
                        MO_STEARNS, 
                        MO_MICHLMAYR, 
                        MO_SCHLOEGL_UNI, 
                        MO_SCHLOEGL_MULTI, 
                        MO_SCHLOEGL_MULTI_OFFSET 
                } ATM_STABILITY;
 
                Meteo(const SnowpackConfig& i_cfg);
                Meteo(const Meteo& mt);
                Meteo& operator=(const Meteo& mt); 
                virtual ~Meteo();
 
                static void projectPrecipitations(const double& SlopeAngle, double& precips, double& hs);
                static bool compHSrate(CurrentMeteo& Mdata, const SnowStation& vecXdata, const double& hs_a3hl6);
                void compMeteo(CurrentMeteo &Mdata, SnowStation &Xdata, const bool runCanopyModel,
                               const bool adjust_height_of_wind_value);
                void compRadiation(const SnowStation &station, mio::SunObject &sun, SnowpackConfig &cfg, CurrentMeteo &Mdata);
                static void radiationOnSlope(const SnowStation &sector, const mio::SunObject &sun, CurrentMeteo &Mdata, SurfaceFluxes &surfFluxes);
                void setStability(const ATM_STABILITY& i_stability);
                static ATM_STABILITY getStability(const std::string& stability_model);
                ATM_STABILITY getStability() const;
 
        private:
                void MicroMet(const SnowStation& Xdata, CurrentMeteo& Mdata, const bool& adjust_VW_height=true) const;
                static double getParameterAverage(mio::IOManager& io, const mio::MeteoData::Parameters& param,
                                                  const mio::Date& current_date, const int& time_span, const int& increment);
                static void RichardsonStability(const double& ta_v, const double& t_surf_v, const double& zref,
                                                const double& vw, const double& z_ratio, double &ustar, double &psi_s);
                static void MOStability(const ATM_STABILITY& use_stability, const double& ta_v, const double& t_surf_v, const double& t_surf,
                                                       const double& zref, const double& vw, const double& z_ratio, double &ustar, double &psi_s, double &psi_m);
 
                Canopy canopy;
                mio::DataGenerator *dataGenerator; //as pointer so we don't have to construct an expensive object if not needed
                double roughness_length, height_of_wind_value;
                std::string variant;
                ATM_STABILITY stability;
                bool research_mode, useCanopyModel;
};
 
#endif //END of Meteo.h