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

Added comments, got rid of unused variable, added check for file existing/readable.

Browse source
This commit is contained in:
Paul Romano 2011-11-10 18:01:51 -05:00
parent 73228f2924
commit 41b1810f72
2 changed files with 49 additions and 35 deletions
Download patch file Download diff file Expand all files Collapse all files

60
src/cross_section.f90
View file

@ -43,10 +43,9 @@ contains
integer :: index_sab ! index in sab_tables
character(10) :: name ! name of isotope, e.g. 92235.03c
character(10) :: alias ! alias of nuclide, e.g. U-235.03c
character(MAX_FILE_LEN) :: path ! path to ACE cross section library
type(Material), pointer :: mat => null()
type(Nuclide), pointer :: nuc => null()
type(SAB_Table), pointer :: sab => null()
type(Material), pointer :: mat => null()
type(Nuclide), pointer :: nuc => null()
type(SAB_Table), pointer :: sab => null()
! ==========================================================================
! COUNT NUMBER OF TABLES AND CREATE DICTIONARIES
@ -172,6 +171,8 @@ contains
end if
end do
call dict_delete(already_read)
! ==========================================================================
! ASSIGN S(A,B) TABLES TO SPECIFIC NUCLIDES WITHIN MATERIALS
@ -210,26 +211,28 @@ contains
subroutine read_ace_table(index_table, index)
integer, intent(in) :: index_table
integer, intent(in) :: index
integer, intent(in) :: index_table ! index in nuclides/sab_tables
integer, intent(in) :: index ! index in xs_listings
integer :: i
integer :: j, j1, j2
integer :: record_length
integer :: location
integer :: entries
integer :: length
integer :: in = 7
integer :: zaids(16)
integer :: filetype
real(8) :: kT
real(8) :: awrs(16)
real(8) :: awr
character(10) :: name
character(10) :: date
character(10) :: mat
character(70) :: comment
character(MAX_FILE_LEN) :: filename ! path to ACE cross section library
integer :: i ! loop index for XSS records
integer :: j, j1, j2 ! indices in XSS
integer :: record_length ! Fortran record length
integer :: location ! location of ACE table
integer :: entries ! number of entries on each record
integer :: length ! length of ACE table
integer :: in = 7 ! file unit
integer :: zaids(16) ! list of ZAIDs (only used for S(a,b))
integer :: filetype ! filetype (ASCII or BINARY)
real(8) :: kT ! temperature of table
real(8) :: awrs(16) ! list of atomic weight ratios (not used)
real(8) :: awr ! atomic weight ratio for table
logical :: file_exists ! does ACE library exist?
character(7) :: readable ! is ACE library readable?
character(10) :: name ! name of ACE table
character(10) :: date ! date ACE library was processed
character(10) :: mat ! material identifier
character(70) :: comment ! comment for ACE table
character(MAX_FILE_LEN) :: filename ! path to ACE cross section library
type(Nuclide), pointer :: nuc => null()
type(SAB_Table), pointer :: sab => null()
type(XsListing), pointer :: listing => null()
@ -242,6 +245,17 @@ contains
entries = listing % entries
filetype = listing % filetype
! Check if ACE library exists and is readable
inquire(FILE=filename, EXIST=file_exists, READ=readable)
if (.not. file_exists) then
message = "ACE library '" // trim(filename) // "' does not exist!"
call fatal_error()
elseif (readable(1:3) == 'NO') then
message = "ACE library '" // trim(filename) // "' is not readable! &
&Change file permissions with chmod command."
call fatal_error()
end if
! display message
message = "Loading ACE cross section table: " // listing % name
call write_message(6)

24
src/cross_section_header.f90
View file

@ -154,18 +154,18 @@ module cross_section_header
!===============================================================================
type XsListing
character(10) :: name
character(10) :: alias
integer :: type
integer :: zaid
integer :: filetype
integer :: location
integer :: recl
integer :: entries
real(8) :: awr
real(8) :: temp
logical :: metastable
character(MAX_FILE_LEN) :: path
character(10) :: name ! table name, e.g. 92235.70c
character(10) :: alias ! table alias, e.g. U-235.70c
integer :: type ! type of table (cont-E neutron, S(A,b), etc)
integer :: zaid ! ZAID identifier = 1000*Z + A
integer :: filetype ! ASCII or BINARY
integer :: location ! location of table within library
integer :: recl ! record length for library
integer :: entries ! number of entries per record
real(8) :: awr ! atomic weight ratio (# of neutron masses)
real(8) :: temp ! temperature
logical :: metastable ! is this nuclide metastable?
character(MAX_FILE_LEN) :: path ! path to library containing table
end type XsListing
!===============================================================================