Merge pull request #73 from njoy/feature/purr-endf8

Feature/purr endf8

Author: jlconlin

docs/testDescription.md

@@ -182,3 +182,9 @@ title: NJOY2016 Test Descriptions
 
   This test is an addition to test 3 using `MATXSR`. While test 3 only uses photons, test 30 uses both neutrons and photons.
 
+## Test Problem 31
+
+[[input](https://raw.githubusercontent.com/njoy/NJOY2016/master/test/31/input)]
+
+  This test runs `PURR` for Pu240 from ENDF/B-VIII.0. This evaluation is one of the few that still gave negative cross sections in the probability tables for ENDF/B-VIII.0. In most cases, this behaviour is due to the fact that the original evaluation has an LSSF flag of 0 (MF3 contains background cross sections) and those background cross sections in the unresolved resonance region are negative (an example would be Na22). In the case of Pu240 however, the LSSF flag is set to 1 (MF3 contains the actual cross sections and the unresolved resonances should only be used for self-shielding). The negative cross section values in the probability table for this nuclide were due to the fact that the total cross section was actually lower than its components.
+

src/groupr.f90

@@ -373,7 +373,7 @@ subroutine groupr
      &8x,f8.1,''s'')') time
    write(nsyse,'(/'' groupr...'',59x,f8.1,''s'')') time
    call ruinb(iaddmt)
-   nwscr=10000
+   nwscr=max(10000,nsigz+ngn+9,nsigz+ngg+8)
    allocate(scr(nwscr))
    if (nendf.lt.0) nend2=-nend2
    if (nendf.lt.0) nend3=-nend3

src/purr.f90

@@ -265,7 +265,16 @@ subroutine purr
       do i=1,nunr
          loc=i
          sb(loc)=sb(loc)-sb(loc+nunr)-sb(loc+2*nunr)-sb(loc+3*nunr)
-         if (abs(sb(loc)).lt.sigmin) sb(loc)=0
+         if (sb(loc).lt.sigmin) then
+            if (sb(loc).lt.-sigmin) then
+               write(strng1,&
+                 '(''total xs less than its components at e=''&
+                 &,1p,e12.4)') eunr(loc)
+               write(strng2,'(''check evaluation file'')')
+               call mess('purr',strng1,strng2)
+            endif
+            sb(loc)=0
+         endif
          sb(loc+nunr)=0
          sb(loc+2*nunr)=0
          sb(loc+3*nunr)=0

src/vers.f90

@@ -3,7 +3,7 @@ module version
    ! These values are updated during the NJOY revision-control process.
    implicit none
    private
-   character(8),public::vers='2016.29'
-   character(8),public::vday='01Mar18'
+   character(8),public::vers='2016.30'
+   character(8),public::vday='21Mar18'
 end module version
 

tests/31/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}/referenceTape30"
+               "${CMAKE_CURRENT_BINARY_DIR}/referenceTape30" COPYONLY )
+	     
+configure_file("${RESOURCES}/n-094_Pu_240-ENDF8.0.endf"
+               "${CMAKE_CURRENT_BINARY_DIR}/tape20" COPYONLY )
+
+add_test( NAME "Test31"
+          COMMAND ${PYTHON_EXECUTABLE} "../execute.py"
+          WORKING_DIRECTORY "${CMAKE_CURRENT_BINARY_DIR}" )

tests/31/input

@@ -0,0 +1,24 @@
+moder
+ 20 -21
+reconr
+ -21 -22/
+ 'reconr/pendf tape for 9440'/
+ 9440 0/
+ .001/
+ 0 /
+broadr
+ -21 -22 -23/
+ 9440 1 0 0/
+ .001/
+ 293.6/
+ 0/
+purr
+ -21 -23 -24/
+ 9440 1 1 20 4/
+ 293.6 /
+ 1.e10 /
+ 0/
+moder
+ -24 30
+stop
+

tests/31/referenceTape30

@@ -40,4 +40,5 @@ add_subdirectory( "27" )
 add_subdirectory( "28" )
 add_subdirectory( "29" )
 add_subdirectory( "30" )
+add_subdirectory( "31" )