Regenerate examples

22 files changed, 1270 insertions, 1266 deletions

diff --git a/Examples/EulerSR/interface.ccl b/Examples/EulerSR/interface.ccl
index 41c6ab6..65853bd 100644
--- a/Examples/EulerSR/interface.ccl
+++ b/Examples/EulerSR/interface.ccl
@@ -186,13 +186,13 @@ CCTK_REAL tau_group type=GF timelevels=3 tags='tensortypealias="Scalar" tensorwe
 } "tau_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 S_grouprhs type=GF timelevels=3 tags='tensortypealias="D" tensorweight=0'
+CCTK_REAL S_grouprhs type=GF timelevels=3 tags='tensortypealias="D" tensorweight=0 Prolongation="None"'
 {
   S1rhs,
   S2rhs,
@@ -200,7 +200,7 @@ CCTK_REAL S_grouprhs type=GF timelevels=3 tags='tensortypealias="D" tensorweight
 } "S_grouprhs"
 
 public:
-CCTK_REAL tau_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0'
+CCTK_REAL tau_grouprhs type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0 Prolongation="None"'
 {
   taurhs
 } "tau_grouprhs"
diff --git a/Examples/EulerSR/schedule.ccl b/Examples/EulerSR/schedule.ccl
index 87edede..8562b14 100644
--- a/Examples/EulerSR/schedule.ccl
+++ b/Examples/EulerSR/schedule.ccl
@@ -202,93 +202,84 @@ if (CCTK_EQUALS(initial_data, "shock"))
   schedule eulersr_initial_shock at CCTK_INITIAL as eulersr_initial
   {
     LANG: C
-    READS:    grid::x
-    WRITES:   EulerSR::epsi
-    WRITES:   EulerSR::rho
-    WRITES:   EulerSR::v1
-    WRITES:   EulerSR::v2
-    WRITES:   EulerSR::v3
+    READS:    grid::x(Everywhere)
+    WRITES:   EulerSR::epsi(Everywhere)
+    WRITES:   EulerSR::rho(Everywhere)
+    WRITES:   EulerSR::v1(Everywhere)
+    WRITES:   EulerSR::v2(Everywhere)
+    WRITES:   EulerSR::v3(Everywhere)
   } "eulersr_initial_shock"
 }
 
 schedule eulersr_cons_calc_zero_rhs in MoL_CalcRHS
 {
   LANG: C
-  WRITES:   EulerSR::Denrhs
-  WRITES:   EulerSR::S1rhs
-  WRITES:   EulerSR::S2rhs
-  WRITES:   EulerSR::S3rhs
-  WRITES:   EulerSR::taurhs
+  WRITES:   EulerSR::Denrhs(Everywhere)
+  WRITES:   EulerSR::S1rhs(Everywhere)
+  WRITES:   EulerSR::S2rhs(Everywhere)
+  WRITES:   EulerSR::S3rhs(Everywhere)
+  WRITES:   EulerSR::taurhs(Everywhere)
 } "eulersr_cons_calc_zero_rhs"
 
 schedule eulersr_cons_calc_conserved at POSTINITIAL
 {
   LANG: C
-  READS:    EulerSR::Den
-  READS:    EulerSR::epsi
-  READS:    EulerSR::h
-  READS:    EulerSR::p
-  READS:    EulerSR::rho
-  READS:    EulerSR::v1
-  READS:    EulerSR::v2
-  READS:    EulerSR::v3
-  READS:    EulerSR::W
-  WRITES:   EulerSR::Den
-  WRITES:   EulerSR::h
-  WRITES:   EulerSR::p
-  WRITES:   EulerSR::S1
-  WRITES:   EulerSR::S2
-  WRITES:   EulerSR::S3
-  WRITES:   EulerSR::tau
-  WRITES:   EulerSR::W
+  READS:    EulerSR::epsi(Everywhere)
+  READS:    EulerSR::rho(Everywhere)
+  READS:    EulerSR::v1(Everywhere)
+  READS:    EulerSR::v2(Everywhere)
+  READS:    EulerSR::v3(Everywhere)
+  WRITES:   EulerSR::Den(Everywhere)
+  WRITES:   EulerSR::h(Everywhere)
+  WRITES:   EulerSR::p(Everywhere)
+  WRITES:   EulerSR::S1(Everywhere)
+  WRITES:   EulerSR::S2(Everywhere)
+  WRITES:   EulerSR::S3(Everywhere)
+  WRITES:   EulerSR::tau(Everywhere)
+  WRITES:   EulerSR::W(Everywhere)
 } "eulersr_cons_calc_conserved"
 
 schedule eulersr_cons_calc_primitives in MoL_PostStep after EulerSR_ApplyBCs
 {
   LANG: C
-  READS:    EulerSR::Den
-  READS:    EulerSR::epsi
-  READS:    EulerSR::h
-  READS:    EulerSR::p
-  READS:    EulerSR::rho
-  READS:    EulerSR::S1
-  READS:    EulerSR::S2
-  READS:    EulerSR::S3
-  READS:    EulerSR::tau
-  READS:    EulerSR::W
-  WRITES:   EulerSR::epsi
-  WRITES:   EulerSR::h
-  WRITES:   EulerSR::p
-  WRITES:   EulerSR::rho
-  WRITES:   EulerSR::v1
-  WRITES:   EulerSR::v2
-  WRITES:   EulerSR::v3
-  WRITES:   EulerSR::W
+  READS:    EulerSR::Den(Everywhere)
+  READS:    EulerSR::p(Everywhere)
+  READS:    EulerSR::S1(Everywhere)
+  READS:    EulerSR::S2(Everywhere)
+  READS:    EulerSR::S3(Everywhere)
+  READS:    EulerSR::tau(Everywhere)
+  WRITES:   EulerSR::epsi(Everywhere)
+  WRITES:   EulerSR::h(Everywhere)
+  WRITES:   EulerSR::p(Everywhere)
+  WRITES:   EulerSR::rho(Everywhere)
+  WRITES:   EulerSR::v1(Everywhere)
+  WRITES:   EulerSR::v2(Everywhere)
+  WRITES:   EulerSR::v3(Everywhere)
+  WRITES:   EulerSR::W(Everywhere)
 } "eulersr_cons_calc_primitives"
 
 schedule group eulersr_cons_calc_reconstruct_1_group in MoL_CalcRHS after eulersr_cons_calc_zero_rhs
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_1"
 
 schedule eulersr_cons_calc_reconstruct_1 in eulersr_cons_calc_reconstruct_1_group
 {
   LANG: C
-  READS:    EulerSR::epsi
-  READS:    EulerSR::rho
-  READS:    EulerSR::v1
-  READS:    EulerSR::v2
-  READS:    EulerSR::v3
-  WRITES:   EulerSR::epsiLeft
-  WRITES:   EulerSR::epsiRight
-  WRITES:   EulerSR::rhoLeft
-  WRITES:   EulerSR::rhoRight
-  WRITES:   EulerSR::v1Left
-  WRITES:   EulerSR::v1Right
-  WRITES:   EulerSR::v2Left
-  WRITES:   EulerSR::v2Right
-  WRITES:   EulerSR::v3Left
-  WRITES:   EulerSR::v3Right
+  READS:    EulerSR::epsi(Everywhere)
+  READS:    EulerSR::rho(Everywhere)
+  READS:    EulerSR::v1(Everywhere)
+  READS:    EulerSR::v2(Everywhere)
+  READS:    EulerSR::v3(Everywhere)
+  WRITES:   EulerSR::epsiLeft(Interior)
+  WRITES:   EulerSR::epsiRight(Interior)
+  WRITES:   EulerSR::rhoLeft(Interior)
+  WRITES:   EulerSR::rhoRight(Interior)
+  WRITES:   EulerSR::v1Left(Interior)
+  WRITES:   EulerSR::v1Right(Interior)
+  WRITES:   EulerSR::v2Left(Interior)
+  WRITES:   EulerSR::v2Right(Interior)
+  WRITES:   EulerSR::v3Left(Interior)
+  WRITES:   EulerSR::v3Right(Interior)
 } "eulersr_cons_calc_reconstruct_1"
 
 schedule eulersr_cons_calc_reconstruct_1_SelectBCs in eulersr_cons_calc_reconstruct_1_bc_group
