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

Interfacing for subroutines to identify real and integer data

This commit is contained in:
Vitesh Shah 2018-10-09 10:57:06 +02:00
parent a560fff2ac
commit 70c746a8f1
2 changed files with 276 additions and 229 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

200
src/CPFEM2.f90 Normal file → Executable file
View File

@ -133,7 +133,7 @@ subroutine CPFEM_init
implicit none implicit none
integer(pInt) :: k,l,m,ph,homog integer(pInt) :: k,l,m,ph,homog
character(len=1024) :: rankStr character(len=1024) :: rankStr
integer(HID_T) :: fileReadID, groupPlasticID integer(HID_T) :: fileReadID, groupPlasticID, groupHomogID
integer :: hdferr integer :: hdferr
mainProcess: if (worldrank == 0) then mainProcess: if (worldrank == 0) then
@ -154,62 +154,76 @@ subroutine CPFEM_init
fileReadID = HDF5_openFile(trim(getSolverJobName())//trim(rankStr)//'.hdf5') fileReadID = HDF5_openFile(trim(getSolverJobName())//trim(rankStr)//'.hdf5')
call HDF5_read(material_phase,fileReadID,'recordedPhase')
!write(6,*) material_phase
call HDF5_read(crystallite_F0,fileReadID,'convergedF') call HDF5_read(crystallite_F0,fileReadID,'convergedF')
call HDF5_read(crystallite_F0,fileReadID,'convergedFp') !write(6,*) crystallite_F0
call HDF5_read(crystallite_F0,fileReadID,'convergedFi') call HDF5_read(crystallite_Fp0,fileReadID,'convergedFp')
call HDF5_read(crystallite_F0,fileReadID,'convergedLp') call HDF5_read(crystallite_Fi0,fileReadID,'convergedFi')
call HDF5_read(crystallite_F0,fileReadID,'convergedLi') call HDF5_read(crystallite_Lp0,fileReadID,'convergedLp')
! call HDF5_readDataSet(fileReadID,'convergeddPdF') call HDF5_read(crystallite_Li0,fileReadID,'convergedLi')
! call HDF5_readDataSet(fileReadID,'convergedTstar') call HDF5_read(crystallite_dPdF0,fileReadID,'convergeddPdF')
call HDF5_read(crystallite_Tstar0_v,fileReadID,'convergedTstar')
! groupPlasticID = HDF5_openGroup2(fileReadID,'PlasticPhases') groupPlasticID = HDF5_openGroup2(fileReadID,'PlasticPhases')
! call HDF5_readDataSet(groupPlasticID,'convergedStateConst') write(6,*) groupPlasticID
do ph = 1_pInt,size(phase_plasticity)
call HDF5_read(plasticState(ph)%state0,groupPlasticID,'convergedStateConst')
write(6,*) plasticState(ph)%state0
enddo
call IO_read_intFile(777,'recordedPhase'//trim(rankStr),modelName,size(material_phase)) groupHomogID = HDF5_openGroup2(fileReadID,'material_Nhomogenization')
read (777,rec=1) material_phase; close (777) write(6,*) groupHomogID
do homog = 1_pInt, material_Nhomogenization
call HDF5_read(homogState(homog)%state0, groupHomogID,'convergedStateHomog')
write(6,*) homogState(homog)%state0
enddo
call IO_read_realFile(777,'convergedF'//trim(rankStr),modelName,size(crystallite_F0)) ! call IO_read_intFile(777,'recordedPhase'//trim(rankStr),modelName,size(material_phase))
read (777,rec=1) crystallite_F0; close (777) ! read (777,rec=1) material_phase; close (777)
call IO_read_realFile(777,'convergedFp'//trim(rankStr),modelName,size(crystallite_Fp0)) ! call IO_read_realFile(777,'convergedF'//trim(rankStr),modelName,size(crystallite_F0))
read (777,rec=1) crystallite_Fp0; close (777) ! read (777,rec=1) crystallite_F0; close (777)
call IO_read_realFile(777,'convergedFi'//trim(rankStr),modelName,size(crystallite_Fi0)) ! call IO_read_realFile(777,'convergedFp'//trim(rankStr),modelName,size(crystallite_Fp0))
read (777,rec=1) crystallite_Fi0; close (777) ! read (777,rec=1) crystallite_Fp0; close (777)
call IO_read_realFile(777,'convergedLp'//trim(rankStr),modelName,size(crystallite_Lp0)) ! call IO_read_realFile(777,'convergedFi'//trim(rankStr),modelName,size(crystallite_Fi0))
read (777,rec=1) crystallite_Lp0; close (777) ! read (777,rec=1) crystallite_Fi0; close (777)
call IO_read_realFile(777,'convergedLi'//trim(rankStr),modelName,size(crystallite_Li0)) ! call IO_read_realFile(777,'convergedLp'//trim(rankStr),modelName,size(crystallite_Lp0))
read (777,rec=1) crystallite_Li0; close (777) ! read (777,rec=1) crystallite_Lp0; close (777)
call IO_read_realFile(777,'convergeddPdF'//trim(rankStr),modelName,size(crystallite_dPdF0)) ! call IO_read_realFile(777,'convergedLi'//trim(rankStr),modelName,size(crystallite_Li0))
read (777,rec=1) crystallite_dPdF0; close (777) ! read (777,rec=1) crystallite_Li0; close (777)
call IO_read_realFile(777,'convergedTstar'//trim(rankStr),modelName,size(crystallite_Tstar0_v)) ! call IO_read_realFile(777,'convergeddPdF'//trim(rankStr),modelName,size(crystallite_dPdF0))
read (777,rec=1) crystallite_Tstar0_v; close (777) ! read (777,rec=1) crystallite_dPdF0; close (777)
call IO_read_realFile(777,'convergedStateConst'//trim(rankStr),modelName) ! call IO_read_realFile(777,'convergedTstar'//trim(rankStr),modelName,size(crystallite_Tstar0_v))
m = 0_pInt ! read (777,rec=1) crystallite_Tstar0_v; close (777)
readPlasticityInstances: do ph = 1_pInt, size(phase_plasticity)
do k = 1_pInt, plasticState(ph)%sizeState
do l = 1, size(plasticState(ph)%state0(1,:))
m = m+1_pInt
read(777,rec=m) plasticState(ph)%state0(k,l)
enddo; enddo
enddo readPlasticityInstances
close (777)
call IO_read_realFile(777,'convergedStateHomog'//trim(rankStr),modelName) ! call IO_read_realFile(777,'convergedStateConst'//trim(rankStr),modelName)
m = 0_pInt ! m = 0_pInt
readHomogInstances: do homog = 1_pInt, material_Nhomogenization ! readPlasticityInstances: do ph = 1_pInt, size(phase_plasticity)
do k = 1_pInt, homogState(homog)%sizeState ! do k = 1_pInt, plasticState(ph)%sizeState
do l = 1, size(homogState(homog)%state0(1,:)) ! do l = 1, size(plasticState(ph)%state0(1,:))
m = m+1_pInt ! m = m+1_pInt
read(777,rec=m) homogState(homog)%state0(k,l) ! read(777,rec=m) plasticState(ph)%state0(k,l)
enddo; enddo ! enddo; enddo
enddo readHomogInstances ! enddo readPlasticityInstances
close (777) ! close (777)
! call IO_read_realFile(777,'convergedStateHomog'//trim(rankStr),modelName)
! m = 0_pInt
! readHomogInstances: do homog = 1_pInt, material_Nhomogenization
! do k = 1_pInt, homogState(homog)%sizeState
! do l = 1, size(homogState(homog)%state0(1,:))
! m = m+1_pInt
! read(777,rec=m) homogState(homog)%state0(k,l)
! enddo; enddo
! enddo readHomogInstances
! close (777)
restartRead = .false. restartRead = .false.
endif endif
@ -270,7 +284,8 @@ subroutine CPFEM_age()
HDF5_closeFile, & HDF5_closeFile, &
HDF5_closeGroup, & HDF5_closeGroup, &
HDF5_addGroup2, & HDF5_addGroup2, &
HDF5_writeScalarDataset3 !HDF5_writeScalarDataset3, &
HDF5_write
!HDF5_addScalarDataset2 !HDF5_addScalarDataset2
use hdf5 use hdf5
use DAMASK_interface, only: & use DAMASK_interface, only: &
@ -321,73 +336,72 @@ if (restartWrite) then
fileHandle = HDF5_createFile(trim(getSolverJobName())//trim(rankStr)//'.hdf5') fileHandle = HDF5_createFile(trim(getSolverJobName())//trim(rankStr)//'.hdf5')
!call HDF5_writeScalarDataset3(fileHandle,testdata,'test',shape(testdata)) call HDF5_write(material_phase,fileHandle,'recordedPhase')
call HDF5_writeScalarDataset3(fileHandle,real(material_phase,pReal),'recordedPhase',shape(material_phase)) call HDF5_write(crystallite_F0,fileHandle,'convergedF')
call HDF5_writeScalarDataset3(fileHandle,crystallite_F0,'convergedF',shape(crystallite_F0)) call HDF5_write(crystallite_Fp0,fileHandle,'convergedFp')
call HDF5_writeScalarDataset3(fileHandle,crystallite_Fp0,'convergedFp',shape(crystallite_Fp0)) call HDF5_write(crystallite_Fi0,fileHandle,'convergedFi')
call HDF5_writeScalarDataset3(fileHandle,crystallite_Fi0,'convergedFi',shape(crystallite_Fi0)) call HDF5_write(crystallite_Lp0,fileHandle,'convergedLp')
call HDF5_writeScalarDataset3(fileHandle,crystallite_Lp0,'convergedLp',shape(crystallite_Lp0)) call HDF5_write(crystallite_Li0,fileHandle,'convergedLi')
call HDF5_writeScalarDataset3(fileHandle,crystallite_Li0,'convergedLi',shape(crystallite_Li0)) call HDF5_write(crystallite_dPdF0,fileHandle,'convergeddPdF')
call HDF5_writeScalarDataset3(fileHandle,crystallite_dPdF0,'convergeddPdF',shape(crystallite_dPdF0)) call HDF5_write(crystallite_Tstar0_v,fileHandle,'convergedTstar')
call HDF5_writeScalarDataset3(fileHandle,crystallite_Tstar0_v,'convergedTstar',shape(crystallite_Tstar0_v))
groupPlastic = HDF5_addGroup2(fileHandle,'PlasticPhases') groupPlastic = HDF5_addGroup2(fileHandle,'PlasticPhases')
do ph = 1_pInt,size(phase_plasticity) do ph = 1_pInt,size(phase_plasticity)
call HDF5_writeScalarDataset3(groupPlastic,plasticState(ph)%state0,'convergedStateConst',shape(plasticState(ph)%state0)) call HDF5_write(plasticState(ph)%state0,groupPlastic,'convergedStateConst')
enddo enddo
groupHomog = HDF5_addGroup2(fileHandle,'material_Nhomogenization') groupHomog = HDF5_addGroup2(fileHandle,'material_Nhomogenization')
do homog = 1_pInt, material_Nhomogenization do homog = 1_pInt, material_Nhomogenization
call HDF5_writeScalarDataset3(groupHomog,homogState(homog)%state0,'convergedStateHomog',shape(homogState(homog)%state0)) call HDF5_write(homogState(homog)%state0,groupHomog,'convergedStateHomog')
enddo enddo
call HDF5_closeFile(fileHandle) call HDF5_closeFile(fileHandle)
call IO_write_jobRealFile(777,'recordedPhase'//trim(rankStr),size(material_phase)) ! call IO_write_jobRealFile(777,'recordedPhase'//trim(rankStr),size(material_phase))
write (777,rec=1) material_phase; close (777) ! write (777,rec=1) material_phase; close (777)
call IO_write_jobRealFile(777,'convergedF'//trim(rankStr),size(crystallite_F0)) ! call IO_write_jobRealFile(777,'convergedF'//trim(rankStr),size(crystallite_F0))
write (777,rec=1) crystallite_F0; close (777) ! write (777,rec=1) crystallite_F0; close (777)
call IO_write_jobRealFile(777,'convergedFp'//trim(rankStr),size(crystallite_Fp0)) ! call IO_write_jobRealFile(777,'convergedFp'//trim(rankStr),size(crystallite_Fp0))
write (777,rec=1) crystallite_Fp0; close (777) ! write (777,rec=1) crystallite_Fp0; close (777)
call IO_write_jobRealFile(777,'convergedFi'//trim(rankStr),size(crystallite_Fi0)) ! call IO_write_jobRealFile(777,'convergedFi'//trim(rankStr),size(crystallite_Fi0))
write (777,rec=1) crystallite_Fi0; close (777) ! write (777,rec=1) crystallite_Fi0; close (777)
call IO_write_jobRealFile(777,'convergedLp'//trim(rankStr),size(crystallite_Lp0)) ! call IO_write_jobRealFile(777,'convergedLp'//trim(rankStr),size(crystallite_Lp0))
write (777,rec=1) crystallite_Lp0; close (777) ! write (777,rec=1) crystallite_Lp0; close (777)
call IO_write_jobRealFile(777,'convergedLi'//trim(rankStr),size(crystallite_Li0)) ! call IO_write_jobRealFile(777,'convergedLi'//trim(rankStr),size(crystallite_Li0))
write (777,rec=1) crystallite_Li0; close (777) ! write (777,rec=1) crystallite_Li0; close (777)
call IO_write_jobRealFile(777,'convergeddPdF'//trim(rankStr),size(crystallite_dPdF0)) ! call IO_write_jobRealFile(777,'convergeddPdF'//trim(rankStr),size(crystallite_dPdF0))
write (777,rec=1) crystallite_dPdF0; close (777) ! write (777,rec=1) crystallite_dPdF0; close (777)
call IO_write_jobRealFile(777,'convergedTstar'//trim(rankStr),size(crystallite_Tstar0_v)) ! call IO_write_jobRealFile(777,'convergedTstar'//trim(rankStr),size(crystallite_Tstar0_v))
write (777,rec=1) crystallite_Tstar0_v; close (777) ! write (777,rec=1) crystallite_Tstar0_v; close (777)
call IO_write_jobRealFile(777,'convergedStateConst'//trim(rankStr)) ! call IO_write_jobRealFile(777,'convergedStateConst'//trim(rankStr))
m = 0_pInt ! m = 0_pInt
writePlasticityInstances: do ph = 1_pInt, size(phase_plasticity) ! writePlasticityInstances: do ph = 1_pInt, size(phase_plasticity)
do k = 1_pInt, plasticState(ph)%sizeState ! do k = 1_pInt, plasticState(ph)%sizeState
do l = 1, size(plasticState(ph)%state0(1,:)) ! do l = 1, size(plasticState(ph)%state0(1,:))
m = m+1_pInt ! m = m+1_pInt
write(777,rec=m) plasticState(ph)%state0(k,l) ! write(777,rec=m) plasticState(ph)%state0(k,l)
enddo; enddo ! enddo; enddo
enddo writePlasticityInstances ! enddo writePlasticityInstances
close (777) ! close (777)
call IO_write_jobRealFile(777,'convergedStateHomog'//trim(rankStr)) ! call IO_write_jobRealFile(777,'convergedStateHomog'//trim(rankStr))
m = 0_pInt ! m = 0_pInt
writeHomogInstances: do homog = 1_pInt, material_Nhomogenization ! writeHomogInstances: do homog = 1_pInt, material_Nhomogenization
do k = 1_pInt, homogState(homog)%sizeState ! do k = 1_pInt, homogState(homog)%sizeState
do l = 1, size(homogState(homog)%state0(1,:)) ! do l = 1, size(homogState(homog)%state0(1,:))
m = m+1_pInt ! m = m+1_pInt
write(777,rec=m) homogState(homog)%state0(k,l) ! write(777,rec=m) homogState(homog)%state0(k,l)
enddo; enddo ! enddo; enddo
enddo writeHomogInstances ! enddo writeHomogInstances
close (777) ! close (777)
endif endif

313
src/HDF5_utilities.f90 Normal file → Executable file
View File

@ -11,10 +11,15 @@ module HDF5_Utilities
integer(pInt), private :: currentInc integer(pInt), private :: currentInc
interface HDF5_read interface HDF5_read
module procedure HDF5_read_double_1 module procedure HDF5_read_pReal
module procedure HDF5_read_double_5 module procedure HDF5_read_pInt
end interface HDF5_read end interface HDF5_read
interface HDF5_write
module procedure HDF5_write_pReal
module procedure HDF5_write_pInt
end interface HDF5_write
public :: & public :: &
HDF5_Utilities_init, & HDF5_Utilities_init, &
HDF5_mappingPhase, & HDF5_mappingPhase, &
@ -36,9 +41,10 @@ module HDF5_Utilities
HDF5_removeLink, & HDF5_removeLink, &
HDF5_createFile, & HDF5_createFile, &
HDF5_closeFile, & HDF5_closeFile, &
HDF5_writeScalarDataset3, & !HDF5_writeScalarDataset3, &
HDF5_addGroup2, HDF5_read, & HDF5_addGroup2, &
HDF5_openFile HDF5_openFile, &
HDF5_read, HDF5_write
contains contains
subroutine HDF5_Utilities_init subroutine HDF5_Utilities_init
@ -75,7 +81,8 @@ subroutine HDF5_createJobFile
#endif #endif
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! initialize HDF5 library and check if integer and float type size match !initialize HDF5 library and check if integer and float type size match
call h5open_f(hdferr) call h5open_f(hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_Utilities_init: h5open_f') if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_Utilities_init: h5open_f')
call h5tget_size_f(H5T_NATIVE_INTEGER,typeSize, hdferr) call h5tget_size_f(H5T_NATIVE_INTEGER,typeSize, hdferr)
@ -129,7 +136,7 @@ end subroutine HDF5_createJobFile
#endif #endif
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! open file ! create a file
!call h5fcreate_f(path,H5F_ACC_TRUNC_F,resultsFile,hdferr) !call h5fcreate_f(path,H5F_ACC_TRUNC_F,resultsFile,hdferr)
call h5fcreate_f(path,H5F_ACC_TRUNC_F,HDF5_createFile,hdferr,access_prp = plist_id) call h5fcreate_f(path,H5F_ACC_TRUNC_F,HDF5_createFile,hdferr,access_prp = plist_id)
if (hdferr < 0) call IO_error(100_pInt,ext_msg=path) if (hdferr < 0) call IO_error(100_pInt,ext_msg=path)
@ -140,7 +147,7 @@ end function HDF5_createFile
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates and initializes HDF5 output file !> @brief close the opened HDF5 output file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_closeJobFile() subroutine HDF5_closeJobFile()
use hdf5 use hdf5
@ -172,7 +179,7 @@ integer(HID_T) function HDF5_openFile(filePath)
end function HDF5_openFile end function HDF5_openFile
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates and initializes HDF5 output file !> @brief close the opened HDF5 output file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_closeFile(fileHandle) subroutine HDF5_closeFile(fileHandle)
use hdf5 use hdf5
@ -188,7 +195,7 @@ subroutine HDF5_closeFile(fileHandle)
end subroutine HDF5_closeFile end subroutine HDF5_closeFile
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file, or if loc is present at the given location !> @brief adds a new group to the results file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
integer(HID_T) function HDF5_addGroup(path) integer(HID_T) function HDF5_addGroup(path)
use hdf5 use hdf5
@ -204,7 +211,7 @@ end function HDF5_addGroup
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file, or if loc is present at the given location !> @brief adds a new group to the fileHandle (additional to addGroup2)
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
integer(HID_T) function HDF5_addGroup2(fileHandle,groupName) integer(HID_T) function HDF5_addGroup2(fileHandle,groupName)
use hdf5 use hdf5
@ -220,7 +227,7 @@ end function HDF5_addGroup2
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file !> @brief open a group from the results file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
integer(HID_T) function HDF5_openGroup(path) integer(HID_T) function HDF5_openGroup(path)
use hdf5 use hdf5
@ -235,7 +242,7 @@ integer(HID_T) function HDF5_openGroup(path)
end function HDF5_openGroup end function HDF5_openGroup
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief open a existing group of the results file !> @brief open an existing group of the results file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
integer(HID_T) function HDF5_openGroup2(FileReadID,path) integer(HID_T) function HDF5_openGroup2(FileReadID,path)
use hdf5 use hdf5
@ -244,7 +251,7 @@ integer(HID_T) function HDF5_openGroup2(FileReadID,path)
character(len=*), intent(in) :: path character(len=*), intent(in) :: path
integer :: hdferr integer :: hdferr
integer(HID_T), intent(in) :: FileReadID integer(HID_T), intent(in) :: FileReadID
write(6,*) FileReadID,'hello';flush(6) !write(6,*) FileReadID,'hello';flush(6)
call h5gopen_f(FileReadID, trim(path), HDF5_openGroup2, hdferr) call h5gopen_f(FileReadID, trim(path), HDF5_openGroup2, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg = 'HDF5_openGroup2: h5gopen_f ('//trim(path)//')') if (hdferr < 0) call IO_error(1_pInt,ext_msg = 'HDF5_openGroup2: h5gopen_f ('//trim(path)//')')
@ -252,7 +259,7 @@ end function HDF5_openGroup2
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file !> @brief ???
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_setLink(path,link) subroutine HDF5_setLink(path,link)
use hdf5 use hdf5
@ -274,7 +281,7 @@ subroutine HDF5_setLink(path,link)
end subroutine HDF5_setLink end subroutine HDF5_setLink
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file !> @brief ???
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_removeLink(link) subroutine HDF5_removeLink(link)
use hdf5 use hdf5
@ -307,7 +314,7 @@ end subroutine HDF5_closeGroup
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief adds a new group to the results file !> @brief adds a StringAttribute to the results file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_addStringAttribute(entity,attrLabel,attrValue) subroutine HDF5_addStringAttribute(entity,attrLabel,attrValue)
use hdf5 use hdf5
@ -1131,7 +1138,7 @@ end subroutine HDF5_mappingCells
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief creates a new 3D Tensor dataset in the given group location !!!TODO: really necessary?
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_addTensor3DDataset(group,Nnodes,tensorSize,label,SIunit) subroutine HDF5_addTensor3DDataset(group,Nnodes,tensorSize,label,SIunit)
use hdf5 use hdf5
@ -1167,7 +1174,7 @@ subroutine HDF5_addTensor3DDataset(group,Nnodes,tensorSize,label,SIunit)
end subroutine HDF5_addTensor3DDataset end subroutine HDF5_addTensor3DDataset
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief creates a new vector dataset in the given group location !!!TODO: really necessary?
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_writeVectorDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset) subroutine HDF5_writeVectorDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset)
use hdf5 use hdf5
@ -1226,8 +1233,8 @@ subroutine HDF5_writeVectorDataset(group,dataset,label,SIunit,dataspace_size,mpi
end subroutine HDF5_writeVectorDataset end subroutine HDF5_writeVectorDataset
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief creates a new tensor dataset in the given group location
! by default, a 3x3 tensor is assumed ! by default, a 3x3 tensor is assumed !!!TODO: really necessary?
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_writeTensorDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset) subroutine HDF5_writeTensorDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset)
use hdf5 use hdf5
@ -1289,7 +1296,7 @@ subroutine HDF5_writeTensorDataset(group,dataset,label,SIunit,dataspace_size,mpi
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief adds a new vector dataset to the given group location
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_addVectorDataset(group,nnodes,vectorSize,label,SIunit) subroutine HDF5_addVectorDataset(group,nnodes,vectorSize,label,SIunit)
use hdf5 use hdf5
@ -1324,7 +1331,7 @@ end subroutine HDF5_addVectorDataset
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief writes to a new scalar dataset in the given group location
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_writeScalarDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset) subroutine HDF5_writeScalarDataset(group,dataset,label,SIunit,dataspace_size,mpiOffset)
use hdf5 use hdf5
@ -1380,142 +1387,170 @@ subroutine HDF5_writeScalarDataset(group,dataset,label,SIunit,dataspace_size,mpi
end subroutine HDF5_writeScalarDataset end subroutine HDF5_writeScalarDataset
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief interfaced subroutine for reading dataset of the type pReal
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_writeScalarDataset3(group,dataset,label,myShape) subroutine HDF5_read_pReal(D,ID1,label)
use hdf5 use hdf5
implicit none implicit none
integer(HID_T), intent(in) :: group real(pReal), dimension(..), intent(out) :: D
integer(pInt), intent(in), dimension(:) :: myShape
integer(HID_T), intent(in) :: ID1
character(len=*), intent(in) :: label character(len=*), intent(in) :: label
real(pReal), intent(in), dimension(*) :: dataset
integer :: hdferr
integer(HID_T) :: dset_id, space_id
!-------------------------------------------------------------------------------------------------- call HDF5_read_pReal2(D,ID1,label,shape(D))
! create dataspace contains
call h5screate_simple_f(size(myShape), int(myShape,HSIZE_T), space_id, hdferr, &
int(myShape,HSIZE_T))
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5screate_simple_f')
!-------------------------------------------------------------------------------------------------- subroutine HDF5_read_pReal2(D,ID2,label,myShape)
! create Dataset
call h5dcreate_f(group, trim(label), H5T_NATIVE_DOUBLE, space_id, dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5dcreate_f')
CALL h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE,dataset,int(myShape,HSIZE_T), hdferr)
!--------------------------------------------------------------------------------------------------
!close types, dataspaces
call h5dclose_f(dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5dclose_f')
call h5sclose_f(space_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5sclose_f')
end subroutine HDF5_writeScalarDataset3
!--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location
!--------------------------------------------------------------------------------------------------
! subroutine HDF5_writeScalarDatasetLoop(group,dataset,label,myShape)
! use hdf5
! implicit none
! integer(HID_T), intent(in) :: group
! integer(pInt), intent(in), dimension(:) :: myShape
! character(len=*), intent(in) :: label
! real(pReal), intent(in), dimension(*) :: dataset
! integer :: hdferr
! integer(HID_T) :: dset_id, space_id
! !--------------------------------------------------------------------------------------------------
! ! create dataspace
! call h5screate_simple_f(size(myShape), int(myShape,HSIZE_T), space_id, hdferr, &
! int(myShape,HSIZE_T))
! if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5screate_simple_f')
! !--------------------------------------------------------------------------------------------------
! ! create Dataset
! call h5dcreate_f(group, trim(label), H5T_NATIVE_DOUBLE, space_id, dset_id, hdferr)
! if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5dcreate_f')
! CALL h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE,dataset,int(myShape,HSIZE_T), hdferr)
! !--------------------------------------------------------------------------------------------------
! !close types, dataspaces
! call h5dclose_f(dset_id, hdferr)
! if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5dclose_f')
! call h5sclose_f(space_id, hdferr)
! if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_writeScalarDataset3: h5sclose_f')
! end subroutine HDF5_writeScalarDatasetLoop
subroutine HDF5_read_double_1(D,ID,label)
implicit none
real(pReal), dimension(:), intent(out) :: D
integer(HID_T), intent(in) :: ID
character(len=*), intent(in) :: label
call HDF5_read_double_generic(D,ID,label,shape(D))
end subroutine HDF5_read_double_1
subroutine HDF5_read_double_5(D,ID,label)
implicit none
real(pReal), dimension(:,:,:,:,:), intent(out) :: D
integer(HID_T), intent(in) :: ID
character(len=*), intent(in) :: label
call HDF5_read_double_generic(D,ID,label,shape(D))
end subroutine HDF5_read_double_5
subroutine HDF5_read_double_generic(D,ID,label,myShape)
use hdf5 use hdf5
implicit none implicit none
real(pReal), dimension(*), intent(out) :: D real(pReal), dimension(*), intent(out) :: D
integer, dimension(:),intent(in) :: myShape integer(HID_T), intent(in) :: ID2
! real, dimension(:), allocatable :: D1 integer(pInt), dimension(:), intent(in) :: myShape
! real, dimension(:,:), allocatable :: D2
! real, dimension(:,:,:), allocatable :: D3
integer(HID_T), intent(in) :: ID
character(len=*), intent(in) :: label character(len=*), intent(in) :: label
integer(HID_T) :: dset_id integer(HID_T) :: dset_id
integer :: hdferr integer :: hdferr
call h5dopen_f(ID,label,dset_id,hdferr) call h5dopen_f(ID2,label,dset_id,hdferr)
call h5dread_f(dset_id,H5T_NATIVE_DOUBLE,D,int(myshape,HID_T),hdferr) call h5dread_f(dset_id,H5T_NATIVE_DOUBLE,D,int(myShape,HSIZE_T),hdferr)
end subroutine
end subroutine end subroutine
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief read datasets from a hdf5 file !> @brief interfaced subroutine for reading dataset of the type pInt
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! subroutine HDF5_read_pInt(D,ID1,label)
!subroutine HDF5_read_double_2(DataSet,FileReadID,NameDataSet) !,myshape) use hdf5
! use hdf5 implicit none
! integer(pInt), dimension(..), intent(out) :: D
! implicit none
! integer(HID_T), intent(in) :: FileReadID integer(HID_T), intent(in) :: ID1
! character(len=*), intent(in) :: NameDataSet character(len=*), intent(in) :: label
! real(pReal), intent(out), dimension(:,:,:,:,:) :: DataSet
! !integer(HSIZE_T), intent(in), dimension(:) :: myshape
! integer(HSIZE_T), dimension(:),allocatable :: myshape call HDF5_read_pInt2(D,ID1,label,shape(D))
! integer(HID_T) :: dset_id contains
! integer :: hdferr
! subroutine HDF5_read_pInt2(D,ID2,label,myShape)
! myShape = int(shape(DataSet),HSIZE_T) use hdf5
! call h5dopen_f(FileReadID,NameDataSet,dset_id,hdferr) implicit none
! call h5dread_f(dset_id,H5T_NATIVE_DOUBLE,DataSet,myshape,hdferr) integer(pInt), dimension(*), intent(out) :: D
! integer(HID_T), intent(in) :: ID2
!end subroutine integer(pInt), dimension(:), intent(in) :: myShape
character(len=*), intent(in) :: label
integer(HID_T) :: dset_id
integer :: hdferr
call h5dopen_f(ID2,label,dset_id,hdferr)
call h5dread_f(dset_id,H5T_NATIVE_INTEGER,D,int(myShape,HSIZE_T),hdferr)
end subroutine
end subroutine
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief interfaced subroutine for writing dataset of the type pReal
!--------------------------------------------------------------------------------------------------
subroutine HDF5_write_pReal(D,ID1,label)
use hdf5
implicit none
real(pReal), dimension(..), intent(out) :: D
integer(HID_T), intent(in) :: ID1
character(len=*), intent(in) :: label
call HDF5_write_pReal2(D,ID1,label,shape(D))
contains
subroutine HDF5_write_pReal2(D,ID2,label,myShape)
use hdf5
implicit none
real(pReal), dimension(*), intent(out) :: D
integer(HID_T), intent(in) :: ID2
integer(pInt), intent(in), dimension(:) :: myShape
character(len=*), intent(in) :: label
integer :: hdferr
integer(HID_T) :: dset_id, space_id
!--------------------------------------------------------------------------------------------------
! create dataspace
call h5screate_simple_f(size(myShape), int(myShape,HSIZE_T), space_id, hdferr, &
int(myShape,HSIZE_T))
if (hdferr < 0) call IO_error(1_pInt,ext_msg='write_pReal: h5screate_simple_f')
!--------------------------------------------------------------------------------------------------
! create Dataset
call h5dcreate_f(ID2, trim(label), H5T_NATIVE_DOUBLE, space_id, dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='write_pReal: h5dcreate_f')
CALL h5dwrite_f(dset_id, H5T_NATIVE_DOUBLE,D,int(myShape,HSIZE_T), hdferr)
!--------------------------------------------------------------------------------------------------
!close types, dataspaces
call h5dclose_f(dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_write_pReal: h5dclose_f')
call h5sclose_f(space_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_write_pReal: h5sclose_f')
end subroutine
end subroutine
!--------------------------------------------------------------------------------------------------
!> @brief interfaced subroutine for writing dataset of the type pInt
!--------------------------------------------------------------------------------------------------
subroutine HDF5_write_pInt(D,ID1,label)
use hdf5
implicit none
integer(pInt), dimension(..), intent(out) :: D
integer(HID_T), intent(in) :: ID1
character(len=*), intent(in) :: label
call HDF5_write_pInt2(D,ID1,label,shape(D))
contains
subroutine HDF5_write_pInt2(D,ID2,label,myShape)
use hdf5
implicit none
integer(pInt), dimension(*), intent(out) :: D
integer(HID_T), intent(in) :: ID2
integer(pInt), intent(in), dimension(:) :: myShape
character(len=*), intent(in) :: label
integer :: hdferr
integer(HID_T) :: dset_id, space_id
!--------------------------------------------------------------------------------------------------
! create dataspace
call h5screate_simple_f(size(myShape), int(myShape,HSIZE_T), space_id, hdferr, &
int(myShape,HSIZE_T))
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_write_pInt: h5screate_simple_f')
!--------------------------------------------------------------------------------------------------
! create Dataset
call h5dcreate_f(ID2, trim(label), H5T_NATIVE_INTEGER, space_id, dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_write_pInt: h5dcreate_f')
CALL h5dwrite_f(dset_id, H5T_NATIVE_INTEGER,D,int(myShape,HSIZE_T), hdferr)
!--------------------------------------------------------------------------------------------------
!close types, dataspaces
call h5dclose_f(dset_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_HDF5_write_pInt: h5dclose_f')
call h5sclose_f(space_id, hdferr)
if (hdferr < 0) call IO_error(1_pInt,ext_msg='HDF5_HDF5_write_pInt: h5sclose_f')
end subroutine
end subroutine
!--------------------------------------------------------------------------------------------------
!> @brief adds a new scalar dataset to the given group location
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_addScalarDataset(group,nnodes,label,SIunit) subroutine HDF5_addScalarDataset(group,nnodes,label,SIunit)
use hdf5 use hdf5
@ -1549,7 +1584,7 @@ subroutine HDF5_addScalarDataset(group,nnodes,label,SIunit)
end subroutine HDF5_addScalarDataset end subroutine HDF5_addScalarDataset
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief creates a new scalar dataset in the given group location !> @brief adds a new scalar dataset in the given group location !!!TODO: to be romoved
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
subroutine HDF5_addScalarDataset2(fileHandle,dataShape,label) subroutine HDF5_addScalarDataset2(fileHandle,dataShape,label)
use hdf5 use hdf5
@ -1581,8 +1616,6 @@ subroutine HDF5_addScalarDataset2(fileHandle,dataShape,label)
end subroutine HDF5_addScalarDataset2 end subroutine HDF5_addScalarDataset2
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
!> @brief copies the current temp results to the actual results file !> @brief copies the current temp results to the actual results file
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------