C The Robertson-Walker metric doesn't work ==> move it to ../../archive/

C
C The argument  rama  is the R(t) in the
C Robertson-Walker metric, and if this is passed correctly then this
C subroutine computes the correct metric.  But the rest of this thorn
C doesn't know to pass this value. :( :(  See Mitica's Cosmo thorn
C for a better way to get this metric.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@177 e296648e-0e4f-0410-bd07-d597d9acff87

9 files changed, 135 insertions, 134 deletions

diff --git a/doc/documentation.tex b/doc/documentation.tex
index ee36522..3467456 100644
--- a/doc/documentation.tex
+++ b/doc/documentation.tex
@@ -244,8 +244,8 @@ Model Name
 \multicolumn{3}{l}{\bf Cosmological spacetimes}				\\
 {\tt "Lemaitre"}
 	& Yes	& Lemaitre-type spacetime				\\
-{\tt "Robertson-Walker"}
-	& Yes	& Robertson-Walker spacetime				\\
+%%{\tt "Robertson-Walker"}
+%%	& Yes	& Robertson-Walker spacetime				\\
 {\tt "de Sitter"}
 	& Yes	& de~Sitter spacetime					\\
 {\tt "de Sitter+Lambda"}
@@ -1000,36 +1000,36 @@ and the scale factor (radius) of the universe at time $t = 0$,
 $R_0 = \verb|Lemaitre__R0|$.
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
-\subsection{Robertson-Walker spacetime}
-\label{AEIThorns/Exact/sect-Robertson-Walker}
-
-\verb|Exact::exact_model = "Robertson-Walker"| specifies a 
-Robertson-Walker spacetime  as described in Hawking and Ellis section~5.3
-and MTW section~27.11 (see also gr-qc/0110031),
-transformed to the usual Cactus $(t,x,y,z)$ Cartesian-topology coordinates.
-The general Robertson-Walker line element in $(t,r,\theta,\phi)$ coordinates
-is
-\begin{equation}
-ds^2 = -dt^2 + R(t)^2 \left[ \frac{dr^2}{1 - kr^2} + r^2 \, d\Omega^2 \right]
-\end{equation}
-
-The physics parameters are
-the scale factor $R(t)$ at time $t = 0$, $R_0 = \verb|Robertson_Walker__R0|$,
-a parameter $\rho = \verb|Robertson_Walker__rho|$ which is related to
-the actual value of the matter density in the Universe,
-the geometry curvature parameter $k = \verb|Robertson_Walker__k|$,
-which can take (only) the values $k=-1$, $0$, or $+1$, corresponding
-to open, flat, or closed 3-geometries, and finally
-the Boolean parameter \verb|Robertson_Walker__pressure| to select
-whether or not to include pressure terms in the model.  If pressure
-is included we have a radiation-dominated universe $p = \frac{1}{3} \rho$;
-if pressure is not included we have a matter-dominated universe $p=0$.
-
-For a good simulation it is necessary to give good numerical values
-for the above parameters (they are very strictly related, through the
-Einstein equations).  See gr-qc/0110031 for some examples.
-
+%%
+%%\subsection{Robertson-Walker spacetime}
+%%\label{AEIThorns/Exact/sect-Robertson-Walker}
+%%
+%%\verb|Exact::exact_model = "Robertson-Walker"| specifies a 
+%%Robertson-Walker spacetime  as described in Hawking and Ellis section~5.3
+%%and MTW section~27.11 (see also gr-qc/0110031),
+%%transformed to the usual Cactus $(t,x,y,z)$ Cartesian-topology coordinates.
+%%The general Robertson-Walker line element in $(t,r,\theta,\phi)$ coordinates
+%%is
+%%\begin{equation}
+%%ds^2 = -dt^2 + R(t)^2 \left[ \frac{dr^2}{1 - kr^2} + r^2 \, d\Omega^2 \right]
+%%\end{equation}
+%%
+%%The physics parameters are
+%%the scale factor $R(t)$ at time $t = 0$, $R_0 = \verb|Robertson_Walker__R0|$,
+%%a parameter $\rho = \verb|Robertson_Walker__rho|$ which is related to
+%%the actual value of the matter density in the Universe,
+%%the geometry curvature parameter $k = \verb|Robertson_Walker__k|$,
+%%which can take (only) the values $k=-1$, $0$, or $+1$, corresponding
+%%to open, flat, or closed 3-geometries, and finally
+%%the Boolean parameter \verb|Robertson_Walker__pressure| to select
+%%whether or not to include pressure terms in the model.  If pressure
+%%is included we have a radiation-dominated universe $p = \frac{1}{3} \rho$;
+%%if pressure is not included we have a matter-dominated universe $p=0$.
+%%
+%%For a good simulation it is necessary to give good numerical values
+%%for the above parameters (they are very strictly related, through the
+%%Einstein equations).  See gr-qc/0110031 for some examples.
+%%
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 \subsection{de~Sitter spacetime}
diff --git a/interface.ccl b/interface.ccl
index d7db872..64ad537 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -26,9 +26,10 @@ decoded_exact_model
 
 # ***** parameters for models which have nonzero stress-energy tensor terms
 
