diff options
author | ryoji <ryoji@971fb155-194f-0410-9daf-e2eca44e59f5> | 2004-11-01 19:41:21 +0000 |
---|---|---|
committer | ryoji <ryoji@971fb155-194f-0410-9daf-e2eca44e59f5> | 2004-11-01 19:41:21 +0000 |
commit | b970f9b97da03a29dc904ba48579ccf9775d941d (patch) | |
tree | b6401c32533cff8a2d5ba6be9d2f2da3d7516f51 | |
parent | ff09ff17f7918ef0d1aa5a201002bd6b331cd0c8 (diff) |
it's same changes as rotating ... stuff.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/DistortedBHIVP/trunk@34 971fb155-194f-0410-9daf-e2eca44e59f5
-rw-r--r-- | src/DistortedBHIVP.F | 2 | ||||
-rw-r--r-- | src/Stab3d.F | 9 | ||||
-rwxr-xr-x | test/test_dbh.par | 2 |
3 files changed, 9 insertions, 4 deletions
diff --git a/src/DistortedBHIVP.F b/src/DistortedBHIVP.F index b76b184..4913781 100644 --- a/src/DistortedBHIVP.F +++ b/src/DistortedBHIVP.F @@ -620,7 +620,7 @@ c Compute eta,q,phi at the each points of cartesian grid interp_coords(1) = CCTK_PointerTo(abseta) interp_coords(2) = CCTK_PointerTo(q) - interp_coords(2) = CCTK_PointerTo(phi) + interp_coords(3) = CCTK_PointerTo(phi) in_array_dims(1) = ne in_array_dims(2) = nq diff --git a/src/Stab3d.F b/src/Stab3d.F index f8ca5e0..aad61b6 100644 --- a/src/Stab3d.F +++ b/src/Stab3d.F @@ -97,12 +97,14 @@ c real*8 x(im*jm*km),r(im*jm*km) c Local variables integer :: i,j,k,kk - real*8 :: p(im*jm*km),Ap(im*jm*km),w(im*jm*km),As(im*jm*km) + real*8, allocatable :: p(:),Ap(:),w(:),As(:) real*8 :: omega, chi,chi1,chi2, delta, deltap, pp * *********************************************************************** * - +c + allocate(p(im*jm*km),Ap(im*jm*km),w(im*jm*km),As(im*jm*km)) + do i = 1,im*jm*km p(i) = 0. Ap(i) = 0. @@ -237,6 +239,9 @@ c rnorm = sum(r**2) rnorm=sqrt(rnorm) if (rnorm .gt. tol) goto 1 +c + deallocate(p,Ap,w,As) +c return end c diff --git a/test/test_dbh.par b/test/test_dbh.par index cdd0a8a..acb04f0 100755 --- a/test/test_dbh.par +++ b/test/test_dbh.par @@ -2,7 +2,7 @@ # single rotating black hole #============================================================================== -ActiveThorns = " DistortedBHIVP adm coordbase coordgauge staticconformal spacemask admbase admcoupling admmacros admanalysis ellbase cartgrid3d time pugh pughreduce pughslab ioascii ioutil iobasic adm admconstraints boundary pughinterp LocalInterp flexio nanchecker " +ActiveThorns = " DistortedBHIVP adm coordbase coordgauge staticconformal spacemask admbase admcoupling admmacros admanalysis ellbase cartgrid3d time pugh pughreduce pughslab ioascii ioutil iobasic adm admconstraints symbase boundary pughinterp LocalInterp flexio nanchecker " IOBasic::outScalar_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift admconstraints::momy admconstraints::ham " IOBasic::outInfo_vars = " admconstraints::ham admbase::alp " IOASCII::out1D_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift admconstraints::momy admconstraints::ham " |