Regenerate code (Kranc commit 519f2cd182c4828e74429e883139fbb8b1a11d3c)

Have tested convergence to exact solution for 3D shifted gauge wave with conformalMethod = 0.

4 files changed, 126 insertions, 12 deletions

diff --git a/ML_BSSN/src/ML_BSSN_RHS1.cc b/ML_BSSN/src/ML_BSSN_RHS1.cc
index 80feba6..ae5fa93 100644
--- a/ML_BSSN/src/ML_BSSN_RHS1.cc
+++ b/ML_BSSN/src/ML_BSSN_RHS1.cc
@@ -1358,6 +1358,15 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED theta = 
       kfmin(ToReal(1.),kexp(kmadd(kmul(rL,ToReal(-1.)),ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1.))));
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt1 = 
+      kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu22,JacPDstandardNth1gt22,kmadd(gtu33,JacPDstandardNth1gt33,kmul(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmul(gtu23,JacPDstandardNth1gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt2 = 
+      kmadd(gtu11,JacPDstandardNth2gt11,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmul(kmadd(gtu12,JacPDstandardNth2gt12,kmadd(gtu13,JacPDstandardNth2gt13,kmul(gtu23,JacPDstandardNth2gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt3 = 
+      kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2.)))));
+    
     CCTK_REAL_VEC beta1rhsL CCTK_ATTRIBUTE_UNUSED ;
     CCTK_REAL_VEC beta2rhsL CCTK_ATTRIBUTE_UNUSED ;
     CCTK_REAL_VEC beta3rhsL CCTK_ATTRIBUTE_UNUSED ;
@@ -1365,13 +1374,31 @@ static void ML_BSSN_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
     if (harmonicShift)
     {
       beta1rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu11,JacPDstandardNth1alpha,kmadd(gtu12,JacPDstandardNth2alpha,kmul(gtu13,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth1gt11,kmul(gtu11,gtu11),kmul(JacPDstandardNth1gt22,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu13,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmul(kmadd(gtu13,JacPDstandardNth1gt33,kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu11,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu12,JacPDstandardNth2gt11,kmadd(gtu13,JacPDstandardNth3gt11,kmadd(kmadd(gtu22,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth1gt33)),ToReal(-1.),kmadd(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu23,JacPDstandardNth3gt12,kmul(gtu33,JacPDstandardNth3gt13))))),ToReal(2.),kmul(gtu23,kmadd(JacPDstandardNth1gt23,ToReal(-2.),kmul(JacPDstandardNth2gt13,ToReal(2.)))))))))),kmul(gtu12,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu23,kmul(JacPDstandardNth3gt22,ToReal(2.)),kmadd(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))),kmul(gtu13,kmul(JacPDstandardNth1gt23,ToReal(4.)))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
       
       beta2rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu12,JacPDstandardNth1alpha,kmadd(gtu22,JacPDstandardNth2alpha,kmul(gtu23,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth2gt22,kmul(gtu22,gtu22),kmul(JacPDstandardNth2gt11,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu23,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu11,kmul(JacPDstandardNth3gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu13,JacPDstandardNth1gt33,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu22,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu23,JacPDstandardNth3gt22,kmadd(kmadd(gtu23,JacPDstandardNth2gt23,kmul(gtu13,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth2gt13,ToReal(-1.),JacPDstandardNth3gt12)))),ToReal(2.),kmadd(gtu11,kmadd(JacPDstandardNth2gt11,ToReal(-1.),kmul(JacPDstandardNth1gt12,ToReal(2.))),kmul(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))))))))),kmul(gtu12,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu13,kmul(JacPDstandardNth3gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth2gt12,ToReal(2.),JacPDstandardNth1gt22),kmadd(gtu33,kmadd(JacPDstandardNth1gt33,ToReal(-1.),kmul(JacPDstandardNth3gt13,ToReal(2.))),kmul(gtu23,kmul(JacPDstandardNth2gt13,ToReal(4.))))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu22,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
       
       beta3rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu13,JacPDstandardNth1alpha,kmadd(gtu23,JacPDstandardNth2alpha,kmul(gtu33,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmul(kmadd(JacPDstandardNth3gt11,kmul(gtu13,gtu13),kmul(JacPDstandardNth3gt22,kmul(gtu23,gtu23))),ToReal(2.)),kmadd(gtu23,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmadd(gtu11,kmul(JacPDstandardNth2gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth3gt23))),ToReal(2.))))))),kmadd(gtu33,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmadd(kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu12,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth3gt12,ToReal(-1.),JacPDstandardNth2gt13)))),ToReal(2.),kmul(gtu11,kmadd(JacPDstandardNth3gt11,ToReal(-1.),kmul(JacPDstandardNth1gt13,ToReal(2.))))))))),kmul(gtu13,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kmul(JacPDstandardNth2gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth1gt22,ToReal(-1.),kmul(JacPDstandardNth2gt12,ToReal(2.))),kmadd(gtu33,kmadd(JacPDstandardNth3gt13,ToReal(2.),JacPDstandardNth1gt33),kmul(gtu23,kmul(JacPDstandardNth3gt12,ToReal(4.))))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu33,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
     }
     else
     {
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
index fa11e9b..3b8d02f 100644
--- a/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
+++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseDtLapseShift.cc
@@ -835,25 +835,55 @@ static void ML_BSSN_convertToADMBaseDtLapseShift_Body(cGH const * restrict const
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu33 = 
       kdiv(kmadd(gt11L,gt22L,kmul(kmul(gt12L,gt12L),ToReal(-1.))),detgt);
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED em4phi = IfThen(conformalMethod == 
+      1.,kmul(phiL,phiL),kexp(kmul(phiL,ToReal(-4.))));
+    
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED eta = 
       kdiv(ToReal(SpatialBetaDriverRadius),kfmax(rL,ToReal(SpatialBetaDriverRadius)));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED theta = 
       kfmin(ToReal(1.),kexp(kmadd(kmul(rL,ToReal(-1.)),ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1.))));
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt1 = 
+      kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu22,JacPDstandardNth1gt22,kmadd(gtu33,JacPDstandardNth1gt33,kmul(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmul(gtu23,JacPDstandardNth1gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt2 = 
+      kmadd(gtu11,JacPDstandardNth2gt11,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmul(kmadd(gtu12,JacPDstandardNth2gt12,kmadd(gtu13,JacPDstandardNth2gt13,kmul(gtu23,JacPDstandardNth2gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt3 = 
+      kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2.)))));
+    
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtalpL = 
       kmadd(kpow(alphaL,harmonicN),kmul(ToReal(-1.),kmul(ToReal(harmonicF),kmadd(kmadd(trKL,ToReal(-1.),AL),ToReal(LapseACoeff),trKL))),kmul(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff)));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetaxL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu11,JacPDstandardNth1alpha,kmadd(gtu12,JacPDstandardNth2alpha,kmul(gtu13,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth1gt11,kmul(gtu11,gtu11),kmul(JacPDstandardNth1gt22,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu13,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmul(kmadd(gtu13,JacPDstandardNth1gt33,kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu11,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu12,JacPDstandardNth2gt11,kmadd(gtu13,JacPDstandardNth3gt11,kmadd(kmadd(gtu22,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth1gt33)),ToReal(-1.),kmadd(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu23,JacPDstandardNth3gt12,kmul(gtu33,JacPDstandardNth3gt13))))),ToReal(2.),kmul(gtu23,kmadd(JacPDstandardNth1gt23,ToReal(-2.),kmul(JacPDstandardNth2gt13,ToReal(2.)))))))))),kmul(gtu12,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu23,kmul(JacPDstandardNth3gt22,ToReal(2.)),kmadd(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))),kmul(gtu13,kmul(JacPDstandardNth1gt23,ToReal(4.))))))))))))))))),kmul(theta,kmul(kadd(Xt1L,kmadd(kmadd(Xt1L,ToReal(-1.),B1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt1L,kmadd(kmadd(Xt1L,ToReal(-1.),B1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetayL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta2,kmadd(beta2L,JacPDupwindNthAnti2beta2,kmadd(beta3L,JacPDupwindNthAnti3beta2,kmadd(JacPDupwindNthSymm1beta2,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta2,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta2,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu12,JacPDstandardNth1alpha,kmadd(gtu22,JacPDstandardNth2alpha,kmul(gtu23,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth2gt22,kmul(gtu22,gtu22),kmul(JacPDstandardNth2gt11,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu23,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu11,kmul(JacPDstandardNth3gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu13,JacPDstandardNth1gt33,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu22,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu23,JacPDstandardNth3gt22,kmadd(kmadd(gtu23,JacPDstandardNth2gt23,kmul(gtu13,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth2gt13,ToReal(-1.),JacPDstandardNth3gt12)))),ToReal(2.),kmadd(gtu11,kmadd(JacPDstandardNth2gt11,ToReal(-1.),kmul(JacPDstandardNth1gt12,ToReal(2.))),kmul(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))))))))),kmul(gtu12,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu13,kmul(JacPDstandardNth3gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth2gt12,ToReal(2.),JacPDstandardNth1gt22),kmadd(gtu33,kmadd(JacPDstandardNth1gt33,ToReal(-1.),kmul(JacPDstandardNth3gt13,ToReal(2.))),kmul(gtu23,kmul(JacPDstandardNth2gt13,ToReal(4.)))))))))))))))))),kmul(theta,kmul(kadd(Xt2L,kmadd(kmadd(Xt2L,ToReal(-1.),B2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta2,kmadd(beta2L,JacPDupwindNthAnti2beta2,kmadd(beta3L,JacPDupwindNthAnti3beta2,kmadd(JacPDupwindNthSymm1beta2,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta2,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta2,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu22,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt2L,kmadd(kmadd(Xt2L,ToReal(-1.),B2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetazL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta3,kmadd(beta2L,JacPDupwindNthAnti2beta3,kmadd(beta3L,JacPDupwindNthAnti3beta3,kmadd(JacPDupwindNthSymm1beta3,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta3,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta3,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu13,JacPDstandardNth1alpha,kmadd(gtu23,JacPDstandardNth2alpha,kmul(gtu33,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmul(kmadd(JacPDstandardNth3gt11,kmul(gtu13,gtu13),kmul(JacPDstandardNth3gt22,kmul(gtu23,gtu23))),ToReal(2.)),kmadd(gtu23,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmadd(gtu11,kmul(JacPDstandardNth2gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth3gt23))),ToReal(2.))))))),kmadd(gtu33,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmadd(kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu12,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth3gt12,ToReal(-1.),JacPDstandardNth2gt13)))),ToReal(2.),kmul(gtu11,kmadd(JacPDstandardNth3gt11,ToReal(-1.),kmul(JacPDstandardNth1gt13,ToReal(2.))))))))),kmul(gtu13,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kmul(JacPDstandardNth2gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth1gt22,ToReal(-1.),kmul(JacPDstandardNth2gt12,ToReal(2.))),kmadd(gtu33,kmadd(JacPDstandardNth3gt13,ToReal(2.),JacPDstandardNth1gt33),kmul(gtu23,kmul(JacPDstandardNth3gt12,ToReal(4.)))))))))))))))))),kmul(theta,kmul(kadd(Xt3L,kmadd(kmadd(Xt3L,ToReal(-1.),B3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta3,kmadd(beta2L,JacPDupwindNthAnti2beta3,kmadd(beta3L,JacPDupwindNthAnti3beta3,kmadd(JacPDupwindNthSymm1beta3,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta3,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta3,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu33,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt3L,kmadd(kmadd(Xt3L,ToReal(-1.),B3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     /* Copy local copies back to grid functions */
diff --git a/ML_CCZ4/src/ML_CCZ4_RHS1.cc b/ML_CCZ4/src/ML_CCZ4_RHS1.cc
index 0aeca46..bccf8f7 100644
--- a/ML_CCZ4/src/ML_CCZ4_RHS1.cc
+++ b/ML_CCZ4/src/ML_CCZ4_RHS1.cc
@@ -2272,6 +2272,15 @@ static void ML_CCZ4_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED theta = 
       kfmin(ToReal(1.),kexp(kmadd(kmul(rL,ToReal(-1.)),ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1.))));
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt1 = 
+      kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu22,JacPDstandardNth1gt22,kmadd(gtu33,JacPDstandardNth1gt33,kmul(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmul(gtu23,JacPDstandardNth1gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt2 = 
+      kmadd(gtu11,JacPDstandardNth2gt11,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmul(kmadd(gtu12,JacPDstandardNth2gt12,kmadd(gtu13,JacPDstandardNth2gt13,kmul(gtu23,JacPDstandardNth2gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt3 = 
+      kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2.)))));
+    
     CCTK_REAL_VEC beta1rhsL CCTK_ATTRIBUTE_UNUSED ;
     CCTK_REAL_VEC beta2rhsL CCTK_ATTRIBUTE_UNUSED ;
     CCTK_REAL_VEC beta3rhsL CCTK_ATTRIBUTE_UNUSED ;
@@ -2279,13 +2288,31 @@ static void ML_CCZ4_RHS1_Body(cGH const * restrict const cctkGH, int const dir,
     if (harmonicShift)
     {
       beta1rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu11,JacPDstandardNth1alpha,kmadd(gtu12,JacPDstandardNth2alpha,kmul(gtu13,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth1gt11,kmul(gtu11,gtu11),kmul(JacPDstandardNth1gt22,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu13,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmul(kmadd(gtu13,JacPDstandardNth1gt33,kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu11,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu12,JacPDstandardNth2gt11,kmadd(gtu13,JacPDstandardNth3gt11,kmadd(kmadd(gtu22,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth1gt33)),ToReal(-1.),kmadd(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu23,JacPDstandardNth3gt12,kmul(gtu33,JacPDstandardNth3gt13))))),ToReal(2.),kmul(gtu23,kmadd(JacPDstandardNth1gt23,ToReal(-2.),kmul(JacPDstandardNth2gt13,ToReal(2.)))))))))),kmul(gtu12,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu23,kmul(JacPDstandardNth3gt22,ToReal(2.)),kmadd(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))),kmul(gtu13,kmul(JacPDstandardNth1gt23,ToReal(4.)))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
       
       beta2rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu12,JacPDstandardNth1alpha,kmadd(gtu22,JacPDstandardNth2alpha,kmul(gtu23,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth2gt22,kmul(gtu22,gtu22),kmul(JacPDstandardNth2gt11,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu23,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu11,kmul(JacPDstandardNth3gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu13,JacPDstandardNth1gt33,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu22,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu23,JacPDstandardNth3gt22,kmadd(kmadd(gtu23,JacPDstandardNth2gt23,kmul(gtu13,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth2gt13,ToReal(-1.),JacPDstandardNth3gt12)))),ToReal(2.),kmadd(gtu11,kmadd(JacPDstandardNth2gt11,ToReal(-1.),kmul(JacPDstandardNth1gt12,ToReal(2.))),kmul(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))))))))),kmul(gtu12,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu13,kmul(JacPDstandardNth3gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth2gt12,ToReal(2.),JacPDstandardNth1gt22),kmadd(gtu33,kmadd(JacPDstandardNth1gt33,ToReal(-1.),kmul(JacPDstandardNth3gt13,ToReal(2.))),kmul(gtu23,kmul(JacPDstandardNth2gt13,ToReal(4.))))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu22,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
       
       beta3rhsL = 
-        kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu13,JacPDstandardNth1alpha,kmadd(gtu23,JacPDstandardNth2alpha,kmul(gtu33,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmul(kmadd(JacPDstandardNth3gt11,kmul(gtu13,gtu13),kmul(JacPDstandardNth3gt22,kmul(gtu23,gtu23))),ToReal(2.)),kmadd(gtu23,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmadd(gtu11,kmul(JacPDstandardNth2gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth3gt23))),ToReal(2.))))))),kmadd(gtu33,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmadd(kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu12,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth3gt12,ToReal(-1.),JacPDstandardNth2gt13)))),ToReal(2.),kmul(gtu11,kmadd(JacPDstandardNth3gt11,ToReal(-1.),kmul(JacPDstandardNth1gt13,ToReal(2.))))))))),kmul(gtu13,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kmul(JacPDstandardNth2gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth1gt22,ToReal(-1.),kmul(JacPDstandardNth2gt12,ToReal(2.))),kmadd(gtu33,kmadd(JacPDstandardNth3gt13,ToReal(2.),JacPDstandardNth1gt33),kmul(gtu23,kmul(JacPDstandardNth3gt12,ToReal(4.))))))))))))))))));
+        kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu33,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+        == 
+        1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.)))))))));
     }
     else
     {
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc
index db5b4d8..aa24a0c 100644
--- a/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc
+++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBaseDtLapseShift.cc
@@ -836,27 +836,57 @@ static void ML_CCZ4_convertToADMBaseDtLapseShift_Body(cGH const * restrict const
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED gtu33 = 
       kdiv(kmadd(gt11L,gt22L,kmul(kmul(gt12L,gt12L),ToReal(-1.))),detgt);
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED em4phi = IfThen(conformalMethod == 
+      1.,kmul(phiL,phiL),kexp(kmul(phiL,ToReal(-4.))));
+    
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED eta = 
       kdiv(ToReal(SpatialBetaDriverRadius),kfmax(rL,ToReal(SpatialBetaDriverRadius)));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED theta = 
       kfmin(ToReal(1.),kexp(kmadd(kmul(rL,ToReal(-1.)),ToReal(ScalarINV(SpatialShiftGammaCoeffRadius)),ToReal(1.))));
     
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt1 = 
+      kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu22,JacPDstandardNth1gt22,kmadd(gtu33,JacPDstandardNth1gt33,kmul(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmul(gtu23,JacPDstandardNth1gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt2 = 
+      kmadd(gtu11,JacPDstandardNth2gt11,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmul(kmadd(gtu12,JacPDstandardNth2gt12,kmadd(gtu13,JacPDstandardNth2gt13,kmul(gtu23,JacPDstandardNth2gt23))),ToReal(2.)))));
+    
+    CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED Ddetgt3 = 
+      kmadd(gtu11,JacPDstandardNth3gt11,kmadd(gtu22,JacPDstandardNth3gt22,kmadd(gtu33,JacPDstandardNth3gt33,kmul(kmadd(gtu12,JacPDstandardNth3gt12,kmadd(gtu13,JacPDstandardNth3gt13,kmul(gtu23,JacPDstandardNth3gt23))),ToReal(2.)))));
+    
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtalpL = 
       kmadd(kpow(alphaL,harmonicN),kmul(ToReal(-1.),kmul(ToReal(harmonicF),kmadd(AL,ToReal(LapseACoeff),kmul(kmadd(trKL,ToReal(-1.),kmul(ThetaL,ToReal(2.))),ToReal(-1. 
       + 
       LapseACoeff))))),kmul(kmadd(beta1L,JacPDupwindNthAnti1alpha,kmadd(beta2L,JacPDupwindNthAnti2alpha,kmadd(beta3L,JacPDupwindNthAnti3alpha,kmadd(JacPDupwindNthSymm1alpha,kfabs(beta1L),kmadd(JacPDupwindNthSymm2alpha,kfabs(beta2L),kmul(JacPDupwindNthSymm3alpha,kfabs(beta3L))))))),ToReal(LapseAdvectionCoeff)));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetaxL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu11,JacPDstandardNth1alpha,kmadd(gtu12,JacPDstandardNth2alpha,kmul(gtu13,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth1gt11,kmul(gtu11,gtu11),kmul(JacPDstandardNth1gt22,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu13,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmul(kmadd(gtu13,JacPDstandardNth1gt33,kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu11,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu12,JacPDstandardNth2gt11,kmadd(gtu13,JacPDstandardNth3gt11,kmadd(kmadd(gtu22,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth1gt33)),ToReal(-1.),kmadd(kmadd(gtu12,JacPDstandardNth1gt12,kmadd(gtu13,JacPDstandardNth1gt13,kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu23,JacPDstandardNth3gt12,kmul(gtu33,JacPDstandardNth3gt13))))),ToReal(2.),kmul(gtu23,kmadd(JacPDstandardNth1gt23,ToReal(-2.),kmul(JacPDstandardNth2gt13,ToReal(2.)))))))))),kmul(gtu12,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu23,kmul(JacPDstandardNth3gt22,ToReal(2.)),kmadd(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))),kmul(gtu13,kmul(JacPDstandardNth1gt23,ToReal(4.))))))))))))))))),kmul(theta,kmul(kadd(Xt1L,kmadd(kmadd(Xt1L,ToReal(-1.),B1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta1,kmadd(beta2L,JacPDupwindNthAnti2beta1,kmadd(beta3L,JacPDupwindNthAnti3beta1,kmadd(JacPDupwindNthSymm1beta1,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta1,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta1,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu11,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt1L,kmadd(kmadd(Xt1L,ToReal(-1.),B1L),ToReal(ShiftBCoeff),kmul(beta1L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetayL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta2,kmadd(beta2L,JacPDupwindNthAnti2beta2,kmadd(beta3L,JacPDupwindNthAnti3beta2,kmadd(JacPDupwindNthSymm1beta2,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta2,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta2,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu12,JacPDstandardNth1alpha,kmadd(gtu22,JacPDstandardNth2alpha,kmul(gtu23,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmadd(JacPDstandardNth2gt22,kmul(gtu22,gtu22),kmul(JacPDstandardNth2gt11,kmul(kmul(gtu12,gtu12),ToReal(2.)))),kmadd(gtu23,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu11,kmul(JacPDstandardNth3gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu13,JacPDstandardNth1gt33,kmul(gtu23,JacPDstandardNth2gt33))),ToReal(2.)))))),kmadd(gtu22,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu23,JacPDstandardNth3gt22,kmadd(kmadd(gtu23,JacPDstandardNth2gt23,kmul(gtu13,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth2gt13,ToReal(-1.),JacPDstandardNth3gt12)))),ToReal(2.),kmadd(gtu11,kmadd(JacPDstandardNth2gt11,ToReal(-1.),kmul(JacPDstandardNth1gt12,ToReal(2.))),kmul(gtu33,kmadd(JacPDstandardNth2gt33,ToReal(-1.),kmul(JacPDstandardNth3gt23,ToReal(2.))))))))),kmul(gtu12,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu13,kmul(JacPDstandardNth3gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth2gt12,ToReal(2.),JacPDstandardNth1gt22),kmadd(gtu33,kmadd(JacPDstandardNth1gt33,ToReal(-1.),kmul(JacPDstandardNth3gt13,ToReal(2.))),kmul(gtu23,kmul(JacPDstandardNth2gt13,ToReal(4.)))))))))))))))))),kmul(theta,kmul(kadd(Xt2L,kmadd(kmadd(Xt2L,ToReal(-1.),B2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta2,kmadd(beta2L,JacPDupwindNthAnti2beta2,kmadd(beta3L,JacPDupwindNthAnti3beta2,kmadd(JacPDupwindNthSymm1beta2,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta2,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta2,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu12,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu22,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt2L,kmadd(kmadd(Xt2L,ToReal(-1.),B2L),ToReal(ShiftBCoeff),kmul(beta2L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     CCTK_REAL_VEC CCTK_ATTRIBUTE_UNUSED dtbetazL = 
-      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta3,kmadd(beta2L,JacPDupwindNthAnti2beta3,kmadd(beta3L,JacPDupwindNthAnti3beta3,kmadd(JacPDupwindNthSymm1beta3,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta3,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta3,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(phiL,kmul(ToReal(0.5),kmadd(phiL,kmul(kmadd(gtu13,JacPDstandardNth1alpha,kmadd(gtu23,JacPDstandardNth2alpha,kmul(gtu33,JacPDstandardNth3alpha))),ToReal(-2.)),kmul(alphaL,kmadd(phiL,kmul(kmadd(JacPDstandardNth3gt11,kmul(gtu13,gtu13),kmul(JacPDstandardNth3gt22,kmul(gtu23,gtu23))),ToReal(2.)),kmadd(gtu23,kmadd(JacPDstandardNth2phi,ToReal(2.),kmul(phiL,kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu33,JacPDstandardNth2gt33,kmadd(gtu11,kmul(JacPDstandardNth2gt11,ToReal(-1.)),kmul(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,JacPDstandardNth1gt22,kmul(gtu33,JacPDstandardNth3gt23))),ToReal(2.))))))),kmadd(gtu33,kmadd(JacPDstandardNth3phi,ToReal(2.),kmul(phiL,kmadd(gtu33,JacPDstandardNth3gt33,kmadd(gtu22,kmul(JacPDstandardNth3gt22,ToReal(-1.)),kmadd(kmadd(gtu22,JacPDstandardNth2gt23,kmul(gtu12,kadd(JacPDstandardNth1gt23,kmadd(JacPDstandardNth3gt12,ToReal(-1.),JacPDstandardNth2gt13)))),ToReal(2.),kmul(gtu11,kmadd(JacPDstandardNth3gt11,ToReal(-1.),kmul(JacPDstandardNth1gt13,ToReal(2.))))))))),kmul(gtu13,kmadd(JacPDstandardNth1phi,ToReal(2.),kmul(phiL,kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kmul(JacPDstandardNth2gt11,ToReal(2.)),kmadd(gtu22,kmadd(JacPDstandardNth1gt22,ToReal(-1.),kmul(JacPDstandardNth2gt12,ToReal(2.))),kmadd(gtu33,kmadd(JacPDstandardNth3gt13,ToReal(2.),JacPDstandardNth1gt33),kmul(gtu23,kmul(JacPDstandardNth3gt12,ToReal(4.)))))))))))))))))),kmul(theta,kmul(kadd(Xt3L,kmadd(kmadd(Xt3L,ToReal(-1.),B3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,ToReal(BetaDriver*(-1. 
+      kmadd(kmadd(beta1L,JacPDupwindNthAnti1beta3,kmadd(beta2L,JacPDupwindNthAnti2beta3,kmadd(beta3L,JacPDupwindNthAnti3beta3,kmadd(JacPDupwindNthSymm1beta3,kfabs(beta1L),kmadd(JacPDupwindNthSymm2beta3,kfabs(beta2L),kmul(JacPDupwindNthSymm3beta3,kfabs(beta3L))))))),ToReal(ShiftAdvectionCoeff),IfThen(harmonicShift,kmul(alphaL,kmul(em4phi,kmul(ToReal(-0.5),kmadd(gtu13,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt11,kmadd(gtu12,kadd(JacPDstandardNth1gt12,JacPDstandardNth2gt11),kmadd(gtu22,JacPDstandardNth2gt12,kmadd(gtu13,kadd(JacPDstandardNth1gt13,JacPDstandardNth3gt11),kmadd(gtu23,kadd(JacPDstandardNth2gt13,JacPDstandardNth3gt12),kmadd(gtu33,JacPDstandardNth3gt13,kmul(JacPDstandardNth1phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt1),kmul(JacPDstandardNth1alpha,ToReal(2.))),kmadd(gtu23,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt12,kmadd(gtu12,kadd(JacPDstandardNth1gt22,JacPDstandardNth2gt12),kmadd(gtu22,JacPDstandardNth2gt22,kmadd(gtu13,kadd(JacPDstandardNth1gt23,JacPDstandardNth3gt12),kmadd(gtu23,kadd(JacPDstandardNth2gt23,JacPDstandardNth3gt22),kmadd(gtu33,JacPDstandardNth3gt23,kmul(JacPDstandardNth2phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt2),kmul(JacPDstandardNth2alpha,ToReal(2.))),kmul(gtu33,kmadd(alphaL,kmadd(kmadd(gtu11,JacPDstandardNth1gt13,kmadd(gtu12,kadd(JacPDstandardNth1gt23,JacPDstandardNth2gt13),kmadd(gtu22,JacPDstandardNth2gt23,kmadd(gtu13,kadd(JacPDstandardNth1gt33,JacPDstandardNth3gt13),kmadd(gtu23,kadd(JacPDstandardNth2gt33,JacPDstandardNth3gt23),kmadd(gtu33,JacPDstandardNth3gt33,kmul(JacPDstandardNth3phi,IfThen(conformalMethod 
+      == 
+      1.,kdiv(ToReal(1),phiL),ToReal(-2.))))))))),ToReal(-2.),Ddetgt3),kmul(JacPDstandardNth3alpha,ToReal(2.))))))))),kmul(theta,kmul(kadd(Xt3L,kmadd(kmadd(Xt3L,ToReal(-1.),B3L),ToReal(ShiftBCoeff),kmul(beta3L,kmul(eta,ToReal(BetaDriver*(-1. 
       + ShiftBCoeff)))))),ToReal(ShiftGammaCoeff)))));
     
     /* Copy local copies back to grid functions */