GRHydro_InitData: initialize c2p2c values as parameters and update routine to test inversions in the gr context.

From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@192 ac85fae7-cede-4708-beff-ae01c7fa1c26

1 files changed, 91 insertions, 0 deletions

diff --git a/param.ccl b/param.ccl
index 67e51e1..bb9dad6 100644
--- a/param.ccl
+++ b/param.ccl
@@ -126,6 +126,8 @@ BOOLEAN attenuate_atmosphere "Attenuate the velocity in the atmosphere"
 {
 } "no"
 
+# Initial magnetic field used by different tests
+
 REAL Bx_init "Initial B-field in the x-dir"
 {
   *:*     :: "Anything"
@@ -141,6 +143,95 @@ REAL Bz_init "Initial B-field in the z-dir"
   *:*     :: "Anything"
 } 0.0
 
+# Initial primitive values:
+REAL rho_init "Initial rest mass density"
+{
+ (0:*    :: "Anything positive."
+} 1.0d-6
+
+REAL velx_init "Initial x velocity" 
+{
+ *:*    :: "Anything."
+} 1.0d-1
+
+REAL vely_init "Initial y velocity" 
+{
+ *:*    :: "Anything."
+} 1.0d-1
+
+REAL velz_init "Initial z velocity" 
+{
+ *:*    :: "Anything."
+} 1.0d-1
+
+REAL eps_init "Initial specific internal energy"
+{
+ (0:*    :: "Anything positive."
+} 1.0d-6
+
+REAL press_init "Initial pressure"
+{
+ (0:*    :: "Anything positive."
+} 6.666666666666667d-7
+
+# Initial conservative values used by c2p tests:
+REAL dens_init "Initial conserved mass density"
+{
+ (0:*    :: "Anything positive."
+} 1.29047362
+
+REAL sx_init "Initial x component of conserved momentum density"
+{
+ *:*    :: "Anything."
+} 0.166666658
+
+REAL sy_init "Initial y component of conserved momentum density"
+{
+ *:*    :: "Anything."
+} 0.166666658
+
+REAL sz_init "Initial z component of conserved momentum density"
+{
+ *:*    :: "Anything."
+} 0.166666658
+
+REAL tau_init "Initial conserved total energy density"
+{
+ (0:*    :: "Anything positive."
+} 0.484123939
+
+# Initial values for 3-metric components. Default to euclidian 3-metric
+REAL gxx_init "Initial xx metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 1.0
+
+REAL gxy_init "Initial xy metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 0.0
+
+REAL gxz_init "Initial xz metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 0.0
+
+REAL gyy_init "Initial yy metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 1.0
+
+REAL gyz_init "Initial yz metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 0.0
+
+REAL gzz_init "Initial zz metric componenent"
+{
+ *:*    :: "Anything, but be carefull to set a positive definite 3-metric!"
+} 1.0
+
+
 KEYWORD monopole_type "Which kind of monopole?"
 {
   "Point" :: "Single point with Bx /= 0"