diff --git a/KiLCA/interface/wave_code_data_64bit.f90 b/KiLCA/interface/wave_code_data_64bit.f90
index 453665a1..beecae6a 100644
--- a/KiLCA/interface/wave_code_data_64bit.f90
+++ b/KiLCA/interface/wave_code_data_64bit.f90
@@ -802,6 +802,7 @@ subroutine interp_profile(dim_old, r_old, q_old, dim_new, r_new, q_new)
 
 integer :: deg = 15, Dmin = 0, Dmax = 0, l, ind;
 
+ind = 0;
 do l = 1,dim_new
 
     call eval_neville_polynom(dim_old, r_old, q_old, deg, r_new(l), Dmin, Dmax, ind, q_new(l));
diff --git a/QL-Balance/src/base/wave_code_data_64bit.f90 b/QL-Balance/src/base/wave_code_data_64bit.f90
index ca11f2f7..1e7562e8 100644
--- a/QL-Balance/src/base/wave_code_data_64bit.f90
+++ b/QL-Balance/src/base/wave_code_data_64bit.f90
@@ -738,6 +738,7 @@ subroutine interp_profile(dim_old, r_old, q_old, dim_new, r_new, q_new)
     real(8), dimension(dim_new) :: r_new, q_new;
     integer :: deg = 9, Dmin = 0, Dmax = 0, l, ind;
 
+    ind = 0;
     do l = 1, dim_new
         call eval_neville_polynom(dim_old, r_old, q_old, deg, r_new(l), Dmin, Dmax, ind, q_new(l));
     end do