Пример #1
0
void 
UZRectConst::tick (const GeometryRect& geo, const std::vector<size_t>&,
                   const Soil&, 
                   SoilWater& soil_water, const SoilHeat&,
                   const Surface&, const Groundwater&, 
                   const double, Treelog&)

{
  const size_t edge_size = geo.edge_size (); // number of edges 

  for (size_t edge = 0; edge != edge_size; ++edge) 
    {
      const double sin_angle = geo.edge_sin_angle (edge);
      const double cos_angle = geo.edge_cos_angle (edge);
      const double q = q_z * sin_angle + q_x * cos_angle;
      soil_water.set_flux (edge, q);
    }
}