Regenerate examples

22 files changed, 1177 insertions, 1183 deletions

diff --git a/Examples/EulerAuto/interface.ccl b/Examples/EulerAuto/interface.ccl
index 0323809..6148db5 100644
--- a/Examples/EulerAuto/interface.ccl
+++ b/Examples/EulerAuto/interface.ccl
@@ -168,19 +168,19 @@ CCTK_REAL S_group type=GF timelevels=3 tags='tensortypealias="U" tensorweight=0'
 } "S_group"
 
 public:
-CCTK_REAL Den_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0'
+CCTK_REAL Den_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0 Prolongation="None"'
 {
   Denrhs
 } "Den_grouprhs"
 
 public:
-CCTK_REAL En_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0'
+CCTK_REAL En_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0 Prolongation="None"'
 {
   Enrhs
 } "En_grouprhs"
 
 public:
-CCTK_REAL S_grouprhs type=GF timelevels=3 tags='tensortypealias="U" tensorweight=0'
+CCTK_REAL S_grouprhs type=GF timelevels=3 tags='tensortypealias="U" tensorweight=0 Prolongation="None"'
 {
   S1rhs,
   S2rhs,
diff --git a/Examples/EulerAuto/schedule.ccl b/Examples/EulerAuto/schedule.ccl
index 1b1eeda..3dde361 100644
--- a/Examples/EulerAuto/schedule.ccl
+++ b/Examples/EulerAuto/schedule.ccl
@@ -187,78 +187,77 @@ if (CCTK_EQUALS(initial_data, "shock"))
   schedule eulerauto_initial_shock at CCTK_INITIAL as eulerauto_initial
   {
     LANG: C
-    READS:    grid::x
-    WRITES:   EulerAuto::p
-    WRITES:   EulerAuto::rho
-    WRITES:   EulerAuto::v1
-    WRITES:   EulerAuto::v2
-    WRITES:   EulerAuto::v3
+    READS:    grid::x(Everywhere)
+    WRITES:   EulerAuto::p(Everywhere)
+    WRITES:   EulerAuto::rho(Everywhere)
+    WRITES:   EulerAuto::v1(Everywhere)
+    WRITES:   EulerAuto::v2(Everywhere)
+    WRITES:   EulerAuto::v3(Everywhere)
   } "eulerauto_initial_shock"
 }
 
 schedule eulerauto_cons_calc_zero_rhs in MoL_CalcRHS
 {
   LANG: C
-  WRITES:   EulerAuto::Denrhs
-  WRITES:   EulerAuto::Enrhs
-  WRITES:   EulerAuto::S1rhs
-  WRITES:   EulerAuto::S2rhs
-  WRITES:   EulerAuto::S3rhs
+  WRITES:   EulerAuto::Denrhs(Everywhere)
+  WRITES:   EulerAuto::Enrhs(Everywhere)
+  WRITES:   EulerAuto::S1rhs(Everywhere)
+  WRITES:   EulerAuto::S2rhs(Everywhere)
+  WRITES:   EulerAuto::S3rhs(Everywhere)
 } "eulerauto_cons_calc_zero_rhs"
 
 schedule eulerauto_cons_calc_conserved at POSTINITIAL
 {
   LANG: C
-  READS:    EulerAuto::p
-  READS:    EulerAuto::rho
-  READS:    EulerAuto::v1
-  READS:    EulerAuto::v2
-  READS:    EulerAuto::v3
-  WRITES:   EulerAuto::Den
-  WRITES:   EulerAuto::En
-  WRITES:   EulerAuto::S1
-  WRITES:   EulerAuto::S2
-  WRITES:   EulerAuto::S3
+  READS:    EulerAuto::p(Everywhere)
+  READS:    EulerAuto::rho(Everywhere)
+  READS:    EulerAuto::v1(Everywhere)
+  READS:    EulerAuto::v2(Everywhere)
+  READS:    EulerAuto::v3(Everywhere)
+  WRITES:   EulerAuto::Den(Everywhere)
+  WRITES:   EulerAuto::En(Everywhere)
+  WRITES:   EulerAuto::S1(Everywhere)
+  WRITES:   EulerAuto::S2(Everywhere)
+  WRITES:   EulerAuto::S3(Everywhere)
 } "eulerauto_cons_calc_conserved"
 
 schedule eulerauto_cons_calc_primitives in MoL_PostStep after EulerAuto_ApplyBCs
 {
   LANG: C
-  READS:    EulerAuto::Den
-  READS:    EulerAuto::En
-  READS:    EulerAuto::S1
-  READS:    EulerAuto::S2
-  READS:    EulerAuto::S3
-  WRITES:   EulerAuto::p
-  WRITES:   EulerAuto::rho
-  WRITES:   EulerAuto::v1
-  WRITES:   EulerAuto::v2
-  WRITES:   EulerAuto::v3
+  READS:    EulerAuto::Den(Everywhere)
+  READS:    EulerAuto::En(Everywhere)
+  READS:    EulerAuto::S1(Everywhere)
+  READS:    EulerAuto::S2(Everywhere)
+  READS:    EulerAuto::S3(Everywhere)
+  WRITES:   EulerAuto::p(Everywhere)
+  WRITES:   EulerAuto::rho(Everywhere)
+  WRITES:   EulerAuto::v1(Everywhere)
+  WRITES:   EulerAuto::v2(Everywhere)
+  WRITES:   EulerAuto::v3(Everywhere)
 } "eulerauto_cons_calc_primitives"
 
 schedule group eulerauto_cons_calc_reconstruct_1_group in MoL_CalcRHS after eulerauto_cons_calc_zero_rhs
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_1"
 
 schedule eulerauto_cons_calc_reconstruct_1 in eulerauto_cons_calc_reconstruct_1_group
 {
   LANG: C
-  READS:    EulerAuto::p
-  READS:    EulerAuto::rho
-  READS:    EulerAuto::v1
-  READS:    EulerAuto::v2
-  READS:    EulerAuto::v3
-  WRITES:   EulerAuto::pLeft
-  WRITES:   EulerAuto::pRight
-  WRITES:   EulerAuto::rhoLeft
-  WRITES:   EulerAuto::rhoRight
-  WRITES:   EulerAuto::v1Left
-  WRITES:   EulerAuto::v1Right
-  WRITES:   EulerAuto::v2Left
-  WRITES:   EulerAuto::v2Right
-  WRITES:   EulerAuto::v3Left
-  WRITES:   EulerAuto::v3Right
+  READS:    EulerAuto::p(Everywhere)
+  READS:    EulerAuto::rho(Everywhere)
+  READS:    EulerAuto::v1(Everywhere)
+  READS:    EulerAuto::v2(Everywhere)
+  READS:    EulerAuto::v3(Everywhere)
+  WRITES:   EulerAuto::pLeft(Interior)
+  WRITES:   EulerAuto::pRight(Interior)
+  WRITES:   EulerAuto::rhoLeft(Interior)
+  WRITES:   EulerAuto::rhoRight(Interior)
+  WRITES:   EulerAuto::v1Left(Interior)
+  WRITES:   EulerAuto::v1Right(Interior)
+  WRITES:   EulerAuto::v2Left(Interior)
+  WRITES:   EulerAuto::v2Right(Interior)
+  WRITES:   EulerAuto::v3Left(Interior)
+  WRITES:   EulerAuto::v3Right(Interior)
 } "eulerauto_cons_calc_reconstruct_1"
 
 schedule eulerauto_cons_calc_reconstruct_1_SelectBCs in eulerauto_cons_calc_reconstruct_1_bc_group
