Skip to content

Commit

Permalink
Fixed a division by zero for IFENG=2 files
Browse files Browse the repository at this point in the history
  • Loading branch information
whaeck committed Jul 13, 2022
1 parent a9695cf commit 78466e5
Show file tree
Hide file tree
Showing 2 changed files with 4 additions and 2 deletions.
4 changes: 3 additions & 1 deletion ReleaseNotes.md
Original file line number Diff line number Diff line change
Expand Up @@ -2,7 +2,9 @@
Given here are some release notes for NJOY2016. Each release is made through a formal [Pull Request](https://github.com/njoy/NJOY2016/pulls) made on GitHub. There are links in this document that point to each of those Pull Requests, where you can see in great details the changes that were made. Often the Pull Requests are made in response to an [issue](https://github.com/njoy/NJOY2016/issues). In such cases, links to those issues are also given.

## [NJOY2016.68](https://github.com/njoy/NJOY2016/pull/xxx)
This updates resolves indexing errors in the calculation of inelastic thermal scattering mubar in ACER module for IFENG=2 ACE files. Test 74 was added to track this issue.
This update fixes a number of minor issues:
- there was an indexing error in the calculation of inelastic thermal scattering mubar in ACER for IFENG=2 ACE files. Test 74 was added to track this issue.
- there was a division by zero that caused a segmentation fault in Debug mode, also for IFENG=2 ACE files.

## [NJOY2016.67](https://github.com/njoy/NJOY2016/pull/243)
This update resolves a number of issues encountered when processing the newly released JENDL5 nuclear data library and TENDL-2021 library for evaluations with sub-actinide fission, and provides fixes for a few minor issues.
Expand Down
2 changes: 1 addition & 1 deletion src/aceth.f90
Original file line number Diff line number Diff line change
Expand Up @@ -2068,7 +2068,7 @@ subroutine tplots(nout,hk)
do j=2,nbini-1
ep=(xss(loc+1)+epl)/2
if (xss(loc+1).gt.epl) then
if (j.eq.2.or.ep/epl.gt.skip) then
if (j.eq.2.or.ep.gt.skip*epl) then
write(nout,'(1p,e14.6,''/'')') ep
ul=-1
write(nout,'(1p,2e14.6,''/'')') -1.,zmin
Expand Down

0 comments on commit 78466e5

Please sign in to comment.