one implicit none is enough

src/DAMASK_interface.f90

@@ -15,6 +15,7 @@
 #define PETSC_MINOR_MIN 10
 #define PETSC_MINOR_MAX 11
 module DAMASK_interface
+
   implicit none
   private
   logical, public, protected :: &
@@ -39,6 +40,7 @@ module DAMASK_interface
     getLoadCaseFile, &
     rectifyPath, &
     makeRelativePath
+
 contains
 
 !--------------------------------------------------------------------------------------------------
@@ -96,7 +98,6 @@ subroutine DAMASK_interface_init
 ===================================================================================================
 #endif
 
-  implicit none
   character(len=1024) :: &
     commandLine, &                                                                                  !< command line call as string
     arg, &                                                                                          !< individual argument
@@ -308,7 +309,6 @@ subroutine setWorkingDirectory(workingDirectoryArg)
     getCWD, &
     setCWD
 
-  implicit none
   character(len=*),  intent(in) :: workingDirectoryArg                                              !< working directory argument
   character(len=1024)           :: workingDirectory                                                 !< working directory argument
   logical                       :: error
@@ -336,7 +336,6 @@ end subroutine setWorkingDirectory
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function getSolverJobName()
 
-  implicit none
   integer :: posExt,posSep
   character(len=1024) :: tempString
 
@@ -363,7 +362,6 @@ character(len=1024) function getGeometryFile(geometryParameter)
   use system_routines, only: &
     getCWD
 
-  implicit none
   character(len=1024), intent(in) :: geometryParameter
   logical                         :: file_exists
   external                        :: quit
@@ -388,7 +386,6 @@ character(len=1024) function getLoadCaseFile(loadCaseParameter)
   use system_routines, only: &
     getCWD
 
-  implicit none
   character(len=1024), intent(in) :: loadCaseParameter
   logical                         :: file_exists
   external                        :: quit
@@ -412,7 +409,6 @@ end function getLoadCaseFile
 !--------------------------------------------------------------------------------------------------
 function rectifyPath(path)
 
-  implicit none
   character(len=*)    :: path
   character(len=1024) :: rectifyPath
   integer :: i,j,k,l
@@ -457,7 +453,6 @@ end function rectifyPath
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function makeRelativePath(a,b)
 
-  implicit none
   character (len=*), intent(in) :: a,b
   character (len=1024)          :: a_cleaned,b_cleaned
   integer :: i,posLastCommonSlash,remainingSlashes
@@ -486,7 +481,6 @@ end function makeRelativePath
 subroutine catchSIGTERM(signal) bind(C)
   use :: iso_c_binding
 
-  implicit none
   integer(C_INT), value :: signal
   SIGTERM = .true.
 
@@ -500,7 +494,6 @@ end subroutine catchSIGTERM
 !--------------------------------------------------------------------------------------------------
 subroutine setSIGTERM(state)
 
-  implicit none
   logical, intent(in) :: state
   SIGTERM = state
 
@@ -513,7 +506,6 @@ end subroutine setSIGTERM
 subroutine catchSIGUSR1(signal) bind(C)
   use :: iso_c_binding
 
-  implicit none
   integer(C_INT), value :: signal
   SIGUSR1 = .true.
 
@@ -527,7 +519,6 @@ end subroutine catchSIGUSR1
 !--------------------------------------------------------------------------------------------------
 subroutine setSIGUSR1(state)
 
-  implicit none
   logical, intent(in) :: state
   SIGUSR1 = state
 
@@ -540,7 +531,6 @@ end subroutine setSIGUSR1
 subroutine catchSIGUSR2(signal) bind(C)
   use :: iso_c_binding
 
-  implicit none
   integer(C_INT), value :: signal
   SIGUSR2 = .true.
 
@@ -554,7 +544,6 @@ end subroutine catchSIGUSR2
 !--------------------------------------------------------------------------------------------------
 subroutine setSIGUSR2(state)
 
-  implicit none
   logical, intent(in) :: state
   SIGUSR2 = state
 

src/prec.f90

@@ -102,7 +102,6 @@ contains
 !--------------------------------------------------------------------------------------------------
 subroutine prec_init
 
-  implicit none
   integer, allocatable, dimension(:) :: realloc_lhs_test
 
   external :: &
@@ -131,7 +130,6 @@ end subroutine prec_init
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function dEq(a,b,tol)
 
-  implicit none
   real(pReal), intent(in)           :: a,b
   real(pReal), intent(in), optional :: tol
   real(pReal)                       :: eps
@@ -155,7 +153,6 @@ end function dEq
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function dNeq(a,b,tol)
 
-  implicit none
   real(pReal), intent(in)           :: a,b
   real(pReal), intent(in), optional :: tol
   real(pReal)                       :: eps
@@ -179,7 +176,6 @@ end function dNeq
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function dEq0(a,tol)
 
-  implicit none
   real(pReal), intent(in)           :: a
   real(pReal), intent(in), optional :: tol
   real(pReal)                       :: eps
@@ -203,7 +199,6 @@ end function dEq0
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function dNeq0(a,tol)
 
