diff --git a/applications/lofar2/model/rtdsp/fourier.py b/applications/lofar2/model/rtdsp/fourier.py
index c2aad788ce41fd176d56a18e522870fec25be625..bfdc16b160c6efc26fd93d772b6aed031de23a8a 100644
--- a/applications/lofar2/model/rtdsp/fourier.py
+++ b/applications/lofar2/model/rtdsp/fourier.py
@@ -34,7 +34,7 @@ from .utilities import c_rtol, c_atol, ceil_pow2, is_even
 # Time domain interpolation using DFT
 ###############################################################################
 
-def fourier_interpolate(HFfilter, Ncoefs):
+def fourier_interpolate(HFfilter, Ncoefs, verbosity=1):
     """Use Fourier interpolation to create final FIR filter coefs.
 
     HF contains filter transfer function for N points, in order 0 to fs. The
@@ -82,7 +82,8 @@ def fourier_interpolate(HFfilter, Ncoefs):
         # K + 1 + K = N values, because N is odd and K = N // 2
     hInterpolated = np.fft.ifft(HFextended)
     if np.allclose(hInterpolated.imag, np.zeros(Ncoefs), rtol=c_rtol, atol=c_atol):
-        print('hInterpolated.imag ~= 0')
+        if verbosity:
+            print('hInterpolated.imag ~= 0')
     else:
         print('WARNING: hInterpolated.imag != 0 (max(abs) = %e)' % np.max(np.abs(hInterpolated.imag)))
     return hInterpolated.real