Skip to content

Commit

Permalink
Increasing nu array size in groupr
Browse files Browse the repository at this point in the history
  • Loading branch information
whaeck committed Aug 23, 2023
1 parent 4d987aa commit ba34c1a
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 2 deletions.
1 change: 1 addition & 0 deletions ReleaseNotes.md
Original file line number Diff line number Diff line change
Expand Up @@ -4,6 +4,7 @@ Given here are some release notes for NJOY2016. Each release is made through a f
## [NJOY2016.72](https://github.com/njoy/NJOY2016/pull/xxx)
This update fixes the following issues:
- Fixed an issue in GROUPR related to an error coming up in production matrix calculations. Depending on when a user asks for a production matrix associated to a reaction, it is possible that the reference frame of the previous reaction is used instead (caused by erronously defining an already declared global variable as local with a "save" attribute). In some circumstances, this causes NJOY2016 to error out (with a message related to unsupported reference frames). No test results had to be updated due to this change.
- Increased allocation of arrays to accommodate ENDF/B-VIII.1 evaluations.

## [NJOY2016.71](https://github.com/njoy/NJOY2016/pull/301)
This update adds the new MF7 MT451 (thermal scattering general information) ENDF format to MODER so that this module will be able to interpret the new MF7 section. No other capability in NJOY2016 currently uses the information in this section.
Expand Down
4 changes: 2 additions & 2 deletions src/groupr.f90
Original file line number Diff line number Diff line change
Expand Up @@ -11385,15 +11385,15 @@ subroutine conver(nin,nout,nscr)
call contio(nin,nout,nscr,scr,nb,nw)
if (mfh.eq.0) go to 110
if (mth.ne.452) go to 595
nnu=8000
nnu=10000
allocate(nu(nnu))
l=1
lnu=l2h
if (lnu.eq.1) call listio(nin,nout,nscr,nu(l),nb,nw)
if (lnu.eq.2) call tab1io(nin,nout,nscr,nu(l),nb,nw)
do while (nb.ne.0)
if (l+nw.gt.nnu) call error('conver',&
'storage for fission nu exceeded',' ')
'storage for fission nu exceeded (increase the size of the nu array)',' ')
l=l+nw
call moreio(nin,nout,nscr,nu(l),nb,nw)
enddo
Expand Down

0 comments on commit ba34c1a

Please sign in to comment.