[Paraview] Wrong fetched data

[deq11009]

Hi,

I can share the following: http://we.tl/3CPwe5y9vA
An easy way to see the problem might be creating a line source in the 
axis of the cylinder, doing Resample and see where interface is located.

Thanks,
Filipe Direito


On 13.07.2015 04:33, Cory Quammen wrote:
> Interesting. Is is possible for you to share the data file(s)?
>
> Thanks,
> Cory
>
> On Mon, Jul 6, 2015 at 10:23 AM, deq11009 <deq11009 at fe.up.pt [12]>
> wrote:
>
>> Hi,
>>
>> Thank you for your reply.
>>
>> No, the mistake is not there. In the meantime I did realize what is
>> exactly happening.
>>
>> I have, say, 200 timesteps.
>> From 0 to 6, data obtained from my code matches the data when I
>> just use the GUI.
>> Starting in timestep 7 (always) till the end, data (for all
>> timesteps) is the data from the last timestep.
>> So, when in python shell I go to timestep 6, I see timestep 6. But
>> when I go to timestep 7, or 100, or 133, I see timestep 199. At the
>> same time, if I go to the GUI, everything is ok.
>>
>> By coincidence, the ".case" time step data is organized in seven
>> columns. Is this related to the problem?
>>
>> Thank you
>> Filipe Direito
>>
>> On 03.07.2015 15:05, Cory Quammen wrote:
>>
>>> Hi Filipe,
>>>
>>> Perhaps the array at index 7 is not the array you are expecting?
>>> Could
>>> you access the array by name instead, .e.g.,
>>>
>>> dAxis.GetPointData().GetArray(variable_name).GetValue(150)
>>>
>>> and see if that works better?
>>>
>>> Thanks,
>>> Cory
>>>
>>> On Wed, Jul 1, 2015 at 1:44 PM, deq11009 <deq11009 at fe.up.pt [6]
>>> [6]>
>>> wrote:
>>>
>>>> Hi,
>>>>
>>>> I am currently using Paraview 4.2, 64, in Windows. (but tried
>>>> 4.3
>>>> and same problem happened)
>>>>
>>>> My case is Ensight transient. Cylinder with 3.6 million cells,
>>>> two-phase flow.
>>>>
>>>> I created a python code, which I run from the shell in the GUI.
>>>>
>>>> The main parts of the code are:
>>>>
>>>> enscase = EnSightReader(CaseFileName = direcp + ensfinal.case)
>>>> ...
>>>> axis = Line()
>>>> ...
>>>> rAxis = ResampleWithDataset(Input = enscase , Source = axis)
>>>> ...
>>>> rAxis.SMProxy.UpdatePipeline(7) # enscase.SMProxy yields same
>>>> problem
>>>> rAxis.UpdatePipelineInformation()
>>>> dAxis = servermanager.Fetch(rAxis)
>>>>
>>>> Now, for example:
>>>> print dAxis.GetPointData().GetArray(7).GetValue(150)
>>>> which gives me the value 150 of the phase array.
>>>>
>>>> This data is wrong. On the other hand, if I go to the GUI
>>>> Pipeline
>>>> and create a Spreadsheet, the values are the real ones
>>>> (confirmed in
>>>> Fluent).
>>>>
>>>> So, the fetched data in the python shell is, for some reason,
>>>> wrong.
>>>> What is my mistake?
>>>>
>>>> Thanks
>>>> Filipe Direito
>>>>
>>>> _______________________________________________
>>>> Powered by www.kitware.com [1] [1]
>>>>
>>>> Visit other Kitware open-source projects at
>>>> http://www.kitware.com/opensource/opensource.html [2] [2]
>>>>
>>>> Please keep messages on-topic and check the ParaView Wiki at:
>>>> http://paraview.org/Wiki/ParaView [3] [3]
>>>>
>>>> Search the list archives at:
>>>> http://markmail.org/search/?q=ParaView [4]
>>>> [4]
>>>>
>>>> Follow this link to subscribe/unsubscribe:
>>>> http://public.kitware.com/mailman/listinfo/paraview [5] [5]
>>
>> _______________________________________________
>> Powered by www.kitware.com [7]
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html [8]
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView [9]
>>
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>> [10]
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview [11]
>
> --
>
> Cory Quammen
> R&D Engineer
> Kitware, Inc.
>
> Links:
> ------
> [1] http://www.kitware.com
> [2] http://www.kitware.com/opensource/opensource.html
> [3] http://paraview.org/Wiki/ParaView
> [4] http://markmail.org/search/?q=ParaView
> [5] http://public.kitware.com/mailman/listinfo/paraview
> [6] mailto:deq11009 at fe.up.pt
> [7] http://www.kitware.com
> [8] http://www.kitware.com/opensource/opensource.html
> [9] http://paraview.org/Wiki/ParaView
> [10] http://markmail.org/search/?q=ParaView
> [11] http://public.kitware.com/mailman/listinfo/paraview
> [12] mailto:deq11009 at fe.up.pt