diff --git a/src/msspec/spec/fortran/eig/common/diagonalize_matrix.f b/src/msspec/spec/fortran/eig/common/diagonalize_matrix.f
index da75446..38dc37d 100644
--- a/src/msspec/spec/fortran/eig/common/diagonalize_matrix.f
+++ b/src/msspec/spec/fortran/eig/common/diagonalize_matrix.f
@@ -17,7 +17,7 @@ c
 c Local variables
 c
       integer    :: lwork
-      complex*16 :: wquery
+      complex*16 :: wquery(1)
       complex*16 :: vl(1,1), vr(1,1)
 c
       real*8,     allocatable :: rwork(:)
@@ -39,7 +39,7 @@ c
          write(iuo1,*) '                        '
       end if
 c
-      lwork = int(wquery)
+      lwork = int(wquery(1))
       allocate(work(lwork))
 c
       call zgeev('n','n',n,a,lda,w,vl,1,vr,1,work,lwork,rwork,info)