diff options
author | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-11-17 14:55:50 +0000 |
---|---|---|
committer | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-11-17 14:55:50 +0000 |
commit | ac099068da91d9528088af849caa7aa2e6e1758a (patch) | |
tree | f8837bb58d475eb0492f4a6a604a0d51bba94423 /src/metric.F | |
parent | d72eb352c6182603bde73125f71c4313f22cea28 (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.F | 293 |
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 |