-# parameters for Robertson-Walker spacetime
-Robertson_Walker___k
-Robertson_Walker___pressure
+### this metric doesn't work and has been moved to ../archive/
+### parameters for Robertson-Walker spacetime
+##Robertson_Walker___k
+##Robertson_Walker___pressure
 } "parameters copied to grid scalars so CalcTmunu code sees them in evolution thorns"
 
 real Exact_pars_real type=SCALAR
@@ -42,9 +43,10 @@ Lemaitre___kappa
 Lemaitre___Lambda
 Lemaitre___epsilon0
 Lemaitre___R0
-# parameters for Robertson-Walker spacetime
-Robertson_Walker___R0
-Robertson_Walker___rho
+### this metric doesn't work and has been moved to ../archive/
+### parameters for Robertson-Walker spacetime
+##Robertson_Walker___R0
+##Robertson_Walker___rho
 # parameters for de Sitter spacetime
 de_Sitter___scale
 # parameters for de Sitter spacetime with cosmological constant
@@ -64,16 +66,17 @@ constant_density_star___radius
 } "parameters copied to grid scalars so CalcTmunu code sees them in evolution thorns"
 
 ################################################################################
-
-#
-# temp variable for some of the stress-energy tensor computations
-#
-protected:
-real Exact_robwal_var TYPE=GF
-{
-   raza
-} "Raza universului Robertson Walker"
-
+##
+### Robertson-Walker metric doesn't work and has been moved to ../archive/
+###
+### temp variable for some of the stress-energy tensor computations
+###
+##protected:
+##real Exact_robwal_var TYPE=GF
+##{
+##   raza
+##} "Raza universului Robertson Walker"
+##
 ################################################################################
 
 private:
diff --git a/param.ccl b/param.ccl
index a8702a7..b0f2a9e 100644
--- a/param.ccl
+++ b/param.ccl
@@ -202,7 +202,7 @@ KEYWORD exact_model "The exact solution/coordinates used in thorn exact"
 # cosmological spacetimes
 #
 "Lemaitre"             :: "Lemaitre-type spacetime"
-"Robertson-Walker"     :: "Robertson-Walker spacetime"
+##"Robertson-Walker"     :: "Robertson-Walker spacetime"
 "de Sitter"            :: "de Sitter spacetime (R-W cosmology, near t=0, p=0)"
 "de Sitter+Lambda"     :: "de Sitter spacetime with cosmological constant"
 "anti-de Sitter+Lambda":: "anti-de Sitter spacetime with cosmological constant"
@@ -622,34 +622,34 @@ REAL Lemaitre__R0 "Lemaitre: scale factor (radius) of the universe at time t=0"
 } 1.0
 
 ################################################################################
