diff --git a/lib/four2a.f90 b/lib/four2a.f90 index e7e0e8018..a63f31a4d 100644 --- a/lib/four2a.f90 +++ b/lib/four2a.f90 @@ -67,7 +67,7 @@ subroutine four2a(a,nfft,ndim,isign,iform) if(nfft.le.NSMALL) then jz=nfft - if(iform.eq.0) jz=nfft/2 + if(iform.le.0) jz=nfft/2 aa(1:jz)=a(1:jz) endif @@ -87,7 +87,7 @@ subroutine four2a(a,nfft,ndim,isign,iform) if(nfft.le.NSMALL) then jz=nfft - if(iform.eq.0) jz=nfft/2 + if(iform.le.0) jz=nfft/2 a(1:jz)=aa(1:jz) endif end if