Skip to content

Commit 6fd1aef

Browse files
committed
fixed index in sammy background array for LRF=3, and the overwriting of number of channels when there is a 2nd fission channel width at a lower energy than a 1st fission channel width
1 parent cec72ac commit 6fd1aef

File tree

2 files changed

+9
-5
lines changed

2 files changed

+9
-5
lines changed

src/reconr.f90

+2-2
Original file line numberDiff line numberDiff line change
@@ -147,9 +147,9 @@ subroutine reconr
147147

148148
!--set samrml options
149149
Want_Partial_Derivs=.false.
150-
Want_Angular_Dist=.false.
150+
Want_Angular_Dist=.true.
151151
Want_SAMRML_BW=.false.
152-
Want_SAMRML_RM=.false.
152+
Want_SAMRML_RM=.true.
153153

154154
!--define local i/o units
155155
nscr1=-10

src/samm.f90

+7-3
Original file line numberDiff line numberDiff line change
@@ -886,7 +886,7 @@ subroutine rdsammy(nin,ier,jnow,nro,naps,mode,el,eh,&
886886
lspin(i,j,ier)=0
887887
chspin(i,j,ier)=0
888888
backgr(i,j,ier)=0
889-
backgrdata(i,igroup,ier,:)=0
889+
backgrdata(i,j,ier,:)=0
890890
enddo
891891
enddo
892892
pari=1
@@ -1023,8 +1023,12 @@ subroutine rdsammy(nin,ier,jnow,nro,naps,mode,el,eh,&
10231023
gamf2=res(jj+5)
10241024
if (gamf.ne.zero) then
10251025
gamma(2,ires,ier)=gamf
1026-
if (nchan(ig,ier).lt.2) then
1027-
nchan(ig,ier)=2
1026+
if (nchan(ig,ier).ne.2) then
1027+
!-- if num channels is already set to 3,
1028+
!-- don't overwrite it with 2
1029+
if (nchan(ig,ier).ne.3) then
1030+
nchan(ig,ier)=2
1031+
endif
10281032
ipp(2,ig,ier)=3
10291033
endif
10301034
endif

0 commit comments

Comments
 (0)