diff options
| author | ryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602> | 2000-02-14 15:26:38 +0000 |
|---|---|---|
| committer | ryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602> | 2000-02-14 15:26:38 +0000 |
| commit | 84bbb7b952d76e567eff4b764fe0e16672418268 (patch) | |
| tree | 41cd95f954080f46229cffa555d762eb4fa087ab /src | |
| parent | f3c90d16c930c2fd714f70968c7e8a3e58a0fd78 (diff) | |
it passed testsuite on t3e.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@11 535fb057-194f-0410-b5a5-c63992f15602
Diffstat (limited to 'src')
| -rw-r--r-- | src/RotatingDBHIVP.F | 6 | ||||
| -rw-r--r-- | src/Stab3d.F | 6 |
2 files changed, 9 insertions, 3 deletions
diff --git a/src/RotatingDBHIVP.F b/src/RotatingDBHIVP.F index e1dfc2b..b65fffc 100644 --- a/src/RotatingDBHIVP.F +++ b/src/RotatingDBHIVP.F @@ -69,6 +69,7 @@ c Need include file from Einstein real*8,parameter :: rbh_tol = 1.0d-7,rbh_eps = 1.0d-10 integer,parameter :: itmax = 30 real*8 pi + integer :: ne,nq,np integer :: nx,ny,nz integer i,j,k,it,ier,nquads,nocts,order,ntries integer npoints,handle,ierror @@ -77,6 +78,11 @@ c Need include file from Einstein pi = 4.0d0*atan(1.0d0) +c DON'T use integer*4 + ne = neta + nq = ntheta + np = nphi + c Set up the grid spacings nx = cctk_lsh(1) ny = cctk_lsh(2) diff --git a/src/Stab3d.F b/src/Stab3d.F index 9cb01ff..0413c60 100644 --- a/src/Stab3d.F +++ b/src/Stab3d.F @@ -6,7 +6,7 @@ c been made. c implicit none c - integer*4,intent(in) :: im,jm,km + integer,intent(in) :: im,jm,km real*8,intent(inout) :: cc(im,jm,km),cn(im,jm,km),cs(im,jm,km), $ ce(im,jm,km),cw(im,jm,km),ct(im,jm,km),cb(im,jm,km) real*8,intent(out) :: eps @@ -86,7 +86,7 @@ c been made. c implicit none c - integer*4,intent(in) :: im,jm,km + integer,intent(in) :: im,jm,km integer,intent(in) :: ncyc real*8,intent(out) :: cc(im*jm*km),cn(im*jm*km) real*8,intent(out) :: cs(im*jm*km),ce(im*jm*km) @@ -251,7 +251,7 @@ c Be careful that the cs are zero on their outer boundary!! c implicit none c - integer*4,intent(in) :: im,jm,km + integer,intent(in) :: im,jm,km real*8,intent(out) :: cc(im*jm*km),cn(im*jm*km) real*8,intent(out) :: cs(im*jm*km),ce(im*jm*km) real*8,intent(out) :: cw(im*jm*km),ct(im*jm*km) |