Transect.hpp

Go to the documentation of this file.

 
#ifndef OPENSWMM_TRANSECT_HPP
#define OPENSWMM_TRANSECT_HPP
 
#include <vector>
#include <string>
 
namespace openswmm {
 
namespace transect {
 
constexpr int N_TRANSECT_TBL = 51;
 
struct TransectData {
    std::string name;
    double n_left   = 0.0;   
    double n_right  = 0.0;   
    double n_channel = 0.0;  
    double length_factor = 1.0;  
 
    // Station-elevation pairs
    std::vector<double> stations;    
    std::vector<double> elevations;  
 
    double x_left_bank  = 0.0;  
    double x_right_bank = 0.0;  
 
    // Computed geometry tables (normalized)
    double y_full = 0.0;
    double w_max  = 0.0;
    double a_full = 0.0;
    double r_full = 0.0;
    // CUSTOM shapes only (legacy TShape.sMax/aMax, shape.c getSmax): max
    // section factor and its area for a shape of UNIT height, computed on the
    // unnormalized tables. Scaled by yFull in the caller (xsect.c:685-686).
    double s_max  = 0.0;
    double a_max  = 0.0;
 
    double area_tbl[N_TRANSECT_TBL]  = {};
    double hrad_tbl[N_TRANSECT_TBL]  = {};
    double width_tbl[N_TRANSECT_TBL] = {};
};

void buildTables(TransectData& td);

void buildCustomTables(TransectData& td, double y_full,
                       const double* curve_x, const double* curve_y, int n_pts);
 
} // namespace transect
} // namespace openswmm
 
#endif // OPENSWMM_TRANSECT_HPP