diff --git a/src/grid/DAMASK_grid.f90 b/src/grid/DAMASK_grid.f90 index f875d5941..573cba13e 100644 --- a/src/grid/DAMASK_grid.f90 +++ b/src/grid/DAMASK_grid.f90 @@ -76,7 +76,8 @@ program DAMASK_spectral myStat, & statUnit = 0, & !< file unit for statistics output lastRestartWritten = 0, & !< total increment # at which last restart information was written - stagIter + stagIter, & + nActiveFields = 0 character(len=6) :: loadcase_string character(len=1024) :: & incInfo diff --git a/src/grid/spectral_utilities.f90 b/src/grid/spectral_utilities.f90 index bb1c70f90..8f47f03ac 100644 --- a/src/grid/spectral_utilities.f90 +++ b/src/grid/spectral_utilities.f90 @@ -24,7 +24,6 @@ module spectral_utilities include 'fftw3-mpi.f03' logical, public :: cutBack = .false. !< cut back of BVP solver in case convergence is not achieved or a material point is terminally ill - integer, public :: nActiveFields = 0 !-------------------------------------------------------------------------------------------------- ! field labels information @@ -38,11 +37,13 @@ module spectral_utilities !-------------------------------------------------------------------------------------------------- ! grid related information information - real(pReal), public :: wgt !< weighting factor 1/Nelems + real(pReal), protected, public :: wgt !< weighting factor 1/Nelems + integer, protected, public :: grid1Red !< grid(1)/2 + real(pReal), protected, public, dimension(3) :: scaledGeomSize !< scaled geometry size for calculation of divergence !-------------------------------------------------------------------------------------------------- ! variables storing information for spectral method and FFTW - integer, public :: grid1Red !< grid(1)/2 + real (C_DOUBLE), public, dimension(:,:,:,:,:), pointer :: tensorField_real !< real representation (some stress or deformation) of field_fourier complex(C_DOUBLE_COMPLEX),public, dimension(:,:,:,:,:), pointer :: tensorField_fourier !< field on which the Fourier transform operates real(C_DOUBLE), public, dimension(:,:,:,:), pointer :: vectorField_real !< vector field real representation for fftw @@ -53,7 +54,7 @@ module spectral_utilities complex(pReal), private, dimension(:,:,:,:), allocatable :: xi1st !< wave vector field for first derivatives complex(pReal), private, dimension(:,:,:,:), allocatable :: xi2nd !< wave vector field for second derivatives real(pReal), private, dimension(3,3,3,3) :: C_ref !< mechanic reference stiffness - real(pReal), protected, public, dimension(3) :: scaledGeomSize !< scaled geometry size for calculation of divergence (Basic, Basic PETSc) + !-------------------------------------------------------------------------------------------------- ! plans for FFTW diff --git a/src/mesh/FEM_mech.f90 b/src/mesh/FEM_mech.f90 index eca81ab36..3a19f67c7 100644 --- a/src/mesh/FEM_mech.f90 +++ b/src/mesh/FEM_mech.f90 @@ -62,7 +62,7 @@ module FEM_mech contains !-------------------------------------------------------------------------------------------------- -!> @brief allocates all neccessary fields and fills them with data, potentially from restart info +!> @brief allocates all neccessary fields and fills them with data !-------------------------------------------------------------------------------------------------- subroutine FEM_mech_init(fieldBC) diff --git a/src/mesh/FEM_utilities.f90 b/src/mesh/FEM_utilities.f90 index 70d347b2f..1303f0df8 100644 --- a/src/mesh/FEM_utilities.f90 +++ b/src/mesh/FEM_utilities.f90 @@ -77,7 +77,6 @@ module FEM_utilities real(pReal) :: time = 0.0_pReal !< length of increment integer :: incs = 0, & !< number of increments outputfrequency = 1, & !< frequency of result writes - restartfrequency = 0, & !< frequency of restart writes logscale = 0 !< linear/logarithmic time inc flag logical :: followFormerTrajectory = .true. !< follow trajectory of former loadcase integer, allocatable :: faceID(:) @@ -145,7 +144,6 @@ subroutine utilities_constitutiveResponse(timeinc,P_av,forwardData) call materialpoint_stressAndItsTangent(.true.,timeinc) ! calculate P field - restartWrite = .false. ! reset restartWrite status cutBack = .false. ! reset cutBack status P_av = sum(sum(materialpoint_P,dim=4),dim=3) * wgt ! average of P