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

6 files changed, 55 insertions, 57 deletions

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	\