{x}(1) -> {x}(*), where {x} = r, c, Wa4.

Author: zoziha

src/dfftb.f90

@@ -3,7 +3,7 @@ subroutine dfftb(n,r,Wsave)
       implicit none
       integer :: n
       real(rk) :: r , Wsave
-      dimension r(1) , Wsave(*)
+      dimension r(*) , Wsave(*)
       if ( n==1 ) return
       call rfftb1(n,r,Wsave,Wsave(n+1),Wsave(2*n+1))
       end subroutine dfftb

src/dfftf.f90

@@ -3,7 +3,7 @@ subroutine dfftf(n,r,Wsave)
       implicit none
       integer :: n
       real(rk) :: r , Wsave
-      dimension r(1) , Wsave(*)
+      dimension r(*) , Wsave(*)
       if ( n==1 ) return
       call rfftf1(n,r,Wsave,Wsave(n+1),Wsave(2*n+1))
       end subroutine dfftf

src/passb5.f90

@@ -8,7 +8,7 @@ subroutine passb5(Ido,l1,Cc,Ch,Wa1,Wa2,Wa3,Wa4)
                   tr4 , tr5 , Wa1 , Wa2 , Wa3 , Wa4
       integer :: i , Ido , k , l1
       dimension Cc(Ido,5,l1) , Ch(Ido,l1,5) , Wa1(*) , Wa2(*) , Wa3(*), &
-                Wa4(1)
+                Wa4(*)
       real(rk),parameter :: pi = acos(-1.0_rk)
       real(rk),parameter :: tr11 = cos(2.0_rk * pi / 5.0_rk)
       real(rk),parameter :: ti11 = sin(2.0_rk * pi / 5.0_rk)

src/passf5.f90

@@ -8,7 +8,7 @@ subroutine passf5(Ido,l1,Cc,Ch,Wa1,Wa2,Wa3,Wa4)
                   tr4 , tr5 , Wa1 , Wa2 , Wa3 , Wa4
       integer :: i , Ido , k , l1
       dimension Cc(Ido,5,l1) , Ch(Ido,l1,5) , Wa1(*) , Wa2(*) , Wa3(*), &
-                Wa4(1)
+                Wa4(*)
       real(rk),parameter :: pi = acos(-1.0_rk)
       real(rk),parameter :: tr11 =  cos(2.0_rk * pi / 5.0_rk)
       real(rk),parameter :: ti11 = -sin(2.0_rk * pi / 5.0_rk)

src/radb5.f90

@@ -8,7 +8,7 @@ subroutine radb5(Ido,l1,Cc,Ch,Wa1,Wa2,Wa3,Wa4)
                   tr4 , tr5 , Wa1 , Wa2 , Wa3 , Wa4
       integer :: i , ic , Ido , idp2 , k , l1
       dimension Cc(Ido,5,l1) , Ch(Ido,l1,5) , Wa1(*) , Wa2(*) , Wa3(*), &
-                Wa4(1)
+                Wa4(*)
       real(rk),parameter :: pi = acos(-1.0_rk)
       real(rk),parameter :: tr11 = cos(2.0_rk * pi / 5.0_rk)
       real(rk),parameter :: ti11 = sin(2.0_rk * pi / 5.0_rk)

src/radf5.f90

@@ -8,7 +8,7 @@ subroutine radf5(Ido,l1,Cc,Ch,Wa1,Wa2,Wa3,Wa4)
                   tr4 , tr5 , Wa1 , Wa2 , Wa3 , Wa4
       integer :: i , ic , Ido , idp2 , k , l1
       dimension Cc(Ido,l1,5) , Ch(Ido,5,l1) , Wa1(*) , Wa2(*) , Wa3(*), &
-                Wa4(1)
+                Wa4(*)
       real(rk),parameter :: pi = acos(-1.0_rk)
       real(rk),parameter :: tr11 = cos(2.0_rk * pi / 5.0_rk)
       real(rk),parameter :: ti11 = sin(2.0_rk * pi / 5.0_rk)

src/zfftb.f90

@@ -3,7 +3,7 @@ subroutine zfftb(n,c,Wsave)
       implicit none
       real(rk) :: c , Wsave
       integer :: iw1 , iw2 , n
-      dimension c(1) , Wsave(*)
+      dimension c(*) , Wsave(*)
       if ( n==1 ) return
       iw1 = n + n + 1
       iw2 = iw1 + n + n

src/zfftf.f90

@@ -3,7 +3,7 @@ subroutine zfftf(n,c,Wsave)
       implicit none
       real(rk) :: c , Wsave
       integer :: iw1 , iw2 , n
-      dimension c(1) , Wsave(*)
+      dimension c(*) , Wsave(*)
       if ( n==1 ) return
       iw1 = n + n + 1
       iw2 = iw1 + n + n