-  implicit none
   real(pReal), intent(in)           :: a
   real(pReal), intent(in), optional :: tol
   real(pReal)                       :: eps
@@ -228,7 +223,6 @@ end function dNeq0
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function cEq(a,b,tol)
 
-  implicit none
   complex(pReal), intent(in)           :: a,b
   real(pReal),    intent(in), optional :: tol
   real(pReal)                          :: eps
@@ -253,7 +247,6 @@ end function cEq
 !--------------------------------------------------------------------------------------------------
 logical elemental pure function cNeq(a,b,tol)
 
-  implicit none
   complex(pReal), intent(in)           :: a,b
   real(pReal),    intent(in), optional :: tol
   real(pReal)                          :: eps

src/source_damage_anisoBrittle.f90

@@ -91,8 +91,6 @@ subroutine source_damage_anisoBrittle_init
    lattice_SchmidMatrix_cleavage, &
    lattice_maxNcleavageFamily
 
- implicit none
-
  integer(pInt) :: Ninstance,phase,instance,source,sourceOffset
  integer(pInt) :: NofMyPhase,p   ,i
  integer(pInt),          dimension(0), parameter :: emptyIntArray    = [integer(pInt)::]
@@ -219,7 +217,6 @@ subroutine source_damage_anisoBrittle_dotState(S, ipc, ip, el)
    lattice_maxNcleavageFamily, &
    lattice_NcleavageSystem
 
- implicit none
  integer(pInt), intent(in) :: &
    ipc, &                                                                                           !< component-ID of integration point
    ip, &                                                                                            !< integration point
