diff options
| author | allen <allen@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2002-05-17 15:36:33 +0000 |
|---|---|---|
| committer | allen <allen@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2002-05-17 15:36:33 +0000 |
| commit | 8ab40af4cb476c6e5060989fa15619a92657ec7b (patch) | |
| tree | 2d1ab29bbf554a56b50d5ae0f9aa859e18d6aa93 /src/setupbrilldata3D.F | |
| parent | 8b4534e2a895baa38abfa539b6e34a498e6fe741 (diff) | |
Tidying so that the conformal solver isn't used to solve for the conformal
factor, in fact the conformal factor is now only set at the end of
the thorn.
Now if the conformal metric is used, psi will be set to the conformal factor
which was actually solved for. This means that right now that you can't use IDBrillData with conformal derivatives, although it would presumably be
easy enough to numerically calculate them, or to return 1 and 0's instead
just for this case.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@76 a678b1cf-93e1-4b43-a69d-d43939e66649
Diffstat (limited to 'src/setupbrilldata3D.F')
| -rw-r--r-- | src/setupbrilldata3D.F | 32 |
1 files changed, 0 insertions, 32 deletions
diff --git a/src/setupbrilldata3D.F b/src/setupbrilldata3D.F index 9daa65f..ecd0fd3 100644 --- a/src/setupbrilldata3D.F +++ b/src/setupbrilldata3D.F @@ -78,38 +78,6 @@ c Initialize psi. brillpsi = one -c Set up conformal metric. - - if (CCTK_EQUALS(metric_type,"static conformal")) then - - conformal_state = 1 - do k=1,nz - do j=1,ny - do i=1,nx - psi(i,j,k) = one - if (CCTK_EQUALS(conformal_storage,"factor+derivs") .or. - & CCTK_EQUALS(conformal_storage,"factor+derivs+2nd derivs")) then - psix(i,j,k) = zero - psiy(i,j,k) = zero - psiz(i,j,k) = zero - conformal_state = 2 - end if - if (CCTK_EQUALS(conformal_storage,"factor+derivs+2nd derivs")) then - 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 - conformal_state = 3 - end if - end do - end do - end do - - end if - - do k=1,nz do j=1,ny do i=1,nx |