Modified Files:

	make.code.defn
Added Files:
	Alvi.F77 Bertotti.F77 Bianchi_I.F77 Goedel.F77
	Kasner_axisymmetric.F77 Kasner_generalized.F77 Kasner_like.F77
	Kerr_BoyerLindquist.F77 Kerr_KerrSchild.F77 Lemaitre.F77
	Milne.F77 Minkowski.F77 Minkowski_funny.F77
	Minkowski_gauge_wave.F77 Minkowski_shift.F77
	Robertson_Walker.F77 Schwarzschild_EF.F77
	Schwarzschild_Lemaitre.F77 Schwarzschild_Novikov.F77
	Schwarzschild_PG.F77 Thorne_fakebinary.F77
	anti_de_Sitter_Lambda.F77 boost_rotation_symmetric.F77
	bowl.F77 constant_density_star.F77 de_Sitter.F77
	de_Sitter_Lambda.F77 multi_BH.F77
Removed Files:
	Alvi.F Bertotti.F Bianchi_I.F Goedel.F Kasner_axisymmetric.F
	Kasner_generalized.F Kasner_like.F Kerr_BoyerLindquist.F
	Kerr_KerrSchild.F Lemaitre.F Milne.F Minkowski.F
	Minkowski_funny.F Minkowski_gauge_wave.F Minkowski_shift.F
	Robertson_Walker.F Schwarzschild_EF.F Schwarzschild_Lemaitre.F
	Schwarzschild_Novikov.F Schwarzschild_PG.F Thorne_fakebinary.F
	anti_de_Sitter_Lambda.F boost_rotation_symmetric.F bowl.F
	constant_density_star.F de_Sitter.F de_Sitter_Lambda.F
	multi_BH.F

Rename files which really are Fortran 77
from foo.F (= Cactus Fortran 90 fixed form)
to foo.F77 (= Fortran 77)

This means they're now compiled with a Fortran 77 compiler.
This should make no difference to the semantics (they were already
Fortran 77 code), but makes it easier to compile this thorn on
platforms which have a Fortran 77 compiler but no Fortran 90 compiler.


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

1 files changed, 72 insertions, 0 deletions

diff --git a/src/metrics/Schwarzschild_Lemaitre.F77 b/src/metrics/Schwarzschild_Lemaitre.F77
new file mode 100644
index 0000000..9ff7835
--- /dev/null
+++ b/src/metrics/Schwarzschild_Lemaitre.F77
@@ -0,0 +1,72 @@
+C  Schwarzschild-Lemaitre metric
+c  (Schwarzschild black hole with cosmological constant)
+C  Proposed by Lemaitre in 1932 
+C Author : D. Vulcanov (Timisoara, Romania)
+C $Header$
+
+#include "cctk.h"
+#include "cctk_Parameters.h"
+
+      subroutine Exact__Schwarzschild_Lemaitre(
+     $     x, y, z, t,
+     $     gdtt, gdtx, gdty, gdtz, 
+     $     gdxx, gdyy, gdzz, gdxy, gdyz, gdzx,
+     $     gutt, gutx, guty, gutz, 
+     $     guxx, guyy, guzz, guxy, guyz, guzx)
+
+      implicit none
+
+      DECLARE_CCTK_PARAMETERS
+
+      CCTK_REAL x, y, z, t
+      CCTK_REAL gdtt, gdtx, gdty, gdtz, 
+     $       gdxx, gdyy, gdzz, gdxy, gdyz, gdzx,
+     $       gutt, gutx, guty, gutz, 
+     $       guxx, guyy, guzz, guxy, guyz, guzx
+
+      logical firstcall
+    
+      CCTK_REAL lam, mas, r2, ppp, unu, doi
+
+      data firstcall /.true./
+      save firstcall
+      if (firstcall) then
+          lam = Schwarzschild_Lemaitre__Lambda
+          mas = Schwarzschild_Lemaitre__mass
+          firstcall = .false.
+      end if
+ 
+      r2 =x*x+y*y+z*z
+
+      ppp=1.0D0 -2.0D0*mas/dsqrt(r2) -r2*lam/3.0D0
+
+      unu=(1.0D0-ppp)/ppp/r2      
+ 
+      doi=(ppp - 1.0D0)/r2
+
+      gdtt = -ppp
+      gdtx = 0.0D0 
+      gdty = 0.0D0
+      gdtz = 0.0D0
+      gdxx = 1.0D0 + x*x*unu
+      gdyy = 1.0D0 + y*y*unu
+      gdzz = 1.0D0 + z*z*unu
+      gdxy =         x*y*unu
+      gdyz =         y*z*unu
+      gdzx =         z*x*unu
+
+
+      gutt = -1.0D0/ppp
+      gutx = 0.0D0
+      guty = 0.0D0
+      gutz = 0.0D0
+      guxx = 1.0D0 + x*x*doi
+      guyy = 1.0D0 + y*y*doi
+      guzz = 1.0D0 + z*z*doi
+      guxy =         x*y*doi
+      guyz =         y*z*doi
+      guzx =         x*z*doi
+
+
+      return
+      end