Program Listing for File kunkel_rh.cpp

↰ Return to documentation for file (src/modules/interp_met/kunkel_rh.cpp)

//
// Canadian Hydrological Model - The Canadian Hydrological Model (CHM) is a novel
// modular unstructured mesh based approach for hydrological modelling
// Copyright (C) 2018 Christopher Marsh
//
// This file is part of Canadian Hydrological Model.
//
// Canadian Hydrological Model 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.
//
// Canadian Hydrological Model 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 Canadian Hydrological Model.  If not, see
// <http://www.gnu.org/licenses/>.
//

#include "kunkel_rh.hpp"
REGISTER_MODULE_CPP(kunkel_rh);

kunkel_rh::kunkel_rh(config_file cfg)
        : module_base("kunkel_rh", parallel::data, cfg)
{
    provides("rh");
    depends_from_met("rh");


    SPDLOG_DEBUG("Successfully instantiated module {}",this->ID);
}

kunkel_rh::~kunkel_rh()
{

}
void kunkel_rh::init(mesh& domain)
{

#pragma omp parallel for
    for (size_t i = 0; i < domain->size_local_faces(); i++)
    {

           auto face = domain->face(i);
           auto& d = face->make_module_data<data>(ID);
           d.interp.init(global_param->interp_algorithm,face->stations().size() );

    }

}
void kunkel_rh::run(mesh_elem &face)
{
    // 1/km
    double lapse_rates[] =
            {-0.09,
             0.0,
             0.09,
             0.11,
             0.11,
             0.12,
             0.14,
             0.15,
             0.11,
             0.07,
             -0.02,
             -0.07
            };

    double lapse = lapse_rates[global_param->month() - 1] / 1000.0; // -> 1/m
    std::vector<boost::tuple<double, double, double> > lowered_values;
    for (auto &s : face->stations())
    {
        if( is_nan((*s)["rh"_s]))
            continue;
        double rh = (*s)["rh"_s];

        double rh_z = rh * exp(lapse * (0.0 - s->z()));

        lowered_values.push_back(boost::make_tuple(s->x(), s->y(), rh_z));
    }


    auto query = boost::make_tuple(face->get_x(), face->get_y(), face->get_z());
    double value = face->get_module_data<data>(ID).interp(lowered_values, query);//C

    double rh = value * exp(lapse * (face->get_z() - 0.0));

    rh = std::min(rh, 100.0);
    rh = std::max(10.0, rh);
    (*face)["rh"_s]= rh;

}