Extend MATXSR to handle dilution dependent spec and prod arrays of GROUPR

ReleaseNotes.md

@@ -1,6 +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.69](https://github.com/njoy/NJOY2016/pull/XXX)
+This update fixes a number of issues:
+  - MATXSR was modified to handle the changes of GROUPR introduced by [NJOY2016.49](https://github.com/njoy/NJOY2016/pull/119). The 10d record now reproduces the spec and prod arrays for the dilutions given in UNRESR,PURR,GROUPR.
+
 ## [NJOY2016.68](https://github.com/njoy/NJOY2016/pull/264)
 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.

src/matxsr.f90

@@ -2214,12 +2214,12 @@ subroutine shufl(nscr,iref,jg1lo,jg1hi,nout3,irec3,irite)
    go to 145
   135 do ik=1,ng2
       jg2=noutg-ig2lo-ik+2
-      a(icdat-1+jg2)=b(irinp+lz+ik-1)
+      a(icdat-1+jg2)=b(irinp+lz+(ik-1)*nz+iz-1)
    enddo
    go to 240
   140 jconst=ig
    jg1=ning+1-ig
-   a(icdat-1+noutg+jg1)=b(irinp+lz+1)
+   a(icdat-1+noutg+jg1)=b(irinp+lz+nz+iz-1)
    go to 240
 
    !--store matrix sub-blocks

tests/77/CMakeLists.txt

@@ -0,0 +1,12 @@
+configure_file("${CMAKE_CURRENT_SOURCE_DIR}/input"
+               "${CMAKE_CURRENT_BINARY_DIR}/input" COPYONLY )
+
+configure_file("${CMAKE_CURRENT_SOURCE_DIR}/referenceTape26"
+               "${CMAKE_CURRENT_BINARY_DIR}/referenceTape26" COPYONLY )
+
+configure_file("${RESOURCES}/n-092_U_235-ENDF8.0.endf"
+               "${CMAKE_CURRENT_BINARY_DIR}/tape20" COPYONLY )
+
+add_test( NAME "Test77"
+          COMMAND ${Python3_EXECUTABLE} "../execute.py"
+          WORKING_DIRECTORY "${CMAKE_CURRENT_BINARY_DIR}" )

tests/77/input

@@ -0,0 +1,50 @@
+ 'moder'
+ 20  -21
+ 'reconr'
+ -21 -22
+ 'pendf tape for  U235'/
+ 9228 3 0
+ .0005 /
+ 'U 235'/
+ 'processed by the njoy nuclear data processing system'/
+ 'see original file for details of evaluation'/
+ 0/
+ 'broadr'
+ -21 -22 -23
+ 9228  2 0 0 0
+ .0005 1.e6 /
+ 300. 900.
+ 0/
+ 'unresr'
+ -21 -23 -24
+ 9228    2 8 1
+ 300. 900.
+ 1.e+10  1.e8 1.e6 1.e4 1.e3 1.e2 1.e1 1.
+ 0/
+'groupr'
+ -21 -24 0 25
+ 9228  9 0 5 5 2 8 1 1
+ 'U 235'/
+ 300. 900.
+ 1.e10  1.e8 1.e6 1.e4 1.e3 1.e2 1.e1 1. /
+ 3 /
+ 5 18 'chi ' /
+ 5 455 'chid' /
+ 6 /
+ 0/
+ 3 /
+ 6 /
+ 0/
+ 0/
+'matxsr'
+ 25  0 +26/
+ 1 'NJOY 2016 CI'/
+ 1 1  1 1
+ '69  group neutron interaction library'/
+ 'n'
+ 69
+ 'nscat'
+ 1
+ 1
+ 'U 235'  9228  /
+'stop'