Regenerate code after Kranc change

18 files changed, 36 insertions, 36 deletions

diff --git a/ML_BSSN/src/ML_BSSN_Advect.cc b/ML_BSSN/src/ML_BSSN_Advect.cc
index 94ea8e6..d74a205 100644
--- a/ML_BSSN/src/ML_BSSN_Advect.cc
+++ b/ML_BSSN/src/ML_BSSN_Advect.cc
@@ -233,7 +233,7 @@ static void ML_BSSN_Advect_Body(cGH const * restrict const cctkGH, int const dir
   CCTK_LOOP3STR(ML_BSSN_Advect,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -2068,7 +2068,7 @@ static void ML_BSSN_Advect_Body(cGH const * restrict const cctkGH, int const dir
       + ShiftAdvectionCoeff)))),ToReal(ShiftBCoeff),B3rhsL);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpharhs[index],alpharhsL);
     vec_store_nta_partial(Arhs[index],ArhsL);
     vec_store_nta_partial(At11rhs[index],At11rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_Dissipation.cc b/ML_BSSN/src/ML_BSSN_Dissipation.cc
index a2645cb..d2a0715 100644
--- a/ML_BSSN/src/ML_BSSN_Dissipation.cc
+++ b/ML_BSSN/src/ML_BSSN_Dissipation.cc
@@ -233,7 +233,7 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
   CCTK_LOOP3STR(ML_BSSN_Dissipation,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -1239,7 +1239,7 @@ static void ML_BSSN_Dissipation_Body(cGH const * restrict const cctkGH, int cons
       kadd(B3rhsL,kmadd(epsdiss1,JacPDdissipationNth1B3,kmadd(epsdiss2,JacPDdissipationNth2B3,kmul(epsdiss3,JacPDdissipationNth3B3))));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpharhs[index],alpharhsL);
     vec_store_nta_partial(Arhs[index],ArhsL);
     vec_store_nta_partial(At11rhs[index],At11rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_InitGamma.cc b/ML_BSSN/src/ML_BSSN_InitGamma.cc
index 968ac6b..c37ebc8 100644
--- a/ML_BSSN/src/ML_BSSN_InitGamma.cc
+++ b/ML_BSSN/src/ML_BSSN_InitGamma.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_InitGamma_Body(cGH const * restrict const cctkGH, int const
   CCTK_LOOP3STR(ML_BSSN_InitGamma,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -241,7 +241,7 @@ static void ML_BSSN_InitGamma_Body(cGH const * restrict const cctkGH, int const
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = ToReal(0.);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(A[index],AL);
     vec_store_nta_partial(B1[index],B1L);
     vec_store_nta_partial(B2[index],B2L);
diff --git a/ML_BSSN/src/ML_BSSN_InitRHS.cc b/ML_BSSN/src/ML_BSSN_InitRHS.cc
index dde5f60..90a93bc 100644
--- a/ML_BSSN/src/ML_BSSN_InitRHS.cc
+++ b/ML_BSSN/src/ML_BSSN_InitRHS.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_InitRHS_Body(cGH const * restrict const cctkGH, int const di
   CCTK_LOOP3STR(ML_BSSN_InitRHS,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -277,7 +277,7 @@ static void ML_BSSN_InitRHS_Body(cGH const * restrict const cctkGH, int const di
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3rhsL = ToReal(0.);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpharhs[index],alpharhsL);
     vec_store_nta_partial(Arhs[index],ArhsL);
     vec_store_nta_partial(At11rhs[index],At11rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_Minkowski.cc b/ML_BSSN/src/ML_BSSN_Minkowski.cc
index d5ac436..812f9e1 100644
--- a/ML_BSSN/src/ML_BSSN_Minkowski.cc
+++ b/ML_BSSN/src/ML_BSSN_Minkowski.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const
   CCTK_LOOP3STR(ML_BSSN_Minkowski,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -278,7 +278,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = ToReal(0.);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(A[index],AL);
     vec_store_nta_partial(alpha[index],alphaL);
     vec_store_nta_partial(At11[index],At11L);
diff --git a/ML_BSSN/src/ML_BSSN_RHS1.cc b/ML_BSSN/src/ML_BSSN_RHS1.cc
index 31b0601..66a02ac 100644
--- a/ML_BSSN/src/ML_BSSN_RHS1.cc
+++ b/ML_BSSN/src/ML_BSSN_RHS1.cc
@@ -230,7 +230,7 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
   CCTK_LOOP3STR(ML_BSSN_RHS1,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -1398,7 +1398,7 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
       kmul(kmadd(B3L,kmul(eta,ToReal(-1.*BetaDriver)),dotXt3),ToReal(ShiftBCoeff));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpharhs[index],alpharhsL);
     vec_store_nta_partial(Arhs[index],ArhsL);
     vec_store_nta_partial(B1rhs[index],B1rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_RHS2.cc b/ML_BSSN/src/ML_BSSN_RHS2.cc
index 02b35fe..d27c0b7 100644
--- a/ML_BSSN/src/ML_BSSN_RHS2.cc
+++ b/ML_BSSN/src/ML_BSSN_RHS2.cc
@@ -209,7 +209,7 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir,
   CCTK_LOOP3STR(ML_BSSN_RHS2,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -1819,7 +1819,7 @@ static void ML_BSSN_RHS2_Body(cGH const * restrict const cctkGH, int const dir,
       kmadd(At33L,kmul(kadd(JacPDstandardNth1beta1,kadd(JacPDstandardNth2beta2,JacPDstandardNth3beta3)),ToReal(-0.666666666666666666666666666667)),kmadd(em4phi,kmadd(g33,kmul(trAts,ToReal(-0.333333333333333333333333333333)),Ats33),kmadd(alphaL,kmadd(At33L,trKL,kmadd(kmadd(At13L,Atm13,kmadd(At23L,Atm23,kmul(At33L,Atm33))),ToReal(-2.),kmul(em4phi,kmul(ToReal(-25.132741228718345),kmadd(g33,kmul(trS,ToReal(-0.333333333333333333333333333333)),eTzzL))))),kmul(kmadd(At13L,JacPDstandardNth3beta1,kmadd(At23L,JacPDstandardNth3beta2,kmul(At33L,JacPDstandardNth3beta3))),ToReal(2.)))));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(At11rhs[index],At11rhsL);
     vec_store_nta_partial(At12rhs[index],At12rhsL);
     vec_store_nta_partial(At13rhs[index],At13rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc
index c69d486..08fc3dd 100644
--- a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc
+++ b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc
@@ -233,7 +233,7 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in
   CCTK_LOOP3STR(ML_BSSN_RHSStaticBoundary,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -313,7 +313,7 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3rhsL = ToReal(0.);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpharhs[index],alpharhsL);
     vec_store_nta_partial(Arhs[index],ArhsL);
     vec_store_nta_partial(At11rhs[index],At11rhsL);
diff --git a/ML_BSSN/src/ML_BSSN_boundary.cc b/ML_BSSN/src/ML_BSSN_boundary.cc
index e838f14..e8e87b5 100644
--- a/ML_BSSN/src/ML_BSSN_boundary.cc
+++ b/ML_BSSN/src/ML_BSSN_boundary.cc
@@ -233,7 +233,7 @@ static void ML_BSSN_boundary_Body(cGH const * restrict const cctkGH, int const d
   CCTK_LOOP3STR(ML_BSSN_boundary,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -314,7 +314,7 @@ static void ML_BSSN_boundary_Body(cGH const * restrict const cctkGH, int const d
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED B3L = ToReal(0.);
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(A[index],AL);
     vec_store_nta_partial(alpha[index],alphaL);
     vec_store_nta_partial(At11[index],At11L);
diff --git a/ML_BSSN/src/ML_BSSN_constraints1.cc b/ML_BSSN/src/ML_BSSN_constraints1.cc
index 096e4d4..5c880fd 100644
--- a/ML_BSSN/src/ML_BSSN_constraints1.cc
+++ b/ML_BSSN/src/ML_BSSN_constraints1.cc
@@ -209,7 +209,7 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con
   CCTK_LOOP3STR(ML_BSSN_constraints1,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -1578,7 +1578,7 @@ static void ML_BSSN_constraints1_Body(cGH const * restrict const cctkGH, int con
       kadd(trR,kmadd(rho,ToReal(-50.26548245743669),kmadd(kmadd(Atm12,Atm21,kmadd(Atm13,Atm31,kmul(Atm23,Atm32))),ToReal(-2.),kmadd(kmadd(Atm11,Atm11,kmadd(Atm22,Atm22,kmul(Atm33,Atm33))),ToReal(-1.),kmul(kmul(trKL,trKL),ToReal(0.666666666666666666666666666667))))));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(H[index],HL);
   }
   CCTK_ENDLOOP3STR(ML_BSSN_constraints1);
diff --git a/ML_BSSN/src/ML_BSSN_constraints2.cc b/ML_BSSN/src/ML_BSSN_constraints2.cc
index a132d2c..362f533 100644
--- a/ML_BSSN/src/ML_BSSN_constraints2.cc
+++ b/ML_BSSN/src/ML_BSSN_constraints2.cc
@@ -218,7 +218,7 @@ static void ML_BSSN_constraints2_Body(cGH const * restrict const cctkGH, int con
   CCTK_LOOP3STR(ML_BSSN_constraints2,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -949,7 +949,7 @@ static void ML_BSSN_constraints2_Body(cGH const * restrict const cctkGH, int con
       kmadd(At11L,gtu11,kmadd(At22L,gtu22,kmadd(At33L,gtu33,kmul(kmadd(At12L,gtu12,kmadd(At13L,gtu13,kmul(At23L,gtu23))),ToReal(2.)))));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(cA[index],cAL);
     vec_store_nta_partial(cS[index],cSL);
     vec_store_nta_partial(cXt1[index],cXt1L);
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc
index 49f4382..8b5be6d 100644
--- a/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc
+++ b/ML_BSSN/src/ML_BSSN_convertFromADMBase.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_convertFromADMBase_Body(cGH const * restrict const cctkGH, i
   CCTK_LOOP3STR(ML_BSSN_convertFromADMBase,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -334,7 +334,7 @@ static void ML_BSSN_convertFromADMBase_Body(cGH const * restrict const cctkGH, i
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED beta3L = betazL;
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpha[index],alphaL);
     vec_store_nta_partial(At11[index],At11L);
     vec_store_nta_partial(At12[index],At12L);
diff --git a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc
index 253b6a0..e9db745 100644
--- a/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc
+++ b/ML_BSSN/src/ML_BSSN_convertFromADMBaseGamma.cc
@@ -215,7 +215,7 @@ static void ML_BSSN_convertFromADMBaseGamma_Body(cGH const * restrict const cctk
   CCTK_LOOP3STR(ML_BSSN_convertFromADMBaseGamma,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -861,7 +861,7 @@ static void ML_BSSN_convertFromADMBaseGamma_Body(cGH const * restrict const cctk
     }
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(A[index],AL);
     vec_store_nta_partial(B1[index],B1L);
     vec_store_nta_partial(B2[index],B2L);
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBase.cc b/ML_BSSN/src/ML_BSSN_convertToADMBase.cc
index 24c249e..36d0c7d 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBase.cc
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBase.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_convertToADMBase_Body(cGH const * restrict const cctkGH, int
   CCTK_LOOP3STR(ML_BSSN_convertToADMBase,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -292,7 +292,7 @@ static void ML_BSSN_convertToADMBase_Body(cGH const * restrict const cctkGH, int
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED betazL = beta3L;
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alp[index],alpL);
     vec_store_nta_partial(betax[index],betaxL);
     vec_store_nta_partial(betay[index],betayL);
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
index 3f243e2..d10b5d6 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
@@ -212,7 +212,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(cGH const * restrict const
   CCTK_LOOP3STR(ML_BSSN_convertToADMBaseDtLapseShift,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -857,7 +857,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(cGH const * restrict const
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(dtalp[index],dtalpL);
     vec_store_nta_partial(dtbetax[index],dtbetaxL);
     vec_store_nta_partial(dtbetay[index],dtbetayL);
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc
index 2f00e65..f50e8c2 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShiftBoundary.cc
@@ -212,7 +212,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
   CCTK_LOOP3STR(ML_BSSN_convertToADMBaseDtLapseShiftBoundary,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -291,7 +291,7 @@ static void ML_BSSN_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * restri
     }
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(dtalp[index],dtalpL);
     vec_store_nta_partial(dtbetax[index],dtbetaxL);
     vec_store_nta_partial(dtbetay[index],dtbetayL);
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc
index 6de84e5..dd877ee 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c
   CCTK_LOOP3STR(ML_BSSN_convertToADMBaseFakeDtLapseShift,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -276,7 +276,7 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c
     }
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(dtalp[index],dtalpL);
     vec_store_nta_partial(dtbetax[index],dtbetaxL);
     vec_store_nta_partial(dtbetay[index],dtbetayL);
diff --git a/ML_BSSN/src/ML_BSSN_enforce.cc b/ML_BSSN/src/ML_BSSN_enforce.cc
index 71b2032..62f0a87 100644
--- a/ML_BSSN/src/ML_BSSN_enforce.cc
+++ b/ML_BSSN/src/ML_BSSN_enforce.cc
@@ -197,7 +197,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
   CCTK_LOOP3STR(ML_BSSN_enforce,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2],
-    imin,imax, CCTK_REAL_VEC_SIZE)
+    kimin,kimax, CCTK_REAL_VEC_SIZE)
   {
     ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
@@ -283,7 +283,7 @@ static void ML_BSSN_enforce_Body(cGH const * restrict const cctkGH, int const di
     alphaL = kfmax(alphaL,ToReal(MinimumLapse));
     
     /* Copy local copies back to grid functions */
-    vec_store_partial_prepare(i,lc_imin,lc_imax);
+    vec_store_partial_prepare(i,kimin,kimax);
     vec_store_nta_partial(alpha[index],alphaL);
     vec_store_nta_partial(At11[index],At11L);
     vec_store_nta_partial(At12[index],At12L);