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

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

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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20101213/51f87ff1/attachment.html>