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);
 
                void compMeteo(CurrentMeteo &Mdata, SnowStation &Xdata, const bool runCanopyModel);
                void compMeteo(CurrentMeteo &Mdata, SnowStation &Xdata, const bool runCanopyModel, const bool i_adjust_height_of_wind_value);
                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 double& roughness_length, 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;
                const double roughness_length_dflt; 
                double height_of_wind_value; 
                ATM_STABILITY stability; 
                bool research_mode; 
                bool useCanopyModel; 
                bool adjust_height_of_wind_value; 
};
 
#endif //END of Meteo.h