[Paraview] ParaView not reading time in netcdf file

Dan Lipsa dan.lipsa at kitware.com
Thu Nov 9 12:01:17 EST 2017


David,
I filed the following bug report

https://gitlab.kitware.com/paraview/paraview/issues/17827

Thanks for reporting this.

Best,
Dan


On Thu, Nov 9, 2017 at 11:36 AM, David Deepwell <ddeepwell at uwaterloo.ca>
wrote:

> Yay, we found a bug!
> Yes, you can attach one (or both) of the files in a bug report.
>
> Thanks for your help,
> David
>
>
> On Nov 9, 2017, at 10:57 AM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>
> I was just looking at that. Indeed, your variable did not depend on time.
> So it seems that the reader does report the correct information:
>
> The following script:
> import vtk
> reader = vtk.vtkNetCDFCFReader()
> reader.SetFileName('output_20.nc')
> reader.UpdateInformation()
> print reader.GetOutputInformation(0).Get(vtk.
> vtkStreamingDemandDrivenPipeline.TIME_RANGE())
> reader.Update()
> output = reader.GetOutput()
> print output.GetClassName()
>
> prints:
> (20.0, 20.0)
> vtkImageData
>
> This is a ParaView bug. Is it OK to attach one of your files to the bug
> report?
> If not, I can file the report without it.
>
> Thanks,
> Dan
>
>
>
> On Thu, Nov 9, 2017 at 10:36 AM, David Deepwell <ddeepwell at uwaterloo.ca>
> wrote:
>
>> I think I was missing an important declaration in the creation of the
>> netcdf file I gave you, so I’m not surprised that the time didn’t work for
>> that file. I’ve fixed that, but the time is still not shown for a single
>> file. Use these instead of the previous file.
>>
>> Cheers,
>> David
>>
>>
>> On Nov 9, 2017, at 9:50 AM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> David, I was thinking this problem comes from the way ParaView loads time
>> dependent datasets. However, I created a copy of the file with _51.nc
>> postfix and ParaView still does not show me the correct time value.
>> Can you send me two files where the time value is shown correctly?
>>
>> Thanks,
>> Dan
>>
>>
>>
>> On Fri, Nov 3, 2017 at 2:39 PM, David Deepwell <ddeepwell at uwaterloo.ca>
>> wrote:
>>
>> I’ve made another discovery.
>> The correct time will be displayed if paraview loads multiple files (or a
>> single file containing multiple times). Not sure why it displays zero for a
>> single file containing a single time.
>>
>> David
>>
>>
>> On Oct 31, 2017, at 3:50 PM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> Hi David,
>> Thanks for your patience. I did look at the vtk reader a little. Did not
>> get to the bottom of it yet. I am trying to determine if the reader has
>> this info and ParaView does not use it, or the reader is the problem.
>> Do check back with me.
>>
>> Thanks,
>> Dan
>>
>>
>> On Tue, Oct 31, 2017 at 2:55 PM, David Deepwell <ddeepwell at uwaterloo.ca>
>> wrote:
>>
>> I’m sure you are busy, but just checking in to see if there is any
>> progress on reading the time.
>>
>> Cheers,
>> David
>>
>>
>> On Oct 23, 2017, at 9:54 AM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> Hi David,
>> I'll try to take a closer look at this today, at least to see if the
>> problem is in the reader or in ParaView.
>>
>> Thanks for your patience,
>> Dan
>>
>>
>> On Mon, Oct 23, 2017 at 9:48 AM, Deepwell, David <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Sorry to bother you again, but have you found anything?
>>
>> David
>>
>>
>> On Oct 17, 2017, at 4:34 PM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> Actually, you are right, it does not work. I'll try to look closer at
>> this in the next couple of days.
>>
>> On Tue, Oct 17, 2017 at 4:11 PM, Deepwell, David <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Hi Dan,
>>
>> That’s strange. I’m also using 5.4.1. Did you get a time of 5 s?
>>
>> Cheers,
>> David
>>
>>
>> On Oct 17, 2017, at 3:59 PM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> David,
>> What version of ParaView are you using? I tried your file with the latest
>> release 5.4.1 and it works fine.
>>
>> Thanks,
>> Dan
>>
>>
>> On Mon, Oct 16, 2017 at 10:18 AM, Deepwell, David <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Ok, no problem.  Just checking in.
>>
>> Thanks,
>> David
>>
>>
>> On Oct 16, 2017, at 9:29 AM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> Sorry, no, I have not had a chance to look at this. I will do that in the
>> next couple of days.
>>
>> Thanks,
>> Dan
>>
>>
>> On Mon, Oct 16, 2017 at 9:05 AM, Deepwell, David <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Hi Dan,
>> Have you been able to see my problem, and find a solution?
>>
>> David
>>
>>
>> On Oct 10, 2017, at 5:46 PM, David Deepwell <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Hi Dan!
>>
>> Thank you for taking a look at my file.
>>
>> David
>> <output_50.nc>
>>
>> On Oct 6, 2017, at 1:38 PM, Dan Lipsa <dan.lipsa at kitware.com> wrote:
>>
>> Hi David,
>> If you send me the file I can take a look.
>>
>> Thanks,
>> Dan
>>
>>
>> On Fri, Oct 6, 2017 at 10:11 AM, Deepwell, David <ddeepwel at uwaterloo.ca>
>> wrote:
>>
>> Hi,
>>
>> I have a netcdf file which I’d like to read in ParaView. The reading and
>> rendering works fine, but the time is not displayed. I’ve declared the
>> variable to be time through the axis attribute, and I’ve given it units
>> since a certain day. Strangely enough, VisIt understands the time value
>> (t=5s) while ParaView does not.
>>
>> Here’s the output from ncdump:
>>
>> netcdf output_50 {
>> dimensions:
>>         x = 1024 ;
>>         z = 128 ;
>>         time = 1 ;
>> variables:
>>         double x(x) ;
>>                 x:units = "m" ;
>>         double z(z) ;
>>                 z:units = "m" ;
>>         double time(time) ;
>>                 time:units = "seconds since 2001-1-1" ;
>>                 time:axis = "T" ;
>>                 time:long_name = "time" ;
>>         double rho(time, x, z) ;
>>                 rho:units = "kg/L" ;
>>
>> // global attributes:
>>                 :description = "SPINS Output 50" ;
>> }
>>
>> Any help will be appreciated. I also have a 1 MB file which I can pass on
>> if needed.
>>
>> Thank you,
>> David
>> _______________________________________________
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Please keep messages on-topic and check the ParaView Wiki at:
>> http://paraview.org/Wiki/ParaView
>>
>> Search the list archives at: http://markmail.org/search/?q=ParaView
>>
>> Follow this link to subscribe/unsubscribe:
>> http://public.kitware.com/mailman/listinfo/paraview
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20171109/435ba10c/attachment.html>


More information about the ParaView mailing list