diff options
| -rwxr-xr-x | interface.ccl | 2 | ||||
| -rwxr-xr-x | param.ccl | 4 | ||||
| -rwxr-xr-x | schedule.ccl | 12 |
3 files changed, 18 insertions, 0 deletions
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" @@ -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 |