# Interface definition for thorn SteadyStateReactor
# $Header$

implements: reactor

inherits: grid

private:

CCTK_REAL fields type=GF dim=3 timelevels=1
{
   theta, fA, fB, fC, fD
} "The dimensionless temperature and concentrations"

CCTK_REAL deltas type=GF dim=3 timelevels=1
{
   dtheta, dfA, dfB, dfC, dfD
} "The linearized quantities"

CCTK_REAL stencil type=SCALAR
{
   ctc, cte, ctw, cfc, cfe, cfw
} "The stencil coefficients"

CCTK_REAL params type=SCALAR
{
   f0, MWave, rho0, G, Dez, cA0
} "Some derived parameters"

CCTK_REAL jacobian type=GF dim=3 timelevels=1
{
   Jttc, Jtte, Jttw, JtAc, JtBc, JtCc, JtDc,
   JAtc, JAte, JAtw, JAAc, JAAe, JAAw, JABc, JABe, JABw, JACc, JACe, JACw,
      JADc, JADe, JADw,
   JBtc, JBte, JBtw, JBAc, JBAe, JBAw, JBBc, JBBe, JBBw, JBCc, JBCe, JBCw,
      JBDc, JBDe, JBDw,
   JCtc, JCte, JCtw, JCAc, JCAe, JCAw, JCBc, JCBe, JCBw, JCCc, JCCe, JCCw,
      JCDc, JCDe, JCDw,
   JDtc, JDte, JDtw, JDAc, JDAe, JDAw, JDBc, JDBe, JDBw, JDCc, JDCe, JDCw,
      JDDc, JDDe, JDDw
} "Representation of the Jacobian matrix"

CCTK_REAL genterms type=GF dim=3 timelevels=1
{
   Hstar, RAstar, RBstar, RCstar, RDstar
} "Generation terms from kinetics"

CCTK_REAL rhss type=GF dim=3 timelevels=1
{
   Rt, RA, RB, RC, RD
} "Right-hand-sides for the linear solves"

CCTK_REAL nonlinterms type=GF dim=3 timelevels=1
{
   Nt, NA, NB, NC, ND
} "Nonlinear terms"

CCTK_REAL auxiliaries type=GF dim=3 timelevels=1
{
   uz
} "Some auxiliary variables"

CCTK_REAL cgworkspace type=GF dim=3 timelevels=1
{
   Ast, AsA, AsB, AsC, AsD,
   Apt, ApA, ApB, ApC, ApD,
   sumGF
} "Some arrays for the CG solver"

CCTK_REAL cgw type=GF dim=3 timelevels=1
{
   wt, wA, wB, wC, wD
} "Some more arrays for the CG solver"

CCTK_REAL cgp type=GF dim=3 timelevels=1
{
   pt, pA, pB, pC, pD
} "Even more arrays for the CG solver"