diff --git a/src/voronoiproj.f90 b/src/voronoiproj.f90
index 1f43bf7..b57f347 100644
--- a/src/voronoiproj.f90
+++ b/src/voronoiproj.f90
@@ -158,14 +158,14 @@ subroutine voronoiproj(leniw,lenrw,colg,nrow,rw,dres,goxd,dvxd,gozd,dvzd,depz,&
       ! using lsmr to solve for the projection coefficients
       !print*, 'LSMR beginning ...'
 
-      !nout = -1
-      nout = 36
-      open(nout,file='lsmrout_sub.txt')
+      nout = -1
+      !nout = 36
+      !open(nout,file='lsmrout_sub.txt')
 
       call LSMR(nd, ncells, leniwgp, lenrwgp,iwgp,rwgp,dres,damp,&
       atol, btol, conlim, itnlim, localSize,nout,&
       xunknown, istop, itn, anorm, acond,rnorm, arnorm, xnorm)
-      close(nout)
+      !close(nout)
       do ii = 1,ncells
         xunknown(ii) = xunknown(ii)/norm(ii)
       enddo