From f627f1645b296226c2b23858f33e28fdfd70f67b Mon Sep 17 00:00:00 2001 From: Bryan Herman Date: Tue, 10 Jan 2012 17:40:12 -0500 Subject: [PATCH] fixed bug in assigning the front-faced outgoing current to cmfd object --- src/cmfd_execute.F90 | 2 +- src/cmfd_utils.F90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/cmfd_execute.F90 b/src/cmfd_execute.F90 index f3a0959bdb..b57929ebaf 100644 --- a/src/cmfd_execute.F90 +++ b/src/cmfd_execute.F90 @@ -197,7 +197,7 @@ contains score_index = sum((bins - 1) * t % stride) + 1 ! incoming cmfd % current(7,h,i,j,k) = t % scores(score_index,1) % val bins(TS_SURFACE) = OUT_FRONT - score_index = sum(t % stride(1:3) * ijk) + OUT_FRONT ! outgoing + score_index = sum((bins - 1) * t % stride) + 1 ! outgoing cmfd % current(8,h,i,j,k) = t % scores(score_index,1) % val ! bottom surface diff --git a/src/cmfd_utils.F90 b/src/cmfd_utils.F90 index f40b944c21..370bc62bff 100644 --- a/src/cmfd_utils.F90 +++ b/src/cmfd_utils.F90 @@ -134,7 +134,7 @@ contains ! get idx from core map nidx = ng*(cmfd % coremap(i,j,k)) - (ng - g) - write(200,*) g,i,j,k,nidx + else ! compute index