SurfaceSpeciesKernel.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: Apache-2.0
 
#ifndef OPENSWMM_ENGINE_QUALITY_SURFACE_SPECIES_KERNEL_HPP
#define OPENSWMM_ENGINE_QUALITY_SURFACE_SPECIES_KERNEL_HPP
 
#include <cstddef>
#include <cstdint>
#include <vector>
 
#include "Landuse.hpp"
 
namespace openswmm {
struct SimulationContext;
namespace surfacespecies {

struct Block {
    int n = 0;                                                      
    int nlu = 0;                                                    
    const std::vector<landuse::BuildupParams>* bu_params = nullptr; 
    const std::vector<landuse::WashoffParams>* wo_params = nullptr; 
    std::vector<double>* buildup = nullptr; 
    std::vector<double>* conc = nullptr;    
    const std::vector<double>* ponded_reacted_mass = nullptr;
    const std::vector<uint8_t>* ponded_reactor_mask = nullptr;
    const std::vector<uint8_t>* lid_reactor_mask = nullptr;
    std::vector<double>* ponded_qual = nullptr;     
    std::vector<double>* runon_qual_rate = nullptr; 
    std::vector<double>* total_load = nullptr;      
    const std::vector<double>* mcf = nullptr;       
    const std::vector<double>* c_rain = nullptr;    
    const std::vector<bool>* snow_only = nullptr;   
    const std::vector<int>* co_pollut = nullptr;    
    const std::vector<double>* co_frac = nullptr;   
    // Ledgers [s], user mass.
    std::vector<double>* led_surface_buildup = nullptr;
    std::vector<double>* led_wet_deposition = nullptr;
    std::vector<double>* led_bmp_removal = nullptr;
    std::vector<double>* led_infil_loss = nullptr;
    std::vector<double>* led_runoff_load = nullptr;
    std::vector<double>* led_sweeping = nullptr;
    std::vector<double>* led_ponded_lost = nullptr; 
    const std::vector<unsigned char>* sweep_now = nullptr;
    const std::vector<double>* sweep_removal = nullptr;

    std::size_t bu_idx(int sc, int lu, int s) const noexcept {
        return static_cast<std::size_t>(sc) * static_cast<std::size_t>(nlu * n) +
               static_cast<std::size_t>(lu) * static_cast<std::size_t>(n) +
               static_cast<std::size_t>(s);
    }
};

void step(SimulationContext& ctx, Block& b, double dt_runoff, double old_runoff_time) noexcept;

void bookPostLid(SimulationContext& ctx, Block& b, double dt_runoff,
                 const std::vector<double>& lid_infil_vol) noexcept;
 
} // namespace surfacespecies
} // namespace openswmm
 
#endif