diff --git a/slycot/src/SB04OD.f b/slycot/src/SB04OD.f
index 929b1d65..d4ccfb0d 100644
--- a/slycot/src/SB04OD.f
+++ b/slycot/src/SB04OD.f
@@ -567,7 +567,8 @@ SUBROUTINE SB04OD( REDUCE, TRANS, JOBD, M, N, A, LDA, B, LDB, C,
 C     $                  LDD, DWORK, DWORK(M+1), DWORK(2*M+1), P, LDP, Q,
 C     $                  LDQ, DWORK(3*M+1), LDWORK-3*M, INFO )
          CALL DGGES( 'Vectors left', 'Vectors right', 'N', 0, N, A, LDA,
-     $               D, LDD, SDIM, DWORK, DWORK(M+1), DWORK(2*M+1), P, LDP, Q,
+     $               D, LDD, SDIM, DWORK, DWORK(M+1), DWORK(2*M+1), P,
+     $               LDP, Q,
      $               LDQ, DWORK(3*M+1), LDWORK-3*M, 0, INFO )
 
 C
@@ -626,8 +627,10 @@ SUBROUTINE SB04OD( REDUCE, TRANS, JOBD, M, N, A, LDA, B, LDB, C,
 C            CALL DGEGS( 'Vectors left', 'Vectors right', N, B, LDB, E,
 C     $                  LDE, DWORK, DWORK(N+1), DWORK(2*N+1), U, LDU, V,
 C     $                  LDV, DWORK(3*N+1), LDWORK-3*N, INFO )
-            CALL DGGES( 'Vectors left', 'Vectors right', 'N', 0, N, B, LDB, E,
-     $                  LDE, SDIM, DWORK, DWORK(N+1), DWORK(2*N+1), U, LDU, V,
+            CALL DGGES( 'Vectors left', 'Vectors right', 'N',
+     $                  0, N, B, LDB, E,
+     $                  LDE, SDIM, DWORK, DWORK(N+1), DWORK(2*N+1),
+     $                  U, LDU, V,
      $                  LDV, DWORK(3*N+1), LDWORK-3*N, 0, INFO )
 C
 C           Undo scaling