Skip to content

Commit

Permalink
Updating pn zaid values with zaoption
Browse files Browse the repository at this point in the history
  • Loading branch information
whaeck committed Feb 7, 2024
1 parent ee89a97 commit de6d9bd
Show file tree
Hide file tree
Showing 2 changed files with 13 additions and 3 deletions.
14 changes: 12 additions & 2 deletions src/acepn.f90
Original file line number Diff line number Diff line change
Expand Up @@ -25,7 +25,7 @@ module acepn
contains

subroutine acephn(nendf,npend,nace,ndir,matd,tempd,iprint,mcnpx,&
ityp,suff,hk,izn,awn)
ityp,suff,hk,izn,awn,izaoption)
!-------------------------------------------------------------------
! Prepare ACE photo-nuclear files.
!-------------------------------------------------------------------
Expand All @@ -36,7 +36,7 @@ subroutine acephn(nendf,npend,nace,ndir,matd,tempd,iprint,mcnpx,&
use endf ! provides endf routines and variables
use acecm ! provides bachaa,eavl,ptleg2,pttab2
! externals
integer::nendf,npend,nace,ndir,matd,iprint,mcnpx,ityp
integer::nendf,npend,nace,ndir,matd,iprint,mcnpx,ityp,izaoption
real(kr)::tempd,suff
integer::izn(16)
real(kr)::awn(16)
Expand Down Expand Up @@ -1825,6 +1825,16 @@ subroutine acephn(nendf,npend,nace,ndir,matd,tempd,iprint,mcnpx,&

!--print and write the photonuclear file
zaid=za+suff
if (izaoption.eq.1.and.is.gt.0) then
zaid=za+300+100*is+suff
endif
if (izaoption.eq.1.and.za.eq.95242) then
if (is.eq.0) then
zaid=za+400+suff
elseif (is.eq.1) then
zaid=za+suff
endif
endif
if (mcnpx.eq.0) then
write(hz,'(f9.2,''u'')') zaid
else
Expand Down
2 changes: 1 addition & 1 deletion src/acer.f90
Original file line number Diff line number Diff line change
Expand Up @@ -479,7 +479,7 @@ subroutine acer
!--prepare photo-nuclear data
else if (iopt.eq.5) then
call acephn(nendf,npend,nace,ndir,matd,tempd,iprint,&
mcnpx,itype,suff,hk,izn,awn)
mcnpx,itype,suff,hk,izn,awn,izaoption)

!--print or edit ace files
else if (iopt.ge.7.and.iopt.le.8) then
Expand Down

0 comments on commit de6d9bd

Please sign in to comment.