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

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

How big is ldimt in SIZE ?

For the given usage, it looks like it should be 5?

Paul



  Sun, 12 Dec 2010, nek5000-users at lists.mcs.anl.gov wrote:

> 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
>