[Nek5000-users] Problem with intp_rstd

[nek5000-users at lists.mcs.anl.gov]

Dear all,


I've been playing around with the 'intp_rstd' subroutine to do calculations on the dealiasing grid (lxd*lxy*lzd). I can map the velocity field onto the fine grid (vx --> vxd) and when I print out values it seems fine. Then I try to calculate vxd**2, and still seems fine. However, when I try to map back to the original grid, the results are complete garbage. Here is a snippet of my code.


        do e=1,nelv    ! mapping to fine grid
          call intp_rstd(vxd(1,1,1,e),vx(1,1,1,e),nx1,nxd,if3d,0) ! 0 --> forward
        enddo

          call vsq(vxd,ntotd)  ! calculating vxd**2

        do e=1,nelv    ! mapping back to original grid
          call intp_rstd(vx(1,1,1,e),vxd(1,1,1,e),nx1,nxd,if3d,1) ! 0 --> backwards
        enddo


Am I doing something wrong?


Thanks,


Juan Diego
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170409/995f9cba/attachment.html>