-
-#
-# parameters for Robertson-Walker spacetime
-#
-
-REAL Robertson_Walker__R0 \
-  "Robertson-Walker: scale factor (radius) of the universe at time t=0"
-{
-(0.0:* :: "any positive real number"
-} 0.1
-
-REAL Robertson_Walker__rho "Robertson-Walker: density parameter"
-{
-0.0:* :: "any real number >= 0"
-} 0.1
-
-INT Robertson_Walker__k "Robertson-Walker: geometry type parameter"
-{
--1:1:1 :: "one of the values -1, 0, 1"
-} 0
-
-# true ==> include pressure terms (radiation-dominated universe)
-# false ==> don't include pressure terms (matter-dominated universe)
-BOOLEAN Robertson_Walker__pressure \
-  "Robertson-Walker: are pressure terms included?"
-{
-} "true"
-
+##
+###
+### parameters for Robertson-Walker spacetime
+###
+##
+##REAL Robertson_Walker__R0 \
+##  "Robertson-Walker: scale factor (radius) of the universe at time t=0"
+##{
+##(0.0:* :: "any positive real number"
+##} 0.1
+##
+##REAL Robertson_Walker__rho "Robertson-Walker: density parameter"
+##{
+##0.0:* :: "any real number >= 0"
+##} 0.1
+##
+##INT Robertson_Walker__k "Robertson-Walker: geometry type parameter"
+##{
+##-1:1:1 :: "one of the values -1, 0, 1"
+##} 0
+##
+### true ==> include pressure terms (radiation-dominated universe)
+### false ==> don't include pressure terms (matter-dominated universe)
+##BOOLEAN Robertson_Walker__pressure \
+##  "Robertson-Walker: are pressure terms included?"
+##{
+##} "true"
+##
 ################################################################################
 
 #
diff --git a/src/decode_pars.F77 b/src/decode_pars.F77
index 7790171..abe55d7 100644
--- a/src/decode_pars.F77
+++ b/src/decode_pars.F77
@@ -83,8 +83,9 @@ c black hole spacetimes
 c cosmological spacetimes
 	elseif (CCTK_Equals(exact_model, "Lemaitre") .ne. 0) then
 	     decoded_exact_model = EXACT__Lemaitre
-	elseif (CCTK_Equals(exact_model, "Robertson-Walker") .ne. 0) then
-	     decoded_exact_model = EXACT__Robertson_Walker
+C this metric doesn't work and has been moved to ../archive/
+CC	elseif (CCTK_Equals(exact_model, "Robertson-Walker") .ne. 0) then
+CC	     decoded_exact_model = EXACT__Robertson_Walker
 	elseif (CCTK_Equals(exact_model, "de Sitter") .ne. 0) then
 	     decoded_exact_model = EXACT__de_Sitter
 	elseif (CCTK_Equals(exact_model, "de Sitter+Lambda") .ne. 0) then
@@ -139,15 +140,16 @@ c
 	Lemaitre___R0       = Lemaitre__R0
 
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c parameters for Robertson-Walker spacetime
-c
-	Robertson_Walker___R0       = Robertson_Walker__R0
-	Robertson_Walker___rho      = Robertson_Walker__rho
-	Robertson_Walker___k        = Robertson_Walker__k
-	Robertson_Walker___pressure = Robertson_Walker__pressure
-
+CCC
+CCC this metric doesn't work and has been moved to ../archive/
+CCc
+CCc parameters for Robertson-Walker spacetime
+CCc
+CC	Robertson_Walker___R0       = Robertson_Walker__R0
+CC	Robertson_Walker___rho      = Robertson_Walker__rho
+CC	Robertson_Walker___k        = Robertson_Walker__k
+CC	Robertson_Walker___pressure = Robertson_Walker__pressure
+CC
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
 c
diff --git a/src/include/Scalar_CalcTmunu.inc b/src/include/Scalar_CalcTmunu.inc
index 50020cd..4a06a41 100644
--- a/src/include/Scalar_CalcTmunu.inc
+++ b/src/include/Scalar_CalcTmunu.inc
@@ -121,36 +121,37 @@ c
           Txz = Txz
 
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c Robertson-Walker spacetime
-C
-        elseif (decoded_exact_model .eq. EXACT__Robertson_Walker) then
-          rr2  = x(i,j,k)*x(i,j,k)+y(i,j,k)*y(i,j,k)+z(i,j,k)*z(i,j,k)
-
-          if (Robertson_Walker___pressure .gt. 0) then
-            aha1  = Robertson_Walker___k * (Robertson_Walker___R0**2)
-     &                                  / (8.0D0*pi*(raza(i,j,k)**2))
-            aha2 = Robertson_Walker___k/(1.0D0 - Robertson_Walker___k*rr2) 
-
-            Ttt =  Ttt + 3.0D0*aha1/(raza(i,j,k)*raza(i,j,k))
-            Txx =  Txx + aha1*(1.0D0 + aha2*x(i,j,k)*x(i,j,k))
-            Tyy =  Tyy + aha1*(1.0D0 + aha2*y(i,j,k)*y(i,j,k))
-            Tzz =  Tzz + aha1*(1.0D0 + aha2*z(i,j,k)*z(i,j,k))
-            Txy =  Txy + aha1*aha2*x(i,j,k)*y(i,j,k)
-            Txz =  Txz + aha1*aha2*x(i,j,k)*z(i,j,k)
-            Tyz =  Tyz + aha1*aha2*y(i,j,k)*y(i,j,k)
-          else
-            Ttt = Ttt+Robertson_Walker___rho * (Robertson_Walker___R0**3)
-     &                                      / (raza(i,j,k)**3)
-            Txx = Txx
-            Tyy = Tyy
-            Tzz = Tzz
-            Txy = Txy
-            Txz = Txz
-            Tyz = Tyz
-          endif
-
+ccc
+ccc this metric doesn't work and has been moved to ../../archive/
+ccc
+ccc Robertson-Walker spacetime
+ccc
+cc        elseif (decoded_exact_model .eq. EXACT__Robertson_Walker) then
+cc          rr2  = x(i,j,k)*x(i,j,k)+y(i,j,k)*y(i,j,k)+z(i,j,k)*z(i,j,k)
+cc
+cc          if (Robertson_Walker___pressure .gt. 0) then
+cc            aha1  = Robertson_Walker___k * (Robertson_Walker___R0**2)
+cc     &                                  / (8.0D0*pi*(raza(i,j,k)**2))
+cc            aha2 = Robertson_Walker___k/(1.0D0 - Robertson_Walker___k*rr2) 
+cc
+cc            Ttt =  Ttt + 3.0D0*aha1/(raza(i,j,k)*raza(i,j,k))
+cc            Txx =  Txx + aha1*(1.0D0 + aha2*x(i,j,k)*x(i,j,k))
+cc            Tyy =  Tyy + aha1*(1.0D0 + aha2*y(i,j,k)*y(i,j,k))
+cc            Tzz =  Tzz + aha1*(1.0D0 + aha2*z(i,j,k)*z(i,j,k))
+cc            Txy =  Txy + aha1*aha2*x(i,j,k)*y(i,j,k)
+cc            Txz =  Txz + aha1*aha2*x(i,j,k)*z(i,j,k)
+cc            Tyz =  Tyz + aha1*aha2*y(i,j,k)*y(i,j,k)
+cc          else
+cc            Ttt = Ttt+Robertson_Walker___rho * (Robertson_Walker___R0**3)
+cc     &                                      / (raza(i,j,k)**3)
+cc            Txx = Txx
+cc            Tyy = Tyy
+cc            Tzz = Tzz
+cc            Txy = Txy
+cc            Txz = Txz
+cc            Tyz = Tyz
+cc          endif
+cc
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
 c
diff --git a/src/include/param_defs.inc b/src/include/param_defs.inc
index 45642ba..bc3b168 100644
--- a/src/include/param_defs.inc
+++ b/src/include/param_defs.inc
@@ -49,7 +49,8 @@ c black hole spacetimes
 
 c cosmological spacetimes
 #define EXACT__Lemaitre			50
-#define EXACT__Robertson_Walker		51
+CC this metric doesn't work, and has been moved to ./archive/
+CC#define EXACT__Robertson_Walker	51
 #define EXACT__de_Sitter		52
 #define EXACT__de_Sitter_Lambda		53
 #define EXACT__anti_de_Sitter_Lambda	54
diff --git a/src/metric.F77 b/src/metric.F77
index 86b4cbe..106b8bb 100644
--- a/src/metric.F77
+++ b/src/metric.F77
@@ -191,14 +191,15 @@ c
      $        guxx, guyy, guzz, guxy, guyz, guxz,
      $        psi, Tmunu_flag)
 
