diff --git a/base/serial/psb_crwextd.f90 b/base/serial/psb_crwextd.f90
index d801428a..7673a935 100644
--- a/base/serial/psb_crwextd.f90
+++ b/base/serial/psb_crwextd.f90
@@ -149,7 +149,7 @@ subroutine psb_cbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -159,8 +159,8 @@ subroutine psb_cbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
@@ -354,7 +354,7 @@ subroutine psb_lcbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -364,8 +364,8 @@ subroutine psb_lcbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
diff --git a/base/serial/psb_drwextd.f90 b/base/serial/psb_drwextd.f90
index 29898578..5a5efae3 100644
--- a/base/serial/psb_drwextd.f90
+++ b/base/serial/psb_drwextd.f90
@@ -149,7 +149,7 @@ subroutine psb_dbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -159,8 +159,8 @@ subroutine psb_dbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
@@ -354,7 +354,7 @@ subroutine psb_ldbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -364,8 +364,8 @@ subroutine psb_ldbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
diff --git a/base/serial/psb_srwextd.f90 b/base/serial/psb_srwextd.f90
index 6fc53334..f92f3970 100644
--- a/base/serial/psb_srwextd.f90
+++ b/base/serial/psb_srwextd.f90
@@ -149,7 +149,7 @@ subroutine psb_sbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -159,8 +159,8 @@ subroutine psb_sbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
@@ -354,7 +354,7 @@ subroutine psb_lsbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -364,8 +364,8 @@ subroutine psb_lsbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
diff --git a/base/serial/psb_zrwextd.f90 b/base/serial/psb_zrwextd.f90
index 949017bc..0ea0bf3d 100644
--- a/base/serial/psb_zrwextd.f90
+++ b/base/serial/psb_zrwextd.f90
@@ -149,7 +149,7 @@ subroutine psb_zbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -159,8 +159,8 @@ subroutine psb_zbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if
@@ -354,7 +354,7 @@ subroutine psb_lzbase_rwextd(nr,a,info,b,rowscale)
           end do
         end do
         do j=i,nr-ma
-          a%irp(ma+i+1) = a%irp(ma+i)
+          a%irp(ma+j+1) = a%irp(ma+j)
         end do
       class default 
 
@@ -364,8 +364,8 @@ subroutine psb_lzbase_rwextd(nr,a,info,b,rowscale)
 
     else
 
-      do i=ma+2,nr+1
-        a%irp(i) = a%irp(i-1)
+      do j=ma+2,nr+1
+        a%irp(j) = a%irp(j-1)
       end do
 
     end if