Merge branch '407-discussion-about-restart-structure' into 'development'

remember at which increment restart happened

Closes #407

See merge request damask/DAMASK!1026

Author: hyharry

src/grid/DAMASK_grid.f90

@@ -17,6 +17,7 @@ program DAMASK_grid
   use parallelization
   use signal
   use CLI
+  use HDF5_utilities
   use IO
   use config
   use math
@@ -116,6 +117,7 @@ program DAMASK_grid
     solver
   character(len=:), allocatable :: &
     fileContent, fname
+  integer(HID_T) :: fileHandle
 
 
 !--------------------------------------------------------------------------------------------------
@@ -386,6 +388,9 @@ program DAMASK_grid
         call MPI_Allreduce(signal_SIGUSR2,sig,1_MPI_INTEGER_KIND,MPI_LOGICAL,MPI_LOR,MPI_COMM_WORLD,err_MPI)
         call parallelization_chkerr(err_MPI)
         if (mod(inc,loadCases(l)%f_restart) == 0 .or. sig) then
+          fileHandle = HDF5_openFile(getSolverJobName()//'_restart.hdf5','w')
+          call HDF5_addAttribute(fileHandle,'increment',totalIncsCounter)
+          call HDF5_closeFile(fileHandle)
           do field = 1, nActiveFields
             select case (ID(field))
               case(FIELD_MECH_ID)

src/grid/grid_mech_FEM.f90

@@ -435,7 +435,7 @@ subroutine grid_mechanical_FEM_restartWrite()
 
   print'(1x,a)', 'saving solver data required for restart'; flush(IO_STDOUT)
 
-  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','w')
+  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','a')
   groupHandle = HDF5_addGroup(fileHandle,'solver')
   call HDF5_write(reshape(F,[3,3,product(cells(1:2))*cells3]),groupHandle,'F')
   call HDF5_write(reshape(F_lastInc,[3,3,product(cells(1:2))*cells3]),groupHandle,'F_lastInc')

src/grid/grid_mech_spectral_Galerkin.f90

@@ -478,12 +478,13 @@ subroutine grid_mechanical_spectral_Galerkin_restartWrite()
   integer(HID_T) :: fileHandle, groupHandle
   real(pREAL), dimension(:,:,:,:), pointer :: F
 
+
   call DMDAVecGetArrayReadF90(DM_mech,F_PETSc,F,err_PETSc)
   CHKERRQ(err_PETSc)
 
   print'(1x,a)', 'saving solver data required for restart'; flush(IO_STDOUT)
 
-  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','w')
+  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','a')
   groupHandle = HDF5_addGroup(fileHandle,'solver')
   call HDF5_write(reshape(F,[3,3,product(cells(1:2))*cells3]),groupHandle,'F')
   call HDF5_write(reshape(F_lastInc,[3,3,product(cells(1:2))*cells3]),groupHandle,'F_lastInc')

src/grid/grid_mech_spectral_basic.f90

@@ -394,14 +394,15 @@ subroutine grid_mechanical_spectral_basic_restartWrite()
   integer(HID_T) :: fileHandle, groupHandle
   real(pREAL), dimension(:,:,:,:), pointer :: F
 
+
   call DMDAVecGetArrayReadF90(DM_mech,F_PETSc,F,err_PETSc)
   CHKERRQ(err_PETSc)
 
   if (num%update_gamma) C_minMaxAvgRestart = C_minMaxAvg
 
   print'(1x,a)', 'saving solver data required for restart'; flush(IO_STDOUT)
 
-  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','w')
+  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','a')
   groupHandle = HDF5_addGroup(fileHandle,'solver')
   call HDF5_write(reshape(F,[3,3,product(cells(1:2))*cells3]),groupHandle,'F')
   call HDF5_write(reshape(F_lastInc,[3,3,product(cells(1:2))*cells3]),groupHandle,'F_lastInc')

src/grid/grid_mech_spectral_polarization.f90

@@ -450,6 +450,7 @@ subroutine grid_mech_spectral_polarization_restartWrite()
   integer(HID_T) :: fileHandle, groupHandle
   real(pREAL), dimension(:,:,:,:), pointer :: FandF_tau, F, F_tau
 
+
   call DMDAVecGetArrayReadF90(DM_mech,FandF_tau_PETSc,FandF_tau,err_PETSc)
   CHKERRQ(err_PETSc)
   F     => FandF_tau(0: 8,:,:,:)
@@ -459,7 +460,7 @@ subroutine grid_mech_spectral_polarization_restartWrite()
 
   print '(1x,a)', 'saving solver data required for restart'; flush(IO_STDOUT)
 
-  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','w')
+  fileHandle  = HDF5_openFile(getSolverJobName()//'_restart.hdf5','a')
   groupHandle = HDF5_addGroup(fileHandle,'solver')
   call HDF5_write(reshape(F,[3,3,product(cells(1:2))*cells3]),groupHandle,'F')
   call HDF5_write(reshape(F_lastInc,[3,3,product(cells(1:2))*cells3]),groupHandle,'F_lastInc')