diff --git a/src/DAMASK_interface.f90 b/src/DAMASK_interface.f90 index ebe625ef6..d7e29865b 100644 --- a/src/DAMASK_interface.f90 +++ b/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 diff --git a/src/prec.f90 b/src/prec.f90 index cba8a68ef..ff05c24e2 100644 --- a/src/prec.f90 +++ b/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 diff --git a/src/source_damage_anisoBrittle.f90 b/src/source_damage_anisoBrittle.f90 index 717ea9fb4..da00db2b2 100644 --- a/src/source_damage_anisoBrittle.f90 +++ b/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 diff --git a/src/source_damage_anisoDuctile.f90 b/src/source_damage_anisoDuctile.f90 index f0774df62..8f6d68a88 100644 --- a/src/source_damage_anisoDuctile.f90 +++ b/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 diff --git a/src/source_damage_isoBrittle.f90 b/src/source_damage_isoBrittle.f90 index 2d13277c7..b60218458 100644 --- a/src/source_damage_isoBrittle.f90 +++ b/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 diff --git a/src/source_damage_isoDuctile.f90 b/src/source_damage_isoDuctile.f90 index 279d436c9..149803693 100644 --- a/src/source_damage_isoDuctile.f90 +++ b/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 diff --git a/src/source_thermal_dissipation.f90 b/src/source_thermal_dissipation.f90 index 23d6e8b58..94452eb47 100644 --- a/src/source_thermal_dissipation.f90 +++ b/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) :: & diff --git a/src/source_thermal_externalheat.f90 b/src/source_thermal_externalheat.f90 index 377513c57..699902ad3 100644 --- a/src/source_thermal_externalheat.f90 +++ b/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 diff --git a/src/system_routines.f90 b/src/system_routines.f90 index d7a27a4f9..066610111 100644 --- a/src/system_routines.f90 +++ b/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 diff --git a/src/thermal_adiabatic.f90 b/src/thermal_adiabatic.f90 index 8fdf55c97..bfc34d1c4 100644 --- a/src/thermal_adiabatic.f90 +++ b/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, & diff --git a/src/thermal_conduction.f90 b/src/thermal_conduction.f90 index 8dd352357..292d8a3f6 100644 --- a/src/thermal_conduction.f90 +++ b/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, & diff --git a/src/thermal_isothermal.f90 b/src/thermal_isothermal.f90 index a161094e7..c15dfa8f1 100644 --- a/src/thermal_isothermal.f90 +++ b/src/thermal_isothermal.f90 @@ -22,7 +22,6 @@ subroutine thermal_isothermal_init() material_Nhomogenization use material - implicit none integer :: & homog, & NofMyHomog