Fix compiler warnings.

Most of them could be fixed by renaming .F77 files to .F
Some had to be fixed by explicitly declaring some variables using
CCTK_DECLARE() (which also only works for .F, not for .F77)



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

1 files changed, 80 insertions, 0 deletions

diff --git a/src/metrics/Schwarzschild_Lemaitre.F b/src/metrics/Schwarzschild_Lemaitre.F
new file mode 100644
index 0000000..aff86a9
--- /dev/null
+++ b/src/metrics/Schwarzschild_Lemaitre.F
@@ -0,0 +1,80 @@
+C  Schwarzschild-Lemaitre metric
+c  (Schwarzschild black hole with cosmological constant)
+C  Proposed by Lemaitre in 1932 
+C
+C Author : D. Vulcanov (Timisoara, Romania)
+C see ../../README for copyright & licensing info
+C
+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,
+     $     psi, Tmunu_flag)
+
+      implicit none
+
+      DECLARE_CCTK_PARAMETERS
+
+c input arguments
+      CCTK_REAL x, y, z
+      CCTK_DECLARE(CCTK_REAL, t,)
+
+c output arguments
+      CCTK_REAL gdtt, gdtx, gdty, gdtz, 
+     $       gdxx, gdyy, gdzz, gdxy, gdyz, gdzx,
+     $       gutt, gutx, guty, gutz, 
+     $       guxx, guyy, guzz, guxy, guyz, guzx
+      CCTK_DECLARE(CCTK_REAL, psi,)
+      LOGICAL   Tmunu_flag
+
+c local variables
+      CCTK_REAL lam, mas
+      CCTK_REAL r2, ppp, unu, doi
+
+C This model has a cosmological constant
+      Tmunu_flag = .true.
+
+      lam = Schwarzschild_Lemaitre__Lambda
+      mas = Schwarzschild_Lemaitre__mass
+ 
+      r2 =x*x+y*y+z*z
+
+      ppp=1.0D0 -2.0D0*mas/sqrt(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