diff --git a/sources/interpolations.F90 b/sources/interpolations.F90
index 0db9c30..08e637f 100644
--- a/sources/interpolations.F90
+++ b/sources/interpolations.F90
@@ -651,8 +651,14 @@ module interpolations
 ! local variables
 !
     logical             :: flag
-    integer             :: i, j, i1, j1, k1, i2, j2, k2
-    real(kind=max_prec) :: sig, fc, fx, fy, fz, xl, xr, yl, yr, zl, zr
+    integer             :: i, j, i1, j1, i2, j2
+#if NDIMS == 3
+    integer             :: k1, k2
+#endif /* NDIMS == 3 */
+    real(kind=max_prec) :: sig, fc, fx, fy, xl, xr, yl, yr
+#if NDIMS == 3
+    real(kind=max_prec) :: fz, zl, zr
+#endif /* NDIMS == 3 */
 
 ! local arrays for derivatives
 !
@@ -824,7 +830,10 @@ module interpolations
 ! local variables
 !
     integer :: i  , j  , k = 1
-    integer :: im1, jm1, km1, ip1, jp1, kp1
+    integer :: im1, jm1, ip1, jp1
+#if NDIMS == 3
+    integer :: km1, kp1
+#endif /* NDIMS == 3 */
 
 ! local vectors
 !
@@ -1067,10 +1076,13 @@ module interpolations
 
 ! local variables
 !
-    logical       :: flag
-    integer       :: i  , j  , k = 1
-    integer       :: il , jl , kl , iu , ju , ku
-    integer       :: im1, jm1, km1, ip1, jp1, kp1
+    logical :: flag
+    integer :: i, il, iu, im1, ip1
+    integer :: j, jl, ju, jm1, jp1
+    integer :: k = 1
+#if NDIMS == 3
+    integer :: kl, ku, km1, kp1
+#endif /* NDIMS == 3 */
 
 ! local arrays for derivatives
 !
@@ -1259,7 +1271,10 @@ module interpolations
 !
     integer      :: i, im1, ip1
     integer      :: j, jm1, jp1
-    integer      :: k = 1, km1, kp1
+    integer      :: k = 1
+#if NDIMS == 3
+    integer      :: km1, kp1
+#endif /* NDIMS == 3 */
     integer      :: m
 #if NDIMS == 3
     integer      :: n, np1, np2