From 652c040e8e063461b8ae1684b7b05e13fd7084d8 Mon Sep 17 00:00:00 2001 From: Bryan Herman Date: Fri, 6 Apr 2012 06:39:54 -0700 Subject: [PATCH] correct diffusion coefficient is now used --- src/cmfd_data.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/cmfd_data.F90 b/src/cmfd_data.F90 index 4eff5abe7b..ce98bbc870 100644 --- a/src/cmfd_data.F90 +++ b/src/cmfd_data.F90 @@ -154,10 +154,10 @@ contains cmfd % diffusion(h,i,j,k) = t % scores(score_index,4) % sum / flux ! calculate diffusion coefficient - cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*cmfd%totalxs(h,i,j,k)) +! cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*cmfd%totalxs(h,i,j,k)) ! cmfd % diffcof(h,i,j,k) = 1.0_8/(3.0_8*(cmfd % totalxs(h,i,j,k) -& ! & cmfd % p1scattxs(h,i,j,k))) -! cmfd % diffcof(h,i,j,k) = cmfd % diffusion(h,i,j,k) + cmfd % diffcof(h,i,j,k) = cmfd % diffusion(h,i,j,k) else if (ital == 2) then @@ -380,7 +380,7 @@ contains ! compute dtilde dtilde = (2*cell_dc*(1-ref_albedo))/(4*cell_dc*(1+ & & ref_albedo)+(1-ref_albedo)*cell_hxyz(xyz_idx)) - dtilde = 0.0_8 + ! dtilde = 0.0_8 else ! not next to a reflector or no core map ! compute dtilde