diff options
| author | allen <allen@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2002-01-06 18:16:09 +0000 |
|---|---|---|
| committer | allen <allen@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2002-01-06 18:16:09 +0000 |
| commit | 71b5bbb5148076bd8a9b39305763136c707ecfe4 (patch) | |
| tree | 222722c382bb15adbbb2feff647b73484c86b889 | |
| parent | e3663f19fdc0c19e3739e2ad1d65d9fa3f1b30fa (diff) | |
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
| -rw-r--r-- | src/setupbrilldata3D.F | 27 |
1 files changed, 26 insertions, 1 deletions
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 |