From b970f9b97da03a29dc904ba48579ccf9775d941d Mon Sep 17 00:00:00 2001 From: ryoji Date: Mon, 1 Nov 2004 19:41:21 +0000 Subject: it's same changes as rotating ... stuff. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/DistortedBHIVP/trunk@34 971fb155-194f-0410-9daf-e2eca44e59f5 --- src/DistortedBHIVP.F | 2 +- src/Stab3d.F | 9 +++++++-- 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 " -- cgit v1.2.3