diff --git a/src/LIB/tensorproductmultigrid_Source/discretisation.f90 b/src/LIB/tensorproductmultigrid_Source/discretisation.f90
index 06352e12b753d5000bd8ca36016856d39b8708c2..0952a99a3082c1c46b763288ab44a1dd510f7cf9 100644
--- a/src/LIB/tensorproductmultigrid_Source/discretisation.f90
+++ b/src/LIB/tensorproductmultigrid_Source/discretisation.f90
@@ -503,27 +503,27 @@ contains
     sigma_j = 2.0_rl*(1.0_rl*iy-0.5_rl)/grid_param%n-1.0_rl
     ! alpha_{i+1,j}
     if (ix == grid_param%n) then
-      alpha_T(1) = 2.0_rl*DSQRT((1.0_rl+(rho_i+0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
+      alpha_T(1) = 2.0_rl*SQRT((1.0_rl+(rho_i+0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
     else
-      alpha_T(1) = DSQRT((1.0_rl+(rho_i+0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
+      alpha_T(1) = SQRT((1.0_rl+(rho_i+0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
     end if
     ! alpha_{i-1,j}
     if (ix == 1) then
-      alpha_T(2) = 2.0_rl*DSQRT((1.0_rl+(rho_i-0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
+      alpha_T(2) = 2.0_rl*SQRT((1.0_rl+(rho_i-0.25_rl*h)**2)/(1.0_rl+sigma_j**2))
     else
-      alpha_T(2) = DSQRT((1.0_rl+(rho_i-0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
+      alpha_T(2) = SQRT((1.0_rl+(rho_i-0.5_rl*h)**2)/(1.0_rl+sigma_j**2))
     end if
     ! alpha_{i,j+1}
     if (iy == grid_param%n) then
-      alpha_T(3) = 2.0_rl*DSQRT((1.0_rl+(sigma_j+0.25_rl*h)**2)/(1.0_rl+rho_i**2))
+      alpha_T(3) = 2.0_rl*SQRT((1.0_rl+(sigma_j+0.25_rl*h)**2)/(1.0_rl+rho_i**2))
     else
-      alpha_T(3) = DSQRT((1.0_rl+(sigma_j+0.5_rl*h)**2)/(1.0_rl+rho_i**2))
+      alpha_T(3) = SQRT((1.0_rl+(sigma_j+0.5_rl*h)**2)/(1.0_rl+rho_i**2))
     end if
     ! alpha_{i,j-1}
     if (iy == 1) then
-      alpha_T(4) = 2.0_rl*DSQRT((1.0_rl+(sigma_j-0.25_rl*h)**2)/(1.0_rl+rho_i**2))
+      alpha_T(4) = 2.0_rl*SQRT((1.0_rl+(sigma_j-0.25_rl*h)**2)/(1.0_rl+rho_i**2))
     else
-      alpha_T(4) = DSQRT((1.0_rl+(sigma_j-0.5_rl*h)**2)/(1.0_rl+rho_i**2))
+      alpha_T(4) = SQRT((1.0_rl+(sigma_j-0.5_rl*h)**2)/(1.0_rl+rho_i**2))
     end if
 #endif
     alpha_T(5) = alpha_T(1) + alpha_T(2) + alpha_T(3) + alpha_T(4)