223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

commit spectral_interface

This commit is contained in:
Jaeyong Jung 2018-06-26 16:09:37 +02:00
parent 8aeb4a115e
commit a8c69dfcad
1 changed files with 15 additions and 150 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

165
src/spectral_interface.f90 Normal file → Executable file
View File

@ -195,17 +195,9 @@ subroutine DAMASK_interface_init()
call quit(1_pInt) call quit(1_pInt)
endif endif
! workingDirectory = trim(storeWorkingDirectory(trim(workingDirArg),trim(geometryArg))) workingDirectory = trim(storeWorkingDirectory(trim(workingDirArg)))
! geometryFile = getGeometryFile(geometryArg) geometryFile = getGeometryFile(geometryArg)
! loadCaseFile = getLoadCaseFile(loadCaseArg) loadCaseFile = getLoadCaseFile(loadCaseArg)
workingDirectory = trim(storeWorkingDirectory2(trim(workingDirArg),trim(geometryArg)))
geometryFile = getGeometryFile2(geometryArg,workingDirectory)
loadCaseFile = getLoadCaseFile2(loadCaseArg,workingDirectory)
! write(*,*) trim(workingDirectory)
! write(*,*) trim(workingDirectory)//'/' ! put '/' next to workingDirectory
call get_environment_variable('USER',userName) call get_environment_variable('USER',userName)
error = getHostName(hostName) error = getHostName(hostName)
@ -224,148 +216,21 @@ subroutine DAMASK_interface_init()
write(6,'(a,i6.6)') ' Restart from increment: ', spectralRestartInc write(6,'(a,i6.6)') ' Restart from increment: ', spectralRestartInc
write(6,'(a,l1,/)') ' Append to result file: ', appendToOutFile write(6,'(a,l1,/)') ' Append to result file: ', appendToOutFile
read(*,*)
end subroutine DAMASK_interface_init end subroutine DAMASK_interface_init
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief extract working directory from given argument or from location of geometry file, !> @brief extract working directory from given argument or from location of geometry file,
!! possibly converting relative arguments to absolut path !! possibly converting relative arguments to absolut path
!> @todo change working directory with call chdir(storeWorkingDirectory)? !> @todo change working directory with call chdir(storeWorkingDirectory)?
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
character(len=1024) function storeWorkingDirectory2(workingDirectoryArg,geometryArg) character(len=1024) function storeWorkingDirectory(workingDirectoryArg)
use system_routines, only: & use system_routines, only: &
isDirectory, & isDirectory, &
getCWD getCWD
implicit none implicit none
character(len=*), intent(in) :: workingDirectoryArg !< working directory argument character(len=*), intent(in) :: workingDirectoryArg !< working directory argument
character(len=*), intent(in) :: geometryArg !< geometry argument
character(len=1024) :: cwd
logical :: error
external :: quit
wdGiven: if (len(workingDirectoryArg)>0) then !< -d is given
absolutePath: if (workingDirectoryArg(1:1) == '/') then !< absolute path is given to workingDirectoryArg
storeWorkingDirectory2 = workingDirectoryArg
else absolutePath !< relative path is given
error = getCWD(cwd)
if (error) call quit(1_pInt)
storeWorkingDirectory2 = trim(cwd)//'/'//workingDirectoryArg !< add relative path to cwd
endif absolutePath
if (storeWorkingDirectory2(len(trim(storeWorkingDirectory2)):len(trim(storeWorkingDirectory2))) /= '/') &
storeWorkingDirectory2 = trim(storeWorkingDirectory2)//'/' ! if path seperator is not given, append it
else wdGiven !< -d is not given
error = getCWD(cwd)
if (error) call quit(1_pInt)
storeWorkingDirectory2 = trim(cwd)//'/'
! if (geometryArg(1:1) == '/') then ! absolute path given as command line argument
! storeWorkingDirectory2 = geometryArg(1:scan(geometryArg,'/',back=.true.))
! else
! error = getCWD(cwd) ! relative path given as command line argument
! if (error) call quit(1_pInt)
! storeWorkingDirectory2 = trim(cwd)//'/'//geometryArg(1:scan(geometryArg,'/',back=.true.)) !< workingDirectory should not depend on geometryArg
! endif
endif wdGiven
storeWorkingDirectory2 = trim(rectifyPath(storeWorkingDirectory2))
if(.not. isDirectory(trim(storeWorkingDirectory2))) then ! check if the directory exists
write(6,'(a20,a,a16)') ' working directory "',trim(storeWorkingDirectory2),'" does not exist'
call quit(1_pInt)
endif
end function storeWorkingDirectory2
!--------------------------------------------------------------------------------------------------
!> @brief basename of geometry file with extension from command line arguments
!--------------------------------------------------------------------------------------------------
character(len=1024) function getGeometryFile2(geometryParameter,workingDirectory)
use system_routines, only: &
getCWD
implicit none
character(len=1024), intent(in) :: &
geometryParameter
character(len=*), intent(in) :: workingDirectory !< working directory
! character(len=1024) :: &
! cwd
integer :: posExt, posSep
logical :: error
external :: quit
getGeometryFile2 = geometryParameter
posExt = scan(getGeometryFile2,'.',back=.true.)
posSep = scan(getGeometryFile2,'/',back=.true.)
if (posExt <= posSep) getGeometryFile2 = trim(getGeometryFile2)//('.geom') ! no extension present
if (scan(getGeometryFile2,'/') /= 1) then ! relative path given as command line argument
! error = getcwd(cwd) ! no more cwd
! cwd = workingDirectory
if (error) call quit(1_pInt)
getGeometryFile2 = rectifyPath(trim(workingDirectory)//'/'//getGeometryFile2)
else
getGeometryFile2 = rectifyPath(getGeometryFile2)
endif
getGeometryFile2 = makeRelativePath(getSolverWorkingDirectoryName(), getGeometryFile2)
end function getGeometryFile2
!--------------------------------------------------------------------------------------------------
!> @brief relative path of loadcase from command line arguments
!--------------------------------------------------------------------------------------------------
character(len=1024) function getLoadCaseFile2(loadCaseParameter,workingDirectory)
use system_routines, only: &
getCWD
implicit none
character(len=1024), intent(in) :: &
loadCaseParameter
character(len=*), intent(in) :: workingDirectory !< working directory
! character(len=1024) :: &
! cwd
integer :: posExt, posSep
logical :: error
external :: quit
getLoadCaseFile2 = loadcaseParameter
posExt = scan(getLoadCaseFile2,'.',back=.true.)
posSep = scan(getLoadCaseFile2,'/',back=.true.)
if (posExt <= posSep) getLoadCaseFile2 = trim(getLoadCaseFile2)//('.load') ! no extension present
if (scan(getLoadCaseFile2,'/') /= 1) then ! relative path given as command line argument
! error = getcwd(cwd)
! cwd = workingDirectory
if (error) call quit(1_pInt)
getLoadCaseFile2 = rectifyPath(trim(workingDirectory)//'/'//getLoadCaseFile2)
else
getLoadCaseFile2 = rectifyPath(getLoadCaseFile2)
endif
getLoadCaseFile2 = makeRelativePath(getSolverWorkingDirectoryName(), getLoadCaseFile2)
end function getLoadCaseFile2
!--------------------------------------------------------------------------------------------------
!> @brief extract working directory from given argument or from location of geometry file,
!! possibly converting relative arguments to absolut path
!> @todo change working directory with call chdir(storeWorkingDirectory)?
!--------------------------------------------------------------------------------------------------
character(len=1024) function storeWorkingDirectory(workingDirectoryArg,geometryArg)
use system_routines, only: &
isDirectory, &
getCWD
implicit none
character(len=*), intent(in) :: workingDirectoryArg !< working directory argument
character(len=*), intent(in) :: geometryArg !< geometry argument
character(len=1024) :: cwd character(len=1024) :: cwd
logical :: error logical :: error
external :: quit external :: quit
@ -381,13 +246,9 @@ character(len=1024) function storeWorkingDirectory(workingDirectoryArg,geometryA
if (storeWorkingDirectory(len(trim(storeWorkingDirectory)):len(trim(storeWorkingDirectory))) /= '/') & if (storeWorkingDirectory(len(trim(storeWorkingDirectory)):len(trim(storeWorkingDirectory))) /= '/') &
storeWorkingDirectory = trim(storeWorkingDirectory)//'/' ! if path seperator is not given, append it storeWorkingDirectory = trim(storeWorkingDirectory)//'/' ! if path seperator is not given, append it
else wdGiven else wdGiven
if (geometryArg(1:1) == '/') then ! absolute path given as command line argument
storeWorkingDirectory = geometryArg(1:scan(geometryArg,'/',back=.true.))
else
error = getCWD(cwd) ! relative path given as command line argument error = getCWD(cwd) ! relative path given as command line argument
if (error) call quit(1_pInt) if (error) call quit(1_pInt)
storeWorkingDirectory = trim(cwd)//'/'//geometryArg(1:scan(geometryArg,'/',back=.true.)) storeWorkingDirectory = trim(cwd)//'/'
endif
endif wdGiven endif wdGiven
storeWorkingDirectory = trim(rectifyPath(storeWorkingDirectory)) storeWorkingDirectory = trim(rectifyPath(storeWorkingDirectory))
@ -458,12 +319,15 @@ character(len=1024) function getGeometryFile(geometryParameter)
if (posExt <= posSep) getGeometryFile = trim(getGeometryFile)//('.geom') ! no extension present if (posExt <= posSep) getGeometryFile = trim(getGeometryFile)//('.geom') ! no extension present
if (scan(getGeometryFile,'/') /= 1) then ! relative path given as command line argument if (scan(getGeometryFile,'/') /= 1) then ! relative path given as command line argument
error = getcwd(cwd) ! error = getcwd(cwd)
if (error) call quit(1_pInt) ! if (error) call quit(1_pInt)
getGeometryFile = rectifyPath(trim(cwd)//'/'//getGeometryFile) ! getGeometryFile = rectifyPath(trim(workingDirectory)//'/'//getGeometryFile)
getGeometryFile = rectifyPath(trim(getSolverWorkingDirectoryName())//trim(getGeometryFile))
else else
getGeometryFile = rectifyPath(getGeometryFile) getGeometryFile = rectifyPath(getGeometryFile)
endif endif
write(*,*) 'getsolv.. ', (getSolverWorkingDirectoryName())
write(*,*) 'getGeometryFile.. ', (getGeometryFile)
getGeometryFile = makeRelativePath(getSolverWorkingDirectoryName(), getGeometryFile) getGeometryFile = makeRelativePath(getSolverWorkingDirectoryName(), getGeometryFile)
@ -492,9 +356,10 @@ character(len=1024) function getLoadCaseFile(loadCaseParameter)
if (posExt <= posSep) getLoadCaseFile = trim(getLoadCaseFile)//('.load') ! no extension present if (posExt <= posSep) getLoadCaseFile = trim(getLoadCaseFile)//('.load') ! no extension present
if (scan(getLoadCaseFile,'/') /= 1) then ! relative path given as command line argument if (scan(getLoadCaseFile,'/') /= 1) then ! relative path given as command line argument
error = getcwd(cwd) ! error = getcwd(cwd)
if (error) call quit(1_pInt) ! if (error) call quit(1_pInt)
getLoadCaseFile = rectifyPath(trim(cwd)//'/'//getLoadCaseFile) ! getLoadCaseFile = rectifyPath(trim(workingDirectory)//'/'//getLoadCaseFile)
getLoadCaseFile = rectifyPath(trim(getSolverWorkingDirectoryName())//trim(getLoadCaseFile))
else else
getLoadCaseFile = rectifyPath(getLoadCaseFile) getLoadCaseFile = rectifyPath(getLoadCaseFile)
endif endif