From 4d10c7d660916c9acecf050533b003f3fbc492d0 Mon Sep 17 00:00:00 2001 From: guzman Date: Fri, 3 May 2002 14:10:24 +0000 Subject: Converted version for Einstein 2 git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@40 b6f3ac56-194f-0410-8878-cdf6079d7f1b --- src/IDAxiOddBrillBH.F | 31 +++++++++++++++++++++++++++---- 1 file changed, 27 insertions(+), 4 deletions(-) diff --git a/src/IDAxiOddBrillBH.F b/src/IDAxiOddBrillBH.F index c49bf48..02865e3 100644 --- a/src/IDAxiOddBrillBH.F +++ b/src/IDAxiOddBrillBH.F @@ -25,9 +25,6 @@ c c @endhistory c@@*/ -c Need include file from Einstein -#include "CactusEinstein/Einstein/src/Einstein.h" - subroutine IDAxiOddBrillBH(CCTK_ARGUMENTS) implicit none @@ -62,6 +59,7 @@ c Perhaps this and others should go into cctk.h integer :: nx,ny,nz integer :: i,j,k,it,ier,nquads,ntries integer :: npoints,handle,ierror + integer make_conformal_derivs pi = 4.0d0*atan(1.0d0) c DO NOT use integer*4 @@ -85,7 +83,32 @@ c -------------- write(*,*) 'Bowen & York Momenta: :',byJ write(*,*) 'Brill Wave: sin^n theta :',n - conformal_state = CONFORMAL_METRIC +c Check if we should create and store conformal factor stuff + + if (CCTK_EQUALS(metric_type, "static conformal")) then + + conformal_state = 1 + + if(CCTK_EQUALS(conformal_storage,"factor+derivs")) then + + conformal_state = 2 + make_conformal_derivs = 1 + + else if (CCTK_EQUALS(conformal_storage,"factor+derivs+2nd derivs")) then + + conformal_state = 3 + make_conformal_derivs = 1 + + endif + + else + + make_conformal_derivs = 0 + + endif + +c conformal_state = CONFORMAL_METRIC + c Sovle on this sized cartesian grid c 3D grid size NE x NT x NP c Add 2 zones for eta coordinate and 4 for theta -- cgit v1.2.3