diff --git a/src/DAMASK_interface.f90 b/src/DAMASK_interface.f90
index 1eed9aab5..b76813fe6 100644
--- a/src/DAMASK_interface.f90
+++ b/src/DAMASK_interface.f90
@@ -15,7 +15,9 @@
 #define PETSC_MINOR_MIN 10
 #define PETSC_MINOR_MAX 11
 module DAMASK_interface
-
+  use prec
+  use system_routines
+  
   implicit none
   private
   logical, public, protected :: &
@@ -48,18 +50,8 @@ contains
 !! information on computation to screen
 !--------------------------------------------------------------------------------------------------
 subroutine DAMASK_interface_init
-  use, intrinsic :: &
-    iso_fortran_env
-  use, intrinsic :: &
-    iso_c_binding
+  use, intrinsic :: iso_fortran_env
   use PETScSys
-  use prec, only: &
-    pReal
-  use system_routines, only: &
-    signalusr1_C, &
-    signalusr2_C, &
-    getHostName, &
-    getCWD
 
 #include <petsc/finclude/petscsys.h>
 #if defined(__GFORTRAN__) &&  __GNUC__<GCC_MIN
@@ -305,9 +297,6 @@ end subroutine DAMASK_interface_init
 !!        possibly converting relative arguments to absolut path
 !--------------------------------------------------------------------------------------------------
 subroutine setWorkingDirectory(workingDirectoryArg)
-  use system_routines, only: &
-    getCWD, &
-    setCWD
 
   character(len=*),  intent(in) :: workingDirectoryArg                                              !< working directory argument
   character(len=1024)           :: workingDirectory                                                 !< working directory argument
@@ -355,8 +344,6 @@ end function getSolverJobName
 !> @brief basename of geometry file with extension from command line arguments
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function getGeometryFile(geometryParameter)
-  use system_routines, only: &
-    getCWD
 
   character(len=1024), intent(in) :: geometryParameter
   logical                         :: file_exists
@@ -379,8 +366,6 @@ end function getGeometryFile
 !> @brief relative path of loadcase from command line arguments
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function getLoadCaseFile(loadCaseParameter)
-  use system_routines, only: &
-    getCWD
 
   character(len=1024), intent(in) :: loadCaseParameter
   logical                         :: file_exists
@@ -475,7 +460,6 @@ end function makeRelativePath
 !> @brief sets global variable SIGTERM to .true.
 !--------------------------------------------------------------------------------------------------
 subroutine catchSIGTERM(signal) bind(C)
-  use :: iso_c_binding
 
   integer(C_INT), value :: signal
   SIGTERM = .true.
@@ -500,7 +484,6 @@ end subroutine setSIGTERM
 !> @brief sets global variable SIGUSR1 to .true.
 !--------------------------------------------------------------------------------------------------
 subroutine catchSIGUSR1(signal) bind(C)
-  use :: iso_c_binding
 
   integer(C_INT), value :: signal
   SIGUSR1 = .true.
@@ -525,7 +508,6 @@ end subroutine setSIGUSR1
 !> @brief sets global variable SIGUSR2 to .true. if program receives SIGUSR2
 !--------------------------------------------------------------------------------------------------
 subroutine catchSIGUSR2(signal) bind(C)
-  use :: iso_c_binding
 
   integer(C_INT), value :: signal
   SIGUSR2 = .true.
diff --git a/src/Lambert.f90 b/src/Lambert.f90
index da0c0b3aa..c7b2c0d49 100644
--- a/src/Lambert.f90
+++ b/src/Lambert.f90
@@ -38,11 +38,7 @@
 !> Modeling and Simulations in Materials Science and Engineering 22, 075013 (2014).
 !--------------------------------------------------------------------------
 module Lambert
-  use prec, only: &
-    pReal
-  use math, only: &
-    PI
-  use future
+  use math
 
   implicit none
   private
@@ -73,9 +69,6 @@ contains
 !> @brief map from 3D cubic grid to 3D ball 
 !--------------------------------------------------------------------------
 function LambertCubeToBall(cube) result(ball)
-  use, intrinsic :: IEEE_ARITHMETIC
-  use prec, only: &
-    dEq0
  
   real(pReal), intent(in), dimension(3) :: cube
   real(pReal),             dimension(3) :: ball, LamXYZ, XYZ
@@ -132,13 +125,6 @@ end function LambertCubeToBall
 !> @brief map from 3D ball to 3D cubic grid  
 !--------------------------------------------------------------------------
 pure function LambertBallToCube(xyz) result(cube)
-  use, intrinsic :: IEEE_ARITHMETIC, only:&
-    IEEE_positive_inf, &
-    IEEE_value
-  use prec, only: &
-    dEq0
-   use math, only: &
-    math_clip
  
   real(pReal), intent(in), dimension(3) :: xyz
   real(pReal),             dimension(3) :: cube, xyz1, xyz3
diff --git a/src/system_routines.f90 b/src/system_routines.f90
index 170e78dee..0611c96db 100644
--- a/src/system_routines.f90
+++ b/src/system_routines.f90
@@ -6,7 +6,6 @@ module system_routines
   use, intrinsic :: ISO_C_Binding
  
   implicit none
-  private
   
   public :: &
     signalterm_C, &