From f6290b709fba33db2fc3b4b3e2a73cc7ec9e8432 Mon Sep 17 00:00:00 2001
From: jthorn <jthorn@0a4070d5-58f5-498f-b6c0-2693e757fa0f>
Date: Fri, 8 Jul 2005 10:41:34 +0000
Subject: Add an option to control whether we generate a conformal metric
 (StaticConformal::conformal_state = 3) or a physical metric ( = 0), as per  
 http://www.cactuscode.org/pipermail/developers/2005-July/001178.html and
 subsequent discussion

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@79 0a4070d5-58f5-498f-b6c0-2693e757fa0f
---
 src/IDAxiBrillBH.F | 25 +++++++++++++++++++++----
 1 file changed, 21 insertions(+), 4 deletions(-)

(limited to 'src')

diff --git a/src/IDAxiBrillBH.F b/src/IDAxiBrillBH.F
index f6553cb..93f6db8 100644
--- a/src/IDAxiBrillBH.F
+++ b/src/IDAxiBrillBH.F
@@ -521,6 +521,21 @@ c       dxgxx = 1/2 \partial gxx / \partial x
          enddo
       enddo
 
+c convert to physical metric if StaticConformal isn't wanted
+      if (generate_StaticConformal_metric .eq. 0) then
+         call CCTK_INFO("converting to physical metric")
+         call ConfToPhysInPlace(nx, ny, nz,
+     $				psi,
+     $				gxx, gxy, gxz,
+     $				     gyy, gyz,
+     $				          gzz)
+c        record that we now have a physical metric
+         conformal_state = 0
+      else
+c        record that we computed psi and its 1st and 2nd derivatives
+         conformal_state = 3
+      end if
+
 c Extrinsic Curvature is identically zero
       kxx = 0.0d0
       kxy = 0.0d0
@@ -531,7 +546,12 @@ c Extrinsic Curvature is identically zero
 
       if (debug .ge. 6) then
          print *, '### final results (again at this 3-D grid point) ###'
-	 print *, '   psi =', psi(debug_i,debug_j,debug_k)
+         if (conformal_state .gt. 0) then
+            print *, '### ... conformal metric with'
+	    print *, '   psi =', psi(debug_i,debug_j,debug_k)
+         else
+            print *, '### ... physical metric with'
+         end if
 	 print *, '   gxx =', gxx(debug_i,debug_j,debug_k)
 	 print *, '   gxy =', gxy(debug_i,debug_j,debug_k)
 	 print *, '   gxz =', gxz(debug_i,debug_j,debug_k)
@@ -562,9 +582,6 @@ c ADM mass
          alp = (2.0d0*r - adm)/(2.0d0*r+adm)
       endif
 
-c conformal_state = 'all' derivatives were calculated
-      conformal_state = 3
-
       deallocate(cc,ce,cw,cn,cs,rhs,psi2d,detapsi2d,dqpsi2d,
      $     detaetapsi2d,detaqpsi2d,dqqpsi2d,
      $     etagrd,qgrd,
-- 
cgit v1.2.3