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

read WMP and write message in load_nuclide capi function

Browse source
This commit is contained in:
liangjg 2019-01-20 17:01:58 -05:00
parent 56e6c4088d
commit 9b93c809eb
1 changed files with 43 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

43
src/nuclide_header.F90
View file

@ -11,6 +11,9 @@ module nuclide_header
use endf_header, only: Function1D, Polynomial, Tabulated1D
use error
use hdf5_interface
use math, only: faddeeva, w_derivative, &
broaden_wmp_polynomials
use multipole_header
use message_passing
use random_lcg, only: prn, future_prn, prn_set_stream
use reaction_header, only: Reaction
@ -728,6 +731,8 @@ contains
integer(C_INT) :: err
integer :: n
logical :: file_exists
character(7) :: readable
integer(HID_T) :: file_id
integer(HID_T) :: group_id
character(:), allocatable :: name_
@ -751,6 +756,9 @@ contains
filename = library_path(LIBRARY_NEUTRON, to_lower(name_))
call write_message('Reading ' // trim(name_) // ' from ' // &
trim(filename), 6)
! Open file and make sure version is sufficient
file_id = file_open(filename, 'r')
call check_data_version(file_id)
@ -769,6 +777,41 @@ contains
! Initialize nuclide grid
call nuclides(n) % init_grid()
! Read WMP data if necessary
nuclides(n) % mp_present = .false.
if (temperature_multipole .and. &
library_present(LIBRARY_WMP, to_lower(name_))) then
filename = library_path(LIBRARY_WMP, to_lower(name_))
! Check if Multipole library exists and is readable
inquire(FILE=filename, EXIST=file_exists, READ=readable)
if (file_exists) then
if (readable(1:3) == 'NO') then
call fatal_error("Multipole library '" // trim(filename) // "' &
&is not readable! Change file permissions.")
end if
! Display message
call write_message("Reading " // trim(name_) // " WMP data from " &
// filename, 6)
! Open file and make sure version is sufficient
file_id = file_open(filename, 'r')
call check_wmp_version(file_id)
! Read nuclide data from HDF5
group_id = open_group(file_id, name_)
allocate(nuclides(n) % multipole)
call nuclides(n) % multipole % from_hdf5(group_id)
nuclides(n) % mp_present = .true.
! Close the group and file.
call close_group(group_id)
call file_close(file_id)
end if
end if
else
err = E_DATA
call set_errmsg("Nuclide '" // trim(name_) // "' is not present &