From 32ff5057a5aeb499616aef3abcb00d93397ff13b Mon Sep 17 00:00:00 2001
From: knarf <knarf@57fe0bb3-ccba-405f-9b23-de0201f165b7>
Date: Mon, 28 Sep 2009 21:03:51 +0000
Subject: define Y_e plus parameter use_Y_e to use it

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinBase/HydroBase/trunk@10 57fe0bb3-ccba-405f-9b23-de0201f165b7
---
 interface.ccl |  2 ++
 param.ccl     |  4 ++++
 schedule.ccl  | 12 ++++++++++++
 3 files changed, 18 insertions(+)

diff --git a/interface.ccl b/interface.ccl
index 742b784..870f94f 100755
--- a/interface.ccl
+++ b/interface.ccl
@@ -11,3 +11,5 @@ CCTK_REAL press type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase:
 CCTK_REAL eps   type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "internal energy"
 
 CCTK_REAL vel[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="U" interpolator="matter"' "velocity"
+
+CCTK_REAL Y_e[3] type = GF Timelevels = 3 tags='ProlongationParameter="HydroBase::prolongation_type" tensortypealias="Scalar" interpolator="matter"' "electron fraction"
diff --git a/param.ccl b/param.ccl
index 7a801d9..573f46c 100755
--- a/param.ccl
+++ b/param.ccl
@@ -18,3 +18,7 @@ STRING prolongation_type "The prolongation operator used by Carpet for HydroBase
   "WENO" :: "Fifth order WENO operators; only fifth order is implemented"
   ".*"   :: "Anything else"
 } "ENO"
+
+BOOLEAN use_Y_e "Enable storage for Y_e"
+{
+} "no"
diff --git a/schedule.ccl b/schedule.ccl
index 2b3b5ac..312aef2 100755
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -11,6 +11,10 @@ if (timelevels == 3)
   STORAGE:press[3]
   STORAGE:eps[3]
   STORAGE:vel[3]
+  if (use_Y_e)
+  {
+    STORAGE:Y_e[3]
+  }
 }
 else if (timelevels == 2)
 {
@@ -18,6 +22,10 @@ else if (timelevels == 2)
   STORAGE:press[2]
   STORAGE:eps[2]
   STORAGE:vel[2]
+  if (use_Y_e)
+  {
+    STORAGE:Y_e[2]
+  }
 }
 else if (timelevels == 1)
 {
@@ -25,6 +33,10 @@ else if (timelevels == 1)
   STORAGE:press[1]
   STORAGE:eps[1]
   STORAGE:vel[1]
+  if (use_Y_e)
+  {
+    STORAGE:Y_e[1]
+  }
 }
 
 schedule group HydroBase_Initial AT Initial after ADMBase_InitialData after ADMBase_InitialGauge after IOUtil_RecoverIDFromDatafiles before ADMBase_PostInitial before SetTmunu
-- 
cgit v1.2.3