Merge pull request #358 from njoy/fix/arrays

Fix/arrays

ReleaseNotes.md

@@ -1,9 +1,10 @@
 # Release Notes—NJOY2016
 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.78
+## [NJOY2016.78](https://github.com/njoy/NJOY2016/pull/xxx)
 This update fixes the following issues:
   - Tape numbers are now allowed to go up to 999 instead of 99 previously. This change was made to accommodate processing of the light water evaluation in ENDF/B-VIII.1 that has 94 temperatures. Test 84 was added to detect this issue in the future.
+  - A few updates were made related to array sizes and array allocation.
 
 ## [NJOY2016.77](https://github.com/njoy/NJOY2016/pull/347)
 This update fixes the following issues:

src/groupr.f90

@@ -6514,6 +6514,7 @@ subroutine getyld(e,enext,idis,yld,mat,mf,mt,lfs,itape)
    enext=tmp(7+2*nr)
   190 continue
    idis=0
+   if (allocated(yield)) deallocate(yield)
    allocate(yield(na))
    do i=1,na
       yield(i)=tmp(i)

src/moder.f90

@@ -54,12 +54,13 @@ subroutine moder
    integer::nin,nout,no,loop,i,nz,inout
    integer::matd,mend,ig,nk,ik,nb,nw
    integer::nscr,ninl,matl
+   integer,parameter::nbuf=1000000
    real(kr)::time
    character(105)::strng
    character(4)::hb(17)
    real(kr)::rb(17)
    equivalence(hb(1),rb(1))
-   real(kr)::a(5200)
+   real(kr),allocatable::a(:)
    real(kr)::z(20)
    character(4)::zc(20)
    equivalence(zc(1),z(1))
@@ -68,6 +69,8 @@ subroutine moder
    ninl=0
    matl=-2
 
+   allocate(a(nbuf)) ! deallocated automatically at end of subroutine
+
    !--read user input, initialize,
    !--and write output header.
    call timer(time)