get_drift_flux_total Function

public function get_drift_flux_total(comm_handler, polygon, drift_flux, sum_over_planes, iseg) result(res)

Computes total flux through polygon segment from drift_flux array, which has to be computed previously via compute_drift_flux. $omp parallel default(none) private(is, ks, ks_glob) & $omp shared(polygon, iseg_start, iseg_end, nexclude_last, drift_flux, res) $omp do reduction(+ : res) $omp end do $omp end parallel

Arguments

Type IntentOptional Attributes Name
type(comm_handler_t), intent(in) :: comm_handler

Communicators

type(polygon_in_mesh_t), intent(in) :: polygon

Polygon

real(kind=GP), intent(in), dimension(polygon%get_np_total()) :: drift_flux

Drift flux through boundary in between polygon points, obtained via compute_drift_flux routine

logical, intent(in) :: sum_over_planes

If true summation over planes is performed. This assumes that polygons have same structure over planes

integer, intent(in), optional :: iseg

Segment for which flux is computed. If not present then flux summed over all segments is returned

Return Value real(kind=GP)


Contents