subroutine dealias1 (u,u_big)
! puts an array into the 'big' size
! doesn't dealias anything
! note: this sort of trashes u
use types,only:rprec
use param,only:ld,ld_big,nx,ny,nz,ny2,dx,dy
use fft
implicit none
integer::jz
real(kind=rprec),dimension(ld,ny,nz),intent(in):: u
real(kind=rprec),dimension(ld_big,ny2,nz),intent(inout)::u_big
real(kind=rprec)::ignore_me,const
! ahh get rid of this
real(kind=rprec),dimension(ld,ny,nz)::temp    
! be careful using u after calling this subroutine!
const=1._rprec/(nx*ny)
temp=const*u

do jz=1,nz
! still need to normalize
   call rfftwnd_f77_one_real_to_complex(forw,temp(:,:,jz),ignore_me)
! check padd syntax
   call padd(u_big(:,:,jz),temp(:,:,jz))
   call rfftwnd_f77_one_complex_to_real(back_big,u_big(:,:,jz),ignore_me)
end do
end subroutine dealias1
