[Nek5000-users] Data Overwritten while using comp_vort3() and hpts()

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

Hello,

I am having an issue with hpts data being overwritten while using the
comp_vort3() and the hpts() routine. However the data dumps (blah0.f0000)
are uncorrupted. Eddy viscosity is computed and stored in temp array,
lambda2 and vorticity components are stored in 4 passive scalar arrays. It
seems like the last two PS values of previous history points are
over-written on vx and vy. I have attached an image of the hpts.out file and
also pasted below is a piece of the usr file :

      common /gaaa/    wo1(lx1,ly1,lz1,lelv)
     &              ,  wo2(lx1,ly1,lz1,lelv)
     &              ,  wo3(lx1,ly1,lz1,lelv)

      if(mod(istep,iostep).eq.0) then
        ! compute lambda2 vortex
        call lambda2(t(1,1,1,1,2))
        ifpsco(1)=.true.
        ! compute vorticity
        if(ldimt.ge.5) call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)
        call hpts()
       ifpsco(2)=.true.
       ifpsco(3)=.true.
       ifpsco(4)=.true.
      endif

t(1,1,1,1,1) has eddy viscosity and ldimt is set to 5. Also, I am wondering
what ldimt1 and ldimt3 means in the SIZE file ? Thanks.

Regards
Shriram
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101212/98816247/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hpts.png
Type: image/png
Size: 58034 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101212/98816247/attachment.png>