diff --git a/src/grid/grid_mech_FEM.f90 b/src/grid/grid_mech_FEM.f90
index a34d880f7..5133d577a 100644
--- a/src/grid/grid_mech_FEM.f90
+++ b/src/grid/grid_mech_FEM.f90
@@ -12,8 +12,6 @@ module grid_mech_FEM
   use PETScdmda
   use PETScsnes
   use prec
-  use CPFEM2
-  use IO
   use debug
   use FEsolving
   use numerics
diff --git a/src/grid/grid_mech_spectral_basic.f90 b/src/grid/grid_mech_spectral_basic.f90
index f05f9bc93..b34034fa3 100644
--- a/src/grid/grid_mech_spectral_basic.f90
+++ b/src/grid/grid_mech_spectral_basic.f90
@@ -15,13 +15,11 @@ module grid_mech_spectral_basic
   use HDF5_utilities
   use math
   use spectral_utilities
-  use IO
   use FEsolving
   use config
   use numerics
   use homogenization
   use mesh_grid
-  use CPFEM2
   use debug
 
   implicit none
diff --git a/src/grid/grid_mech_spectral_polarisation.f90 b/src/grid/grid_mech_spectral_polarisation.f90
index 33c3e4e72..4986a3457 100644
--- a/src/grid/grid_mech_spectral_polarisation.f90
+++ b/src/grid/grid_mech_spectral_polarisation.f90
@@ -16,13 +16,11 @@ module grid_mech_spectral_polarisation
   use math
   use rotations
   use spectral_utilities
-  use IO
   use FEsolving
   use config
   use numerics
   use homogenization
   use mesh_grid
-  use CPFEM2
   use debug
  
   implicit none