Fix inconsistent index positions

No change to generated code

3 files changed, 7 insertions, 7 deletions

diff --git a/m/McLachlan_ADMConstraints.m b/m/McLachlan_ADMConstraints.m
index 9d26a2b..19346e7 100644
--- a/m/McLachlan_ADMConstraints.m
+++ b/m/McLachlan_ADMConstraints.m
@@ -136,7 +136,7 @@ ADMConstraintsCalc =
   Where -> Interior,
   After -> "MoL_PostStep",
   Shorthands -> {detg, gu[ua,ub], G[ua,lb,lc],
-                 R[la,lb], trR, Km[la,lb], trK,
+                 R[la,lb], trR, Km[ua,lb], trK,
                  rho, S[la]},
   Equations -> 
   {
diff --git a/m/McLachlan_BSSN.m b/m/McLachlan_BSSN.m
index 8e7dc6d..8b97c74 100644
--- a/m/McLachlan_BSSN.m
+++ b/m/McLachlan_BSSN.m
@@ -799,7 +799,7 @@ RHSRadiativeBoundaryCalc =
     em4phi     -> IfThen [conformalMethod, phi^2, Exp[-4 phi]],
     gu[ua,ub]  -> em4phi gtu[ua,ub],
     
-    nn[la] -> normal[la],
+    nn[la] -> Euc[la,lb] normal[ub],
     nu[ua] -> gu[ua,ub] nn[lb],
     nlen2  -> nu[ua] nn[la],
     nlen   -> Sqrt[nlen2],
@@ -884,7 +884,7 @@ constraintsCalc =
                  Gt[ua,lb,lc], Gtl[la,lb,lc], Gtlu[la,lb,uc], Xtn[ua],
                  e4phi, em4phi,
                  g[la,lb], detg, gu[ua,ub], ddetg[la], G[ua,lb,lc],
-                 Rt[la,lb], Rphi[la,lb], R[la,lb], trR, Atm[la,lb],
+                 Rt[la,lb], Rphi[la,lb], R[la,lb], trR, Atm[ua,lb],
                  gK[la,lb,lc], cdphi[la], cdphi2[la,lb],
                  rho, S[la], fac1, fac2},
   Equations -> 
diff --git a/m/hydro.m b/m/hydro.m
index 85b9e45..a69194f 100644
--- a/m/hydro.m
+++ b/m/hydro.m
@@ -113,8 +113,8 @@ prim2conCalc =
   {
     vol     -> h^3,
     mass    -> vol rho,
-    mom[la] -> mass vel[la],
-    ene     -> (1/2) mass vel[ua] vel[la] + mass eps
+    mom[la] -> mass Euc[la,lb] vel[ub],
+    ene     -> (1/2) mass Euc[la,lb] vel[ub] vel[ua] + mass eps
   }
 };
 
@@ -129,8 +129,8 @@ con2primCalc =
   Equations -> 
   {
     rho     -> mass / vol,
-    vel[ua] -> mom[ua] / mass,
-    eps     -> ene / mass - (1/2) vel[ua] vel[la],
+    vel[ua] -> Euc[ua,ub] mom[lb] / mass,
+    eps     -> ene / mass - (1/2) Euc[la,lc] vel[ua] vel[uc],
     
     press   -> Gamma rho eps
   }