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

flux is not output by group in vtk file

Browse source
This commit is contained in:
Bryan Herman 2012-01-14 16:25:06 -05:00
parent 820fe07212
commit 9d3689c650
1 changed files with 22 additions and 10 deletions
Download patch file Download diff file Expand all files Collapse all files

32
src/cmfd_utils.F90
View file

@ -876,6 +876,7 @@ contains
integer :: i ! x loop counter
integer :: j ! y loop counter
integer :: k ! z loop counter
integer :: g ! group counter
integer :: nx ! number of mesh cells in x direction
integer :: ny ! number of mesh cells in y direction
integer :: nz ! number of mesh cells in z direction
@ -893,13 +894,15 @@ contains
type(StructuredMesh), pointer :: m => null() ! pointer to mesh
! vtk specific variables
integer :: E_IO ! error code
integer :: nn ! number of nodes
integer :: nc ! number of cells
integer :: con(8) ! connectivity vector
integer :: off(1:1) ! offset, number of nodes in cell
integer(1) :: cell_id(1:1) ! cell type
real(8) :: real_buffer(1:1) ! real data buffer 8-byte
integer :: E_IO ! error code
integer :: nn ! number of nodes
integer :: nc ! number of cells
integer :: con(8) ! connectivity vector
integer :: off(1:1) ! offset, number of nodes in cell
integer(1) :: cell_id(1:1) ! cell type
real(8) :: real_buffer(1:1) ! real data buffer 8-byte
character(len=40) :: varname ! name of output variable
character(len=3) :: str_g ! string for energy group #
! extract spatial and energy indices from object
nx = cmfd % indices(1)
@ -948,9 +951,18 @@ contains
! open data block in vtk file
E_IO = VTK_DAT_XML('cell','open')
! write out flux
real_buffer = (/cmfd%flux(1,i,j,k)/)
E_IO = VTK_VAR_XML(nc,'flux',real_buffer)
! loop around energy
GROUP: do g = 1,ng
! convert group int to str
write(str_g,'(I3)') g
! write out flux
real_buffer = (/cmfd%flux(g,i,j,k)/)
varname = 'flux_'//trim(adjustl(str_g))
E_IO = VTK_VAR_XML(nc,varname,real_buffer)
end do GROUP
! close data block in vtk file
E_IO = VTK_DAT_XML('cell','close')