diff --git a/src/grid/spectral_utilities.f90 b/src/grid/spectral_utilities.f90
index 487b0816f..d4a4fbc09 100644
--- a/src/grid/spectral_utilities.f90
+++ b/src/grid/spectral_utilities.f90
@@ -398,18 +398,6 @@ subroutine utilities_updateGamma(C)
 end subroutine utilities_updateGamma
 
 
-!--------------------------------------------------------------------------------------------------
-!> @brief backward FFT of data in field_fourier to field_real
-!> @details Does an weighted inverse FFT transform from complex to real
-!--------------------------------------------------------------------------------------------------
-subroutine utilities_FFTvectorBackward()
-
-  call fftw_mpi_execute_dft_c2r(planVectorBack,vectorField_fourier,vectorField_real)
-  vectorField_real = vectorField_real * wgt                                                         ! normalize the result by number of elements
-
-end subroutine utilities_FFTvectorBackward
-
-
 !--------------------------------------------------------------------------------------------------
 !> @brief doing convolution gamma_hat * field_real, ensuring that average value = fieldAim
 !--------------------------------------------------------------------------------------------------