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

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

Hi Shriram,

I can't recall if I checked into this before - this might
be one of those problems that doesn't show up on our platform,
but I'm willing to give it a shot if you can tar up the
requisite files and put on our ftp site.   How many processors
were you using?

Paul


On Mon, 13 Dec 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hi Paul,
>
> Yes. ldimt is set to 5.
>
> Regards
> Shriram Jagannathan
>
>
> On 13 December 2010 02:00, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>>
>>
>> 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
>>>
>>>  _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>
>