@@ -304,81 +295,72 @@ schedule eulersr_cons_calc_reconstruct_1_SelectBCs in eulersr_cons_calc_reconstr
 
 schedule group ApplyBCs as eulersr_cons_calc_reconstruct_1_ApplyBCs in eulersr_cons_calc_reconstruct_1_bc_group after eulersr_cons_calc_reconstruct_1_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_reconstruct_1"
 
 schedule group eulersr_cons_calc_reconstruct_1_bc_group in eulersr_cons_calc_reconstruct_1_group after eulersr_cons_calc_reconstruct_1
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_1"
 
 schedule group eulersr_cons_calc_reconstruct_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_1"
 
 schedule eulersr_cons_calc_intercell_conserved_1 in MoL_CalcRHS after eulersr_cons_calc_reconstruct_1
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::h
-  READS:    EulerSR::p
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::v1Left
-  READS:    EulerSR::v1Right
-  READS:    EulerSR::v2Left
-  READS:    EulerSR::v2Right
-  READS:    EulerSR::v3Left
-  READS:    EulerSR::v3Right
-  READS:    EulerSR::W
-  WRITES:   EulerSR::DenLeft
-  WRITES:   EulerSR::DenRight
-  WRITES:   EulerSR::h
-  WRITES:   EulerSR::p
-  WRITES:   EulerSR::S1Left
-  WRITES:   EulerSR::S1Right
-  WRITES:   EulerSR::S2Left
-  WRITES:   EulerSR::S2Right
-  WRITES:   EulerSR::S3Left
-  WRITES:   EulerSR::S3Right
-  WRITES:   EulerSR::tauLeft
-  WRITES:   EulerSR::tauRight
-  WRITES:   EulerSR::W
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::v1Left(Everywhere)
+  READS:    EulerSR::v1Right(Everywhere)
+  READS:    EulerSR::v2Left(Everywhere)
+  READS:    EulerSR::v2Right(Everywhere)
+  READS:    EulerSR::v3Left(Everywhere)
+  READS:    EulerSR::v3Right(Everywhere)
+  WRITES:   EulerSR::DenLeft(Everywhere)
+  WRITES:   EulerSR::DenRight(Everywhere)
+  WRITES:   EulerSR::h(Everywhere)
+  WRITES:   EulerSR::p(Everywhere)
+  WRITES:   EulerSR::S1Left(Everywhere)
+  WRITES:   EulerSR::S1Right(Everywhere)
+  WRITES:   EulerSR::S2Left(Everywhere)
+  WRITES:   EulerSR::S2Right(Everywhere)
+  WRITES:   EulerSR::S3Left(Everywhere)
+  WRITES:   EulerSR::S3Right(Everywhere)
+  WRITES:   EulerSR::tauLeft(Everywhere)
+  WRITES:   EulerSR::tauRight(Everywhere)
+  WRITES:   EulerSR::W(Everywhere)
 } "eulersr_cons_calc_intercell_conserved_1"
 
 schedule group eulersr_cons_calc_flux_1_group in MoL_CalcRHS after eulersr_cons_calc_intercell_conserved_1
 {
-  # no language specified
 } "eulersr_cons_calc_flux_1"
 
 schedule eulersr_cons_calc_flux_1 in eulersr_cons_calc_flux_1_group
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::S1Left
-  READS:    EulerSR::S1Right
-  READS:    EulerSR::S2Left
-  READS:    EulerSR::S2Right
-  READS:    EulerSR::S3Left
-  READS:    EulerSR::S3Right
-  READS:    EulerSR::tauLeft
-  READS:    EulerSR::tauRight
-  READS:    EulerSR::v1Left
-  READS:    EulerSR::v1Right
-  WRITES:   EulerSR::DenFlux
-  WRITES:   EulerSR::S1Flux
-  WRITES:   EulerSR::S2Flux
-  WRITES:   EulerSR::S3Flux
-  WRITES:   EulerSR::tauFlux
+  READS:    EulerSR::DenLeft(Everywhere)
+  READS:    EulerSR::DenRight(Everywhere)
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::S1Left(Everywhere)
+  READS:    EulerSR::S1Right(Everywhere)
+  READS:    EulerSR::S2Left(Everywhere)
+  READS:    EulerSR::S2Right(Everywhere)
+  READS:    EulerSR::S3Left(Everywhere)
+  READS:    EulerSR::S3Right(Everywhere)
+  READS:    EulerSR::tauLeft(Everywhere)
+  READS:    EulerSR::tauRight(Everywhere)
+  READS:    EulerSR::v1Left(Everywhere)
+  READS:    EulerSR::v1Right(Everywhere)
+  WRITES:   EulerSR::DenFlux(Interior)
+  WRITES:   EulerSR::S1Flux(Interior)
+  WRITES:   EulerSR::S2Flux(Interior)
+  WRITES:   EulerSR::S3Flux(Interior)
+  WRITES:   EulerSR::tauFlux(Interior)
 } "eulersr_cons_calc_flux_1"
 
 schedule eulersr_cons_calc_flux_1_SelectBCs in eulersr_cons_calc_flux_1_bc_group
@@ -394,62 +376,58 @@ schedule eulersr_cons_calc_flux_1_SelectBCs in eulersr_cons_calc_flux_1_bc_group
 
 schedule group ApplyBCs as eulersr_cons_calc_flux_1_ApplyBCs in eulersr_cons_calc_flux_1_bc_group after eulersr_cons_calc_flux_1_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_flux_1"
 
 schedule group eulersr_cons_calc_flux_1_bc_group in eulersr_cons_calc_flux_1_group after eulersr_cons_calc_flux_1
 {
-  # no language specified
 } "eulersr_cons_calc_flux_1"
 
 schedule group eulersr_cons_calc_flux_1_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_flux_1"
 
 schedule eulersr_cons_calc_rhs_1 in MoL_CalcRHS after eulersr_cons_calc_flux_1
 {
   LANG: C
-  READS:    EulerSR::DenFlux
-  READS:    EulerSR::Denrhs
-  READS:    EulerSR::S1Flux
-  READS:    EulerSR::S1rhs
-  READS:    EulerSR::S2Flux
-  READS:    EulerSR::S2rhs
-  READS:    EulerSR::S3Flux
-  READS:    EulerSR::S3rhs
-  READS:    EulerSR::tauFlux
-  READS:    EulerSR::taurhs
-  WRITES:   EulerSR::Denrhs
-  WRITES:   EulerSR::S1rhs
-  WRITES:   EulerSR::S2rhs
-  WRITES:   EulerSR::S3rhs
-  WRITES:   EulerSR::taurhs
+  READS:    EulerSR::DenFlux(Everywhere)
+  READS:    EulerSR::S1Flux(Everywhere)
+  READS:    EulerSR::S2Flux(Everywhere)
+  READS:    EulerSR::S3Flux(Everywhere)
+  READS:    EulerSR::tauFlux(Everywhere)
+  READS:    EulerSR::Denrhs(Everywhere)
+  READS:    EulerSR::S1rhs(Everywhere)
+  READS:    EulerSR::S2rhs(Everywhere)
+  READS:    EulerSR::S3rhs(Everywhere)
+  READS:    EulerSR::taurhs(Everywhere)
+  WRITES:   EulerSR::Denrhs(Interior)
+  WRITES:   EulerSR::S1rhs(Interior)
+  WRITES:   EulerSR::S2rhs(Interior)
+  WRITES:   EulerSR::S3rhs(Interior)
+  WRITES:   EulerSR::taurhs(Interior)
 } "eulersr_cons_calc_rhs_1"
 
 schedule group eulersr_cons_calc_reconstruct_2_group in MoL_CalcRHS after eulersr_cons_calc_rhs_1
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_2"
 
 schedule eulersr_cons_calc_reconstruct_2 in eulersr_cons_calc_reconstruct_2_group
 {
   LANG: C
-  READS:    EulerSR::epsi
-  READS:    EulerSR::rho
-  READS:    EulerSR::v1
-  READS:    EulerSR::v2
-  READS:    EulerSR::v3
-  WRITES:   EulerSR::epsiLeft
-  WRITES:   EulerSR::epsiRight
-  WRITES:   EulerSR::rhoLeft
-  WRITES:   EulerSR::rhoRight
-  WRITES:   EulerSR::v1Left
-  WRITES:   EulerSR::v1Right
-  WRITES:   EulerSR::v2Left
-  WRITES:   EulerSR::v2Right
-  WRITES:   EulerSR::v3Left
-  WRITES:   EulerSR::v3Right
+  READS:    EulerSR::epsi(Everywhere)
+  READS:    EulerSR::rho(Everywhere)
+  READS:    EulerSR::v1(Everywhere)
+  READS:    EulerSR::v2(Everywhere)
+  READS:    EulerSR::v3(Everywhere)
+  WRITES:   EulerSR::epsiLeft(Interior)
+  WRITES:   EulerSR::epsiRight(Interior)
+  WRITES:   EulerSR::rhoLeft(Interior)
+  WRITES:   EulerSR::rhoRight(Interior)
+  WRITES:   EulerSR::v1Left(Interior)
+  WRITES:   EulerSR::v1Right(Interior)
+  WRITES:   EulerSR::v2Left(Interior)
+  WRITES:   EulerSR::v2Right(Interior)
+  WRITES:   EulerSR::v3Left(Interior)
+  WRITES:   EulerSR::v3Right(Interior)
 } "eulersr_cons_calc_reconstruct_2"
 
 schedule eulersr_cons_calc_reconstruct_2_SelectBCs in eulersr_cons_calc_reconstruct_2_bc_group
@@ -465,81 +443,72 @@ schedule eulersr_cons_calc_reconstruct_2_SelectBCs in eulersr_cons_calc_reconstr
 
 schedule group ApplyBCs as eulersr_cons_calc_reconstruct_2_ApplyBCs in eulersr_cons_calc_reconstruct_2_bc_group after eulersr_cons_calc_reconstruct_2_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_reconstruct_2"
 
 schedule group eulersr_cons_calc_reconstruct_2_bc_group in eulersr_cons_calc_reconstruct_2_group after eulersr_cons_calc_reconstruct_2
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_2"
 
 schedule group eulersr_cons_calc_reconstruct_2_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_2"
 
 schedule eulersr_cons_calc_intercell_conserved_2 in MoL_CalcRHS after eulersr_cons_calc_reconstruct_2
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::h
-  READS:    EulerSR::p
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::v1Left
-  READS:    EulerSR::v1Right
-  READS:    EulerSR::v2Left
-  READS:    EulerSR::v2Right
-  READS:    EulerSR::v3Left
-  READS:    EulerSR::v3Right
-  READS:    EulerSR::W
-  WRITES:   EulerSR::DenLeft
-  WRITES:   EulerSR::DenRight
-  WRITES:   EulerSR::h
-  WRITES:   EulerSR::p
-  WRITES:   EulerSR::S1Left
-  WRITES:   EulerSR::S1Right
-  WRITES:   EulerSR::S2Left
-  WRITES:   EulerSR::S2Right
-  WRITES:   EulerSR::S3Left
-  WRITES:   EulerSR::S3Right
-  WRITES:   EulerSR::tauLeft
-  WRITES:   EulerSR::tauRight
-  WRITES:   EulerSR::W
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::v1Left(Everywhere)
+  READS:    EulerSR::v1Right(Everywhere)
+  READS:    EulerSR::v2Left(Everywhere)
+  READS:    EulerSR::v2Right(Everywhere)
+  READS:    EulerSR::v3Left(Everywhere)
+  READS:    EulerSR::v3Right(Everywhere)
+  WRITES:   EulerSR::DenLeft(Everywhere)
+  WRITES:   EulerSR::DenRight(Everywhere)
+  WRITES:   EulerSR::h(Everywhere)
+  WRITES:   EulerSR::p(Everywhere)
+  WRITES:   EulerSR::S1Left(Everywhere)
+  WRITES:   EulerSR::S1Right(Everywhere)
+  WRITES:   EulerSR::S2Left(Everywhere)
+  WRITES:   EulerSR::S2Right(Everywhere)
+  WRITES:   EulerSR::S3Left(Everywhere)
+  WRITES:   EulerSR::S3Right(Everywhere)
+  WRITES:   EulerSR::tauLeft(Everywhere)
+  WRITES:   EulerSR::tauRight(Everywhere)
+  WRITES:   EulerSR::W(Everywhere)
 } "eulersr_cons_calc_intercell_conserved_2"
 
 schedule group eulersr_cons_calc_flux_2_group in MoL_CalcRHS after eulersr_cons_calc_intercell_conserved_2
 {
-  # no language specified
 } "eulersr_cons_calc_flux_2"
 
 schedule eulersr_cons_calc_flux_2 in eulersr_cons_calc_flux_2_group
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::S1Left
-  READS:    EulerSR::S1Right
-  READS:    EulerSR::S2Left
-  READS:    EulerSR::S2Right
-  READS:    EulerSR::S3Left
-  READS:    EulerSR::S3Right
-  READS:    EulerSR::tauLeft
-  READS:    EulerSR::tauRight
-  READS:    EulerSR::v2Left
-  READS:    EulerSR::v2Right
-  WRITES:   EulerSR::DenFlux
-  WRITES:   EulerSR::S1Flux
-  WRITES:   EulerSR::S2Flux
-  WRITES:   EulerSR::S3Flux
-  WRITES:   EulerSR::tauFlux
+  READS:    EulerSR::DenLeft(Everywhere)
+  READS:    EulerSR::DenRight(Everywhere)
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::S1Left(Everywhere)
+  READS:    EulerSR::S1Right(Everywhere)
+  READS:    EulerSR::S2Left(Everywhere)
+  READS:    EulerSR::S2Right(Everywhere)
+  READS:    EulerSR::S3Left(Everywhere)
+  READS:    EulerSR::S3Right(Everywhere)
+  READS:    EulerSR::tauLeft(Everywhere)
+  READS:    EulerSR::tauRight(Everywhere)
+  READS:    EulerSR::v2Left(Everywhere)
+  READS:    EulerSR::v2Right(Everywhere)
+  WRITES:   EulerSR::DenFlux(Interior)
+  WRITES:   EulerSR::S1Flux(Interior)
+  WRITES:   EulerSR::S2Flux(Interior)
+  WRITES:   EulerSR::S3Flux(Interior)
+  WRITES:   EulerSR::tauFlux(Interior)
 } "eulersr_cons_calc_flux_2"
 
 schedule eulersr_cons_calc_flux_2_SelectBCs in eulersr_cons_calc_flux_2_bc_group
