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

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

Hi Paul,

We had discussed this earlier and I had to fix the restart file to get it to
work. Now it doesn't seem to come from the restart file though. I am running
it on 128 procs (~24k elements) on our cluster. I can try running it on a
different machine here to see if I get the same problem. Do you have any
instructions on how to upload to the ftp site ?

Regards
Shriram Jagannathan


On 13 December 2010 08:10, <nek5000-users at lists.mcs.anl.gov> wrote:

>
> 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
>>>
>>>
>>  _______________________________________________
> 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/c757884c/attachment.html>