Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Increasing array size #343

Merged
merged 4 commits into from
Aug 12, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
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.77](https://github.com/njoy/NJOY2016/pull/xxx)
This update fixes the following issues:
- Background cross section values for reactions other than total, elastic, fission and capture were not handled properly in the unresolved resonance region. The background for total, elastic, fission and capture is integrated into the unresolved resonance cross section values in the genunr subroutine. In the emerge subroutine, the background in the unresolved resonance region was therefore zeroed out for any resonance reaction. This has never been an issue but now we have LRF=7 evaluations that can define reactions other than total, elastic, fission and capture. Test 82 was added to detect this issue in the future.
- Increased the size of internal arrays in VIEWR.

## [NJOY2016.76](https://github.com/njoy/NJOY2016/pull/340)
This update fixes the following issues:
Expand Down
9 changes: 6 additions & 3 deletions src/viewr.f90
Original file line number Diff line number Diff line change
Expand Up @@ -340,7 +340,7 @@ subroutine viewr
real(kr)::xx,yy,zz
character(80)::text
integer,parameter::mmax=20000 !same in plotr and viewr
integer,parameter::maxaa=200000
integer,parameter::maxaa=500000
real(kr),dimension(15)::z
real(kr),dimension(maxaa)::aa
real(kr),dimension(mmax)::x,y,b,dxm,dxp,dym,dyp
Expand Down Expand Up @@ -1269,10 +1269,11 @@ subroutine set3d(iplot,xyz,nxyz)
integer::iplot,nxyz
real(kr)::xyz(nxyz)
! internals
integer,parameter::length=2000
integer::i,j,nn,k,ncurv,major,minor,itop,ibot,l,n
real(kr)::wt,xlo,xhi,ylo,yhi,zlo,zhi,xn,yn,zn,yy
real(kr)::top,bot,xop,yop,zop
integer::lll(400)
integer::lll(length)
real(kr)::x(2000),y(2000),z(2000)
real(kr),parameter::big=1.e10_kr
real(kr),parameter::d0=.001e0_kr
Expand Down Expand Up @@ -1322,6 +1323,9 @@ subroutine set3d(iplot,xyz,nxyz)
enddo
i=i+2+2*nn
j=j+1
if (j.gt.length) then
call error('set3d','array overflow, increase the length parameter',' ')
endif
endif
enddo
ncurv=j-1
Expand Down Expand Up @@ -1663,4 +1667,3 @@ subroutine ascalv(m,z1,z2,major,minor)
end subroutine ascalv

end module viewm