aboutsummaryrefslogtreecommitdiff
path: root/src/metric.F
diff options
context:
space:
mode:
authorjthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87>2002-11-17 14:55:50 +0000
committerjthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87>2002-11-17 14:55:50 +0000
commitac099068da91d9528088af849caa7aa2e6e1758a (patch)
treef8837bb58d475eb0492f4a6a604a0d51bba94423 /src/metric.F
parentd72eb352c6182603bde73125f71c4313f22cea28 (diff)
Modified Files:
make.code.defn Added Files: Bona_Masso_data.F77 blended_boundary.F77 boundary.F77 decode_pars.F77 gauge.F77 initialize.F77 metric.F77 xyz_blended_boundary.F77 Removed Files: Bona_Masso_data.F blended_boundary.F boundary.F decode_pars.F gauge.F initialize.F metric.F xyz_blended_boundary.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. Also small bugfix: blended_boundary.F --> .F77 xyz_blended_boundary.F --> .F77 * change declaration of ierr from CCTK_REAL --> integer git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@135 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src/metric.F')
-rw-r--r--src/metric.F293
1 files changed, 0 insertions, 293 deletions
diff --git a/src/metric.F b/src/metric.F
deleted file mode 100644
index 23797f7..0000000
--- a/src/metric.F
+++ /dev/null
@@ -1,293 +0,0 @@
-c This subroutine calculates the 4-metric and its inverse at an event,
-c by decoding decoded_exact_model and calling the appropriate subroutine
-C $Header$
-
-#include "cctk.h"
-#include "cctk_Parameters.h"
-#include "cctk_Arguments.h"
-
-#include "param_defs.inc"
-
- subroutine Exact__metric(
- $ decoded_exact_model,
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz,rama)
-
- implicit none
- DECLARE_CCTK_FUNCTIONS
-
-c arguments
- CCTK_INT decoded_exact_model
- CCTK_REAL x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz, rama
-
-c local variables
- character*100 warn_buffer
-
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c Minkowski spacetime
-c
-
- if (decoded_exact_model .eq. EXACT__Minkowski) then
- call Exact__Minkowski(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Minkowski_shift) then
- call Exact__Minkowski_shift(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Minkowski_funny) then
- call Exact__Minkowski_funny(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Minkowski_gauge_wave) then
- call Exact__Minkowski_gauge_wave(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c black hole spacetimes
-c
-
- elseif (decoded_exact_model .eq. EXACT__Schwarzschild_EF) then
- call Exact__Schwarzschild_EF(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Schwarzschild_PG) then
- call Exact__Schwarzschild_PG(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Schwarzschild_Novikov) then
- call Exact__Schwarzschild_Novikov(x,y,z,t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Kerr_BoyerLindquist) then
- call Exact__Kerr_BoyerLindquist(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Kerr_KerrSchild) then
- call Exact__Kerr_KerrSchild(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Schwarzschild_Lemaitre) then
- call Exact__Schwarzschild_Lemaitre(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__multi_BH) then
- call Exact__multi_BH(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
-c
-c not fully implemented yet -- see Nina Jansen for details
-c
-c elseif (decoded_exact_model .eq. EXACT__Alvi) then
-c call Exact__Alvi(
-c $ x, y, z, t,
-c $ gdtt, gdtx, gdty, gdtz,
-c $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
-c $ gutt, gutx, guty, gutz,
-c $ guxx, guyy, guzz, guxy, guyz, guxz)
-c
-
- elseif (decoded_exact_model .eq. EXACT__Thorne_fakebinary) then
- call Exact__Thorne_fakebinary(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c cosmological spacetimes
-c
-
- elseif (decoded_exact_model .eq. EXACT__Lemaitre) then
- call Exact__Lemaitre(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- 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,rama)
-
- elseif (decoded_exact_model .eq. EXACT__de_Sitter) then
- call Exact__de_Sitter(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__de_Sitter_Lambda) then
- call Exact__de_Sitter_Lambda(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__anti_de_Sitter_Lambda) then
- call Exact__anti_de_Sitter_Lambda(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Bianchi_I) then
- call Exact__Bianchi_I(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Goedel) then
- call Exact__Goedel(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Bertotti) then
- call Exact__Bertotti(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Kasner_like) then
- call Exact__Kasner_like(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Kasner_axisymmetric) then
- call Exact__Kasner_axisymmetric(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Kasner_generalized) then
- call Exact__Kasner_generalized(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__Milne) then
- call Exact__Milne(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-c
-c miscellaneous spacetimes
-c
-
- elseif (decoded_exact_model .eq. EXACT__boost_rotation_symmetric) then
- call Exact__boost_rotation_symmetric(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__bowl) then
- call Exact__bowl(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
- elseif (decoded_exact_model .eq. EXACT__constant_density_star) then
- call Exact__constant_density_star(
- $ x, y, z, t,
- $ gdtt, gdtx, gdty, gdtz,
- $ gdxx, gdyy, gdzz, gdxy, gdyz, gdxz,
- $ gutt, gutx, guty, gutz,
- $ guxx, guyy, guzz, guxy, guyz, guxz)
-
-cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
- else
- write (warn_buffer, '(a,i8)')
- $ 'Unknown decoded_exact_model = ', decoded_exact_model
- call CCTK_WARN(0, warn_buffer)
- endif
-
- return
- end