.. _program_listing_file_src_modules_slope_iswr.cpp:

Program Listing for File slope_iswr.cpp
=======================================

|exhale_lsh| :ref:`Return to documentation for file <file_src_modules_slope_iswr.cpp>` (``src/modules/slope_iswr.cpp``)

.. |exhale_lsh| unicode:: U+021B0 .. UPWARDS ARROW WITH TIP LEFTWARDS

.. code-block:: 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 "slope_iswr.hpp"
   REGISTER_MODULE_CPP(slope_iswr);
   
   slope_iswr::slope_iswr(config_file cfg)
           : module_base("slope_iswr", parallel::data, cfg)
   {
       depends("iswr");
       depends("iswr_diffuse");
       depends("iswr_direct");
   
       depends("solar_el");
       depends("solar_az");
   
   
       provides("iswr");
       provides("iswr_direct");
       provides("solar_angle");
   
       optional("shadow");
   
       SPDLOG_DEBUG("Successfully instantiated module {}",this->ID);
   
       assume_no_slope = cfg.get("no_slope",false);
   }
   void slope_iswr::run(mesh_elem& face)
   {
   
   
       double A = (*face)["solar_az"_s] * mio::Cst::to_rad;
       double E = (*face)["solar_el"_s] * mio::Cst::to_rad;
   
       //radiation data
       //solar vector
       //xyz cartesian
       arma::vec S;
       S << cos(E) * sin(A) << arma::endr
         << cos(E) * cos(A) << arma::endr
         << sin(E) << arma::endr;
   
       arma::vec N;
       if (assume_no_slope)
       {
   
           N << 0 << arma::endr
           << 0 << arma::endr
           << 1 << arma::endr;
       }
       else
       {
           Vector_3 n = face->normal();
   
   
           N << n[0] << arma::endr
           << n[1] << arma::endr
           << n[2] << arma::endr;
       }
   
   
       double angle = acos(arma::dot(S,N));
       angle = cos(angle);
   
       if(angle < 0.0 || E < 0.0523598776) //3deg -> rad
           angle = 0.0;
   
       (*face)["solar_angle"_s]=angle;
   
       //if we have remote shadowing
       if(has_optional("shadow"))
       {
           if((*face)["shadow"_s] == 1)
           {
               angle = 0;
           }
       }
   
       double swr =  angle * (*face)["iswr_direct"_s];
       double diff = (*face)["iswr_diffuse"_s];
   
       (*face)["iswr_direct"_s]=swr ;
       (*face)["iswr"_s]= swr + diff ;
   
   }
   
   slope_iswr::~slope_iswr()
   {
   
   
   }