@@ -555,62 +524,58 @@ schedule eulersr_cons_calc_flux_2_SelectBCs in eulersr_cons_calc_flux_2_bc_group
 
 schedule group ApplyBCs as eulersr_cons_calc_flux_2_ApplyBCs in eulersr_cons_calc_flux_2_bc_group after eulersr_cons_calc_flux_2_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_flux_2"
 
 schedule group eulersr_cons_calc_flux_2_bc_group in eulersr_cons_calc_flux_2_group after eulersr_cons_calc_flux_2
 {
-  # no language specified
 } "eulersr_cons_calc_flux_2"
 
 schedule group eulersr_cons_calc_flux_2_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_flux_2"
 
 schedule eulersr_cons_calc_rhs_2 in MoL_CalcRHS after eulersr_cons_calc_flux_2
 {
   LANG: C
-  READS:    EulerSR::DenFlux
-  READS:    EulerSR::Denrhs
-  READS:    EulerSR::S1Flux
-  READS:    EulerSR::S1rhs
-  READS:    EulerSR::S2Flux
-  READS:    EulerSR::S2rhs
-  READS:    EulerSR::S3Flux
-  READS:    EulerSR::S3rhs
-  READS:    EulerSR::tauFlux
-  READS:    EulerSR::taurhs
-  WRITES:   EulerSR::Denrhs
-  WRITES:   EulerSR::S1rhs
-  WRITES:   EulerSR::S2rhs
-  WRITES:   EulerSR::S3rhs
-  WRITES:   EulerSR::taurhs
+  READS:    EulerSR::DenFlux(Everywhere)
+  READS:    EulerSR::S1Flux(Everywhere)
+  READS:    EulerSR::S2Flux(Everywhere)
+  READS:    EulerSR::S3Flux(Everywhere)
+  READS:    EulerSR::tauFlux(Everywhere)
+  READS:    EulerSR::Denrhs(Everywhere)
+  READS:    EulerSR::S1rhs(Everywhere)
+  READS:    EulerSR::S2rhs(Everywhere)
+  READS:    EulerSR::S3rhs(Everywhere)
+  READS:    EulerSR::taurhs(Everywhere)
+  WRITES:   EulerSR::Denrhs(Interior)
+  WRITES:   EulerSR::S1rhs(Interior)
+  WRITES:   EulerSR::S2rhs(Interior)
+  WRITES:   EulerSR::S3rhs(Interior)
+  WRITES:   EulerSR::taurhs(Interior)
 } "eulersr_cons_calc_rhs_2"
 
 schedule group eulersr_cons_calc_reconstruct_3_group in MoL_CalcRHS after eulersr_cons_calc_rhs_2
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_3"
 
 schedule eulersr_cons_calc_reconstruct_3 in eulersr_cons_calc_reconstruct_3_group
 {
   LANG: C
-  READS:    EulerSR::epsi
-  READS:    EulerSR::rho
-  READS:    EulerSR::v1
-  READS:    EulerSR::v2
-  READS:    EulerSR::v3
-  WRITES:   EulerSR::epsiLeft
-  WRITES:   EulerSR::epsiRight
-  WRITES:   EulerSR::rhoLeft
-  WRITES:   EulerSR::rhoRight
-  WRITES:   EulerSR::v1Left
-  WRITES:   EulerSR::v1Right
-  WRITES:   EulerSR::v2Left
-  WRITES:   EulerSR::v2Right
-  WRITES:   EulerSR::v3Left
-  WRITES:   EulerSR::v3Right
+  READS:    EulerSR::epsi(Everywhere)
+  READS:    EulerSR::rho(Everywhere)
+  READS:    EulerSR::v1(Everywhere)
+  READS:    EulerSR::v2(Everywhere)
+  READS:    EulerSR::v3(Everywhere)
+  WRITES:   EulerSR::epsiLeft(Interior)
+  WRITES:   EulerSR::epsiRight(Interior)
+  WRITES:   EulerSR::rhoLeft(Interior)
+  WRITES:   EulerSR::rhoRight(Interior)
+  WRITES:   EulerSR::v1Left(Interior)
+  WRITES:   EulerSR::v1Right(Interior)
+  WRITES:   EulerSR::v2Left(Interior)
+  WRITES:   EulerSR::v2Right(Interior)
+  WRITES:   EulerSR::v3Left(Interior)
+  WRITES:   EulerSR::v3Right(Interior)
 } "eulersr_cons_calc_reconstruct_3"
 
 schedule eulersr_cons_calc_reconstruct_3_SelectBCs in eulersr_cons_calc_reconstruct_3_bc_group
@@ -626,81 +591,72 @@ schedule eulersr_cons_calc_reconstruct_3_SelectBCs in eulersr_cons_calc_reconstr
 
 schedule group ApplyBCs as eulersr_cons_calc_reconstruct_3_ApplyBCs in eulersr_cons_calc_reconstruct_3_bc_group after eulersr_cons_calc_reconstruct_3_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_reconstruct_3"
 
 schedule group eulersr_cons_calc_reconstruct_3_bc_group in eulersr_cons_calc_reconstruct_3_group after eulersr_cons_calc_reconstruct_3
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_3"
 
 schedule group eulersr_cons_calc_reconstruct_3_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_reconstruct_3"
 
 schedule eulersr_cons_calc_intercell_conserved_3 in MoL_CalcRHS after eulersr_cons_calc_reconstruct_3
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::h
-  READS:    EulerSR::p
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::v1Left
-  READS:    EulerSR::v1Right
-  READS:    EulerSR::v2Left
-  READS:    EulerSR::v2Right
-  READS:    EulerSR::v3Left
-  READS:    EulerSR::v3Right
-  READS:    EulerSR::W
-  WRITES:   EulerSR::DenLeft
-  WRITES:   EulerSR::DenRight
-  WRITES:   EulerSR::h
-  WRITES:   EulerSR::p
-  WRITES:   EulerSR::S1Left
-  WRITES:   EulerSR::S1Right
-  WRITES:   EulerSR::S2Left
-  WRITES:   EulerSR::S2Right
-  WRITES:   EulerSR::S3Left
-  WRITES:   EulerSR::S3Right
-  WRITES:   EulerSR::tauLeft
-  WRITES:   EulerSR::tauRight
-  WRITES:   EulerSR::W
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::v1Left(Everywhere)
+  READS:    EulerSR::v1Right(Everywhere)
+  READS:    EulerSR::v2Left(Everywhere)
+  READS:    EulerSR::v2Right(Everywhere)
+  READS:    EulerSR::v3Left(Everywhere)
+  READS:    EulerSR::v3Right(Everywhere)
+  WRITES:   EulerSR::DenLeft(Everywhere)
+  WRITES:   EulerSR::DenRight(Everywhere)
+  WRITES:   EulerSR::h(Everywhere)
+  WRITES:   EulerSR::p(Everywhere)
+  WRITES:   EulerSR::S1Left(Everywhere)
+  WRITES:   EulerSR::S1Right(Everywhere)
+  WRITES:   EulerSR::S2Left(Everywhere)
+  WRITES:   EulerSR::S2Right(Everywhere)
+  WRITES:   EulerSR::S3Left(Everywhere)
+  WRITES:   EulerSR::S3Right(Everywhere)
+  WRITES:   EulerSR::tauLeft(Everywhere)
+  WRITES:   EulerSR::tauRight(Everywhere)
+  WRITES:   EulerSR::W(Everywhere)
 } "eulersr_cons_calc_intercell_conserved_3"
 
 schedule group eulersr_cons_calc_flux_3_group in MoL_CalcRHS after eulersr_cons_calc_intercell_conserved_3
 {
-  # no language specified
 } "eulersr_cons_calc_flux_3"
 
 schedule eulersr_cons_calc_flux_3 in eulersr_cons_calc_flux_3_group
 {
   LANG: C
-  READS:    EulerSR::DenLeft
-  READS:    EulerSR::DenRight
-  READS:    EulerSR::epsiLeft
-  READS:    EulerSR::epsiRight
-  READS:    EulerSR::rhoLeft
-  READS:    EulerSR::rhoRight
-  READS:    EulerSR::S1Left
-  READS:    EulerSR::S1Right
-  READS:    EulerSR::S2Left
-  READS:    EulerSR::S2Right
-  READS:    EulerSR::S3Left
-  READS:    EulerSR::S3Right
-  READS:    EulerSR::tauLeft
-  READS:    EulerSR::tauRight
-  READS:    EulerSR::v3Left
-  READS:    EulerSR::v3Right
-  WRITES:   EulerSR::DenFlux
-  WRITES:   EulerSR::S1Flux
-  WRITES:   EulerSR::S2Flux
-  WRITES:   EulerSR::S3Flux
-  WRITES:   EulerSR::tauFlux
+  READS:    EulerSR::DenLeft(Everywhere)
+  READS:    EulerSR::DenRight(Everywhere)
+  READS:    EulerSR::epsiLeft(Everywhere)
+  READS:    EulerSR::epsiRight(Everywhere)
+  READS:    EulerSR::rhoLeft(Everywhere)
+  READS:    EulerSR::rhoRight(Everywhere)
+  READS:    EulerSR::S1Left(Everywhere)
+  READS:    EulerSR::S1Right(Everywhere)
+  READS:    EulerSR::S2Left(Everywhere)
+  READS:    EulerSR::S2Right(Everywhere)
+  READS:    EulerSR::S3Left(Everywhere)
+  READS:    EulerSR::S3Right(Everywhere)
+  READS:    EulerSR::tauLeft(Everywhere)
+  READS:    EulerSR::tauRight(Everywhere)
+  READS:    EulerSR::v3Left(Everywhere)
+  READS:    EulerSR::v3Right(Everywhere)
+  WRITES:   EulerSR::DenFlux(Interior)
+  WRITES:   EulerSR::S1Flux(Interior)
+  WRITES:   EulerSR::S2Flux(Interior)
+  WRITES:   EulerSR::S3Flux(Interior)
+  WRITES:   EulerSR::tauFlux(Interior)
 } "eulersr_cons_calc_flux_3"
 
 schedule eulersr_cons_calc_flux_3_SelectBCs in eulersr_cons_calc_flux_3_bc_group