-      elseif (decoded_exact_model .eq. EXACT__Robertson_Walker) then
-         call Exact__Robertson_Walker(
-     $        x, y, z, t,
-     $        gdtt, gdtx, gdty, gdtz,
-     $        gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
-     $        gutt, gutx, guty, gutz,
-     $        guxx, guyy, guzz, guxy, guyz, guxz,
-     $        psi, Tmunu_flag, rama)
+CC this metric doesn't work and has been moved to ../archive/
+CC      elseif (decoded_exact_model .eq. EXACT__Robertson_Walker) then
+CC         call Exact__Robertson_Walker(
+CC     $        x, y, z, t,
+CC     $        gdtt, gdtx, gdty, gdtz,
+CC     $        gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
+CC     $        gutt, gutx, guty, gutz,
+CC     $        guxx, guyy, guzz, guxy, guyz, guxz,
+CC     $        psi, Tmunu_flag, rama)
 
       elseif (decoded_exact_model .eq. EXACT__de_Sitter) then
          call Exact__de_Sitter(
diff --git a/src/metrics/README b/src/metrics/README
index dba8912..f027d9b 100644
--- a/src/metrics/README
+++ b/src/metrics/README
@@ -20,7 +20,6 @@ inverse 4-metric $g^{ab}$:
      $       guxx, guyy, guzz, guxy, guyz, guzx
       CCTK_REAL psi
       LOGICAL   Tmunu_flag
-      CCTK_REAL rama
 	...
       RETURN
       END
@@ -34,11 +33,6 @@ to set the stress-energy tensor (i.e. if it's non-vacuum and/or it has
 a cosmological constant, or .false. if it doesn't (i.e. if it's vacuum
 and has no cosmological constant).
 
-The "Robertson/Walker" model has an extra argument  rama .
-Alas, Mitica neglected to document what it means :( :( :(,
-and so far I (Jonathan) haven't found any higher-level routine which
-uses it.
-
 
 
 FIXME: Wouldn't it be cleaner to just compute $g_{ab}$,
diff --git a/src/metrics/make.code.defn b/src/metrics/make.code.defn
index 184c302..b50ab0f 100644
--- a/src/metrics/make.code.defn
+++ b/src/metrics/make.code.defn
@@ -22,7 +22,6 @@ SRCS = Minkowski.F77			\
        Thorne_fakebinary.F77		\
 					\
        Lemaitre.F77			\
-       Robertson_Walker.F77		\
        de_Sitter.F77			\
        de_Sitter_Lambda.F77		\
        anti_de_Sitter_Lambda.F77	\