@@ -279,7 +276,6 @@ subroutine source_damage_anisobrittle_getRateAndItsTangent(localphiDot, dLocalph
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent
@@ -307,7 +303,6 @@ function source_damage_anisoBrittle_postResults(phase, constituent)
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent

src/source_damage_anisoDuctile.f90

@@ -82,7 +82,6 @@ subroutine source_damage_anisoDuctile_init
    config_phase
 
    
- implicit none
  integer(pInt) :: Ninstance,phase,instance,source,sourceOffset
  integer(pInt) :: NofMyPhase,p ,i
 
@@ -194,7 +193,6 @@ subroutine source_damage_anisoDuctile_dotState(ipc, ip, el)
    damage, &
    damageMapping
 
- implicit none
  integer(pInt), intent(in) :: &
    ipc, &                                                                                           !< component-ID of integration point
    ip, &                                                                                            !< integration point
@@ -231,7 +229,6 @@ subroutine source_damage_anisoDuctile_getRateAndItsTangent(localphiDot, dLocalph
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent
@@ -259,7 +256,6 @@ function source_damage_anisoDuctile_postResults(phase, constituent)
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent

src/source_damage_isoBrittle.f90

@@ -74,7 +74,6 @@ subroutine source_damage_isoBrittle_init
    config_phase, &
    material_Nphase
 
- implicit none
 
  integer(pInt) :: Ninstance,phase,instance,source,sourceOffset
  integer(pInt) :: NofMyPhase,p,i   
@@ -176,7 +175,6 @@ subroutine source_damage_isoBrittle_deltaState(C, Fe, ipc, ip, el)
    math_sym33to6, &
    math_I3
 
- implicit none
  integer(pInt), intent(in) :: &
    ipc, &                                                                                           !< component-ID of integration point
    ip, &                                                                                            !< integration point
@@ -221,7 +219,6 @@ subroutine source_damage_isoBrittle_getRateAndItsTangent(localphiDot, dLocalphiD
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent
@@ -251,7 +248,6 @@ function source_damage_isoBrittle_postResults(phase, constituent)
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent

src/source_damage_isoDuctile.f90

@@ -74,7 +74,6 @@ subroutine source_damage_isoDuctile_init
    config_phase, &
    material_Nphase
 
- implicit none
 
  integer(pInt) :: Ninstance,phase,instance,source,sourceOffset
  integer(pInt) :: NofMyPhase,p,i
@@ -177,7 +176,6 @@ subroutine source_damage_isoDuctile_dotState(ipc, ip, el)
    damage, &
    damageMapping
 
- implicit none
  integer(pInt), intent(in) :: &
    ipc, &                                                                                           !< component-ID of integration point
    ip, &                                                                                            !< integration point
@@ -206,7 +204,6 @@ subroutine source_damage_isoDuctile_getRateAndItsTangent(localphiDot, dLocalphiD
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent
@@ -234,7 +231,6 @@ function source_damage_isoDuctile_postResults(phase, constituent)
  use material, only: &
    sourceState
 
- implicit none
  integer(pInt), intent(in) :: &
    phase, &
    constituent

src/source_thermal_dissipation.f90

@@ -56,7 +56,6 @@ subroutine source_thermal_dissipation_init
     config_phase, &
     material_Nphase
  
-  implicit none
   integer :: Ninstance,instance,source,sourceOffset
   integer :: NofMyPhase,p   
  
@@ -103,7 +102,6 @@ end subroutine source_thermal_dissipation_init
 !--------------------------------------------------------------------------------------------------
 subroutine source_thermal_dissipation_getRateAndItsTangent(TDot, dTDOT_dT, Tstar, Lp, phase)
 
-  implicit none
   integer, intent(in) :: &
     phase
   real(pReal),  intent(in), dimension(3,3) :: &

src/source_thermal_externalheat.f90

@@ -63,7 +63,6 @@ subroutine source_thermal_externalheat_init
     config_phase, &
     material_Nphase
  
-  implicit none
  
   integer :: maxNinstance,instance,source,sourceOffset,NofMyPhase,p   
  
@@ -120,7 +119,6 @@ subroutine source_thermal_externalheat_dotState(phase, of)
   use material, only: &
     sourceState
  
-  implicit none
   integer, intent(in) :: &
     phase, &
     of
@@ -140,7 +138,6 @@ subroutine source_thermal_externalheat_getRateAndItsTangent(TDot, dTDot_dT, phas
   use material, only: &
     sourceState
  
-  implicit none
   integer, intent(in) :: &
     phase, &
     of

src/system_routines.f90

@@ -81,7 +81,6 @@ contains
 !--------------------------------------------------------------------------------------------------
 logical function isDirectory(path)
 
- implicit none
  character(len=*), intent(in) :: path
  character(kind=C_CHAR), dimension(1024) :: strFixedLength                                         ! C string as array
  integer :: i 
@@ -100,7 +99,6 @@ end function isDirectory
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function getCWD()
 
- implicit none
  character(kind=C_CHAR), dimension(1024) :: charArray                                               ! C string is an array
  integer(C_INT) :: stat
  integer :: i
@@ -126,7 +124,7 @@ end function getCWD
 !> @brief gets the current host name
 !--------------------------------------------------------------------------------------------------
 character(len=1024) function getHostName()
- implicit none
+
  character(kind=C_CHAR), dimension(1024) :: charArray                                              ! C string is an array
  integer(C_INT) :: stat
  integer :: i
@@ -152,7 +150,7 @@ end function getHostName
 !> @brief changes the current working directory
 !--------------------------------------------------------------------------------------------------
 logical function setCWD(path)
- implicit none
+
  character(len=*), intent(in) :: path
  character(kind=C_CHAR), dimension(1024) :: strFixedLength                                          ! C string is an array
  integer :: i

src/thermal_adiabatic.f90

@@ -57,7 +57,6 @@ subroutine thermal_adiabatic_init
   use config, only: &
     config_homogenization
  
-  implicit none
   integer :: maxNinstance,section,instance,i,sizeState,NofMyHomog   
   character(len=65536),   dimension(0), parameter :: emptyStringArray = [character(len=65536)::]
   character(len=65536), dimension(:), allocatable :: outputs
@@ -124,7 +123,6 @@ function thermal_adiabatic_updateState(subdt, ip, el)
     temperatureRate, &
     thermalMapping
   
-  implicit none
   integer,     intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -181,7 +179,6 @@ subroutine thermal_adiabatic_getSourceAndItsTangent(Tdot, dTdot_dT, T, ip, el)
     crystallite_S, &
     crystallite_Lp  
  
-  implicit none
   integer,     intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -246,7 +243,6 @@ function thermal_adiabatic_getSpecificHeat(ip,el)
   use mesh, only: &
     mesh_element
  
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -282,7 +278,6 @@ function thermal_adiabatic_getMassDensity(ip,el)
   use mesh, only: &
     mesh_element
     
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -312,7 +307,6 @@ function thermal_adiabatic_postResults(homog,instance,of) result(postResults)
   use material, only: &
     temperature
  
-  implicit none
   integer, intent(in) :: &
     homog, &
     instance, &

src/thermal_conduction.f90

@@ -58,7 +58,6 @@ subroutine thermal_conduction_init
   use config, only: &
     config_homogenization
  
-  implicit none
   integer :: maxNinstance,section,instance,i
   integer :: sizeState
   integer :: NofMyHomog   
@@ -135,7 +134,6 @@ subroutine thermal_conduction_getSourceAndItsTangent(Tdot, dTdot_dT, T, ip, el)
     crystallite_S, &
     crystallite_Lp  
  
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -205,7 +203,6 @@ function thermal_conduction_getConductivity33(ip,el)
   use crystallite, only: &
     crystallite_push33ToRef
  
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -239,7 +236,6 @@ function thermal_conduction_getSpecificHeat(ip,el)
   use mesh, only: &
     mesh_element
  
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -273,7 +269,6 @@ function thermal_conduction_getMassDensity(ip,el)
   use mesh, only: &
     mesh_element
  
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -306,7 +301,6 @@ subroutine thermal_conduction_putTemperatureAndItsRate(T,Tdot,ip,el)
     temperatureRate, &
     thermalMapping
 
-  implicit none
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
@@ -332,7 +326,6 @@ function thermal_conduction_postResults(homog,instance,of) result(postResults)
   use material, only: &
     temperature
  
-  implicit none
   integer,              intent(in) :: &
     homog, &
     instance, &

src/thermal_isothermal.f90

@@ -22,7 +22,6 @@ subroutine thermal_isothermal_init()
    material_Nhomogenization
  use material
  
- implicit none
  integer :: &
    homog, &
    NofMyHomog