@@ -716,37 +672,34 @@ schedule eulersr_cons_calc_flux_3_SelectBCs in eulersr_cons_calc_flux_3_bc_group
 
 schedule group ApplyBCs as eulersr_cons_calc_flux_3_ApplyBCs in eulersr_cons_calc_flux_3_bc_group after eulersr_cons_calc_flux_3_SelectBCs
 {
-  # no language specified
 } "Apply BCs for groups set in eulersr_cons_calc_flux_3"
 
 schedule group eulersr_cons_calc_flux_3_bc_group in eulersr_cons_calc_flux_3_group after eulersr_cons_calc_flux_3
 {
-  # no language specified
 } "eulersr_cons_calc_flux_3"
 
 schedule group eulersr_cons_calc_flux_3_bc_group in MoL_PseudoEvolutionBoundaries after MoL_PostStep
 {
-  # no language specified
 } "eulersr_cons_calc_flux_3"
 
 schedule eulersr_cons_calc_rhs_3 in MoL_CalcRHS after eulersr_cons_calc_flux_3
 {
   LANG: C
-  READS:    EulerSR::DenFlux
-  READS:    EulerSR::Denrhs
-  READS:    EulerSR::S1Flux
-  READS:    EulerSR::S1rhs
-  READS:    EulerSR::S2Flux
-  READS:    EulerSR::S2rhs
-  READS:    EulerSR::S3Flux
-  READS:    EulerSR::S3rhs
-  READS:    EulerSR::tauFlux
-  READS:    EulerSR::taurhs
-  WRITES:   EulerSR::Denrhs
-  WRITES:   EulerSR::S1rhs
-  WRITES:   EulerSR::S2rhs
-  WRITES:   EulerSR::S3rhs
-  WRITES:   EulerSR::taurhs
+  READS:    EulerSR::DenFlux(Everywhere)
+  READS:    EulerSR::S1Flux(Everywhere)
+  READS:    EulerSR::S2Flux(Everywhere)
+  READS:    EulerSR::S3Flux(Everywhere)
+  READS:    EulerSR::tauFlux(Everywhere)
+  READS:    EulerSR::Denrhs(Everywhere)
+  READS:    EulerSR::S1rhs(Everywhere)
+  READS:    EulerSR::S2rhs(Everywhere)
+  READS:    EulerSR::S3rhs(Everywhere)
+  READS:    EulerSR::taurhs(Everywhere)
+  WRITES:   EulerSR::Denrhs(Interior)
+  WRITES:   EulerSR::S1rhs(Interior)
+  WRITES:   EulerSR::S2rhs(Interior)
+  WRITES:   EulerSR::S3rhs(Interior)
+  WRITES:   EulerSR::taurhs(Interior)
 } "eulersr_cons_calc_rhs_3"
 
 schedule EulerSR_SelectBoundConds in MoL_PostStep
@@ -772,5 +725,4 @@ schedule EulerSR_RegisterVars in MoL_Register
 
 schedule group ApplyBCs as EulerSR_ApplyBCs in MoL_PostStep after EulerSR_SelectBoundConds
 {
-  # no language specified
 } "Apply boundary conditions controlled by thorn Boundary"
diff --git a/Examples/EulerSR/src/Boundaries.cc b/Examples/EulerSR/src/Boundaries.cc
index a862753..5d051c3 100644
--- a/Examples/EulerSR/src/Boundaries.cc
+++ b/Examples/EulerSR/src/Boundaries.cc
@@ -30,7 +30,7 @@ extern "C" void EulerSR_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 EulerSR_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 EulerSR_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)
@@ -161,7 +161,7 @@ extern "C" void EulerSR_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(tau_group_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_tau_group_bound = -1;
+    static CCTK_INT handle_tau_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_tau_group_bound < 0) handle_tau_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_tau_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_tau_group_bound , tau_group_bound_limit, "LIMIT") < 0)
@@ -180,7 +180,7 @@ extern "C" void EulerSR_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 EulerSR_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)
@@ -218,7 +218,7 @@ extern "C" void EulerSR_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)
@@ -237,7 +237,7 @@ extern "C" void EulerSR_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)
@@ -256,7 +256,7 @@ extern "C" void EulerSR_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(tau_bound, "radiative"))
   {
    /* select radiation boundary condition */
-    static CCTK_INT handle_tau_bound = -1;
+    static CCTK_INT handle_tau_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_tau_bound < 0) handle_tau_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_tau_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_tau_bound , tau_bound_limit, "LIMIT") < 0)
@@ -275,7 +275,7 @@ extern "C" void EulerSR_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 EulerSR_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)
@@ -309,7 +309,7 @@ extern "C" void EulerSR_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(tau_group_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_tau_group_bound = -1;
+    static CCTK_INT handle_tau_group_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_tau_group_bound < 0) handle_tau_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_tau_group_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_tau_group_bound ,tau_group_bound_scalar, "SCALAR") < 0)
@@ -326,7 +326,7 @@ extern "C" void EulerSR_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 EulerSR_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)
@@ -360,7 +360,7 @@ extern "C" void EulerSR_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)
@@ -377,7 +377,7 @@ extern "C" void EulerSR_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)
@@ -394,7 +394,7 @@ extern "C" void EulerSR_SelectBoundConds(CCTK_ARGUMENTS)
   if (CCTK_EQUALS(tau_bound, "scalar"))
   {
    /* select scalar boundary condition */
-    static CCTK_INT handle_tau_bound = -1;
+    static CCTK_INT handle_tau_bound CCTK_ATTRIBUTE_UNUSED  = -1;
     if (handle_tau_bound < 0) handle_tau_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
     if (handle_tau_bound < 0) CCTK_WARN(0, "could not create table!");
     if (Util_TableSetReal(handle_tau_bound ,tau_bound_scalar, "SCALAR") < 0)
diff --git a/Examples/EulerSR/src/Differencing.h b/Examples/EulerSR/src/Differencing.h
index 913cfa8..323c0fe 100644
--- a/Examples/EulerSR/src/Differencing.h
+++ b/Examples/EulerSR/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/EulerSR/src/RegisterMoL.cc b/Examples/EulerSR/src/RegisterMoL.cc
index 9ef9a4d..8c1e4a1 100644
--- a/Examples/EulerSR/src/RegisterMoL.cc
+++ b/Examples/EulerSR/src/RegisterMoL.cc
@@ -9,7 +9,7 @@ extern "C" void EulerSR_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("EulerSR::Den"),  CCTK_VarIndex("EulerSR::Denrhs"));
diff --git a/Examples/EulerSR/src/Startup.cc b/Examples/EulerSR/src/Startup.cc
index dc65910..b0abe67 100644
--- a/Examples/EulerSR/src/Startup.cc
+++ b/Examples/EulerSR/src/Startup.cc
@@ -4,7 +4,7 @@
 
 extern "C" int EulerSR_Startup(void)
 {
-  const char * banner = "EulerSR";
+  const char * banner CCTK_ATTRIBUTE_UNUSED  = "EulerSR";
   CCTK_RegisterBanner(banner);
   return 0;
 }
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_conserved.cc b/Examples/EulerSR/src/eulersr_cons_calc_conserved.cc
index ba11a6e..1d21aa4 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_conserved.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_conserved.cc
@@ -31,34 +31,34 @@ static void eulersr_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,19 +74,19 @@ static void eulersr_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 DenL = Den[index];
-    CCTK_REAL epsiL = epsi[index];
-    CCTK_REAL hL = h[index];
-    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 WL = W[index];
+    CCTK_REAL DenL CCTK_ATTRIBUTE_UNUSED = Den[index];
+    CCTK_REAL epsiL CCTK_ATTRIBUTE_UNUSED = epsi[index];
+    CCTK_REAL hL CCTK_ATTRIBUTE_UNUSED = h[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];
+    CCTK_REAL WL CCTK_ATTRIBUTE_UNUSED = W[index];
     
     
     /* Include user supplied include files */
