diff --git a/lib/lorentzian.f90 b/lib/lorentzian.f90
index 4ea349d7b..1eec9baf7 100644
--- a/lib/lorentzian.f90
+++ b/lib/lorentzian.f90
@@ -10,9 +10,9 @@ subroutine lorentzian(y,npts,a)
   real y(npts)
   real a(4)
   real deltaa(4)
-  real x(1000)
+  real x(2000)
 
-  if(npts.gt.1000) stop 'Error in lorentzian'
+  if(npts.gt.2000) stop 'Error in lorentzian'
 
   a=0.
   df=12000.0/8192.0                               !df = 1.465 Hz