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

neutron balance residual info now written to hdf5 at the end of every cycle

Browse source
This commit is contained in:
Bryan Herman 2012-02-01 16:38:14 -05:00
parent 0b878b76be
commit 6e2783aa6f
2 changed files with 27 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

3
src/cmfd_header.F90
View file

@ -56,6 +56,9 @@ module cmfd_header
real(8), allocatable :: phi(:)
real(8) :: keff
! residual for neutron balance
real(8), allocatable :: resnb(:,:,:,:)
end type cmfd_obj
contains

32
src/cmfd_output.F90
View file

@ -37,9 +37,9 @@ contains
character(30) :: label
! open cmfd file for output
filename = "cmfd_"//trim(to_str(current_cycle))//".out"
open(FILE=filename, UNIT=UNIT_CMFD, STATUS='replace', ACTION='write', &
IOSTAT=io_error)
! filename = "cmfd_"//trim(to_str(current_cycle))//".out"
! open(FILE=filename, UNIT=UNIT_CMFD, STATUS='replace', ACTION='write', &
! IOSTAT=io_error)
! extract spatial and energy indices from object
nx = cmfd % indices(1)
@ -47,6 +47,9 @@ contains
nz = cmfd % indices(3)
ng = cmfd % indices(4)
! allocate res dataspace
if (.not. allocated(cmfd%resnb)) allocate(cmfd%resnb(ng,nx,ny,nz))
! begin loop around space and energy groups
ZLOOP: do k = 1,nz
@ -89,12 +92,15 @@ contains
! normalize by flux
res = res/cmfd%flux(g,i,j,k)
! bank res in cmfd object
cmfd%resnb(g,i,j,k) = res
! write output
label = "MESH (" // trim(int4_to_str(i)) // ". " // &
& trim(int4_to_str(j)) // ", " // trim(int4_to_str(k)) // &
& ") GROUP " // trim(int4_to_str(g))
write(UNIT=UNIT_CMFD, FMT='(A,T35,A)') label, &
& trim(real_to_str(res))
! label = "MESH (" // trim(int4_to_str(i)) // ". " // &
! & trim(int4_to_str(j)) // ", " // trim(int4_to_str(k)) // &
! & ") GROUP " // trim(int4_to_str(g))
! write(UNIT=UNIT_CMFD, FMT='(A,T35,A)') label, &
! & trim(real_to_str(res))
end do GROUPG
@ -258,6 +264,16 @@ contains
call h5sclose_f(dataspace_id,hdf5_err)
call h5dclose_f(dataset_id,hdf5_err)
! write neutron balance residual information from cmfd object
dim4 = (/ng,nx,ny,nz/)
call h5screate_simple_f(4,dim4,dataspace_id,hdf5_err)
call h5dcreate_f(hdf5_output_file,trim(cycname)//"/resnb",H5T_NATIVE_DOUBLE, &
& dataspace_id,dataset_id,hdf5_err)
call h5dwrite_f(dataset_id,H5T_NATIVE_DOUBLE,cmfd%resnb,dim4,hdf5_err)
call h5sclose_f(dataspace_id,hdf5_err)
call h5dclose_f(dataset_id,hdf5_err)
! only write the following if core map is active
if (cmfd_coremap) then