diff --git a/sources/interpolations.F90 b/sources/interpolations.F90 index 80d0225..f8ed3ea 100644 --- a/sources/interpolations.F90 +++ b/sources/interpolations.F90 @@ -852,7 +852,7 @@ module interpolations real(kind=8), dimension(:,:,:) , intent(in) :: q real(kind=8), dimension(:,:,:,:,:), intent(out) :: qi - integer :: i , j , k = 1 + integer :: i , j , k integer :: im1, jm1, ip1, jp1 #if NDIMS == 3 integer :: km1, kp1 @@ -866,6 +866,10 @@ module interpolations !------------------------------------------------------------------------------- ! +#if NDIMS == 2 + k = 1 +#endif /* NDIMS == 2 */ + ! copy ghost zones ! do j = 1, NDIMS @@ -982,10 +986,14 @@ module interpolations real(kind=8), dimension(:,:,:) , intent(in) :: q real(kind=8), dimension(:,:,:,:,:), intent(out) :: qi - integer :: i, j, k = 1 + integer :: i, j, k !------------------------------------------------------------------------------- ! +#if NDIMS == 2 + k = 1 +#endif /* NDIMS == 2 */ + ! copy ghost zones ! do j = 1, NDIMS @@ -1085,7 +1093,7 @@ module interpolations logical :: flag integer :: i, il, iu, im1, ip1 integer :: j, jl, ju, jm1, jp1 - integer :: k = 1 + integer :: k #if NDIMS == 3 integer :: kl, ku, km1, kp1 #endif /* NDIMS == 3 */ @@ -1099,6 +1107,10 @@ module interpolations !------------------------------------------------------------------------------- ! +#if NDIMS == 2 + k = 1 +#endif /* NDIMS == 2 */ + ! copy ghost zones ! do j = 1, NDIMS @@ -1268,7 +1280,7 @@ module interpolations integer :: i, im1, ip1 integer :: j, jm1, jp1 - integer :: k = 1 + integer :: k #if NDIMS == 3 integer :: km1, kp1 #endif /* NDIMS == 3 */ @@ -1287,6 +1299,10 @@ module interpolations !------------------------------------------------------------------------------- ! +#if NDIMS == 2 + k = 1 +#endif /* NDIMS == 2 */ + #if NDIMS == 3 do k = 1, nn km1 = max( 1, k - 1)