Merge pull request #343 from njoy/fix/viewr-array

Increasing array size

ReleaseNotes.md

@@ -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:

src/viewr.f90

@@ -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
@@ -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
@@ -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
@@ -1663,4 +1667,3 @@ subroutine ascalv(m,z1,z2,major,minor)
    end subroutine ascalv
 
 end module viewm
-