diff --git a/src/grid_mech_spectral_basic.f90 b/src/grid_mech_spectral_basic.f90
index f698fdc6d..2daebefbd 100644
--- a/src/grid_mech_spectral_basic.f90
+++ b/src/grid_mech_spectral_basic.f90
@@ -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.)
diff --git a/src/grid_mech_spectral_polarisation.f90 b/src/grid_mech_spectral_polarisation.f90
index 55df4455f..3bd30a360 100644
--- a/src/grid_mech_spectral_polarisation.f90
+++ b/src/grid_mech_spectral_polarisation.f90
@@ -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.)