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 /src/Stab3d.F | |
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
Diffstat (limited to 'src/Stab3d.F')
-rw-r--r-- | src/Stab3d.F | 9 |
1 files changed, 7 insertions, 2 deletions
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 |