CCTK_REAL chi_floor "Floor for chi" STEERABLE=always
{
  (0:* :: ""
} 1.0e-10

CCTK_REAL alphaG_floor "Floor for alphaG" STEERABLE=always
{
  (0:* :: ""
} 1.0e-10

  READS: chi(everywhere)

  READS: alphaG(everywhere)
  WRITES: chi(everywhere)

  WRITES: alphaG(everywhere)

#include <cctk_Parameters.h>

  DECLARE_CCTK_PARAMETERS;

  const GF3D<CCTK_REAL, 0, 0, 0> gf_chi_(cctkGH, chi);

  const GF3D<CCTK_REAL, 0, 0, 0> gf_alphaG_(cctkGH, alphaG);

        const CCTK_REAL chi_old = gf_chi_(p.I);
        const CCTK_REAL alphaG_old = gf_alphaG_(p.I);

        // Enforce floors

        const CCTK_REAL chi = fmax(chi_floor, chi_old);
        const CCTK_REAL alphaG = fmax(alphaG_floor, alphaG_old);

        const CCTK_REAL chi_old = 1 / cbrt(detgammat_old);

        const CCTK_REAL chi1_old = 1 / cbrt(detgammat_old);

            [&](int a, int b) { return chi_old * gammat_old(a, b); });

            [&](int a, int b) { return chi1_old * gammat_old(a, b); });

        gf_chi_(p.I) = chi;

        gf_alphaG_(p.I) = alphaG;