diff options
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 |