From ef16a6cb575e9f907a2c2babc091aca5d9c7d86d Mon Sep 17 00:00:00 2001 From: Paul Romano Date: Thu, 1 Feb 2018 12:35:29 -0600 Subject: [PATCH] Fix bug in stochastic volume calculation with void cell --- src/volume_calc.F90 | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/src/volume_calc.F90 b/src/volume_calc.F90 index 15c0497699..e374f21064 100644 --- a/src/volume_calc.F90 +++ b/src/volume_calc.F90 @@ -193,11 +193,13 @@ contains if (this % domain_type == FILTER_MATERIAL) then i_material = p % material - do i_domain = 1, size(this % domain_id) - if (materials(i_material) % id == this % domain_id(i_domain)) then - call check_hit(i_domain, i_material, indices, hits, n_mat) - end if - end do + if (i_material /= MATERIAL_VOID) then + do i_domain = 1, size(this % domain_id) + if (materials(i_material) % id == this % domain_id(i_domain)) then + call check_hit(i_domain, i_material, indices, hits, n_mat) + end if + end do + end if elseif (this % domain_type == FILTER_CELL) THEN do level = 1, p % n_coord