From 71b5bbb5148076bd8a9b39305763136c707ecfe4 Mon Sep 17 00:00:00 2001 From: allen Date: Sun, 6 Jan 2002 18:16:09 +0000 Subject: Initialise derivatives for conformal metric and set conformal_state grid variable git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@62 a678b1cf-93e1-4b43-a69d-d43939e66649 --- src/setupbrilldata3D.F | 27 ++++++++++++++++++++++++++- 1 file changed, 26 insertions(+), 1 deletion(-) diff --git a/src/setupbrilldata3D.F b/src/setupbrilldata3D.F index edeb21a..f4550d5 100644 --- a/src/setupbrilldata3D.F +++ b/src/setupbrilldata3D.F @@ -82,7 +82,32 @@ c Initialize psi. c Set up conformal metric. - psi = one + if (use_conformal == 1) then + + conformal_state = CONFORMAL_METRIC + do k=1,nz + do j=1,ny + do i=1,nx + psi(i,j,k) = one + psix(i,j,k) = zero + psiy(i,j,k) = zero + psiz(i,j,k) = zero + psixx(i,j,k) = zero + psiyy(i,j,k) = zero + psizz(i,j,k) = zero + psixy(i,j,k) = zero + psixz(i,j,k) = zero + psiyz(i,j,k) = zero + end do + end do + end do + + else + + conformal_state = NOCONFORMAL_METRIC + + end if + do k=1,nz do j=1,ny -- cgit v1.2.3