need to read reference stiffness

This commit is contained in:
Martin Diehl 2019-04-16 14:53:36 +02:00
parent d0a4cbf8d5
commit ed8270d432
2 changed files with 6 additions and 0 deletions

View File

@ -116,6 +116,7 @@ subroutine grid_mech_spectral_basic_init
F ! pointer to solution data
PetscInt, dimension(worldsize) :: localK
integer(HID_T) :: fileHandle
integer :: fileUnit
character(len=1024) :: rankStr
write(6,'(/,a)') ' <<<+- grid_mech_spectral_basic init -+>>>'
@ -201,6 +202,8 @@ subroutine grid_mech_spectral_basic_init
call HDF5_read(fileHandle,C_volAvgLastInc,'C_volAvgLastInc')
call HDF5_closeFile(fileHandle)
fileUnit = IO_open_jobFile_binary('C_ref')
read(fileUnit) C_minMaxAvg; close(fileUnit)
endif restartRead
call utilities_updateGamma(C_minMaxAvg,.true.)

View File

@ -124,6 +124,7 @@ subroutine grid_mech_spectral_polarisation_init
F_tau ! specific (sub)pointer
PetscInt, dimension(worldsize) :: localK
integer(HID_T) :: fileHandle
integer :: fileUnit
character(len=1024) :: rankStr
write(6,'(/,a)') ' <<<+- grid_mech_spectral_polarisation init -+>>>'
@ -214,6 +215,8 @@ subroutine grid_mech_spectral_polarisation_init
call HDF5_read(fileHandle,C_volAvgLastInc,'C_volAvgLastInc')
call HDF5_closeFile(fileHandle)
fileUnit = IO_open_jobFile_binary('C_ref')
read(fileUnit) C_minMaxAvg; close(fileUnit)
endif restartRead
call utilities_updateGamma(C_minMaxAvg,.true.)