@@ -274,77 +273,73 @@ schedule eulerauto_cons_calc_reconstruct_1_SelectBCs in eulerauto_cons_calc_reco
 
 schedule group ApplyBCs as eulerauto_cons_calc_reconstruct_1_ApplyBCs in eulerauto_cons_calc_reconstruct_1_bc_group after eulerauto_cons_calc_reconstruct_1_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_reconstruct_1"
 
 schedule group eulerauto_cons_calc_reconstruct_1_bc_group in eulerauto_cons_calc_reconstruct_1_group after eulerauto_cons_calc_reconstruct_1
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_1"
 
 schedule group eulerauto_cons_calc_reconstruct_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_1"
 
 schedule eulerauto_cons_calc_intercell_conserved_1 in MoL_CalcRHS after eulerauto_cons_calc_reconstruct_1
 {
   LANG: C
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenLeft
-  WRITES:   EulerAuto::DenRight
-  WRITES:   EulerAuto::EnLeft
-  WRITES:   EulerAuto::EnRight
-  WRITES:   EulerAuto::S1Left
-  WRITES:   EulerAuto::S1Right
-  WRITES:   EulerAuto::S2Left
-  WRITES:   EulerAuto::S2Right
-  WRITES:   EulerAuto::S3Left
-  WRITES:   EulerAuto::S3Right
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenLeft(Everywhere)
+  WRITES:   EulerAuto::DenRight(Everywhere)
+  WRITES:   EulerAuto::EnLeft(Everywhere)
+  WRITES:   EulerAuto::EnRight(Everywhere)
+  WRITES:   EulerAuto::S1Left(Everywhere)
+  WRITES:   EulerAuto::S1Right(Everywhere)
+  WRITES:   EulerAuto::S2Left(Everywhere)
+  WRITES:   EulerAuto::S2Right(Everywhere)
+  WRITES:   EulerAuto::S3Left(Everywhere)
+  WRITES:   EulerAuto::S3Right(Everywhere)
 } "eulerauto_cons_calc_intercell_conserved_1"
 
 schedule group eulerauto_cons_calc_flux_1_group in MoL_CalcRHS after eulerauto_cons_calc_intercell_conserved_1
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_1"
 
 schedule eulerauto_cons_calc_flux_1 in eulerauto_cons_calc_flux_1_group
 {
   LANG: C
-  READS:    EulerAuto::DenLeft
-  READS:    EulerAuto::DenRight
-  READS:    EulerAuto::EnLeft
-  READS:    EulerAuto::EnRight
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::S1Left
-  READS:    EulerAuto::S1Right
-  READS:    EulerAuto::S2Left
-  READS:    EulerAuto::S2Right
-  READS:    EulerAuto::S3Left
-  READS:    EulerAuto::S3Right
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenFlux
-  WRITES:   EulerAuto::EnFlux
-  WRITES:   EulerAuto::S1Flux
-  WRITES:   EulerAuto::S2Flux
-  WRITES:   EulerAuto::S3Flux
+  READS:    EulerAuto::DenLeft(Everywhere)
+  READS:    EulerAuto::DenRight(Everywhere)
+  READS:    EulerAuto::EnLeft(Everywhere)
+  READS:    EulerAuto::EnRight(Everywhere)
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::S1Left(Everywhere)
+  READS:    EulerAuto::S1Right(Everywhere)
+  READS:    EulerAuto::S2Left(Everywhere)
+  READS:    EulerAuto::S2Right(Everywhere)
+  READS:    EulerAuto::S3Left(Everywhere)
+  READS:    EulerAuto::S3Right(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenFlux(Interior)
+  WRITES:   EulerAuto::EnFlux(Interior)
+  WRITES:   EulerAuto::S1Flux(Interior)
+  WRITES:   EulerAuto::S2Flux(Interior)
+  WRITES:   EulerAuto::S3Flux(Interior)
 } "eulerauto_cons_calc_flux_1"
 
 schedule eulerauto_cons_calc_flux_1_SelectBCs in eulerauto_cons_calc_flux_1_bc_group
@@ -360,62 +355,58 @@ schedule eulerauto_cons_calc_flux_1_SelectBCs in eulerauto_cons_calc_flux_1_bc_g
 
 schedule group ApplyBCs as eulerauto_cons_calc_flux_1_ApplyBCs in eulerauto_cons_calc_flux_1_bc_group after eulerauto_cons_calc_flux_1_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_flux_1"
 
 schedule group eulerauto_cons_calc_flux_1_bc_group in eulerauto_cons_calc_flux_1_group after eulerauto_cons_calc_flux_1
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_1"
 
 schedule group eulerauto_cons_calc_flux_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_1"
 
 schedule eulerauto_cons_calc_rhs_1 in MoL_CalcRHS after eulerauto_cons_calc_flux_1
 {
   LANG: C
-  READS:    EulerAuto::DenFlux
-  READS:    EulerAuto::Denrhs
-  READS:    EulerAuto::EnFlux
-  READS:    EulerAuto::Enrhs
-  READS:    EulerAuto::S1Flux
-  READS:    EulerAuto::S1rhs
-  READS:    EulerAuto::S2Flux
-  READS:    EulerAuto::S2rhs
-  READS:    EulerAuto::S3Flux
-  READS:    EulerAuto::S3rhs
-  WRITES:   EulerAuto::Denrhs
-  WRITES:   EulerAuto::Enrhs
-  WRITES:   EulerAuto::S1rhs
-  WRITES:   EulerAuto::S2rhs
-  WRITES:   EulerAuto::S3rhs
+  READS:    EulerAuto::DenFlux(Everywhere)
+  READS:    EulerAuto::EnFlux(Everywhere)
+  READS:    EulerAuto::S1Flux(Everywhere)
+  READS:    EulerAuto::S2Flux(Everywhere)
+  READS:    EulerAuto::S3Flux(Everywhere)
+  READS:    EulerAuto::Denrhs(Everywhere)
+  READS:    EulerAuto::Enrhs(Everywhere)
+  READS:    EulerAuto::S1rhs(Everywhere)
+  READS:    EulerAuto::S2rhs(Everywhere)
+  READS:    EulerAuto::S3rhs(Everywhere)
+  WRITES:   EulerAuto::Denrhs(Interior)
+  WRITES:   EulerAuto::Enrhs(Interior)
+  WRITES:   EulerAuto::S1rhs(Interior)
+  WRITES:   EulerAuto::S2rhs(Interior)
+  WRITES:   EulerAuto::S3rhs(Interior)
 } "eulerauto_cons_calc_rhs_1"
 
 schedule group eulerauto_cons_calc_reconstruct_2_group in MoL_CalcRHS after eulerauto_cons_calc_rhs_1
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_2"
 
 schedule eulerauto_cons_calc_reconstruct_2 in eulerauto_cons_calc_reconstruct_2_group
 {
   LANG: C
-  READS:    EulerAuto::p
-  READS:    EulerAuto::rho
-  READS:    EulerAuto::v1
-  READS:    EulerAuto::v2
-  READS:    EulerAuto::v3
-  WRITES:   EulerAuto::pLeft
-  WRITES:   EulerAuto::pRight
-  WRITES:   EulerAuto::rhoLeft
-  WRITES:   EulerAuto::rhoRight
-  WRITES:   EulerAuto::v1Left
-  WRITES:   EulerAuto::v1Right
-  WRITES:   EulerAuto::v2Left
-  WRITES:   EulerAuto::v2Right
-  WRITES:   EulerAuto::v3Left
-  WRITES:   EulerAuto::v3Right
+  READS:    EulerAuto::p(Everywhere)
+  READS:    EulerAuto::rho(Everywhere)
+  READS:    EulerAuto::v1(Everywhere)
+  READS:    EulerAuto::v2(Everywhere)
+  READS:    EulerAuto::v3(Everywhere)
+  WRITES:   EulerAuto::pLeft(Interior)
+  WRITES:   EulerAuto::pRight(Interior)
+  WRITES:   EulerAuto::rhoLeft(Interior)
+  WRITES:   EulerAuto::rhoRight(Interior)
+  WRITES:   EulerAuto::v1Left(Interior)
+  WRITES:   EulerAuto::v1Right(Interior)
+  WRITES:   EulerAuto::v2Left(Interior)
+  WRITES:   EulerAuto::v2Right(Interior)
+  WRITES:   EulerAuto::v3Left(Interior)
+  WRITES:   EulerAuto::v3Right(Interior)
 } "eulerauto_cons_calc_reconstruct_2"
 
 schedule eulerauto_cons_calc_reconstruct_2_SelectBCs in eulerauto_cons_calc_reconstruct_2_bc_group
@@ -431,77 +422,73 @@ schedule eulerauto_cons_calc_reconstruct_2_SelectBCs in eulerauto_cons_calc_reco
 
 schedule group ApplyBCs as eulerauto_cons_calc_reconstruct_2_ApplyBCs in eulerauto_cons_calc_reconstruct_2_bc_group after eulerauto_cons_calc_reconstruct_2_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_reconstruct_2"
 
 schedule group eulerauto_cons_calc_reconstruct_2_bc_group in eulerauto_cons_calc_reconstruct_2_group after eulerauto_cons_calc_reconstruct_2
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_2"
 
 schedule group eulerauto_cons_calc_reconstruct_2_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_2"
 
 schedule eulerauto_cons_calc_intercell_conserved_2 in MoL_CalcRHS after eulerauto_cons_calc_reconstruct_2
 {
   LANG: C
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenLeft
-  WRITES:   EulerAuto::DenRight
-  WRITES:   EulerAuto::EnLeft
-  WRITES:   EulerAuto::EnRight
-  WRITES:   EulerAuto::S1Left
-  WRITES:   EulerAuto::S1Right
-  WRITES:   EulerAuto::S2Left
-  WRITES:   EulerAuto::S2Right
-  WRITES:   EulerAuto::S3Left
-  WRITES:   EulerAuto::S3Right
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenLeft(Everywhere)
+  WRITES:   EulerAuto::DenRight(Everywhere)
+  WRITES:   EulerAuto::EnLeft(Everywhere)
+  WRITES:   EulerAuto::EnRight(Everywhere)
+  WRITES:   EulerAuto::S1Left(Everywhere)
+  WRITES:   EulerAuto::S1Right(Everywhere)
+  WRITES:   EulerAuto::S2Left(Everywhere)
+  WRITES:   EulerAuto::S2Right(Everywhere)
+  WRITES:   EulerAuto::S3Left(Everywhere)
+  WRITES:   EulerAuto::S3Right(Everywhere)
 } "eulerauto_cons_calc_intercell_conserved_2"
 
 schedule group eulerauto_cons_calc_flux_2_group in MoL_CalcRHS after eulerauto_cons_calc_intercell_conserved_2
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_2"
 
 schedule eulerauto_cons_calc_flux_2 in eulerauto_cons_calc_flux_2_group
 {
   LANG: C
-  READS:    EulerAuto::DenLeft
-  READS:    EulerAuto::DenRight
-  READS:    EulerAuto::EnLeft
-  READS:    EulerAuto::EnRight
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::S1Left
-  READS:    EulerAuto::S1Right
-  READS:    EulerAuto::S2Left
-  READS:    EulerAuto::S2Right
-  READS:    EulerAuto::S3Left
-  READS:    EulerAuto::S3Right
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenFlux
-  WRITES:   EulerAuto::EnFlux
-  WRITES:   EulerAuto::S1Flux
-  WRITES:   EulerAuto::S2Flux
-  WRITES:   EulerAuto::S3Flux
+  READS:    EulerAuto::DenLeft(Everywhere)
+  READS:    EulerAuto::DenRight(Everywhere)
+  READS:    EulerAuto::EnLeft(Everywhere)
+  READS:    EulerAuto::EnRight(Everywhere)
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::S1Left(Everywhere)
+  READS:    EulerAuto::S1Right(Everywhere)
+  READS:    EulerAuto::S2Left(Everywhere)
+  READS:    EulerAuto::S2Right(Everywhere)
+  READS:    EulerAuto::S3Left(Everywhere)
+  READS:    EulerAuto::S3Right(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenFlux(Interior)
+  WRITES:   EulerAuto::EnFlux(Interior)
+  WRITES:   EulerAuto::S1Flux(Interior)
+  WRITES:   EulerAuto::S2Flux(Interior)
+  WRITES:   EulerAuto::S3Flux(Interior)
 } "eulerauto_cons_calc_flux_2"
 
 schedule eulerauto_cons_calc_flux_2_SelectBCs in eulerauto_cons_calc_flux_2_bc_group
@@ -517,62 +504,58 @@ schedule eulerauto_cons_calc_flux_2_SelectBCs in eulerauto_cons_calc_flux_2_bc_g
 
 schedule group ApplyBCs as eulerauto_cons_calc_flux_2_ApplyBCs in eulerauto_cons_calc_flux_2_bc_group after eulerauto_cons_calc_flux_2_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_flux_2"
 
 schedule group eulerauto_cons_calc_flux_2_bc_group in eulerauto_cons_calc_flux_2_group after eulerauto_cons_calc_flux_2
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_2"
 
 schedule group eulerauto_cons_calc_flux_2_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_2"
 
 schedule eulerauto_cons_calc_rhs_2 in MoL_CalcRHS after eulerauto_cons_calc_flux_2
 {
   LANG: C
-  READS:    EulerAuto::DenFlux
-  READS:    EulerAuto::Denrhs
-  READS:    EulerAuto::EnFlux
-  READS:    EulerAuto::Enrhs
-  READS:    EulerAuto::S1Flux
-  READS:    EulerAuto::S1rhs
-  READS:    EulerAuto::S2Flux
-  READS:    EulerAuto::S2rhs
-  READS:    EulerAuto::S3Flux
-  READS:    EulerAuto::S3rhs
-  WRITES:   EulerAuto::Denrhs
-  WRITES:   EulerAuto::Enrhs
-  WRITES:   EulerAuto::S1rhs
-  WRITES:   EulerAuto::S2rhs
-  WRITES:   EulerAuto::S3rhs
+  READS:    EulerAuto::DenFlux(Everywhere)
+  READS:    EulerAuto::EnFlux(Everywhere)
+  READS:    EulerAuto::S1Flux(Everywhere)
+  READS:    EulerAuto::S2Flux(Everywhere)
+  READS:    EulerAuto::S3Flux(Everywhere)
+  READS:    EulerAuto::Denrhs(Everywhere)
+  READS:    EulerAuto::Enrhs(Everywhere)
+  READS:    EulerAuto::S1rhs(Everywhere)
+  READS:    EulerAuto::S2rhs(Everywhere)
+  READS:    EulerAuto::S3rhs(Everywhere)
+  WRITES:   EulerAuto::Denrhs(Interior)
+  WRITES:   EulerAuto::Enrhs(Interior)
+  WRITES:   EulerAuto::S1rhs(Interior)
+  WRITES:   EulerAuto::S2rhs(Interior)
+  WRITES:   EulerAuto::S3rhs(Interior)
 } "eulerauto_cons_calc_rhs_2"
 
 schedule group eulerauto_cons_calc_reconstruct_3_group in MoL_CalcRHS after eulerauto_cons_calc_rhs_2
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_3"
 
 schedule eulerauto_cons_calc_reconstruct_3 in eulerauto_cons_calc_reconstruct_3_group
 {
   LANG: C
-  READS:    EulerAuto::p
-  READS:    EulerAuto::rho
-  READS:    EulerAuto::v1
-  READS:    EulerAuto::v2
-  READS:    EulerAuto::v3
-  WRITES:   EulerAuto::pLeft
-  WRITES:   EulerAuto::pRight
-  WRITES:   EulerAuto::rhoLeft
-  WRITES:   EulerAuto::rhoRight
-  WRITES:   EulerAuto::v1Left
-  WRITES:   EulerAuto::v1Right
-  WRITES:   EulerAuto::v2Left
-  WRITES:   EulerAuto::v2Right
-  WRITES:   EulerAuto::v3Left
-  WRITES:   EulerAuto::v3Right
+  READS:    EulerAuto::p(Everywhere)
+  READS:    EulerAuto::rho(Everywhere)
+  READS:    EulerAuto::v1(Everywhere)
+  READS:    EulerAuto::v2(Everywhere)
+  READS:    EulerAuto::v3(Everywhere)
+  WRITES:   EulerAuto::pLeft(Interior)
+  WRITES:   EulerAuto::pRight(Interior)
+  WRITES:   EulerAuto::rhoLeft(Interior)
+  WRITES:   EulerAuto::rhoRight(Interior)
+  WRITES:   EulerAuto::v1Left(Interior)
+  WRITES:   EulerAuto::v1Right(Interior)
+  WRITES:   EulerAuto::v2Left(Interior)
+  WRITES:   EulerAuto::v2Right(Interior)
+  WRITES:   EulerAuto::v3Left(Interior)
+  WRITES:   EulerAuto::v3Right(Interior)
 } "eulerauto_cons_calc_reconstruct_3"
 
 schedule eulerauto_cons_calc_reconstruct_3_SelectBCs in eulerauto_cons_calc_reconstruct_3_bc_group
@@ -588,77 +571,73 @@ schedule eulerauto_cons_calc_reconstruct_3_SelectBCs in eulerauto_cons_calc_reco
 
 schedule group ApplyBCs as eulerauto_cons_calc_reconstruct_3_ApplyBCs in eulerauto_cons_calc_reconstruct_3_bc_group after eulerauto_cons_calc_reconstruct_3_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_reconstruct_3"
 
 schedule group eulerauto_cons_calc_reconstruct_3_bc_group in eulerauto_cons_calc_reconstruct_3_group after eulerauto_cons_calc_reconstruct_3
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_3"
 
 schedule group eulerauto_cons_calc_reconstruct_3_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_reconstruct_3"
 
 schedule eulerauto_cons_calc_intercell_conserved_3 in MoL_CalcRHS after eulerauto_cons_calc_reconstruct_3
 {
   LANG: C
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenLeft
-  WRITES:   EulerAuto::DenRight
-  WRITES:   EulerAuto::EnLeft
-  WRITES:   EulerAuto::EnRight
-  WRITES:   EulerAuto::S1Left
-  WRITES:   EulerAuto::S1Right
-  WRITES:   EulerAuto::S2Left
-  WRITES:   EulerAuto::S2Right
-  WRITES:   EulerAuto::S3Left
-  WRITES:   EulerAuto::S3Right
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenLeft(Everywhere)
+  WRITES:   EulerAuto::DenRight(Everywhere)
+  WRITES:   EulerAuto::EnLeft(Everywhere)
+  WRITES:   EulerAuto::EnRight(Everywhere)
+  WRITES:   EulerAuto::S1Left(Everywhere)
+  WRITES:   EulerAuto::S1Right(Everywhere)
+  WRITES:   EulerAuto::S2Left(Everywhere)
+  WRITES:   EulerAuto::S2Right(Everywhere)
+  WRITES:   EulerAuto::S3Left(Everywhere)
+  WRITES:   EulerAuto::S3Right(Everywhere)
 } "eulerauto_cons_calc_intercell_conserved_3"
 
 schedule group eulerauto_cons_calc_flux_3_group in MoL_CalcRHS after eulerauto_cons_calc_intercell_conserved_3
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_3"
 
 schedule eulerauto_cons_calc_flux_3 in eulerauto_cons_calc_flux_3_group
 {
   LANG: C
-  READS:    EulerAuto::DenLeft
-  READS:    EulerAuto::DenRight
-  READS:    EulerAuto::EnLeft
-  READS:    EulerAuto::EnRight
-  READS:    EulerAuto::pLeft
-  READS:    EulerAuto::pRight
-  READS:    EulerAuto::rhoLeft
-  READS:    EulerAuto::rhoRight
-  READS:    EulerAuto::S1Left
-  READS:    EulerAuto::S1Right
-  READS:    EulerAuto::S2Left
-  READS:    EulerAuto::S2Right
-  READS:    EulerAuto::S3Left
-  READS:    EulerAuto::S3Right
-  READS:    EulerAuto::v1Left
-  READS:    EulerAuto::v1Right
-  READS:    EulerAuto::v2Left
-  READS:    EulerAuto::v2Right
-  READS:    EulerAuto::v3Left
-  READS:    EulerAuto::v3Right
-  WRITES:   EulerAuto::DenFlux
-  WRITES:   EulerAuto::EnFlux
-  WRITES:   EulerAuto::S1Flux
-  WRITES:   EulerAuto::S2Flux
-  WRITES:   EulerAuto::S3Flux
+  READS:    EulerAuto::DenLeft(Everywhere)
+  READS:    EulerAuto::DenRight(Everywhere)
+  READS:    EulerAuto::EnLeft(Everywhere)
+  READS:    EulerAuto::EnRight(Everywhere)
+  READS:    EulerAuto::pLeft(Everywhere)
+  READS:    EulerAuto::pRight(Everywhere)
+  READS:    EulerAuto::rhoLeft(Everywhere)
+  READS:    EulerAuto::rhoRight(Everywhere)
+  READS:    EulerAuto::S1Left(Everywhere)
+  READS:    EulerAuto::S1Right(Everywhere)
+  READS:    EulerAuto::S2Left(Everywhere)
+  READS:    EulerAuto::S2Right(Everywhere)
+  READS:    EulerAuto::S3Left(Everywhere)
+  READS:    EulerAuto::S3Right(Everywhere)
+  READS:    EulerAuto::v1Left(Everywhere)
+  READS:    EulerAuto::v1Right(Everywhere)
+  READS:    EulerAuto::v2Left(Everywhere)
+  READS:    EulerAuto::v2Right(Everywhere)
+  READS:    EulerAuto::v3Left(Everywhere)
+  READS:    EulerAuto::v3Right(Everywhere)
+  WRITES:   EulerAuto::DenFlux(Interior)
+  WRITES:   EulerAuto::EnFlux(Interior)
+  WRITES:   EulerAuto::S1Flux(Interior)
+  WRITES:   EulerAuto::S2Flux(Interior)
+  WRITES:   EulerAuto::S3Flux(Interior)
 } "eulerauto_cons_calc_flux_3"
 
 schedule eulerauto_cons_calc_flux_3_SelectBCs in eulerauto_cons_calc_flux_3_bc_group
@@ -674,37 +653,34 @@ schedule eulerauto_cons_calc_flux_3_SelectBCs in eulerauto_cons_calc_flux_3_bc_g
 
 schedule group ApplyBCs as eulerauto_cons_calc_flux_3_ApplyBCs in eulerauto_cons_calc_flux_3_bc_group after eulerauto_cons_calc_flux_3_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulerauto_cons_calc_flux_3"
 
 schedule group eulerauto_cons_calc_flux_3_bc_group in eulerauto_cons_calc_flux_3_group after eulerauto_cons_calc_flux_3
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_3"
 
 schedule group eulerauto_cons_calc_flux_3_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulerauto_cons_calc_flux_3"
 
 schedule eulerauto_cons_calc_rhs_3 in MoL_CalcRHS after eulerauto_cons_calc_flux_3
 {
   LANG: C
-  READS:    EulerAuto::DenFlux
-  READS:    EulerAuto::Denrhs
-  READS:    EulerAuto::EnFlux
-  READS:    EulerAuto::Enrhs
-  READS:    EulerAuto::S1Flux
-  READS:    EulerAuto::S1rhs
-  READS:    EulerAuto::S2Flux
-  READS:    EulerAuto::S2rhs
-  READS:    EulerAuto::S3Flux
-  READS:    EulerAuto::S3rhs
-  WRITES:   EulerAuto::Denrhs
-  WRITES:   EulerAuto::Enrhs
-  WRITES:   EulerAuto::S1rhs
-  WRITES:   EulerAuto::S2rhs
-  WRITES:   EulerAuto::S3rhs
+  READS:    EulerAuto::DenFlux(Everywhere)
+  READS:    EulerAuto::EnFlux(Everywhere)
+  READS:    EulerAuto::S1Flux(Everywhere)
+  READS:    EulerAuto::S2Flux(Everywhere)
+  READS:    EulerAuto::S3Flux(Everywhere)
+  READS:    EulerAuto::Denrhs(Everywhere)
+  READS:    EulerAuto::Enrhs(Everywhere)
+  READS:    EulerAuto::S1rhs(Everywhere)
+  READS:    EulerAuto::S2rhs(Everywhere)
+  READS:    EulerAuto::S3rhs(Everywhere)
+  WRITES:   EulerAuto::Denrhs(Interior)
+  WRITES:   EulerAuto::Enrhs(Interior)
+  WRITES:   EulerAuto::S1rhs(Interior)
+  WRITES:   EulerAuto::S2rhs(Interior)
+  WRITES:   EulerAuto::S3rhs(Interior)
 } "eulerauto_cons_calc_rhs_3"
 
 schedule EulerAuto_SelectBoundConds in MoL_PostStep
@@ -730,5 +706,4 @@ schedule EulerAuto_RegisterVars in MoL_Register
 
 schedule group ApplyBCs as EulerAuto_ApplyBCs in MoL_PostStep after EulerAuto_SelectBoundConds
 {
-  # no language specified
 } "Apply boundary conditions controlled by thorn Boundary"
diff --git a/Examples/EulerAuto/src/Boundaries.cc b/Examples/EulerAuto/src/Boundaries.cc
index 8d3fdb0..51485a0 100644
--- a/Examples/EulerAuto/src/Boundaries.cc
+++ b/Examples/EulerAuto/src/Boundaries.cc
@@ -30,7 +30,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   
   if (CCTK_EQUALS(Den_group_bound, "none"  ) ||
       CCTK_EQUALS(Den_group_bound, "static") ||
@@ -123,7 +123,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(Den_group_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_Den_group_bound = -1;
+    static CCTK_INT handle_Den_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_Den_group_bound < 0) handle_Den_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_Den_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_Den_group_bound , Den_group_bound_limit, "LIMIT") < 0)
@@ -142,7 +142,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(En_group_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_En_group_bound = -1;
+    static CCTK_INT handle_En_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_En_group_bound < 0) handle_En_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_En_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_En_group_bound , En_group_bound_limit, "LIMIT") < 0)
@@ -161,7 +161,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S_group_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_S_group_bound = -1;
+    static CCTK_INT handle_S_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S_group_bound < 0) handle_S_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S_group_bound , S_group_bound_limit, "LIMIT") < 0)
@@ -180,7 +180,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(Den_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_Den_bound = -1;
+    static CCTK_INT handle_Den_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_Den_bound < 0) handle_Den_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_Den_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_Den_bound , Den_bound_limit, "LIMIT") < 0)
@@ -199,7 +199,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(En_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_En_bound = -1;
+    static CCTK_INT handle_En_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_En_bound < 0) handle_En_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_En_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_En_bound , En_bound_limit, "LIMIT") < 0)
@@ -218,7 +218,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S1_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_S1_bound = -1;
+    static CCTK_INT handle_S1_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S1_bound < 0) handle_S1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S1_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S1_bound , S1_bound_limit, "LIMIT") < 0)
@@ -237,7 +237,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S2_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_S2_bound = -1;
+    static CCTK_INT handle_S2_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S2_bound < 0) handle_S2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S2_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S2_bound , S2_bound_limit, "LIMIT") < 0)
@@ -256,7 +256,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S3_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_S3_bound = -1;
+    static CCTK_INT handle_S3_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S3_bound < 0) handle_S3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S3_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S3_bound , S3_bound_limit, "LIMIT") < 0)
@@ -275,7 +275,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(Den_group_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_Den_group_bound = -1;
+    static CCTK_INT handle_Den_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_Den_group_bound < 0) handle_Den_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_Den_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_Den_group_bound ,Den_group_bound_scalar, "SCALAR") < 0)
@@ -292,7 +292,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(En_group_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_En_group_bound = -1;
+    static CCTK_INT handle_En_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_En_group_bound < 0) handle_En_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_En_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_En_group_bound ,En_group_bound_scalar, "SCALAR") < 0)
@@ -309,7 +309,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S_group_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_S_group_bound = -1;
+    static CCTK_INT handle_S_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S_group_bound < 0) handle_S_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S_group_bound ,S_group_bound_scalar, "SCALAR") < 0)
@@ -326,7 +326,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(Den_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_Den_bound = -1;
+    static CCTK_INT handle_Den_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_Den_bound < 0) handle_Den_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_Den_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_Den_bound ,Den_bound_scalar, "SCALAR") < 0)
@@ -343,7 +343,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(En_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_En_bound = -1;
+    static CCTK_INT handle_En_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_En_bound < 0) handle_En_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_En_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_En_bound ,En_bound_scalar, "SCALAR") < 0)
@@ -360,7 +360,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S1_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_S1_bound = -1;
+    static CCTK_INT handle_S1_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S1_bound < 0) handle_S1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S1_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S1_bound ,S1_bound_scalar, "SCALAR") < 0)
@@ -377,7 +377,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S2_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_S2_bound = -1;
+    static CCTK_INT handle_S2_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S2_bound < 0) handle_S2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S2_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S2_bound ,S2_bound_scalar, "SCALAR") < 0)
@@ -394,7 +394,7 @@ extern "C" void EulerAuto_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(S3_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_S3_bound = -1;
+    static CCTK_INT handle_S3_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_S3_bound < 0) handle_S3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_S3_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_S3_bound ,S3_bound_scalar, "SCALAR") < 0)
diff --git a/Examples/EulerAuto/src/Differencing.h b/Examples/EulerAuto/src/Differencing.h
index 913cfa8..323c0fe 100644
--- a/Examples/EulerAuto/src/Differencing.h
+++ b/Examples/EulerAuto/src/Differencing.h
@@ -1,72 +1,72 @@
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffPlus1(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1)
+#  define DiffPlus1(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1)
 #else
 #  define DiffPlus1(u) (DiffPlus1_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffPlus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffPlus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffPlus2(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1)
