From ef41db22413aa8deb541c4813b924e0a61c53ca2 Mon Sep 17 00:00:00 2001 From: Wim Haeck Date: Thu, 16 Jan 2020 14:13:01 -0700 Subject: [PATCH] Added an error message to detect zero awp values for non photons in some TENDL evaluations --- src/heatr.f90 | 8 ++++++++ src/vers.f90 | 4 ++-- 2 files changed, 10 insertions(+), 2 deletions(-) diff --git a/src/heatr.f90 b/src/heatr.f90 index 3d0cb7d4..8ab80de3 100644 --- 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 index 86073801..436dec3d 100644 --- 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