diff --git a/fvn_linear/fvn_linear.f90 b/fvn_linear/fvn_linear.f90
index 3ec48fd..8bfa9cf 100644
--- a/fvn_linear/fvn_linear.f90
+++ b/fvn_linear/fvn_linear.f90
@@ -1409,8 +1409,8 @@ bmat=reshape ( (/  (y(i),i=1,np)   /)  ,shape = (/ np,1 /))
 
 ! query workspace size
 call dgels('N',np,deg+1,1,mat,np,bmat,np,twork,-1,info)
-lwork=twork(1)
-allocate(work(int(lwork)))
+lwork=int(twork(1))
+allocate(work(lwork))
 ! real(kind=sp_kind) call
 call dgels('N',np,deg+1,1,mat,np,bmat,np,work,lwork,info)
 
@@ -1465,8 +1465,8 @@ bmat=reshape ( (/  (y(i),i=1,np)   /)  ,shape = (/ np,1 /))
 
 ! query workspace size
 call sgels('N',np,deg+1,1,mat,np,bmat,np,twork,-1,info)
-lwork=twork(1)
-allocate(work(int(lwork)))
+lwork=int(twork(1))
+allocate(work(lwork))
 ! real(kind=sp_kind) call
 call sgels('N',np,deg+1,1,mat,np,bmat,np,work,lwork,info)
 
@@ -1528,8 +1528,8 @@ bmat=reshape ( (/  (y(i),i=1,np)   /)  ,shape = (/ np,1 /))
 
 ! query workspace size
 call dgelss(np,deg+1,1,mat,np,bmat,np,singval,-1.,rank,twork,-1,info)
-lwork=twork(1)
-allocate(work(int(lwork)))
+lwork=int(twork(1))
+allocate(work(lwork))
 ! real(kind=sp_kind) call
 call dgelss(np,deg+1,1,mat,np,bmat,np,singval,-1.,rank,work,lwork,info)
 
@@ -1584,8 +1584,8 @@ bmat=reshape ( (/  (y(i),i=1,np)   /)  ,shape = (/ np,1 /))
 
 ! query workspace size
 call sgelss(np,deg+1,1,mat,np,bmat,np,singval,-1.,rank,twork,-1,info)
-lwork=twork(1)
-allocate(work(int(lwork)))
+lwork=int(twork(1))
+allocate(work(lwork))
 ! real(kind=sp_kind) call
 call sgelss(np,deg+1,1,mat,np,bmat,np,singval,-1.,rank,work,lwork,info)