@@ -94,23 +94,25 @@ static void eulersr_cons_calc_conserved_Body(cGH const * restrict const cctkGH,
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL Wx = 1 - SQR(v1L) - SQR(v2L) - SQR(v3L);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Wx = 1. - 1.*(SQR(v1L) + SQR(v2L) 
+      + SQR(v3L));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiL*rhoL*(-1 + ToReal(gamma));
+    pL = epsiL*rhoL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiL + pL*INV(rhoL);
+    hL = 1. + epsiL + pL*INV(rhoL);
     
     DenL = rhoL*WL;
     
-    CCTK_REAL S1L = hL*rhoL*v1L*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1L = hL*rhoL*v1L*SQR(WL);
     
-    CCTK_REAL S2L = hL*rhoL*v2L*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2L = hL*rhoL*v2L*SQR(WL);
     
-    CCTK_REAL S3L = hL*rhoL*v3L*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3L = hL*rhoL*v3L*SQR(WL);
     
-    CCTK_REAL tauL = -DenL - pL + hL*rhoL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauL = -1.*(DenL + pL) + 
+      hL*rhoL*SQR(WL);
     
     /* Copy local copies back to grid functions */
     Den[index] = DenL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_flux_1.cc b/Examples/EulerSR/src/eulersr_cons_calc_flux_1.cc
index a618b63..4a3ed77 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_flux_1.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_flux_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_flux_1_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 */
   
@@ -98,78 +98,85 @@ static void eulersr_cons_calc_flux_1_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 DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[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 tauLeftL = tauLeft[index];
-    CCTK_REAL tauRightL = tauRight[index];
-    CCTK_REAL v1LeftL = v1Left[index];
-    CCTK_REAL v1RightL = v1Right[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[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 tauLeftL CCTK_ATTRIBUTE_UNUSED = tauLeft[index];
+    CCTK_REAL tauRightL CCTK_ATTRIBUTE_UNUSED = tauRight[index];
+    CCTK_REAL v1LeftL CCTK_ATTRIBUTE_UNUSED = v1Left[index];
+    CCTK_REAL v1RightL CCTK_ATTRIBUTE_UNUSED = v1Right[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const ShiftMinus1DenRight = ShiftMinus1(&DenRight[index]);
-    CCTK_REAL const ShiftMinus1epsiRight = ShiftMinus1(&epsiRight[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 ShiftMinus1tauRight = ShiftMinus1(&tauRight[index]);
-    CCTK_REAL const ShiftMinus1v1Right = ShiftMinus1(&v1Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus1DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1epsiRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&epsiRight[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*/ ShiftMinus1tauRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&tauRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus1v1Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus1(&v1Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = DenLeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = DenLeftL*v1LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus1DenRight*ShiftMinus1v1Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus1DenRight*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 = S1LeftL*v1LeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = S1LeftL*v1LeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    CCTK_REAL S1FluxRight = ShiftMinus1S1Right*ShiftMinus1v1Right + 
-      ShiftMinus1epsiRight*ShiftMinus1rhoRight*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = 
+      ShiftMinus1S1Right*ShiftMinus1v1Right + 
+      ShiftMinus1epsiRight*ShiftMinus1rhoRight*(-1. + ToReal(gamma));
     
-    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 = S2LeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = S2LeftL*v1LeftL;
     
-    CCTK_REAL S2FluxRight = ShiftMinus1S2Right*ShiftMinus1v1Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = 
+      ShiftMinus1S2Right*ShiftMinus1v1Right;
     
-    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 = S3LeftL*v1LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = S3LeftL*v1LeftL;
     
-    CCTK_REAL S3FluxRight = ShiftMinus1S3Right*ShiftMinus1v1Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = 
+      ShiftMinus1S3Right*ShiftMinus1v1Right;
     
-    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 tauFluxLeft = v1LeftL*(tauLeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxLeft = v1LeftL*(tauLeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxRight = ShiftMinus1v1Right*(ShiftMinus1tauRight + 
-      ShiftMinus1epsiRight*ShiftMinus1rhoRight*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxRight = 
+      ShiftMinus1v1Right*(ShiftMinus1tauRight + 
+      ShiftMinus1epsiRight*ShiftMinus1rhoRight*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxL = 0.5*(tauFluxLeft + tauFluxRight + (-tauLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxL = 0.5*(tauFluxLeft + 
+      tauFluxRight + (-1.*tauLeftL + 
       ShiftMinus1tauRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_flux_2.cc b/Examples/EulerSR/src/eulersr_cons_calc_flux_2.cc
index 9c264a5..9365d09 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_flux_2.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_flux_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_flux_2_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 */
   
@@ -98,78 +98,85 @@ static void eulersr_cons_calc_flux_2_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 DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[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 tauLeftL = tauLeft[index];
-    CCTK_REAL tauRightL = tauRight[index];
-    CCTK_REAL v2LeftL = v2Left[index];
-    CCTK_REAL v2RightL = v2Right[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[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 tauLeftL CCTK_ATTRIBUTE_UNUSED = tauLeft[index];
+    CCTK_REAL tauRightL CCTK_ATTRIBUTE_UNUSED = tauRight[index];
+    CCTK_REAL v2LeftL CCTK_ATTRIBUTE_UNUSED = v2Left[index];
+    CCTK_REAL v2RightL CCTK_ATTRIBUTE_UNUSED = v2Right[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const ShiftMinus2DenRight = ShiftMinus2(&DenRight[index]);
-    CCTK_REAL const ShiftMinus2epsiRight = ShiftMinus2(&epsiRight[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 ShiftMinus2tauRight = ShiftMinus2(&tauRight[index]);
-    CCTK_REAL const ShiftMinus2v2Right = ShiftMinus2(&v2Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus2DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2epsiRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&epsiRight[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*/ ShiftMinus2tauRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&tauRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus2v2Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus2(&v2Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = DenLeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = DenLeftL*v2LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus2DenRight*ShiftMinus2v2Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus2DenRight*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 = S1LeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = S1LeftL*v2LeftL;
     
-    CCTK_REAL S1FluxRight = ShiftMinus2S1Right*ShiftMinus2v2Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = 
+      ShiftMinus2S1Right*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 = S2LeftL*v2LeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = S2LeftL*v2LeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    CCTK_REAL S2FluxRight = ShiftMinus2S2Right*ShiftMinus2v2Right + 
-      ShiftMinus2epsiRight*ShiftMinus2rhoRight*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = 
+      ShiftMinus2S2Right*ShiftMinus2v2Right + 
+      ShiftMinus2epsiRight*ShiftMinus2rhoRight*(-1. + ToReal(gamma));
     
-    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 = S3LeftL*v2LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = S3LeftL*v2LeftL;
     
-    CCTK_REAL S3FluxRight = ShiftMinus2S3Right*ShiftMinus2v2Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = 
+      ShiftMinus2S3Right*ShiftMinus2v2Right;
     
-    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 tauFluxLeft = v2LeftL*(tauLeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxLeft = v2LeftL*(tauLeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxRight = ShiftMinus2v2Right*(ShiftMinus2tauRight + 
-      ShiftMinus2epsiRight*ShiftMinus2rhoRight*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxRight = 
+      ShiftMinus2v2Right*(ShiftMinus2tauRight + 
+      ShiftMinus2epsiRight*ShiftMinus2rhoRight*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxL = 0.5*(tauFluxLeft + tauFluxRight + (-tauLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxL = 0.5*(tauFluxLeft + 
+      tauFluxRight + (-1.*tauLeftL + 
       ShiftMinus2tauRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_flux_3.cc b/Examples/EulerSR/src/eulersr_cons_calc_flux_3.cc
index ede78c8..01c4fe2 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_flux_3.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_flux_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_flux_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_flux_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_flux_3_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 */
   
@@ -98,78 +98,85 @@ static void eulersr_cons_calc_flux_3_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 DenLeftL = DenLeft[index];
-    CCTK_REAL DenRightL = DenRight[index];
-    CCTK_REAL epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[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 tauLeftL = tauLeft[index];
-    CCTK_REAL tauRightL = tauRight[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 epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[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 tauLeftL CCTK_ATTRIBUTE_UNUSED = tauLeft[index];
+    CCTK_REAL tauRightL CCTK_ATTRIBUTE_UNUSED = tauRight[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 ShiftMinus3epsiRight = ShiftMinus3(&epsiRight[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 ShiftMinus3tauRight = ShiftMinus3(&tauRight[index]);
-    CCTK_REAL const ShiftMinus3v3Right = ShiftMinus3(&v3Right[index]);
+    CCTK_REAL /*const*/ ShiftMinus3DenRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&DenRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3epsiRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&epsiRight[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*/ ShiftMinus3tauRight CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&tauRight[index]);
+    CCTK_REAL /*const*/ ShiftMinus3v3Right CCTK_ATTRIBUTE_UNUSED  = ShiftMinus3(&v3Right[index]);
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL DenFluxLeft = DenLeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxLeft = DenLeftL*v3LeftL;
     
-    CCTK_REAL DenFluxRight = ShiftMinus3DenRight*ShiftMinus3v3Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED DenFluxRight = 
+      ShiftMinus3DenRight*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 = S1LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxLeft = S1LeftL*v3LeftL;
     
-    CCTK_REAL S1FluxRight = ShiftMinus3S1Right*ShiftMinus3v3Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1FluxRight = 
+      ShiftMinus3S1Right*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 = S2LeftL*v3LeftL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxLeft = S2LeftL*v3LeftL;
     
-    CCTK_REAL S2FluxRight = ShiftMinus3S2Right*ShiftMinus3v3Right;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2FluxRight = 
+      ShiftMinus3S2Right*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 = S3LeftL*v3LeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxLeft = S3LeftL*v3LeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    CCTK_REAL S3FluxRight = ShiftMinus3S3Right*ShiftMinus3v3Right + 
-      ShiftMinus3epsiRight*ShiftMinus3rhoRight*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3FluxRight = 
+      ShiftMinus3S3Right*ShiftMinus3v3Right + 
+      ShiftMinus3epsiRight*ShiftMinus3rhoRight*(-1. + ToReal(gamma));
     
-    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 tauFluxLeft = v3LeftL*(tauLeftL + 
-      epsiLeftL*rhoLeftL*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxLeft = v3LeftL*(tauLeftL + 
+      epsiLeftL*rhoLeftL*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxRight = ShiftMinus3v3Right*(ShiftMinus3tauRight + 
-      ShiftMinus3epsiRight*ShiftMinus3rhoRight*(-1 + ToReal(gamma)));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxRight = 
+      ShiftMinus3v3Right*(ShiftMinus3tauRight + 
+      ShiftMinus3epsiRight*ShiftMinus3rhoRight*(-1. + ToReal(gamma)));
     
-    CCTK_REAL tauFluxL = 0.5*(tauFluxLeft + tauFluxRight + (-tauLeftL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauFluxL = 0.5*(tauFluxLeft + 
+      tauFluxRight + (-1.*tauLeftL + 
       ShiftMinus3tauRight)*ToReal(hlleAlpha));
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_1.cc b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_1.cc
index 773d8f5..4480d3d 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_1.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_1.cc
@@ -31,34 +31,34 @@ static void eulersr_cons_calc_intercell_conserved_1_Body(cGH const * restrict co
   /* 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,25 +74,25 @@ static void eulersr_cons_calc_intercell_conserved_1_Body(cGH const * restrict co
     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 epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[index];
-    CCTK_REAL hL = h[index];
-    CCTK_REAL pL = p[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 WL = W[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[index];
+    CCTK_REAL hL CCTK_ATTRIBUTE_UNUSED = h[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[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];
+    CCTK_REAL WL CCTK_ATTRIBUTE_UNUSED = W[index];
     
     
     /* Include user supplied include files */
@@ -100,41 +100,49 @@ static void eulersr_cons_calc_intercell_conserved_1_Body(cGH const * restrict co
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL Wx = 1 - SQR(v1LeftL) - SQR(v2LeftL) - SQR(v3LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Wx = 1. - 1.*(SQR(v1LeftL) + 
+      SQR(v2LeftL) + SQR(v3LeftL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    pL = epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiLeftL + pL*INV(rhoLeftL);
+    hL = 1. + epsiLeftL + pL*INV(rhoLeftL);
     
     DenLeftL = rhoLeftL*WL;
     
-    CCTK_REAL S1LeftL = hL*rhoLeftL*v1LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = 
+      hL*rhoLeftL*v1LeftL*SQR(WL);
     
-    CCTK_REAL S2LeftL = hL*rhoLeftL*v2LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = 
+      hL*rhoLeftL*v2LeftL*SQR(WL);
     
-    CCTK_REAL S3LeftL = hL*rhoLeftL*v3LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = 
+      hL*rhoLeftL*v3LeftL*SQR(WL);
     
-    CCTK_REAL tauLeftL = -DenLeftL - pL + hL*rhoLeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauLeftL = -1.*(DenLeftL + pL) + 
+      hL*rhoLeftL*SQR(WL);
     
-    Wx = 1 - SQR(v1RightL) - SQR(v2RightL) - SQR(v3RightL);
+    Wx = 1. - 1.*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiRightL*rhoRightL*(-1 + ToReal(gamma));
+    pL = epsiRightL*rhoRightL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiRightL + pL*INV(rhoRightL);
+    hL = 1. + epsiRightL + pL*INV(rhoRightL);
     
     DenRightL = rhoRightL*WL;
     
-    CCTK_REAL S1RightL = hL*rhoRightL*v1RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = 
+      hL*rhoRightL*v1RightL*SQR(WL);
     
-    CCTK_REAL S2RightL = hL*rhoRightL*v2RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = 
+      hL*rhoRightL*v2RightL*SQR(WL);
     
-    CCTK_REAL S3RightL = hL*rhoRightL*v3RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = 
+      hL*rhoRightL*v3RightL*SQR(WL);
     
-    CCTK_REAL tauRightL = -DenRightL - pL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauRightL = -1.*(DenRightL + pL) + 
       hL*rhoRightL*SQR(WL);
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_2.cc b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_2.cc
index a40fee7..5799b05 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_2.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_2.cc
@@ -31,34 +31,34 @@ static void eulersr_cons_calc_intercell_conserved_2_Body(cGH const * restrict co
   /* 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,25 +74,25 @@ static void eulersr_cons_calc_intercell_conserved_2_Body(cGH const * restrict co
     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 epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[index];
-    CCTK_REAL hL = h[index];
-    CCTK_REAL pL = p[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 WL = W[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[index];
+    CCTK_REAL hL CCTK_ATTRIBUTE_UNUSED = h[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[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];
+    CCTK_REAL WL CCTK_ATTRIBUTE_UNUSED = W[index];
     
     
     /* Include user supplied include files */
@@ -100,41 +100,49 @@ static void eulersr_cons_calc_intercell_conserved_2_Body(cGH const * restrict co
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL Wx = 1 - SQR(v1LeftL) - SQR(v2LeftL) - SQR(v3LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Wx = 1. - 1.*(SQR(v1LeftL) + 
+      SQR(v2LeftL) + SQR(v3LeftL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    pL = epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiLeftL + pL*INV(rhoLeftL);
+    hL = 1. + epsiLeftL + pL*INV(rhoLeftL);
     
     DenLeftL = rhoLeftL*WL;
     
-    CCTK_REAL S1LeftL = hL*rhoLeftL*v1LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = 
+      hL*rhoLeftL*v1LeftL*SQR(WL);
     
-    CCTK_REAL S2LeftL = hL*rhoLeftL*v2LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = 
+      hL*rhoLeftL*v2LeftL*SQR(WL);
     
-    CCTK_REAL S3LeftL = hL*rhoLeftL*v3LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = 
+      hL*rhoLeftL*v3LeftL*SQR(WL);
     
-    CCTK_REAL tauLeftL = -DenLeftL - pL + hL*rhoLeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauLeftL = -1.*(DenLeftL + pL) + 
+      hL*rhoLeftL*SQR(WL);
     
-    Wx = 1 - SQR(v1RightL) - SQR(v2RightL) - SQR(v3RightL);
+    Wx = 1. - 1.*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiRightL*rhoRightL*(-1 + ToReal(gamma));
+    pL = epsiRightL*rhoRightL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiRightL + pL*INV(rhoRightL);
+    hL = 1. + epsiRightL + pL*INV(rhoRightL);
     
     DenRightL = rhoRightL*WL;
     
-    CCTK_REAL S1RightL = hL*rhoRightL*v1RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = 
+      hL*rhoRightL*v1RightL*SQR(WL);
     
-    CCTK_REAL S2RightL = hL*rhoRightL*v2RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = 
+      hL*rhoRightL*v2RightL*SQR(WL);
     
-    CCTK_REAL S3RightL = hL*rhoRightL*v3RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = 
+      hL*rhoRightL*v3RightL*SQR(WL);
     
-    CCTK_REAL tauRightL = -DenRightL - pL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauRightL = -1.*(DenRightL + pL) + 
       hL*rhoRightL*SQR(WL);
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_3.cc b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_3.cc
index 1ac3fae..7925290 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_3.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_intercell_conserved_3.cc
@@ -31,34 +31,34 @@ static void eulersr_cons_calc_intercell_conserved_3_Body(cGH const * restrict co
   /* 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,25 +74,25 @@ static void eulersr_cons_calc_intercell_conserved_3_Body(cGH const * restrict co
     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 epsiLeftL = epsiLeft[index];
-    CCTK_REAL epsiRightL = epsiRight[index];
-    CCTK_REAL hL = h[index];
-    CCTK_REAL pL = p[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 WL = W[index];
+    CCTK_REAL DenLeftL CCTK_ATTRIBUTE_UNUSED = DenLeft[index];
+    CCTK_REAL DenRightL CCTK_ATTRIBUTE_UNUSED = DenRight[index];
+    CCTK_REAL epsiLeftL CCTK_ATTRIBUTE_UNUSED = epsiLeft[index];
+    CCTK_REAL epsiRightL CCTK_ATTRIBUTE_UNUSED = epsiRight[index];
+    CCTK_REAL hL CCTK_ATTRIBUTE_UNUSED = h[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[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];
+    CCTK_REAL WL CCTK_ATTRIBUTE_UNUSED = W[index];
     
     
     /* Include user supplied include files */
@@ -100,41 +100,49 @@ static void eulersr_cons_calc_intercell_conserved_3_Body(cGH const * restrict co
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL Wx = 1 - SQR(v1LeftL) - SQR(v2LeftL) - SQR(v3LeftL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Wx = 1. - 1.*(SQR(v1LeftL) + 
+      SQR(v2LeftL) + SQR(v3LeftL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiLeftL*rhoLeftL*(-1 + ToReal(gamma));
+    pL = epsiLeftL*rhoLeftL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiLeftL + pL*INV(rhoLeftL);
+    hL = 1. + epsiLeftL + pL*INV(rhoLeftL);
     
     DenLeftL = rhoLeftL*WL;
     
-    CCTK_REAL S1LeftL = hL*rhoLeftL*v1LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1LeftL = 
+      hL*rhoLeftL*v1LeftL*SQR(WL);
     
-    CCTK_REAL S2LeftL = hL*rhoLeftL*v2LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2LeftL = 
+      hL*rhoLeftL*v2LeftL*SQR(WL);
     
-    CCTK_REAL S3LeftL = hL*rhoLeftL*v3LeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3LeftL = 
+      hL*rhoLeftL*v3LeftL*SQR(WL);
     
-    CCTK_REAL tauLeftL = -DenLeftL - pL + hL*rhoLeftL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauLeftL = -1.*(DenLeftL + pL) + 
+      hL*rhoLeftL*SQR(WL);
     
-    Wx = 1 - SQR(v1RightL) - SQR(v2RightL) - SQR(v3RightL);
+    Wx = 1. - 1.*(SQR(v1RightL) + SQR(v2RightL) + SQR(v3RightL));
     
     WL = INV(sqrt(Wx));
     
-    pL = epsiRightL*rhoRightL*(-1 + ToReal(gamma));
+    pL = epsiRightL*rhoRightL*(-1. + ToReal(gamma));
     
-    hL = 1 + epsiRightL + pL*INV(rhoRightL);
+    hL = 1. + epsiRightL + pL*INV(rhoRightL);
     
     DenRightL = rhoRightL*WL;
     
-    CCTK_REAL S1RightL = hL*rhoRightL*v1RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S1RightL = 
+      hL*rhoRightL*v1RightL*SQR(WL);
     
-    CCTK_REAL S2RightL = hL*rhoRightL*v2RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2RightL = 
+      hL*rhoRightL*v2RightL*SQR(WL);
     
-    CCTK_REAL S3RightL = hL*rhoRightL*v3RightL*SQR(WL);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3RightL = 
+      hL*rhoRightL*v3RightL*SQR(WL);
     
-    CCTK_REAL tauRightL = -DenRightL - pL + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED tauRightL = -1.*(DenRightL + pL) + 
       hL*rhoRightL*SQR(WL);
     
     /* Copy local copies back to grid functions */
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_primitives.cc b/Examples/EulerSR/src/eulersr_cons_calc_primitives.cc
index 8f58a2a..fcd1353 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_primitives.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_primitives.cc
@@ -31,34 +31,34 @@ static void eulersr_cons_calc_primitives_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,20 +74,20 @@ static void eulersr_cons_calc_primitives_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 DenL = Den[index];
-    CCTK_REAL epsiL = epsi[index];
-    CCTK_REAL hL = h[index];
-    CCTK_REAL pL = p[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL S1L = S1[index];
-    CCTK_REAL S2L = S2[index];
-    CCTK_REAL S3L = S3[index];
-    CCTK_REAL tauL = tau[index];
-    CCTK_REAL WL = W[index];
+    CCTK_REAL DenL CCTK_ATTRIBUTE_UNUSED = Den[index];
+    CCTK_REAL epsiL CCTK_ATTRIBUTE_UNUSED = epsi[index];
+    CCTK_REAL hL CCTK_ATTRIBUTE_UNUSED = h[index];
+    CCTK_REAL pL CCTK_ATTRIBUTE_UNUSED = p[index];
+    CCTK_REAL rhoL CCTK_ATTRIBUTE_UNUSED = rho[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];
+    CCTK_REAL tauL CCTK_ATTRIBUTE_UNUSED = tau[index];
+    CCTK_REAL WL CCTK_ATTRIBUTE_UNUSED = W[index];
     
     
     /* Include user supplied include files */
@@ -95,34 +95,36 @@ static void eulersr_cons_calc_primitives_Body(cGH const * restrict const cctkGH,
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL pBar = pL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pBar = pL;
     
-    CCTK_REAL f = 10;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED f = 10.;
     
-    CCTK_REAL Z = DenL + tauL + pBar;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Z = DenL + tauL + pBar;
     
-    CCTK_REAL Ssq = SQR(S1L) + SQR(S2L) + SQR(S3L);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED Ssq = SQR(S1L) + SQR(S2L) + 
+      SQR(S3L);
     
-    CCTK_REAL vsq = Ssq*INV(SQR(Z));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED vsq = Ssq*INV(SQR(Z));
     
-    WL = INV(sqrt(1 - vsq));
+    WL = INV(sqrt(1. - 1.*vsq));
     
     rhoL = DenL*INV(WL);
     
     hL = Z*INV(rhoL*SQR(WL));
     
-    epsiL = hL - (rhoL + pBar)*INV(rhoL);
+    epsiL = hL - 1.*(rhoL + pBar)*INV(rhoL);
     
-    CCTK_REAL pEOS = epsiL*rhoL*(-1 + ToReal(gamma));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED pEOS = epsiL*rhoL*(-1. + 
+      ToReal(gamma));
     
-    f = -pBar + pEOS;
+    f = -1.*pBar + pEOS;
     
-    CCTK_REAL cs = sqrt(epsiL*INV(hL)*(-1 + 
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED cs = sqrt(epsiL*INV(hL)*(-1. + 
       ToReal(gamma))*ToReal(gamma));
     
-    CCTK_REAL df = -1 + vsq*SQR(cs);
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED df = -1. + vsq*SQR(cs);
     
-    pBar = pBar - f*INV(df);
+    pBar = pBar - 1.*f*INV(df);
     
     Z = DenL + tauL + pBar;
     
@@ -130,23 +132,23 @@ static void eulersr_cons_calc_primitives_Body(cGH const * restrict const cctkGH,
     
     vsq = Ssq*INV(SQR(Z));
     
-    WL = INV(sqrt(1 - vsq));
+    WL = INV(sqrt(1. - 1.*vsq));
     
     rhoL = DenL*INV(WL);
     
     hL = Z*INV(rhoL*SQR(WL));
     
-    epsiL = hL - (rhoL + pBar)*INV(rhoL);
+    epsiL = hL - 1.*(rhoL + pBar)*INV(rhoL);
     
-    pEOS = epsiL*rhoL*(-1 + ToReal(gamma));
+    pEOS = epsiL*rhoL*(-1. + ToReal(gamma));
     
-    f = -pBar + pEOS;
+    f = -1.*pBar + pEOS;
     
-    cs = sqrt(epsiL*INV(hL)*(-1 + ToReal(gamma))*ToReal(gamma));
+    cs = sqrt(epsiL*INV(hL)*(-1. + ToReal(gamma))*ToReal(gamma));
     
-    df = -1 + vsq*SQR(cs);
+    df = -1. + vsq*SQR(cs);
     
-    pBar = pBar - f*INV(df);
+    pBar = pBar - 1.*f*INV(df);
     
     Z = DenL + tauL + pBar;
     
@@ -154,23 +156,23 @@ static void eulersr_cons_calc_primitives_Body(cGH const * restrict const cctkGH,
     
     vsq = Ssq*INV(SQR(Z));
     
-    WL = INV(sqrt(1 - vsq));
+    WL = INV(sqrt(1. - 1.*vsq));
     
     rhoL = DenL*INV(WL);
     
     hL = Z*INV(rhoL*SQR(WL));
     
-    epsiL = hL - (rhoL + pBar)*INV(rhoL);
+    epsiL = hL - 1.*(rhoL + pBar)*INV(rhoL);
     
-    pEOS = epsiL*rhoL*(-1 + ToReal(gamma));
+    pEOS = epsiL*rhoL*(-1. + ToReal(gamma));
     
-    f = -pBar + pEOS;
+    f = -1.*pBar + pEOS;
     
-    cs = sqrt(epsiL*INV(hL)*(-1 + ToReal(gamma))*ToReal(gamma));
+    cs = sqrt(epsiL*INV(hL)*(-1. + ToReal(gamma))*ToReal(gamma));
     
-    df = -1 + vsq*SQR(cs);
+    df = -1. + vsq*SQR(cs);
     
-    pBar = pBar - f*INV(df);
+    pBar = pBar - 1.*f*INV(df);
     
     Z = DenL + tauL + pBar;
     
@@ -178,23 +180,23 @@ static void eulersr_cons_calc_primitives_Body(cGH const * restrict const cctkGH,
     
     vsq = Ssq*INV(SQR(Z));
     
-    WL = INV(sqrt(1 - vsq));
+    WL = INV(sqrt(1. - 1.*vsq));
     
     rhoL = DenL*INV(WL);
     
     hL = Z*INV(rhoL*SQR(WL));
     
-    epsiL = hL - (rhoL + pBar)*INV(rhoL);
+    epsiL = hL - 1.*(rhoL + pBar)*INV(rhoL);
     
-    pEOS = epsiL*rhoL*(-1 + ToReal(gamma));
+    pEOS = epsiL*rhoL*(-1. + ToReal(gamma));
     
-    f = -pBar + pEOS;
+    f = -1.*pBar + pEOS;
     
-    cs = sqrt(epsiL*INV(hL)*(-1 + ToReal(gamma))*ToReal(gamma));
+    cs = sqrt(epsiL*INV(hL)*(-1. + ToReal(gamma))*ToReal(gamma));
     
-    df = -1 + vsq*SQR(cs);
+    df = -1. + vsq*SQR(cs);
     
-    pBar = pBar - f*INV(df);
+    pBar = pBar - 1.*f*INV(df);
     
     Z = DenL + tauL + pBar;
     
@@ -202,31 +204,34 @@ static void eulersr_cons_calc_primitives_Body(cGH const * restrict const cctkGH,
     
     vsq = Ssq*INV(SQR(Z));
     
-    WL = INV(sqrt(1 - vsq));
+    WL = INV(sqrt(1. - 1.*vsq));
     
     rhoL = DenL*INV(WL);
     
     hL = Z*INV(rhoL*SQR(WL));
     
-    epsiL = hL - (rhoL + pBar)*INV(rhoL);
+    epsiL = hL - 1.*(rhoL + pBar)*INV(rhoL);
     
-    pEOS = epsiL*rhoL*(-1 + ToReal(gamma));
+    pEOS = epsiL*rhoL*(-1. + ToReal(gamma));
     
-    f = -pBar + pEOS;
+    f = -1.*pBar + pEOS;
     
-    cs = sqrt(epsiL*INV(hL)*(-1 + ToReal(gamma))*ToReal(gamma));
+    cs = sqrt(epsiL*INV(hL)*(-1. + ToReal(gamma))*ToReal(gamma));
     
-    df = -1 + vsq*SQR(cs);
+    df = -1. + vsq*SQR(cs);
     
-    pBar = pBar - f*INV(df);
+    pBar = pBar - 1.*f*INV(df);
     
     pL = pBar;
     
-    CCTK_REAL v1L = S1L*INV(hL*rhoL*SQR(WL));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1L = 
+      S1L*INV(hL*rhoL*SQR(WL));
     
-    CCTK_REAL v2L = S2L*INV(hL*rhoL*SQR(WL));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2L = 
+      S2L*INV(hL*rhoL*SQR(WL));
     
-    CCTK_REAL v3L = S3L*INV(hL*rhoL*SQR(WL));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3L = 
+      S3L*INV(hL*rhoL*SQR(WL));
     
     /* Copy local copies back to grid functions */
     epsi[index] = epsiL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_1.cc b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_1.cc
index 98c55bb..3d45f61 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_1.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::epsi_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::epsi_lr_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
   /* 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 eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
     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 epsiL = epsi[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL epsiL CCTK_ATTRIBUTE_UNUSED = epsi[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 DiffPlus1epsi = DiffPlus1(&epsi[index]);
-    CCTK_REAL const DiffMinus1epsi = DiffMinus1(&epsi[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*/ DiffPlus1epsi CCTK_ATTRIBUTE_UNUSED  = DiffPlus1(&epsi[index]);
+    CCTK_REAL /*const*/ DiffMinus1epsi CCTK_ATTRIBUTE_UNUSED  = DiffMinus1(&epsi[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 eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
     
     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 = DiffMinus1epsi;
     
@@ -170,9 +170,9 @@ static void eulersr_cons_calc_reconstruct_1_Body(cGH const * restrict const cctk
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL epsiLeftL = epsiL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiLeftL = epsiL - 0.5*slope;
     
-    CCTK_REAL epsiRightL = epsiL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiRightL = epsiL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     epsiLeft[index] = epsiLeftL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_2.cc b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_2.cc
index a0b8409..0205248 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_2.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::epsi_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::epsi_lr_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
   /* 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 eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
     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 epsiL = epsi[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL epsiL CCTK_ATTRIBUTE_UNUSED = epsi[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 DiffPlus2epsi = DiffPlus2(&epsi[index]);
-    CCTK_REAL const DiffMinus2epsi = DiffMinus2(&epsi[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*/ DiffPlus2epsi CCTK_ATTRIBUTE_UNUSED  = DiffPlus2(&epsi[index]);
+    CCTK_REAL /*const*/ DiffMinus2epsi CCTK_ATTRIBUTE_UNUSED  = DiffMinus2(&epsi[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 eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
     
     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 = DiffMinus2epsi;
     
@@ -170,9 +170,9 @@ static void eulersr_cons_calc_reconstruct_2_Body(cGH const * restrict const cctk
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL epsiLeftL = epsiL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiLeftL = epsiL - 0.5*slope;
     
-    CCTK_REAL epsiRightL = epsiL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiRightL = epsiL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     epsiLeft[index] = epsiLeftL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_3.cc b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_3.cc
index 6536df9..718cf41 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_3.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_reconstruct_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::epsi_lr_group","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::epsi_lr_group.");
@@ -55,34 +55,34 @@ static void eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
   /* 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 eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
     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 epsiL = epsi[index];
-    CCTK_REAL rhoL = rho[index];
-    CCTK_REAL v1L = v1[index];
-    CCTK_REAL v2L = v2[index];
-    CCTK_REAL v3L = v3[index];
+    CCTK_REAL epsiL CCTK_ATTRIBUTE_UNUSED = epsi[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 DiffPlus3epsi = DiffPlus3(&epsi[index]);
-    CCTK_REAL const DiffMinus3epsi = DiffMinus3(&epsi[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*/ DiffPlus3epsi CCTK_ATTRIBUTE_UNUSED  = DiffPlus3(&epsi[index]);
+    CCTK_REAL /*const*/ DiffMinus3epsi CCTK_ATTRIBUTE_UNUSED  = DiffMinus3(&epsi[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 eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
     
     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 eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
     
     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 = DiffMinus3epsi;
     
@@ -170,9 +170,9 @@ static void eulersr_cons_calc_reconstruct_3_Body(cGH const * restrict const cctk
     
     slope = VanLeer(slopeL,slopeR);
     
-    CCTK_REAL epsiLeftL = epsiL - 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiLeftL = epsiL - 0.5*slope;
     
-    CCTK_REAL epsiRightL = epsiL + 0.5*slope;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiRightL = epsiL + 0.5*slope;
     
     /* Copy local copies back to grid functions */
     epsiLeft[index] = epsiLeftL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_rhs_1.cc b/Examples/EulerSR/src/eulersr_cons_calc_rhs_1.cc
index cfcf0cb..dcb91a1 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_rhs_1.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_rhs_1.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulersr_cons_calc_rhs_1_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 */
   
@@ -92,41 +92,41 @@ static void eulersr_cons_calc_rhs_1_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 DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[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 tauFluxL = tauFlux[index];
-    CCTK_REAL taurhsL = taurhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[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];
+    CCTK_REAL tauFluxL CCTK_ATTRIBUTE_UNUSED = tauFlux[index];
+    CCTK_REAL taurhsL CCTK_ATTRIBUTE_UNUSED = taurhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus1DenFlux = PDplus1(&DenFlux[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 PDplus1tauFlux = PDplus1(&tauFlux[index]);
+    CCTK_REAL /*const*/ PDplus1DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&DenFlux[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]);
+    CCTK_REAL /*const*/ PDplus1tauFlux CCTK_ATTRIBUTE_UNUSED  = PDplus1(&tauFlux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus1DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus1DenFlux;
     
-    S1rhsL = S1rhsL - PDplus1S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus1S1Flux;
     
-    S2rhsL = S2rhsL - PDplus1S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus1S2Flux;
     
-    S3rhsL = S3rhsL - PDplus1S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus1S3Flux;
     
-    taurhsL = taurhsL - PDplus1tauFlux;
+    taurhsL = taurhsL - 1.*PDplus1tauFlux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_rhs_2.cc b/Examples/EulerSR/src/eulersr_cons_calc_rhs_2.cc
index e5ca803..0c379a1 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_rhs_2.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_rhs_2.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulersr_cons_calc_rhs_2_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 */
   
@@ -92,41 +92,41 @@ static void eulersr_cons_calc_rhs_2_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 DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[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 tauFluxL = tauFlux[index];
-    CCTK_REAL taurhsL = taurhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[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];
+    CCTK_REAL tauFluxL CCTK_ATTRIBUTE_UNUSED = tauFlux[index];
+    CCTK_REAL taurhsL CCTK_ATTRIBUTE_UNUSED = taurhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus2DenFlux = PDplus2(&DenFlux[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 PDplus2tauFlux = PDplus2(&tauFlux[index]);
+    CCTK_REAL /*const*/ PDplus2DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&DenFlux[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]);
+    CCTK_REAL /*const*/ PDplus2tauFlux CCTK_ATTRIBUTE_UNUSED  = PDplus2(&tauFlux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus2DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus2DenFlux;
     
-    S1rhsL = S1rhsL - PDplus2S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus2S1Flux;
     
-    S2rhsL = S2rhsL - PDplus2S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus2S2Flux;
     
-    S3rhsL = S3rhsL - PDplus2S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus2S3Flux;
     
-    taurhsL = taurhsL - PDplus2tauFlux;
+    taurhsL = taurhsL - 1.*PDplus2tauFlux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_rhs_3.cc b/Examples/EulerSR/src/eulersr_cons_calc_rhs_3.cc
index b494cdc..f75b204 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_rhs_3.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_rhs_3.cc
@@ -27,7 +27,7 @@ extern "C" void eulersr_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 */, "EulerSR::Den_grouprhs","flat");
   if (ierr < 0)
     CCTK_WARN(1, "Failed to register flat BC for EulerSR::Den_grouprhs.");
@@ -49,34 +49,34 @@ static void eulersr_cons_calc_rhs_3_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 */
   
@@ -92,41 +92,41 @@ static void eulersr_cons_calc_rhs_3_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 DenFluxL = DenFlux[index];
-    CCTK_REAL DenrhsL = Denrhs[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 tauFluxL = tauFlux[index];
-    CCTK_REAL taurhsL = taurhs[index];
+    CCTK_REAL DenFluxL CCTK_ATTRIBUTE_UNUSED = DenFlux[index];
+    CCTK_REAL DenrhsL CCTK_ATTRIBUTE_UNUSED = Denrhs[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];
+    CCTK_REAL tauFluxL CCTK_ATTRIBUTE_UNUSED = tauFlux[index];
+    CCTK_REAL taurhsL CCTK_ATTRIBUTE_UNUSED = taurhs[index];
     
     
     /* Include user supplied include files */
     
     /* Precompute derivatives */
-    CCTK_REAL const PDplus3DenFlux = PDplus3(&DenFlux[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 PDplus3tauFlux = PDplus3(&tauFlux[index]);
+    CCTK_REAL /*const*/ PDplus3DenFlux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&DenFlux[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]);
+    CCTK_REAL /*const*/ PDplus3tauFlux CCTK_ATTRIBUTE_UNUSED  = PDplus3(&tauFlux[index]);
     
     /* Calculate temporaries and grid functions */
-    DenrhsL = DenrhsL - PDplus3DenFlux;
+    DenrhsL = DenrhsL - 1.*PDplus3DenFlux;
     
-    S1rhsL = S1rhsL - PDplus3S1Flux;
+    S1rhsL = S1rhsL - 1.*PDplus3S1Flux;
     
-    S2rhsL = S2rhsL - PDplus3S2Flux;
+    S2rhsL = S2rhsL - 1.*PDplus3S2Flux;
     
-    S3rhsL = S3rhsL - PDplus3S3Flux;
+    S3rhsL = S3rhsL - 1.*PDplus3S3Flux;
     
-    taurhsL = taurhsL - PDplus3tauFlux;
+    taurhsL = taurhsL - 1.*PDplus3tauFlux;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerSR/src/eulersr_cons_calc_zero_rhs.cc b/Examples/EulerSR/src/eulersr_cons_calc_zero_rhs.cc
index 382e07c..4bbcf4f 100644
--- a/Examples/EulerSR/src/eulersr_cons_calc_zero_rhs.cc
+++ b/Examples/EulerSR/src/eulersr_cons_calc_zero_rhs.cc
@@ -31,34 +31,34 @@ static void eulersr_cons_calc_zero_rhs_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 */
   
@@ -74,7 +74,7 @@ static void eulersr_cons_calc_zero_rhs_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 */
     
@@ -85,35 +85,35 @@ static void eulersr_cons_calc_zero_rhs_Body(cGH const * restrict const cctkGH, i
     /* 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.;
     
-    S1rhsL = 0;
+    S1rhsL = 0.;
     
-    S1rhsL = 0;
+    S1rhsL = 0.;
     
-    CCTK_REAL S2rhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S2rhsL = 0.;
     
-    S2rhsL = 0;
+    S2rhsL = 0.;
     
-    S2rhsL = 0;
+    S2rhsL = 0.;
     
-    CCTK_REAL S3rhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED S3rhsL = 0.;
     
-    S3rhsL = 0;
+    S3rhsL = 0.;
     
-    S3rhsL = 0;
+    S3rhsL = 0.;
     
-    CCTK_REAL taurhsL = 0;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED taurhsL = 0.;
     
-    taurhsL = 0;
+    taurhsL = 0.;
     
-    taurhsL = 0;
+    taurhsL = 0.;
     
     /* Copy local copies back to grid functions */
     Denrhs[index] = DenrhsL;
diff --git a/Examples/EulerSR/src/eulersr_initial_shock.cc b/Examples/EulerSR/src/eulersr_initial_shock.cc
index 80a8705..044b04f 100644
--- a/Examples/EulerSR/src/eulersr_initial_shock.cc
+++ b/Examples/EulerSR/src/eulersr_initial_shock.cc
@@ -31,34 +31,34 @@ static void eulersr_initial_shock_Body(cGH const * restrict const cctkGH, int co
   /* 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 eulersr_initial_shock_Body(cGH const * restrict const cctkGH, int co
     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,22 +86,22 @@ static void eulersr_initial_shock_Body(cGH const * restrict const cctkGH, int co
     /* Precompute derivatives */
     
     /* Calculate temporaries and grid functions */
-    CCTK_REAL X = xL;
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED X = xL;
     
-    CCTK_REAL rhoL = ToReal(rhoL0) + StepFunction(X)*(-ToReal(rhoL0) + 
-      ToReal(rhoR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED rhoL = ToReal(rhoL0) + 
+      StepFunction(X)*(-1.*ToReal(rhoL0) + ToReal(rhoR0));
     
-    CCTK_REAL v1L = ToReal(vL0) + StepFunction(X)*(-ToReal(vL0) + 
-      ToReal(vR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v1L = ToReal(vL0) + 
+      StepFunction(X)*(-1.*ToReal(vL0) + ToReal(vR0));
     
-    CCTK_REAL v2L = ToReal(vL0) + StepFunction(X)*(-ToReal(vL0) + 
-      ToReal(vR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v2L = ToReal(vL0) + 
+      StepFunction(X)*(-1.*ToReal(vL0) + ToReal(vR0));
     
-    CCTK_REAL v3L = ToReal(vL0) + StepFunction(X)*(-ToReal(vL0) + 
-      ToReal(vR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED v3L = ToReal(vL0) + 
+      StepFunction(X)*(-1.*ToReal(vL0) + ToReal(vR0));
     
-    CCTK_REAL epsiL = ToReal(epsiL0) + StepFunction(X)*(-ToReal(epsiL0) + 
-      ToReal(epsiR0));
+    CCTK_REAL CCTK_ATTRIBUTE_UNUSED epsiL = ToReal(epsiL0) + 
+      StepFunction(X)*(-1.*ToReal(epsiL0) + ToReal(epsiR0));
     
     /* Copy local copies back to grid functions */
     epsi[index] = epsiL;