* add hybrid EOS

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEOS/EOS_Omni/EOS_Omni@3 8e189c6b-2ab8-4400-aa02-70a9cfce18b9

1 files changed, 33 insertions, 2 deletions

diff --git a/param.ccl b/param.ccl
index 6b8acf7..870d04d 100644
--- a/param.ccl
+++ b/param.ccl
@@ -4,7 +4,12 @@
 restricted:
 
 # poly EOS
-CCTK_REAL poly_gamma "Adiabatic Index for poly EOS"
+REAL poly_gamma "Adiabatic Index for poly EOS"
+{
+ : :: ""
+} 2.0
+
+REAL poly_gamma_ini "Initial Adiabatic Index for poly EOS"
 {
  : :: ""
 } 2.0
@@ -15,7 +20,7 @@ REAL poly_k "Polytropic constant in c=G=Msun=1"
 } 100.0
 
 # gamma-law EOS
-CCTK_REAL gl_gamma "Adiabatic Index for gamma-law EOS"
+REAL gl_gamma "Adiabatic Index for gamma-law EOS"
 {
  : :: ""
 } 2.0
@@ -24,3 +29,29 @@ REAL gl_k "Polytropic constant in c=G=Msun=1 for gamma-law EOS"
 {
  : :: ""
 } 100.0
+
+# hybrid EOS
+REAL hybrid_gamma1 "subnuclear adiabatic Index for hybrid EOS"
+{
+ : :: ""
+} 1.325
+
+REAL hybrid_gamma2 "subnuclear adiabatic Index for hybrid EOS"
+{
+ : :: ""
+} 2.5
+
+REAL hybrid_gamma_th "Thermal gamma for hybrid EOS"
+{
+ : :: ""
+} 1.5
+
+REAL hybrid_k1 "Polytropic constant in c=G=Msun=1 for hybrid EOS" 
+{
+ : :: ""
+} 0.4640517
+
+REAL hybrid_rho_nuc "Density at which to switch betwen gammas; c=G=Msun=1" 
+{
+ : :: ""
+} 3.238607e-4