Merge pull request #145 from njoy/fix/tendl-zero-awr

Added an error message to detect zero awp values for non photons
njoy · Jan 22, 2020 · e7eae75 · e7eae75
2 parents e2d4c84 + ef41db2
commit e7eae75
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 2 deletions.
diff --git a/src/heatr.f90 b/src/heatr.f90
@@ -2791,6 +2791,14 @@ subroutine sixbar(e,ebar,yld,dame,nin,c,ncmax,nscr,b,nbmax,&
    awp=c2h
    law=l2h
    l=l+nw
+
+   !--error for awp=0. for non photons
+   if (zap.ne.0.and.awp.eq.0) then
+     write(strng,'(''awp is zero for for mt'',i3,'' particle '',i5)')&
+       mth,nint(zap)
+     call error('sixbar',strng,' ')
+   endif
+
    do while (nb.ne.0)
       call moreio(nin,0,0,c(l),nb,nw)
       l=l+nw

diff --git a/src/vers.f90 b/src/vers.f90
@@ -3,6 +3,6 @@ module version
    ! These values are updated during the NJOY revision-control process.
    implicit none
    private
-   character(8),public::vers='2016.53'
-   character(8),public::vday='07Nov19'
+   character(8),public::vers='2016.54'
+   character(8),public::vday='16Jan20'
 end module version