diff --git a/src/grid/grid_mech_FEM.f90 b/src/grid/grid_mech_FEM.f90
index baf3e59e1..a479319c3 100644
--- a/src/grid/grid_mech_FEM.f90
+++ b/src/grid/grid_mech_FEM.f90
@@ -201,7 +201,7 @@ subroutine grid_mech_FEM_init
     F         = spread(spread(spread(math_I3,3,grid(1)),4,grid(2)),5,grid3)
   endif restartRead
   materialpoint_F0 = reshape(F_lastInc, [3,3,1,product(grid(1:2))*grid3])                           ! set starting condition for materialpoint_stressAndItsTangent
-  call utilities_updateIPcoords(F)
+  call utilities_updateCoords(F)
   call utilities_constitutiveResponse(P_current,temp33_Real,C_volAvg,devNull, &                     ! stress field, stress avg, global average of stiffness and (min+max)/2
                                       F, &                                                          ! target F
                                       0.0_pReal, &                                                  ! time increment
@@ -326,7 +326,7 @@ subroutine grid_mech_FEM_forward(guess,timeinc,timeinc_old,loadCaseTime,deformat
 
     endif
     call CPFEM_age                                                                                  ! age state and kinematics
-    call utilities_updateIPcoords(F)
+    call utilities_updateCoords(F)
 
     C_volAvgLastInc    = C_volAvg
  
diff --git a/src/grid/grid_mech_spectral_basic.f90 b/src/grid/grid_mech_spectral_basic.f90
index f0c564035..54190cc30 100644
--- a/src/grid/grid_mech_spectral_basic.f90
+++ b/src/grid/grid_mech_spectral_basic.f90
@@ -170,7 +170,7 @@ subroutine grid_mech_spectral_basic_init
   endif restartRead
  
   materialpoint_F0 = reshape(F_lastInc, [3,3,1,product(grid(1:2))*grid3])                           ! set starting condition for materialpoint_stressAndItsTangent
-  call Utilities_updateIPcoords(reshape(F,shape(F_lastInc)))
+  call Utilities_updateCoords(reshape(F,shape(F_lastInc)))
   call Utilities_constitutiveResponse(P,temp33_Real,C_volAvg,C_minMaxAvg, &                         ! stress field, stress avg, global average of stiffness and (min+max)/2
                                       reshape(F,shape(F_lastInc)), &                                ! target F
                                       0.0_pReal, &                                                  ! time increment
@@ -300,7 +300,7 @@ subroutine grid_mech_spectral_basic_forward(guess,timeinc,timeinc_old,loadCaseTi
     endif
 
     call CPFEM_age                                                                                  ! age state and kinematics
-    call utilities_updateIPcoords(F)
+    call utilities_updateCoords(F)
 
     C_volAvgLastInc    = C_volAvg
     C_minMaxAvgLastInc = C_minMaxAvg
diff --git a/src/grid/grid_mech_spectral_polarisation.f90 b/src/grid/grid_mech_spectral_polarisation.f90
index 0ceaac572..e0f12c9ee 100644
--- a/src/grid/grid_mech_spectral_polarisation.f90
+++ b/src/grid/grid_mech_spectral_polarisation.f90
@@ -183,7 +183,7 @@ subroutine grid_mech_spectral_polarisation_init
   endif restartRead
  
   materialpoint_F0 = reshape(F_lastInc, [3,3,1,product(grid(1:2))*grid3])                           ! set starting condition for materialpoint_stressAndItsTangent
-  call Utilities_updateIPcoords(reshape(F,shape(F_lastInc)))
+  call Utilities_updateCoords(reshape(F,shape(F_lastInc)))
   call Utilities_constitutiveResponse(P,temp33_Real,C_volAvg,C_minMaxAvg, &                         ! stress field, stress avg, global average of stiffness and (min+max)/2
                                       reshape(F,shape(F_lastInc)), &                                ! target F
                                       0.0_pReal, &                                                  ! time increment
@@ -324,7 +324,7 @@ subroutine grid_mech_spectral_polarisation_forward(guess,timeinc,timeinc_old,loa
     endif
 
     call CPFEM_age                                                                                  ! age state and kinematics
-    call utilities_updateIPcoords(F)
+    call utilities_updateCoords(F)
 
     C_volAvgLastInc    = C_volAvg
     C_minMaxAvgLastInc = C_minMaxAvg
diff --git a/src/grid/spectral_utilities.f90 b/src/grid/spectral_utilities.f90
index bf462cf1c..f19434415 100644
--- a/src/grid/spectral_utilities.f90
+++ b/src/grid/spectral_utilities.f90
@@ -157,7 +157,7 @@ module spectral_utilities
     utilities_constitutiveResponse, &
     utilities_calculateRate, &
     utilities_forwardField, &
-    utilities_updateIPcoords, &
+    utilities_updateCoords, &
     FIELD_UNDEFINED_ID, &
     FIELD_MECH_ID, &
     FIELD_THERMAL_ID, &
@@ -1024,7 +1024,7 @@ end function utilities_getFreqDerivative
 ! using integration in Fourier space. Similar as in mesh.f90, but using data already defined for
 ! convolution
 !--------------------------------------------------------------------------------------------------
-subroutine utilities_updateIPcoords(F)
+subroutine utilities_updateCoords(F)
  
   real(pReal),   dimension(3,3,grid(1),grid(2),grid3), intent(in) :: F
   real(pReal),   dimension(3,  grid(1),grid(2),grid3)             :: IPcoords
@@ -1069,6 +1069,6 @@ subroutine utilities_updateIPcoords(F)
   
   call discretization_setIPcoords(reshape(IPcoords,[3,grid(1)*grid(2)*grid3]))
 
-end subroutine utilities_updateIPcoords
+end subroutine utilities_updateCoords
 
 end module spectral_utilities