crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 06:05:58 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

tested a surface current extraction in print cmfd

Browse source
This commit is contained in:
Bryan Herman 2011-11-12 19:00:43 -05:00
parent 289e5c15ee
commit eadf4f0945
1 changed files with 14 additions and 6 deletions
Download patch file Download diff file Expand all files Collapse all files

20
src/cmfd_utils.f90
View file

@ -64,26 +64,34 @@ contains
type(StructuredMesh), pointer :: m ! pointer for mesh object
! associate pointers with objects
t => tallies(1)
t => tallies(3)
m => meshes(t % mesh)
! set all bins to 1
bins = 1
! get mesh indices, first we will first force to 1,1,1
ijk = (/ 1, 1, 1 /)
! ijk = (/ 1, 1, 1 /)
! apply filters, here we will just try a mesh filter first
bins(T_MESH) = mesh_indices_to_bin(m,ijk)
! bins(T_MESH) = mesh_indices_to_bin(m,ijk)
! calculate score index from bins
score_index = sum((bins - 1) * t%stride) + 1
! score_index = sum((bins - 1) * t%stride) + 1
! get value from tally object
tally_val = t%scores(score_index,1)%val
! tally_val = t%scores(score_index,2)%val
! write value to file
write(7,*) "Tally value is:",tally_val
! write(7,*) "Tally value is:",tally_val
! Left Surface
ijk = (/ 1-1, 1, 1 /)
score_index = sum(t % stride(1:3) * ijk) + IN_RIGHT
print *, "Outgiong Current from Left", t % scores(score_index,1) % val
score_index = sum(t % stride(1:3) * ijk) + OUT_RIGHT
print *, "Incoming Current from Left", t % scores(score_index,1) % val
end subroutine print_cmfd