+#  define DiffPlus2(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1)
 #else
 #  define DiffPlus2(u) (DiffPlus2_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffPlus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffPlus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffPlus3(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1)
+#  define DiffPlus3(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1)
 #else
 #  define DiffPlus3(u) (DiffPlus3_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffPlus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffPlus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffMinus1(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,-1,0,0))*p1o1)
+#  define DiffMinus1(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,-1,0,0))*p1o1)
 #else
 #  define DiffMinus1(u) (DiffMinus1_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,-1,0,0))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,-1,0,0))*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffMinus2(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,-1,0))*p1o1)
+#  define DiffMinus2(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,-1,0))*p1o1)
 #else
 #  define DiffMinus2(u) (DiffMinus2_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,-1,0))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,-1,0))*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define DiffMinus3(u) ((KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,0,-1))*p1o1)
+#  define DiffMinus3(u) ((KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,0,-1))*p1o1)
 #else
 #  define DiffMinus3(u) (DiffMinus3_impl(u,p1o1,cdj,cdk))
 static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (KRANC_GFOFFSET3D(u,0,0,0) - KRANC_GFOFFSET3D(u,0,0,-1))*p1o1;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (KRANC_GFOFFSET3D(u,0,0,0) - 1.*KRANC_GFOFFSET3D(u,0,0,-1))*p1o1;
 }
 #endif
 
