diff --git a/src/csf/sigma_vector.irp.f b/src/csf/sigma_vector.irp.f
index beb827ba..026dd839 100644
--- a/src/csf/sigma_vector.irp.f
+++ b/src/csf/sigma_vector.irp.f
@@ -51,16 +51,24 @@
       if(cfg_seniority_index(i+2) > ncfgpersomo) then
           ncfgpersomo = cfg_seniority_index(i+2)
       else
-        k = 0
-        if ((i+2+k) < size(cfg_seniority_index,1)) then
-          do while(cfg_seniority_index(i+2+k) < ncfgpersomo)
-            k = k + 2
-            if ((i+2+k) >= size(cfg_seniority_index,1)) then
-              exit
-            endif
-            ncfgpersomo = cfg_seniority_index(i+2+k)
-          enddo
-        endif
+        ! l = i+k+2 
+        ! Loop over l with a constraint to ensure that l <= size(cfg_seniority_index,1)-1
+        ! Old version commented just below
+        do l = min(size(cfg_seniority_index,1)-1, i+2), size(cfg_seniority_index,1)-1, 2
+          if (cfg_seniority_index(l) >= ncfgpersomo) then
+            ncfgpersomo = cfg_seniority_index(l)
+          endif
+        enddo
+        !k = 0
+        !if ((i+2+k) < size(cfg_seniority_index,1)) then
+        !  do while(cfg_seniority_index(i+2+k) < ncfgpersomo)
+        !    k = k + 2
+        !    if ((i+2+k) >= size(cfg_seniority_index,1)) then
+        !      exit
+        !    endif
+        !    ncfgpersomo = cfg_seniority_index(i+2+k)
+        !  enddo
+        !endif
       endif
     endif
     ncfg = ncfgpersomo - ncfgprev
@@ -79,16 +87,26 @@
     if(cfg_seniority_index(i+2) > ncfgprev) then
       ncfgprev = cfg_seniority_index(i+2)
     else
-      k = 0
-      if ((i+2+k) < size(cfg_seniority_index,1)) then
-        do while(cfg_seniority_index(i+2+k) < ncfgprev)
-          k = k + 2
-          if ((i+2+k) >= size(cfg_seniority_index,1)) then
-            exit
-          endif
-          ncfgprev = cfg_seniority_index(i+2+k)
-        enddo
-      endif
+      ! l = i+k+2 
+      ! Loop over l with a constraint to ensure that l <= size(cfg_seniority_index,1)-1
+      ! Old version commented just below
+      do l = min(size(cfg_seniority_index,1)-1, i+2), size(cfg_seniority_index,1)-1, 2
+        print*,'l'
+        if (cfg_seniority_index(l) >= ncfgprev) then
+          ncfgprev = cfg_seniority_index(l)
+        endif
+        print*,'ncfgprev', ncfgprev
+      enddo
+      !k = 0
+      !if ((i+2+k) < size(cfg_seniority_index,1)) then
+      !  do while(cfg_seniority_index(i+2+k) < ncfgprev)
+      !    k = k + 2
+      !    if ((i+2+k) >= size(cfg_seniority_index,1)) then
+      !      exit
+      !    endif
+      !    ncfgprev = cfg_seniority_index(i+2+k)
+      !  enddo
+      !endif
     endif
   enddo
 END_PROVIDER