diff --git a/src/grid/grid_mech_spectral_basic.f90 b/src/grid/grid_mech_spectral_basic.f90
index 70e764549..4bd9c5a96 100644
--- a/src/grid/grid_mech_spectral_basic.f90
+++ b/src/grid/grid_mech_spectral_basic.f90
@@ -375,7 +375,7 @@ end subroutine grid_mechanical_spectral_basic_forward
 
 
 !--------------------------------------------------------------------------------------------------
-!> @brief Update coordinates
+!> @brief Update coordinates.
 !--------------------------------------------------------------------------------------------------
 subroutine grid_mechanical_spectral_basic_updateCoords()
 
@@ -384,7 +384,7 @@ subroutine grid_mechanical_spectral_basic_updateCoords()
 
   call DMDAVecGetArrayReadF90(da,solution_vec,F,err_PETSc)
   CHKERRQ(err_PETSc)
-  call utilities_updateCoords(F)
+  call utilities_updateCoords(reshape(F,[3,3,size(F,2),size(F,3),size(F,4)]))
   call DMDAVecRestoreArrayReadF90(da,solution_vec,F,err_PETSc)
   CHKERRQ(err_PETSc)
 
@@ -392,7 +392,7 @@ end subroutine grid_mechanical_spectral_basic_updateCoords
 
 
 !--------------------------------------------------------------------------------------------------
-!> @brief Write current solver and constitutive data for restart to file
+!> @brief Write current solver and constitutive data for restart to file.
 !--------------------------------------------------------------------------------------------------
 subroutine grid_mechanical_spectral_basic_restartWrite()
 
diff --git a/src/grid/grid_mech_spectral_polarisation.f90 b/src/grid/grid_mech_spectral_polarisation.f90
index 168794316..60e7d676e 100644
--- a/src/grid/grid_mech_spectral_polarisation.f90
+++ b/src/grid/grid_mech_spectral_polarisation.f90
@@ -432,7 +432,7 @@ end subroutine grid_mechanical_spectral_polarisation_forward
 
 
 !--------------------------------------------------------------------------------------------------
-!> @brief Update coordinates
+!> @brief Update coordinates.
 !--------------------------------------------------------------------------------------------------
 subroutine grid_mechanical_spectral_polarisation_updateCoords()
 
@@ -441,7 +441,7 @@ subroutine grid_mechanical_spectral_polarisation_updateCoords()
 
   call DMDAVecGetArrayReadF90(da,solution_vec,FandF_tau,err_PETSc)
   CHKERRQ(err_PETSc)
-  call utilities_updateCoords(FandF_tau(0:8,:,:,:))
+  call utilities_updateCoords(reshape(FandF_tau(0:8,:,:,:),[3,3,size(FandF_tau,2),size(FandF_tau,3),size(FandF_tau,4)]))
   call DMDAVecRestoreArrayReadF90(da,solution_vec,FandF_tau,err_PETSc)
   CHKERRQ(err_PETSc)
 
@@ -449,7 +449,7 @@ end subroutine grid_mechanical_spectral_polarisation_updateCoords
 
 
 !--------------------------------------------------------------------------------------------------
-!> @brief Write current solver and constitutive data for restart to file
+!> @brief Write current solver and constitutive data for restart to file.
 !--------------------------------------------------------------------------------------------------
 subroutine grid_mechanical_spectral_polarisation_restartWrite()