@@ -77,7 +77,7 @@ static CCTK_REAL DiffMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL co
 static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
   return KRANC_GFOFFSET3D(u,-1,0,0)*p1o1;
 }
 #endif
@@ -89,7 +89,7 @@ static CCTK_REAL ShiftMinus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL c
 static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
   return KRANC_GFOFFSET3D(u,0,-1,0)*p1o1;
 }
 #endif
@@ -101,44 +101,44 @@ static CCTK_REAL ShiftMinus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL c
 static CCTK_REAL ShiftMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL ShiftMinus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1o1, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
   return KRANC_GFOFFSET3D(u,0,0,-1)*p1o1;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define PDplus1(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx)
+#  define PDplus1(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx)
 #else
 #  define PDplus1(u) (PDplus1_impl(u,p1odx,cdj,cdk))
 static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL PDplus1_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odx, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,1,0,0))*p1odx;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define PDplus2(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody)
+#  define PDplus2(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody)
 #else
 #  define PDplus2(u) (PDplus2_impl(u,p1ody,cdj,cdk))
 static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL PDplus2_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1ody, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,1,0))*p1ody;
 }
 #endif
 
 #ifndef KRANC_DIFF_FUNCTIONS
-#  define PDplus3(u) ((-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz)
+#  define PDplus3(u) ((-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz)
 #else
 #  define PDplus3(u) (PDplus3_impl(u,p1odz,cdj,cdk))
 static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk) CCTK_ATTRIBUTE_NOINLINE CCTK_ATTRIBUTE_UNUSED;
 static CCTK_REAL PDplus3_impl(CCTK_REAL const* restrict const u, CCTK_REAL const p1odz, ptrdiff_t const cdj, ptrdiff_t const cdk)
 {
-  ptrdiff_t const cdi=sizeof(CCTK_REAL);
-  return (-KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz;
+  ptrdiff_t const cdi CCTK_ATTRIBUTE_UNUSED = sizeof(CCTK_REAL);
+  return (-1.*KRANC_GFOFFSET3D(u,0,0,0) + KRANC_GFOFFSET3D(u,0,0,1))*p1odz;
 }
 #endif
 
diff --git a/Examples/EulerAuto/src/RegisterMoL.cc b/Examples/EulerAuto/src/RegisterMoL.cc
index 76554bb..f68cb49 100644
--- a/Examples/EulerAuto/src/RegisterMoL.cc
+++ b/Examples/EulerAuto/src/RegisterMoL.cc
@@ -9,7 +9,7 @@ extern "C" void EulerAuto_RegisterVars(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   
   /* Register all the evolved grid functions with MoL */
   ierr += MoLRegisterEvolved(CCTK_VarIndex("EulerAuto::Den"),  CCTK_VarIndex("EulerAuto::Denrhs"));
diff --git a/Examples/EulerAuto/src/Startup.cc b/Examples/EulerAuto/src/Startup.cc
index 419637c..5072260 100644
--- a/Examples/EulerAuto/src/Startup.cc
+++ b/Examples/EulerAuto/src/Startup.cc
@@ -4,7 +4,7 @@
 
 extern "C" int EulerAuto_Startup(void)
 {
-  const char * banner = "EulerAuto";
+  const char * banner CCTK_ATTRIBUTE_UNUSED  = "EulerAuto";
   CCTK_RegisterBanner(banner);
   return 0;
 }
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_conserved.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_conserved.cc
index 9696876..dc1af3b 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_conserved.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_conserved.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_conserved_Body(cGH const * restrict const cctkGH
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,15 +74,15 @@ static void eulerauto_cons_calc_conserved_Body(cGH const * restrict const cctkGH
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pL = p[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[index];
+    CCTK_REAL rhoL CCTK_ATTRIBUTE_UNUSED = rho[index];
+    CCTK_REAL v1L CCTK_ATTRIBUTE_UNUSED = v1[index];
+    CCTK_REAL v2L CCTK_ATTRIBUTE_UNUSED = v2[index];
+    CCTK_REAL v3L CCTK_ATTRIBUTE_UNUSED = v3[index];
     
     
     /* Include user supplied include files */
@@ -90,16 +90,16 @@ static void eulerauto_cons_calc_conserved_Body(cGH const * restrict const cctkGH
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenL = rhoL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenL = rhoL;
     
-    CCTK_REAL S1L = rhoL*v1L;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1L = rhoL*v1L;
     
-    CCTK_REAL S2L = rhoL*v2L;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2L = rhoL*v2L;
     
-    CCTK_REAL S3L = rhoL*v3L;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3L = rhoL*v3L;
     
-    CCTK_REAL EnL = 0.5*INV(-1 + ToReal(gamma))*(2*pL + 
-      rhoL*(SQR(v1L) + SQR(v2L) + SQR(v3L))*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnL = pL*INV(-1. + ToReal(gamma)) + 
+      0.5*rhoL*(SQR(v1L) + SQR(v2L) + SQR(v3L));
     
     /* Copy local copies back to grid functions */
     Den[index] = DenL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_1.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_1.cc
index 01c2e08..e2ed7b6 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_1.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_flux_1_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_flux_1_Body(cGH const * restrict const cctkGH, i
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,85 +98,91 @@ static void eulerauto_cons_calc_flux_1_Body(cGH const * restrict const cctkGH, i
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL EnLeftL = EnLeft[index];
-    CCTK_REAL EnRightL = EnRight[index];
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL S1LeftL = S1Left[index];
-    CCTK_REAL S1RightL = S1Right[index];
-    CCTK_REAL S2LeftL = S2Left[index];
-    CCTK_REAL S2RightL = S2Right[index];
-    CCTK_REAL S3LeftL = S3Left[index];
-    CCTK_REAL S3RightL = S3Right[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL EnLeftL CCTK_ATTRIBUTE_UNUSED = EnLeft[index];
+    CCTK_REAL EnRightL CCTK_ATTRIBUTE_UNUSED = EnRight[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL S1LeftL CCTK_ATTRIBUTE_UNUSED = S1Left[index];
+    CCTK_REAL S1RightL CCTK_ATTRIBUTE_UNUSED = S1Right[index];
+    CCTK_REAL S2LeftL CCTK_ATTRIBUTE_UNUSED = S2Left[index];
+    CCTK_REAL S2RightL CCTK_ATTRIBUTE_UNUSED = S2Right[index];
+    CCTK_REAL S3LeftL CCTK_ATTRIBUTE_UNUSED = S3Left[index];
+    CCTK_REAL S3RightL CCTK_ATTRIBUTE_UNUSED = S3Right[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const ShiftMinus1DenRight = ShiftMinus1(&DenRight[index]);
-    CCTK_REAL const ShiftMinus1EnRight = ShiftMinus1(&EnRight[index]);
-    CCTK_REAL const ShiftMinus1pRight = ShiftMinus1(&pRight[index]);
-    CCTK_REAL const ShiftMinus1rhoRight = ShiftMinus1(&rhoRight[index]);
-    CCTK_REAL const ShiftMinus1S1Right = ShiftMinus1(&S1Right[index]);
-    CCTK_REAL const ShiftMinus1S2Right = ShiftMinus1(&S2Right[index]);
-    CCTK_REAL const ShiftMinus1S3Right = ShiftMinus1(&S3Right[index]);
-    CCTK_REAL const ShiftMinus1v1Right = ShiftMinus1(&v1Right[index]);
-    CCTK_REAL const ShiftMinus1v2Right = ShiftMinus1(&v2Right[index]);
-    CCTK_REAL const ShiftMinus1v3Right = ShiftMinus1(&v3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1EnRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&EnRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1pRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&pRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1rhoRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&rhoRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1S1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&S1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1S2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&S2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1S3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&S3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1v1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&v1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1v2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&v2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1v3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&v3Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = rhoLeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = rhoLeftL*v1LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus1rhoRight*ShiftMinus1v1Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus1rhoRight*ShiftMinus1v1Right;
     
-    CCTK_REAL DenFluxL = 0.5*(DenFluxLeft + DenFluxRight + (-DenLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxL = 0.5*(DenFluxLeft + 
+      DenFluxRight + (-1.*DenLeftL + 
       ShiftMinus1DenRight)*ToReal(hlleAlpha));
     
-    CCTK_REAL S1FluxLeft = pLeftL + rhoLeftL*SQR(v1LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = pLeftL + 
+      rhoLeftL*SQR(v1LeftL);
     
-    CCTK_REAL S1FluxRight = ShiftMinus1pRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = ShiftMinus1pRight + 
       ShiftMinus1rhoRight*SQR(ShiftMinus1v1Right);
     
-    CCTK_REAL S1FluxL = 0.5*(S1FluxLeft + S1FluxRight + (-S1LeftL + 
-      ShiftMinus1S1Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxL = 0.5*(S1FluxLeft + 
+      S1FluxRight + (-1.*S1LeftL + ShiftMinus1S1Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S2FluxLeft = rhoLeftL*v1LeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = 
+      rhoLeftL*v1LeftL*v2LeftL;
     
-    CCTK_REAL S2FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = 
       ShiftMinus1rhoRight*ShiftMinus1v1Right*ShiftMinus1v2Right;
     
-    CCTK_REAL S2FluxL = 0.5*(S2FluxLeft + S2FluxRight + (-S2LeftL + 
-      ShiftMinus1S2Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxL = 0.5*(S2FluxLeft + 
+      S2FluxRight + (-1.*S2LeftL + ShiftMinus1S2Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S3FluxLeft = rhoLeftL*v1LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = 
+      rhoLeftL*v1LeftL*v3LeftL;
     
-    CCTK_REAL S3FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = 
       ShiftMinus1rhoRight*ShiftMinus1v1Right*ShiftMinus1v3Right;
     
-    CCTK_REAL S3FluxL = 0.5*(S3FluxLeft + S3FluxRight + (-S3LeftL + 
-      ShiftMinus1S3Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxL = 0.5*(S3FluxLeft + 
+      S3FluxRight + (-1.*S3LeftL + ShiftMinus1S3Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL EnFluxLeft = (EnLeftL + pLeftL)*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxLeft = (EnLeftL + 
+      pLeftL)*v1LeftL;
     
-    CCTK_REAL EnFluxRight = (ShiftMinus1EnRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxRight = (ShiftMinus1EnRight + 
       ShiftMinus1pRight)*ShiftMinus1v1Right;
     
-    CCTK_REAL EnFluxL = 0.5*(EnFluxLeft + EnFluxRight + (-EnLeftL + 
-      ShiftMinus1EnRight)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxL = 0.5*(EnFluxLeft + 
+      EnFluxRight + (-1.*EnLeftL + ShiftMinus1EnRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
     DenFlux[index] = DenFluxL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_2.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_2.cc
index 8ad72ef..e6fed14 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_2.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_flux_2_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_flux_2_Body(cGH const * restrict const cctkGH, i
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,85 +98,91 @@ static void eulerauto_cons_calc_flux_2_Body(cGH const * restrict const cctkGH, i
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL EnLeftL = EnLeft[index];
-    CCTK_REAL EnRightL = EnRight[index];
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL S1LeftL = S1Left[index];
-    CCTK_REAL S1RightL = S1Right[index];
-    CCTK_REAL S2LeftL = S2Left[index];
-    CCTK_REAL S2RightL = S2Right[index];
-    CCTK_REAL S3LeftL = S3Left[index];
-    CCTK_REAL S3RightL = S3Right[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL EnLeftL CCTK_ATTRIBUTE_UNUSED = EnLeft[index];
+    CCTK_REAL EnRightL CCTK_ATTRIBUTE_UNUSED = EnRight[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL S1LeftL CCTK_ATTRIBUTE_UNUSED = S1Left[index];
+    CCTK_REAL S1RightL CCTK_ATTRIBUTE_UNUSED = S1Right[index];
+    CCTK_REAL S2LeftL CCTK_ATTRIBUTE_UNUSED = S2Left[index];
+    CCTK_REAL S2RightL CCTK_ATTRIBUTE_UNUSED = S2Right[index];
+    CCTK_REAL S3LeftL CCTK_ATTRIBUTE_UNUSED = S3Left[index];
+    CCTK_REAL S3RightL CCTK_ATTRIBUTE_UNUSED = S3Right[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const ShiftMinus2DenRight = ShiftMinus2(&DenRight[index]);
-    CCTK_REAL const ShiftMinus2EnRight = ShiftMinus2(&EnRight[index]);
-    CCTK_REAL const ShiftMinus2pRight = ShiftMinus2(&pRight[index]);
-    CCTK_REAL const ShiftMinus2rhoRight = ShiftMinus2(&rhoRight[index]);
-    CCTK_REAL const ShiftMinus2S1Right = ShiftMinus2(&S1Right[index]);
-    CCTK_REAL const ShiftMinus2S2Right = ShiftMinus2(&S2Right[index]);
-    CCTK_REAL const ShiftMinus2S3Right = ShiftMinus2(&S3Right[index]);
-    CCTK_REAL const ShiftMinus2v1Right = ShiftMinus2(&v1Right[index]);
-    CCTK_REAL const ShiftMinus2v2Right = ShiftMinus2(&v2Right[index]);
-    CCTK_REAL const ShiftMinus2v3Right = ShiftMinus2(&v3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2EnRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&EnRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2pRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&pRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2rhoRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&rhoRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2S1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&S1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2S2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&S2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2S3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&S3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2v1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&v1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2v2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&v2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2v3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&v3Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = rhoLeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = rhoLeftL*v2LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus2rhoRight*ShiftMinus2v2Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus2rhoRight*ShiftMinus2v2Right;
     
-    CCTK_REAL DenFluxL = 0.5*(DenFluxLeft + DenFluxRight + (-DenLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxL = 0.5*(DenFluxLeft + 
+      DenFluxRight + (-1.*DenLeftL + 
       ShiftMinus2DenRight)*ToReal(hlleAlpha));
     
-    CCTK_REAL S1FluxLeft = rhoLeftL*v1LeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = 
+      rhoLeftL*v1LeftL*v2LeftL;
     
-    CCTK_REAL S1FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = 
       ShiftMinus2rhoRight*ShiftMinus2v1Right*ShiftMinus2v2Right;
     
-    CCTK_REAL S1FluxL = 0.5*(S1FluxLeft + S1FluxRight + (-S1LeftL + 
-      ShiftMinus2S1Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxL = 0.5*(S1FluxLeft + 
+      S1FluxRight + (-1.*S1LeftL + ShiftMinus2S1Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S2FluxLeft = pLeftL + rhoLeftL*SQR(v2LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = pLeftL + 
+      rhoLeftL*SQR(v2LeftL);
     
-    CCTK_REAL S2FluxRight = ShiftMinus2pRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = ShiftMinus2pRight + 
       ShiftMinus2rhoRight*SQR(ShiftMinus2v2Right);
     
-    CCTK_REAL S2FluxL = 0.5*(S2FluxLeft + S2FluxRight + (-S2LeftL + 
-      ShiftMinus2S2Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxL = 0.5*(S2FluxLeft + 
+      S2FluxRight + (-1.*S2LeftL + ShiftMinus2S2Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S3FluxLeft = rhoLeftL*v2LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = 
+      rhoLeftL*v2LeftL*v3LeftL;
     
-    CCTK_REAL S3FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = 
       ShiftMinus2rhoRight*ShiftMinus2v2Right*ShiftMinus2v3Right;
     
-    CCTK_REAL S3FluxL = 0.5*(S3FluxLeft + S3FluxRight + (-S3LeftL + 
-      ShiftMinus2S3Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxL = 0.5*(S3FluxLeft + 
+      S3FluxRight + (-1.*S3LeftL + ShiftMinus2S3Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL EnFluxLeft = (EnLeftL + pLeftL)*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxLeft = (EnLeftL + 
+      pLeftL)*v2LeftL;
     
-    CCTK_REAL EnFluxRight = (ShiftMinus2EnRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxRight = (ShiftMinus2EnRight + 
       ShiftMinus2pRight)*ShiftMinus2v2Right;
     
-    CCTK_REAL EnFluxL = 0.5*(EnFluxLeft + EnFluxRight + (-EnLeftL + 
-      ShiftMinus2EnRight)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxL = 0.5*(EnFluxLeft + 
+      EnFluxRight + (-1.*EnLeftL + ShiftMinus2EnRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
     DenFlux[index] = DenFluxL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_3.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_3.cc
index 752fec3..7aed459 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_flux_3.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_flux_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_flux_3_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_flux_3_Body(cGH const * restrict const cctkGH, i
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,85 +98,91 @@ static void eulerauto_cons_calc_flux_3_Body(cGH const * restrict const cctkGH, i
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL EnLeftL = EnLeft[index];
-    CCTK_REAL EnRightL = EnRight[index];
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL S1LeftL = S1Left[index];
-    CCTK_REAL S1RightL = S1Right[index];
-    CCTK_REAL S2LeftL = S2Left[index];
-    CCTK_REAL S2RightL = S2Right[index];
-    CCTK_REAL S3LeftL = S3Left[index];
-    CCTK_REAL S3RightL = S3Right[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL EnLeftL CCTK_ATTRIBUTE_UNUSED = EnLeft[index];
+    CCTK_REAL EnRightL CCTK_ATTRIBUTE_UNUSED = EnRight[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL S1LeftL CCTK_ATTRIBUTE_UNUSED = S1Left[index];
+    CCTK_REAL S1RightL CCTK_ATTRIBUTE_UNUSED = S1Right[index];
+    CCTK_REAL S2LeftL CCTK_ATTRIBUTE_UNUSED = S2Left[index];
+    CCTK_REAL S2RightL CCTK_ATTRIBUTE_UNUSED = S2Right[index];
+    CCTK_REAL S3LeftL CCTK_ATTRIBUTE_UNUSED = S3Left[index];
+    CCTK_REAL S3RightL CCTK_ATTRIBUTE_UNUSED = S3Right[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const ShiftMinus3DenRight = ShiftMinus3(&DenRight[index]);
-    CCTK_REAL const ShiftMinus3EnRight = ShiftMinus3(&EnRight[index]);
-    CCTK_REAL const ShiftMinus3pRight = ShiftMinus3(&pRight[index]);
-    CCTK_REAL const ShiftMinus3rhoRight = ShiftMinus3(&rhoRight[index]);
-    CCTK_REAL const ShiftMinus3S1Right = ShiftMinus3(&S1Right[index]);
-    CCTK_REAL const ShiftMinus3S2Right = ShiftMinus3(&S2Right[index]);
-    CCTK_REAL const ShiftMinus3S3Right = ShiftMinus3(&S3Right[index]);
-    CCTK_REAL const ShiftMinus3v1Right = ShiftMinus3(&v1Right[index]);
-    CCTK_REAL const ShiftMinus3v2Right = ShiftMinus3(&v2Right[index]);
-    CCTK_REAL const ShiftMinus3v3Right = ShiftMinus3(&v3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3EnRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&EnRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3pRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&pRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3rhoRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&rhoRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3S1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&S1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3S2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&S2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3S3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&S3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3v1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&v1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3v2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&v2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3v3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&v3Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = rhoLeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = rhoLeftL*v3LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus3rhoRight*ShiftMinus3v3Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus3rhoRight*ShiftMinus3v3Right;
     
-    CCTK_REAL DenFluxL = 0.5*(DenFluxLeft + DenFluxRight + (-DenLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxL = 0.5*(DenFluxLeft + 
+      DenFluxRight + (-1.*DenLeftL + 
       ShiftMinus3DenRight)*ToReal(hlleAlpha));
     
-    CCTK_REAL S1FluxLeft = rhoLeftL*v1LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = 
+      rhoLeftL*v1LeftL*v3LeftL;
     
-    CCTK_REAL S1FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = 
       ShiftMinus3rhoRight*ShiftMinus3v1Right*ShiftMinus3v3Right;
     
-    CCTK_REAL S1FluxL = 0.5*(S1FluxLeft + S1FluxRight + (-S1LeftL + 
-      ShiftMinus3S1Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxL = 0.5*(S1FluxLeft + 
+      S1FluxRight + (-1.*S1LeftL + ShiftMinus3S1Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S2FluxLeft = rhoLeftL*v2LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = 
+      rhoLeftL*v2LeftL*v3LeftL;
     
-    CCTK_REAL S2FluxRight = 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = 
       ShiftMinus3rhoRight*ShiftMinus3v2Right*ShiftMinus3v3Right;
     
-    CCTK_REAL S2FluxL = 0.5*(S2FluxLeft + S2FluxRight + (-S2LeftL + 
-      ShiftMinus3S2Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxL = 0.5*(S2FluxLeft + 
+      S2FluxRight + (-1.*S2LeftL + ShiftMinus3S2Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL S3FluxLeft = pLeftL + rhoLeftL*SQR(v3LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = pLeftL + 
+      rhoLeftL*SQR(v3LeftL);
     
-    CCTK_REAL S3FluxRight = ShiftMinus3pRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = ShiftMinus3pRight + 
       ShiftMinus3rhoRight*SQR(ShiftMinus3v3Right);
     
-    CCTK_REAL S3FluxL = 0.5*(S3FluxLeft + S3FluxRight + (-S3LeftL + 
-      ShiftMinus3S3Right)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxL = 0.5*(S3FluxLeft + 
+      S3FluxRight + (-1.*S3LeftL + ShiftMinus3S3Right)*ToReal(hlleAlpha));
     
-    CCTK_REAL EnFluxLeft = (EnLeftL + pLeftL)*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxLeft = (EnLeftL + 
+      pLeftL)*v3LeftL;
     
-    CCTK_REAL EnFluxRight = (ShiftMinus3EnRight + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxRight = (ShiftMinus3EnRight + 
       ShiftMinus3pRight)*ShiftMinus3v3Right;
     
-    CCTK_REAL EnFluxL = 0.5*(EnFluxLeft + EnFluxRight + (-EnLeftL + 
-      ShiftMinus3EnRight)*ToReal(hlleAlpha));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnFluxL = 0.5*(EnFluxLeft + 
+      EnFluxRight + (-1.*EnLeftL + ShiftMinus3EnRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
     DenFlux[index] = DenFluxL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_1.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_1.cc
index 6f22748..4370d97 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_1.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_1.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_intercell_conserved_1_Body(cGH const * restrict
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,20 +74,20 @@ static void eulerauto_cons_calc_intercell_conserved_1_Body(cGH const * restrict
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
@@ -95,29 +95,29 @@ static void eulerauto_cons_calc_intercell_conserved_1_Body(cGH const * restrict
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenLeftL = rhoLeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenLeftL = rhoLeftL;
     
-    CCTK_REAL S1LeftL = rhoLeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = rhoLeftL*v1LeftL;
     
-    CCTK_REAL S2LeftL = rhoLeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = rhoLeftL*v2LeftL;
     
-    CCTK_REAL S3LeftL = rhoLeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = rhoLeftL*v3LeftL;
     
-    CCTK_REAL EnLeftL = 0.5*INV(-1 + ToReal(gamma))*(2*pLeftL + 
-      rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + SQR(v3LeftL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnLeftL = pLeftL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + 
+      SQR(v3LeftL));
     
-    CCTK_REAL DenRightL = rhoRightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenRightL = rhoRightL;
     
-    CCTK_REAL S1RightL = rhoRightL*v1RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = rhoRightL*v1RightL;
     
-    CCTK_REAL S2RightL = rhoRightL*v2RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = rhoRightL*v2RightL;
     
-    CCTK_REAL S3RightL = rhoRightL*v3RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = rhoRightL*v3RightL;
     
-    CCTK_REAL EnRightL = 0.5*INV(-1 + ToReal(gamma))*(2*pRightL + 
-      rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnRightL = pRightL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + 
+      SQR(v3RightL));
     
     /* Copy local copies back to grid functions */
     DenLeft[index] = DenLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_2.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_2.cc
index 3d9488e..820bdfb 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_2.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_2.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_intercell_conserved_2_Body(cGH const * restrict
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,20 +74,20 @@ static void eulerauto_cons_calc_intercell_conserved_2_Body(cGH const * restrict
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
@@ -95,29 +95,29 @@ static void eulerauto_cons_calc_intercell_conserved_2_Body(cGH const * restrict
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenLeftL = rhoLeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenLeftL = rhoLeftL;
     
-    CCTK_REAL S1LeftL = rhoLeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = rhoLeftL*v1LeftL;
     
-    CCTK_REAL S2LeftL = rhoLeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = rhoLeftL*v2LeftL;
     
-    CCTK_REAL S3LeftL = rhoLeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = rhoLeftL*v3LeftL;
     
-    CCTK_REAL EnLeftL = 0.5*INV(-1 + ToReal(gamma))*(2*pLeftL + 
-      rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + SQR(v3LeftL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnLeftL = pLeftL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + 
+      SQR(v3LeftL));
     
-    CCTK_REAL DenRightL = rhoRightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenRightL = rhoRightL;
     
-    CCTK_REAL S1RightL = rhoRightL*v1RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = rhoRightL*v1RightL;
     
-    CCTK_REAL S2RightL = rhoRightL*v2RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = rhoRightL*v2RightL;
     
-    CCTK_REAL S3RightL = rhoRightL*v3RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = rhoRightL*v3RightL;
     
-    CCTK_REAL EnRightL = 0.5*INV(-1 + ToReal(gamma))*(2*pRightL + 
-      rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnRightL = pRightL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + 
+      SQR(v3RightL));
     
     /* Copy local copies back to grid functions */
     DenLeft[index] = DenLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_3.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_3.cc
index 096f13c..f6e7174 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_3.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_intercell_conserved_3.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_intercell_conserved_3_Body(cGH const * restrict
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,20 +74,20 @@ static void eulerauto_cons_calc_intercell_conserved_3_Body(cGH const * restrict
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pLeftL = pLeft[index];
-    CCTK_REAL pRightL = pRight[index];
-    CCTK_REAL rhoLeftL = rhoLeft[index];
-    CCTK_REAL rhoRightL = rhoRight[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
-    CCTK_REAL v3LeftL = v3Left[index];
-    CCTK_REAL v3RightL = v3Right[index];
+    CCTK_REAL pLeftL CCTK_ATTRIBUTE_UNUSED = pLeft[index];
+    CCTK_REAL pRightL CCTK_ATTRIBUTE_UNUSED = pRight[index];
+    CCTK_REAL rhoLeftL CCTK_ATTRIBUTE_UNUSED = rhoLeft[index];
+    CCTK_REAL rhoRightL CCTK_ATTRIBUTE_UNUSED = rhoRight[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
+    CCTK_REAL v3LeftL CCTK_ATTRIBUTE_UNUSED = v3Left[index];
+    CCTK_REAL v3RightL CCTK_ATTRIBUTE_UNUSED = v3Right[index];
     
     
     /* Include user supplied include files */
@@ -95,29 +95,29 @@ static void eulerauto_cons_calc_intercell_conserved_3_Body(cGH const * restrict
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenLeftL = rhoLeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenLeftL = rhoLeftL;
     
-    CCTK_REAL S1LeftL = rhoLeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = rhoLeftL*v1LeftL;
     
-    CCTK_REAL S2LeftL = rhoLeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = rhoLeftL*v2LeftL;
     
-    CCTK_REAL S3LeftL = rhoLeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = rhoLeftL*v3LeftL;
     
-    CCTK_REAL EnLeftL = 0.5*INV(-1 + ToReal(gamma))*(2*pLeftL + 
-      rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + SQR(v3LeftL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnLeftL = pLeftL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoLeftL*(SQR(v1LeftL) + SQR(v2LeftL) + 
+      SQR(v3LeftL));
     
-    CCTK_REAL DenRightL = rhoRightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenRightL = rhoRightL;
     
-    CCTK_REAL S1RightL = rhoRightL*v1RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = rhoRightL*v1RightL;
     
-    CCTK_REAL S2RightL = rhoRightL*v2RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = rhoRightL*v2RightL;
     
-    CCTK_REAL S3RightL = rhoRightL*v3RightL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = rhoRightL*v3RightL;
     
-    CCTK_REAL EnRightL = 0.5*INV(-1 + ToReal(gamma))*(2*pRightL + 
-      rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL))*(-1 + 
-      ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnRightL = pRightL*INV(-1. + 
+      ToReal(gamma)) + 0.5*rhoRightL*(SQR(v1RightL) + SQR(v2RightL) + 
+      SQR(v3RightL));
     
     /* Copy local copies back to grid functions */
     DenLeft[index] = DenLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_primitives.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_primitives.cc
index b2cafeb..bcd991c 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_primitives.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_primitives.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_primitives_Body(cGH const * restrict const cctkG
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,15 +74,15 @@ static void eulerauto_cons_calc_primitives_Body(cGH const * restrict const cctkG
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenL = Den[index];
-    CCTK_REAL EnL = En[index];
-    CCTK_REAL S1L = S1[index];
-    CCTK_REAL S2L = S2[index];
-    CCTK_REAL S3L = S3[index];
+    CCTK_REAL DenL CCTK_ATTRIBUTE_UNUSED = Den[index];
+    CCTK_REAL EnL CCTK_ATTRIBUTE_UNUSED = En[index];
+    CCTK_REAL S1L CCTK_ATTRIBUTE_UNUSED = S1[index];
+    CCTK_REAL S2L CCTK_ATTRIBUTE_UNUSED = S2[index];
+    CCTK_REAL S3L CCTK_ATTRIBUTE_UNUSED = S3[index];
     
     
     /* Include user supplied include files */
@@ -90,16 +90,17 @@ static void eulerauto_cons_calc_primitives_Body(cGH const * restrict const cctkG
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL rhoL = DenL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoL = DenL;
     
-    CCTK_REAL v1L = S1L*INV(DenL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1L = S1L*INV(DenL);
     
-    CCTK_REAL v2L = S2L*INV(DenL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2L = S2L*INV(DenL);
     
-    CCTK_REAL v3L = S3L*INV(DenL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3L = S3L*INV(DenL);
     
-    CCTK_REAL pL = 0.5*INV(DenL)*(2*DenL*EnL - SQR(S1L) - 
-      SQR(S2L) - SQR(S3L))*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pL = (EnL - 
+      0.5*INV(DenL)*(SQR(S1L) + SQR(S2L) + SQR(S3L)))*(-1. + 
+      ToReal(gamma));
     
     /* Copy local copies back to grid functions */
     p[index] = pL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_1.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_1.cc
index f0dbe5f..4ab6b82 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_1.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_reconstruct_1_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::p_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::p_lr_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,41 +98,41 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pL = p[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[index];
+    CCTK_REAL rhoL CCTK_ATTRIBUTE_UNUSED = rho[index];
+    CCTK_REAL v1L CCTK_ATTRIBUTE_UNUSED = v1[index];
+    CCTK_REAL v2L CCTK_ATTRIBUTE_UNUSED = v2[index];
+    CCTK_REAL v3L CCTK_ATTRIBUTE_UNUSED = v3[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const DiffPlus1p = DiffPlus1(&p[index]);
-    CCTK_REAL const DiffMinus1p = DiffMinus1(&p[index]);
-    CCTK_REAL const DiffPlus1rho = DiffPlus1(&rho[index]);
-    CCTK_REAL const DiffMinus1rho = DiffMinus1(&rho[index]);
-    CCTK_REAL const DiffPlus1v1 = DiffPlus1(&v1[index]);
-    CCTK_REAL const DiffMinus1v1 = DiffMinus1(&v1[index]);
-    CCTK_REAL const DiffPlus1v2 = DiffPlus1(&v2[index]);
-    CCTK_REAL const DiffMinus1v2 = DiffMinus1(&v2[index]);
-    CCTK_REAL const DiffPlus1v3 = DiffPlus1(&v3[index]);
-    CCTK_REAL const DiffMinus1v3 = DiffMinus1(&v3[index]);
+    CCTK_REAL /*const*/ DiffPlus1p CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&p[index]);
+    CCTK_REAL /*const*/ DiffMinus1p CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&p[index]);
+    CCTK_REAL /*const*/ DiffPlus1rho CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&rho[index]);
+    CCTK_REAL /*const*/ DiffMinus1rho CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&rho[index]);
+    CCTK_REAL /*const*/ DiffPlus1v1 CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&v1[index]);
+    CCTK_REAL /*const*/ DiffMinus1v1 CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&v1[index]);
+    CCTK_REAL /*const*/ DiffPlus1v2 CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&v2[index]);
+    CCTK_REAL /*const*/ DiffMinus1v2 CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&v2[index]);
+    CCTK_REAL /*const*/ DiffPlus1v3 CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&v3[index]);
+    CCTK_REAL /*const*/ DiffMinus1v3 CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&v3[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL slopeL = DiffMinus1rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeL = DiffMinus1rho;
     
-    CCTK_REAL slopeR = DiffPlus1rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeR = DiffPlus1rho;
     
-    CCTK_REAL slope = VanLeer(slopeL,slopeR);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL rhoLeftL = rhoL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoLeftL = rhoL - 0.5*slope;
     
-    CCTK_REAL rhoRightL = rhoL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoRightL = rhoL + 0.5*slope;
     
     slopeL = DiffMinus1v1;
     
@@ -140,9 +140,9 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v1LeftL = v1L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1LeftL = v1L - 0.5*slope;
     
-    CCTK_REAL v1RightL = v1L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1RightL = v1L + 0.5*slope;
     
     slopeL = DiffMinus1v2;
     
@@ -150,9 +150,9 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v2LeftL = v2L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2LeftL = v2L - 0.5*slope;
     
-    CCTK_REAL v2RightL = v2L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2RightL = v2L + 0.5*slope;
     
     slopeL = DiffMinus1v3;
     
@@ -160,9 +160,9 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v3LeftL = v3L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3LeftL = v3L - 0.5*slope;
     
-    CCTK_REAL v3RightL = v3L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3RightL = v3L + 0.5*slope;
     
     slopeL = DiffMinus1p;
     
@@ -170,9 +170,9 @@ static void eulerauto_cons_calc_reconstruct_1_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL pLeftL = pL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pLeftL = pL - 0.5*slope;
     
-    CCTK_REAL pRightL = pL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pRightL = pL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     pLeft[index] = pLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_2.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_2.cc
index f47b0d7..e445f17 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_2.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_reconstruct_2_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::p_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::p_lr_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,41 +98,41 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pL = p[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[index];
+    CCTK_REAL rhoL CCTK_ATTRIBUTE_UNUSED = rho[index];
+    CCTK_REAL v1L CCTK_ATTRIBUTE_UNUSED = v1[index];
+    CCTK_REAL v2L CCTK_ATTRIBUTE_UNUSED = v2[index];
+    CCTK_REAL v3L CCTK_ATTRIBUTE_UNUSED = v3[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const DiffPlus2p = DiffPlus2(&p[index]);
-    CCTK_REAL const DiffMinus2p = DiffMinus2(&p[index]);
-    CCTK_REAL const DiffPlus2rho = DiffPlus2(&rho[index]);
-    CCTK_REAL const DiffMinus2rho = DiffMinus2(&rho[index]);
-    CCTK_REAL const DiffPlus2v1 = DiffPlus2(&v1[index]);
-    CCTK_REAL const DiffMinus2v1 = DiffMinus2(&v1[index]);
-    CCTK_REAL const DiffPlus2v2 = DiffPlus2(&v2[index]);
-    CCTK_REAL const DiffMinus2v2 = DiffMinus2(&v2[index]);
-    CCTK_REAL const DiffPlus2v3 = DiffPlus2(&v3[index]);
-    CCTK_REAL const DiffMinus2v3 = DiffMinus2(&v3[index]);
+    CCTK_REAL /*const*/ DiffPlus2p CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&p[index]);
+    CCTK_REAL /*const*/ DiffMinus2p CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&p[index]);
+    CCTK_REAL /*const*/ DiffPlus2rho CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&rho[index]);
+    CCTK_REAL /*const*/ DiffMinus2rho CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&rho[index]);
+    CCTK_REAL /*const*/ DiffPlus2v1 CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&v1[index]);
+    CCTK_REAL /*const*/ DiffMinus2v1 CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&v1[index]);
+    CCTK_REAL /*const*/ DiffPlus2v2 CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&v2[index]);
+    CCTK_REAL /*const*/ DiffMinus2v2 CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&v2[index]);
+    CCTK_REAL /*const*/ DiffPlus2v3 CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&v3[index]);
+    CCTK_REAL /*const*/ DiffMinus2v3 CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&v3[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL slopeL = DiffMinus2rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeL = DiffMinus2rho;
     
-    CCTK_REAL slopeR = DiffPlus2rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeR = DiffPlus2rho;
     
-    CCTK_REAL slope = VanLeer(slopeL,slopeR);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL rhoLeftL = rhoL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoLeftL = rhoL - 0.5*slope;
     
-    CCTK_REAL rhoRightL = rhoL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoRightL = rhoL + 0.5*slope;
     
     slopeL = DiffMinus2v1;
     
@@ -140,9 +140,9 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v1LeftL = v1L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1LeftL = v1L - 0.5*slope;
     
-    CCTK_REAL v1RightL = v1L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1RightL = v1L + 0.5*slope;
     
     slopeL = DiffMinus2v2;
     
@@ -150,9 +150,9 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v2LeftL = v2L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2LeftL = v2L - 0.5*slope;
     
-    CCTK_REAL v2RightL = v2L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2RightL = v2L + 0.5*slope;
     
     slopeL = DiffMinus2v3;
     
@@ -160,9 +160,9 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v3LeftL = v3L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3LeftL = v3L - 0.5*slope;
     
-    CCTK_REAL v3RightL = v3L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3RightL = v3L + 0.5*slope;
     
     slopeL = DiffMinus2p;
     
@@ -170,9 +170,9 @@ static void eulerauto_cons_calc_reconstruct_2_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL pLeftL = pL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pLeftL = pL - 0.5*slope;
     
-    CCTK_REAL pRightL = pL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pRightL = pL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     pLeft[index] = pLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_3.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_3.cc
index 3075063..66e7159 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_3.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_reconstruct_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_reconstruct_3_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::p_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::p_lr_group.");
@@ -55,34 +55,34 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -98,41 +98,41 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL pL = p[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[index];
+    CCTK_REAL rhoL CCTK_ATTRIBUTE_UNUSED = rho[index];
+    CCTK_REAL v1L CCTK_ATTRIBUTE_UNUSED = v1[index];
+    CCTK_REAL v2L CCTK_ATTRIBUTE_UNUSED = v2[index];
+    CCTK_REAL v3L CCTK_ATTRIBUTE_UNUSED = v3[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const DiffPlus3p = DiffPlus3(&p[index]);
-    CCTK_REAL const DiffMinus3p = DiffMinus3(&p[index]);
-    CCTK_REAL const DiffPlus3rho = DiffPlus3(&rho[index]);
-    CCTK_REAL const DiffMinus3rho = DiffMinus3(&rho[index]);
-    CCTK_REAL const DiffPlus3v1 = DiffPlus3(&v1[index]);
-    CCTK_REAL const DiffMinus3v1 = DiffMinus3(&v1[index]);
-    CCTK_REAL const DiffPlus3v2 = DiffPlus3(&v2[index]);
-    CCTK_REAL const DiffMinus3v2 = DiffMinus3(&v2[index]);
-    CCTK_REAL const DiffPlus3v3 = DiffPlus3(&v3[index]);
-    CCTK_REAL const DiffMinus3v3 = DiffMinus3(&v3[index]);
+    CCTK_REAL /*const*/ DiffPlus3p CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&p[index]);
+    CCTK_REAL /*const*/ DiffMinus3p CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&p[index]);
+    CCTK_REAL /*const*/ DiffPlus3rho CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&rho[index]);
+    CCTK_REAL /*const*/ DiffMinus3rho CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&rho[index]);
+    CCTK_REAL /*const*/ DiffPlus3v1 CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&v1[index]);
+    CCTK_REAL /*const*/ DiffMinus3v1 CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&v1[index]);
+    CCTK_REAL /*const*/ DiffPlus3v2 CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&v2[index]);
+    CCTK_REAL /*const*/ DiffMinus3v2 CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&v2[index]);
+    CCTK_REAL /*const*/ DiffPlus3v3 CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&v3[index]);
+    CCTK_REAL /*const*/ DiffMinus3v3 CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&v3[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL slopeL = DiffMinus3rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeL = DiffMinus3rho;
     
-    CCTK_REAL slopeR = DiffPlus3rho;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slopeR = DiffPlus3rho;
     
-    CCTK_REAL slope = VanLeer(slopeL,slopeR);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL rhoLeftL = rhoL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoLeftL = rhoL - 0.5*slope;
     
-    CCTK_REAL rhoRightL = rhoL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoRightL = rhoL + 0.5*slope;
     
     slopeL = DiffMinus3v1;
     
@@ -140,9 +140,9 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v1LeftL = v1L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1LeftL = v1L - 0.5*slope;
     
-    CCTK_REAL v1RightL = v1L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1RightL = v1L + 0.5*slope;
     
     slopeL = DiffMinus3v2;
     
@@ -150,9 +150,9 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v2LeftL = v2L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2LeftL = v2L - 0.5*slope;
     
-    CCTK_REAL v2RightL = v2L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2RightL = v2L + 0.5*slope;
     
     slopeL = DiffMinus3v3;
     
@@ -160,9 +160,9 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL v3LeftL = v3L - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3LeftL = v3L - 0.5*slope;
     
-    CCTK_REAL v3RightL = v3L + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3RightL = v3L + 0.5*slope;
     
     slopeL = DiffMinus3p;
     
@@ -170,9 +170,9 @@ static void eulerauto_cons_calc_reconstruct_3_Body(cGH const * restrict const cc
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL pLeftL = pL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pLeftL = pL - 0.5*slope;
     
-    CCTK_REAL pRightL = pL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pRightL = pL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     pLeft[index] = pLeftL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_1.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_1.cc
index 5f1e9ed..863aaa3 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_1.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_rhs_1_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulerauto_cons_calc_rhs_1_Body(cGH const * restrict const cctkGH, in
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -92,41 +92,41 @@ static void eulerauto_cons_calc_rhs_1_Body(cGH const * restrict const cctkGH, in
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[index];
-    CCTK_REAL EnFluxL = EnFlux[index];
-    CCTK_REAL EnrhsL = Enrhs[index];
-    CCTK_REAL S1FluxL = S1Flux[index];
-    CCTK_REAL S1rhsL = S1rhs[index];
-    CCTK_REAL S2FluxL = S2Flux[index];
-    CCTK_REAL S2rhsL = S2rhs[index];
-    CCTK_REAL S3FluxL = S3Flux[index];
-    CCTK_REAL S3rhsL = S3rhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[index];
+    CCTK_REAL EnFluxL CCTK_ATTRIBUTE_UNUSED = EnFlux[index];
+    CCTK_REAL EnrhsL CCTK_ATTRIBUTE_UNUSED = Enrhs[index];
+    CCTK_REAL S1FluxL CCTK_ATTRIBUTE_UNUSED = S1Flux[index];
+    CCTK_REAL S1rhsL CCTK_ATTRIBUTE_UNUSED = S1rhs[index];
+    CCTK_REAL S2FluxL CCTK_ATTRIBUTE_UNUSED = S2Flux[index];
+    CCTK_REAL S2rhsL CCTK_ATTRIBUTE_UNUSED = S2rhs[index];
+    CCTK_REAL S3FluxL CCTK_ATTRIBUTE_UNUSED = S3Flux[index];
+    CCTK_REAL S3rhsL CCTK_ATTRIBUTE_UNUSED = S3rhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus1DenFlux = PDplus1(&DenFlux[index]);
-    CCTK_REAL const PDplus1EnFlux = PDplus1(&EnFlux[index]);
-    CCTK_REAL const PDplus1S1Flux = PDplus1(&S1Flux[index]);
-    CCTK_REAL const PDplus1S2Flux = PDplus1(&S2Flux[index]);
-    CCTK_REAL const PDplus1S3Flux = PDplus1(&S3Flux[index]);
+    CCTK_REAL /*const*/ PDplus1DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&DenFlux[index]);
+    CCTK_REAL /*const*/ PDplus1EnFlux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&EnFlux[index]);
+    CCTK_REAL /*const*/ PDplus1S1Flux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&S1Flux[index]);
+    CCTK_REAL /*const*/ PDplus1S2Flux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&S2Flux[index]);
+    CCTK_REAL /*const*/ PDplus1S3Flux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&S3Flux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus1DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus1DenFlux;
     
-    EnrhsL = EnrhsL - PDplus1EnFlux;
+    EnrhsL = EnrhsL - 1.*PDplus1EnFlux;
     
-    S1rhsL = S1rhsL - PDplus1S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus1S1Flux;
     
-    S2rhsL = S2rhsL - PDplus1S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus1S2Flux;
     
-    S3rhsL = S3rhsL - PDplus1S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus1S3Flux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_2.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_2.cc
index 796a85b..e185672 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_2.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_rhs_2_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulerauto_cons_calc_rhs_2_Body(cGH const * restrict const cctkGH, in
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -92,41 +92,41 @@ static void eulerauto_cons_calc_rhs_2_Body(cGH const * restrict const cctkGH, in
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[index];
-    CCTK_REAL EnFluxL = EnFlux[index];
-    CCTK_REAL EnrhsL = Enrhs[index];
-    CCTK_REAL S1FluxL = S1Flux[index];
-    CCTK_REAL S1rhsL = S1rhs[index];
-    CCTK_REAL S2FluxL = S2Flux[index];
-    CCTK_REAL S2rhsL = S2rhs[index];
-    CCTK_REAL S3FluxL = S3Flux[index];
-    CCTK_REAL S3rhsL = S3rhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[index];
+    CCTK_REAL EnFluxL CCTK_ATTRIBUTE_UNUSED = EnFlux[index];
+    CCTK_REAL EnrhsL CCTK_ATTRIBUTE_UNUSED = Enrhs[index];
+    CCTK_REAL S1FluxL CCTK_ATTRIBUTE_UNUSED = S1Flux[index];
+    CCTK_REAL S1rhsL CCTK_ATTRIBUTE_UNUSED = S1rhs[index];
+    CCTK_REAL S2FluxL CCTK_ATTRIBUTE_UNUSED = S2Flux[index];
+    CCTK_REAL S2rhsL CCTK_ATTRIBUTE_UNUSED = S2rhs[index];
+    CCTK_REAL S3FluxL CCTK_ATTRIBUTE_UNUSED = S3Flux[index];
+    CCTK_REAL S3rhsL CCTK_ATTRIBUTE_UNUSED = S3rhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus2DenFlux = PDplus2(&DenFlux[index]);
-    CCTK_REAL const PDplus2EnFlux = PDplus2(&EnFlux[index]);
-    CCTK_REAL const PDplus2S1Flux = PDplus2(&S1Flux[index]);
-    CCTK_REAL const PDplus2S2Flux = PDplus2(&S2Flux[index]);
-    CCTK_REAL const PDplus2S3Flux = PDplus2(&S3Flux[index]);
+    CCTK_REAL /*const*/ PDplus2DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&DenFlux[index]);
+    CCTK_REAL /*const*/ PDplus2EnFlux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&EnFlux[index]);
+    CCTK_REAL /*const*/ PDplus2S1Flux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&S1Flux[index]);
+    CCTK_REAL /*const*/ PDplus2S2Flux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&S2Flux[index]);
+    CCTK_REAL /*const*/ PDplus2S3Flux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&S3Flux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus2DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus2DenFlux;
     
-    EnrhsL = EnrhsL - PDplus2EnFlux;
+    EnrhsL = EnrhsL - 1.*PDplus2EnFlux;
     
-    S1rhsL = S1rhsL - PDplus2S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus2S1Flux;
     
-    S2rhsL = S2rhsL - PDplus2S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus2S2Flux;
     
-    S3rhsL = S3rhsL - PDplus2S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus2S3Flux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_3.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_3.cc
index 3aca46c..eab8718 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_3.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_rhs_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulerauto_cons_calc_rhs_3_SelectBCs(CCTK_ARGUMENTS)
   DECLARE_CCTK_ARGUMENTS;
   DECLARE_CCTK_PARAMETERS;
   
-  CCTK_INT ierr = 0;
+  CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED  = 0;
   ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "EulerAuto::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerAuto::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulerauto_cons_calc_rhs_3_Body(cGH const * restrict const cctkGH, in
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -92,41 +92,41 @@ static void eulerauto_cons_calc_rhs_3_Body(cGH const * restrict const cctkGH, in
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[index];
-    CCTK_REAL EnFluxL = EnFlux[index];
-    CCTK_REAL EnrhsL = Enrhs[index];
-    CCTK_REAL S1FluxL = S1Flux[index];
-    CCTK_REAL S1rhsL = S1rhs[index];
-    CCTK_REAL S2FluxL = S2Flux[index];
-    CCTK_REAL S2rhsL = S2rhs[index];
-    CCTK_REAL S3FluxL = S3Flux[index];
-    CCTK_REAL S3rhsL = S3rhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[index];
+    CCTK_REAL EnFluxL CCTK_ATTRIBUTE_UNUSED = EnFlux[index];
+    CCTK_REAL EnrhsL CCTK_ATTRIBUTE_UNUSED = Enrhs[index];
+    CCTK_REAL S1FluxL CCTK_ATTRIBUTE_UNUSED = S1Flux[index];
+    CCTK_REAL S1rhsL CCTK_ATTRIBUTE_UNUSED = S1rhs[index];
+    CCTK_REAL S2FluxL CCTK_ATTRIBUTE_UNUSED = S2Flux[index];
+    CCTK_REAL S2rhsL CCTK_ATTRIBUTE_UNUSED = S2rhs[index];
+    CCTK_REAL S3FluxL CCTK_ATTRIBUTE_UNUSED = S3Flux[index];
+    CCTK_REAL S3rhsL CCTK_ATTRIBUTE_UNUSED = S3rhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus3DenFlux = PDplus3(&DenFlux[index]);
-    CCTK_REAL const PDplus3EnFlux = PDplus3(&EnFlux[index]);
-    CCTK_REAL const PDplus3S1Flux = PDplus3(&S1Flux[index]);
-    CCTK_REAL const PDplus3S2Flux = PDplus3(&S2Flux[index]);
-    CCTK_REAL const PDplus3S3Flux = PDplus3(&S3Flux[index]);
+    CCTK_REAL /*const*/ PDplus3DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&DenFlux[index]);
+    CCTK_REAL /*const*/ PDplus3EnFlux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&EnFlux[index]);
+    CCTK_REAL /*const*/ PDplus3S1Flux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&S1Flux[index]);
+    CCTK_REAL /*const*/ PDplus3S2Flux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&S2Flux[index]);
+    CCTK_REAL /*const*/ PDplus3S3Flux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&S3Flux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus3DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus3DenFlux;
     
-    EnrhsL = EnrhsL - PDplus3EnFlux;
+    EnrhsL = EnrhsL - 1.*PDplus3EnFlux;
     
-    S1rhsL = S1rhsL - PDplus3S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus3S1Flux;
     
-    S2rhsL = S2rhsL - PDplus3S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus3S2Flux;
     
-    S3rhsL = S3rhsL - PDplus3S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus3S3Flux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerAuto/src/eulerauto_cons_calc_zero_rhs.cc b/Examples/EulerAuto/src/eulerauto_cons_calc_zero_rhs.cc
index 267d6c6..fc655a5 100644
--- a/Examples/EulerAuto/src/eulerauto_cons_calc_zero_rhs.cc
+++ b/Examples/EulerAuto/src/eulerauto_cons_calc_zero_rhs.cc
@@ -31,34 +31,34 @@ static void eulerauto_cons_calc_zero_rhs_Body(cGH const * restrict const cctkGH,
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,7 +74,7 @@ static void eulerauto_cons_calc_zero_rhs_Body(cGH const * restrict const cctkGH,
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
@@ -85,35 +85,35 @@ static void eulerauto_cons_calc_zero_rhs_Body(cGH const * restrict const cctkGH,
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenrhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenrhsL = 0.;
     
-    DenrhsL = 0;
+    DenrhsL = 0.;
     
-    DenrhsL = 0;
+    DenrhsL = 0.;
     
-    CCTK_REAL S1rhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1rhsL = 0.;
     
-    CCTK_REAL S2rhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2rhsL = 0.;
     
-    CCTK_REAL S3rhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3rhsL = 0.;
     
-    S1rhsL = 0;
+    S1rhsL = 0.;
     
-    S2rhsL = 0;
+    S2rhsL = 0.;
     
-    S3rhsL = 0;
+    S3rhsL = 0.;
     
-    S1rhsL = 0;
+    S1rhsL = 0.;
     
-    S2rhsL = 0;
+    S2rhsL = 0.;
     
-    S3rhsL = 0;
+    S3rhsL = 0.;
     
-    CCTK_REAL EnrhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED EnrhsL = 0.;
     
-    EnrhsL = 0;
+    EnrhsL = 0.;
     
-    EnrhsL = 0;
+    EnrhsL = 0.;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerAuto/src/eulerauto_initial_shock.cc b/Examples/EulerAuto/src/eulerauto_initial_shock.cc
index 3becd2f..bd086e5 100644
--- a/Examples/EulerAuto/src/eulerauto_initial_shock.cc
+++ b/Examples/EulerAuto/src/eulerauto_initial_shock.cc
@@ -31,34 +31,34 @@ static void eulerauto_initial_shock_Body(cGH const * restrict const cctkGH, int
   /* Include user-supplied include files */
   
   /* Initialise finite differencing variables */
-  ptrdiff_t const di = 1;
-  ptrdiff_t const dj = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const dk = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
-  ptrdiff_t const cdi = sizeof(CCTK_REAL) * di;
-  ptrdiff_t const cdj = sizeof(CCTK_REAL) * dj;
-  ptrdiff_t const cdk = sizeof(CCTK_REAL) * dk;
-  CCTK_REAL const dx = ToReal(CCTK_DELTA_SPACE(0));
-  CCTK_REAL const dy = ToReal(CCTK_DELTA_SPACE(1));
-  CCTK_REAL const dz = ToReal(CCTK_DELTA_SPACE(2));
-  CCTK_REAL const dt = ToReal(CCTK_DELTA_TIME);
-  CCTK_REAL const t = ToReal(cctk_time);
-  CCTK_REAL const dxi = INV(dx);
-  CCTK_REAL const dyi = INV(dy);
-  CCTK_REAL const dzi = INV(dz);
-  CCTK_REAL const khalf = 0.5;
-  CCTK_REAL const kthird = 1/3.0;
-  CCTK_REAL const ktwothird = 2.0/3.0;
-  CCTK_REAL const kfourthird = 4.0/3.0;
-  CCTK_REAL const keightthird = 8.0/3.0;
-  CCTK_REAL const hdxi = 0.5 * dxi;
-  CCTK_REAL const hdyi = 0.5 * dyi;
-  CCTK_REAL const hdzi = 0.5 * dzi;
+  ptrdiff_t /*const*/ di CCTK_ATTRIBUTE_UNUSED  = 1;
+  ptrdiff_t /*const*/ dj CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,1,0) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ dk CCTK_ATTRIBUTE_UNUSED  = CCTK_GFINDEX3D(cctkGH,0,0,1) - CCTK_GFINDEX3D(cctkGH,0,0,0);
+  ptrdiff_t /*const*/ cdi CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * di;
+  ptrdiff_t /*const*/ cdj CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dj;
+  ptrdiff_t /*const*/ cdk CCTK_ATTRIBUTE_UNUSED  = sizeof(CCTK_REAL) * dk;
+  CCTK_REAL /*const*/ dx CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(0));
+  CCTK_REAL /*const*/ dy CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(1));
+  CCTK_REAL /*const*/ dz CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_SPACE(2));
+  CCTK_REAL /*const*/ dt CCTK_ATTRIBUTE_UNUSED  = ToReal(CCTK_DELTA_TIME);
+  CCTK_REAL /*const*/ t CCTK_ATTRIBUTE_UNUSED  = ToReal(cctk_time);
+  CCTK_REAL /*const*/ dxi CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ dyi CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ dzi CCTK_ATTRIBUTE_UNUSED  = INV(dz);
+  CCTK_REAL /*const*/ khalf CCTK_ATTRIBUTE_UNUSED  = 0.5;
+  CCTK_REAL /*const*/ kthird CCTK_ATTRIBUTE_UNUSED  = 1/3.0;
+  CCTK_REAL /*const*/ ktwothird CCTK_ATTRIBUTE_UNUSED  = 2.0/3.0;
+  CCTK_REAL /*const*/ kfourthird CCTK_ATTRIBUTE_UNUSED  = 4.0/3.0;
+  CCTK_REAL /*const*/ keightthird CCTK_ATTRIBUTE_UNUSED  = 8.0/3.0;
+  CCTK_REAL /*const*/ hdxi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dxi;
+  CCTK_REAL /*const*/ hdyi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dyi;
+  CCTK_REAL /*const*/ hdzi CCTK_ATTRIBUTE_UNUSED  = 0.5 * dzi;
   
   /* Initialize predefined quantities */
-  CCTK_REAL const p1o1 = 1;
-  CCTK_REAL const p1odx = INV(dx);
-  CCTK_REAL const p1ody = INV(dy);
-  CCTK_REAL const p1odz = INV(dz);
+  CCTK_REAL /*const*/ p1o1 CCTK_ATTRIBUTE_UNUSED  = 1.;
+  CCTK_REAL /*const*/ p1odx CCTK_ATTRIBUTE_UNUSED  = INV(dx);
+  CCTK_REAL /*const*/ p1ody CCTK_ATTRIBUTE_UNUSED  = INV(dy);
+  CCTK_REAL /*const*/ p1odz CCTK_ATTRIBUTE_UNUSED  = INV(dz);
   
   /* Assign local copies of arrays functions */
   
@@ -74,11 +74,11 @@ static void eulerauto_initial_shock_Body(cGH const * restrict const cctkGH, int
     i,j,k, imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
     cctk_ash[0],cctk_ash[1],cctk_ash[2])
   {
-    ptrdiff_t const index = di*i + dj*j + dk*k;
+    ptrdiff_t /*const*/ index CCTK_ATTRIBUTE_UNUSED  = di*i + dj*j + dk*k;
     
     /* Assign local copies of grid functions */
     
-    CCTK_REAL xL = x[index];
+    CCTK_REAL xL CCTK_ATTRIBUTE_UNUSED = x[index];
     
     
     /* Include user supplied include files */
@@ -86,18 +86,18 @@ static void eulerauto_initial_shock_Body(cGH const * restrict const cctkGH, int
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL rhoL = ToReal(rhoL0) + StepFunction(-0.5 + 
-      xL)*(-ToReal(rhoL0) + ToReal(rhoR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoL = ToReal(rhoL0) + 
+      StepFunction(-0.5 + xL)*(-1.*ToReal(rhoL0) + ToReal(rhoR0));
     
-    CCTK_REAL v1L = ToReal(vL0) + StepFunction(-0.5 + xL)*(-ToReal(vL0) 
-      + ToReal(vR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1L = ToReal(vL0) + StepFunction(-0.5 
+      + xL)*(-1.*ToReal(vL0) + ToReal(vR0));
     
-    CCTK_REAL v2L = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2L = 0.;
     
-    CCTK_REAL v3L = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3L = 0.;
     
-    CCTK_REAL pL = ToReal(pL0) + StepFunction(-0.5 + xL)*(-ToReal(pL0) + 
-      ToReal(pR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pL = ToReal(pL0) + StepFunction(-0.5 + 
+      xL)*(-1.*ToReal(pL0) + ToReal(pR0));
     
     /* Copy local copies back to grid functions */
     p[index] = pL;