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();
 
                void editMeteoData(mio::MeteoData& md) const;
                bool validMeteoData(const mio::MeteoData& md, const std::string& StationName, const unsigned int& nslopes) const;
                void copyMeteoData(const mio::MeteoData& md, CurrentMeteo& Mdata,
                                   const double& prevailing_wind_dir, bool &iswr_is_net) const;
                static double getHS_last3hours(mio::IOManager &io, const mio::Date& current_date);
                
                static void setShortWave(CurrentMeteo& Mdata, const SnowStation& Xdata, const bool& iswr_is_net);
                void dataForCurrentTimeStep(CurrentMeteo& Mdata, SurfaceFluxes& surfFluxes, std::vector<SnowStation>& vecXdata,
                                                                        const Slope& slope, SnowpackConfig& cfg,
                                                                        mio::SunObject &sun,
                                                                        double& precip, const double& lw_in, const double hs_a3hl6,
                                                                        double& tot_mass_in,
                                                                        const bool& iswr_is_net);
                
                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; 
                double roughness_length; 
                double height_of_wind_value; 
                double thresh_rain; 
                double wind_scaling_factor; 
                std::string variant; 
                std::string sw_mode; 
                ATM_STABILITY stability; 
                bool research_mode; 
                bool useCanopyModel; 
                bool enforce_snow_height; 
                bool advective_heat; 
                bool soil_flux; 
                bool force_sw_mode; 
};
 
#endif //END of Meteo.h