updating BB for SLBW/MLBW - fix j index and change it to call cssammy for reconstruction

Author: amlewis1578

src/reconr.f90

@@ -2611,19 +2611,19 @@ subroutine sigma(e,sig,sigd)
             if (mode.eq.0) then
                ! no res parameters
                call csnorp(e,sigp)
-            else if (mode.eq.1) then
+            else if (mode.eq.1.and.nmtres.eq.0) then
                ! single level breit-wigner
                call csslbw(e,sigp)
-            else if (mode.eq.2.and.tempr.eq.zero) then
+            else if (mode.eq.2.and.tempr.eq.zero.and.nmtres.eq.0) then
                ! multilevel breit-wigner
                call csmlbw(e,sigp)
-            else if (mode.eq.2.and.tempr.gt.zero) then
+            else if (mode.eq.2.and.tempr.gt.zero.and.nmtres.eq.0) then
                ! multilevel breit-wigner using g/h method
                call csmlbw2(e,sigp)
             else if (mode.eq.3.and.(isammy.eq.0.or.nmtres.eq.0)) then
                ! multi-level reich-moore (r-matrix)
                call csrmat(e,sigp)
-            else if (mode.eq.3.and.isammy.eq.1.and.nmtres.gt.0) then
+            else if (mode.le.3.and.isammy.eq.1.and.nmtres.gt.0) then
                ! sammy calculation
                call cssammy(e,sigp,siga,sigx,mmtres,nmtres,ncoef,nresp,i)
             else if (mode.eq.4) then

src/samm.f90

@@ -694,7 +694,7 @@ subroutine rdsammy(nin,ier,jnow,nro,naps,mode,el,eh,&
             lspin(i,j,ier)=0
             chspin(i,j,ier)=0
             backgr(i,j,ier)=0
-            backgrdata(i,igroup,ier,:)=0
+            backgrdata(i,j,ier,:)=0
          enddo
       enddo
       pari=1