<div dir="ltr"><font color="#212121" size="2">Hi Cynthia, </font><div><font color="#212121" size="2"><br></font></div><div><font color="#212121" size="2">If you could modify the </font><span style="color:rgb(33,33,33)">vtkNetCDFCFReader (</span><font color="#212121"><a href="https://github.com/Kitware/VTK/blob/master/IO/NetCDF/vtkNetCDFCFReader.cxx">https://github.com/Kitware/VTK/blob/master/IO/NetCDF/vtkNetCDFCFReader.cxx</a>)</font><span style="color:rgb(33,33,33)"> to work with WRF formats, that would be an awesome addition. I know few folks that can use this work. Unfortunately we do not have any direct funding to accomplish this but we will be more than happy to guide or help a little if needed. </span></div><div><span style="color:rgb(33,33,33)"><br></span></div><div><span style="color:rgb(33,33,33)">Let us know and we can discuss more offline. </span></div><div><span style="color:rgb(33,33,33)"><br></span></div><div><span style="color:rgb(33,33,33)">Thanks, </span></div><div><span style="color:rgb(33,33,33)"><br></span></div><br><div class="gmail_quote"><div dir="ltr">On Tue, Aug 8, 2017 at 4:40 PM Dan Lipsa <<a href="mailto:dan.lipsa@kitware.com">dan.lipsa@kitware.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Cynthia,<div>It is possible to change the reader. You'll need to compile ParaView from source. The class you need to change is vtkNetCDFCFReader.</div><div>Probably you'll need an extra parameter to the reader that would be the time unit. The time unit read from the file should overides this parameter.</div><div><br></div><div>You could also submit your changes back to VTK (most readers in ParaView are actually part of VTK). This seems like something other people might need.</div></div><div dir="ltr"><div><br></div><div>Dan</div><div><br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 8, 2017 at 12:43 PM, Hart, Cynthia Joyce - (cyndyb) <span dir="ltr"><<a href="mailto:cyndyb@email.arizona.edu" target="_blank">cyndyb@email.arizona.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="white" lang="EN-US" link="blue" vlink="purple">
<div class="m_7649014869735212880m_-3292657193733308466WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Hello,
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">To make a long story short, this solution you used below, does not work well for the weather forecasting group here on my campus that I am working with. Is there any possibility
of obtaining the source code of the Paraview netCDF reader? If I could edit that source code to accept the non- CF compliant WRF outputs, I could make use of Paraview in an automated script to create animations for our forecasters much, much more readily
than manipulating netCDF file headers. The data I’d sent to you was but a very small subset of the hundreds of netCDF files we’d need to change, every time we do a model run, and this is a process we’d like to actually automate daily. I am unable to force
the output to be compliant as it is created. Can you help? <u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Sincerely,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">Cynthia Hart<u></u><u></u></span></p><span>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"><u></u> <u></u></span></p>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-family:"Calibri",sans-serif;color:black">From:
</span></b><span style="font-family:"Calibri",sans-serif;color:black">Dan Lipsa <<a href="mailto:dan.lipsa@kitware.com" target="_blank">dan.lipsa@kitware.com</a>><br>
<b>Date: </b>Friday, June 23, 2017 at 10:57 AM<br>
<b>To: </b>Aashish Chaudhary <<a href="mailto:aashish.chaudhary@kitware.com" target="_blank">aashish.chaudhary@kitware.com</a>><br>
<b>Cc: </b>"Hart, Cynthia Joyce - (cyndyb)" <<a href="mailto:cyndyb@email.arizona.edu" target="_blank">cyndyb@email.arizona.edu</a>>, "<a href="mailto:paraview@paraview.org" target="_blank">paraview@paraview.org</a>" <<a href="mailto:paraview@paraview.org" target="_blank">paraview@paraview.org</a>><br>
<b>Subject: </b>Re: [Paraview] troubles with netCDF reader ?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</span><div><span class="m_7649014869735212880HOEnZb"><font color="#888888">
<p class="MsoNormal">Cynthia, <u></u><u></u></p></font></span><div><div class="m_7649014869735212880h5">
<div>
<p class="MsoNormal">It seems that the files you send us are not CF compliant.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">The problem is that the Time dimension does not have a unit.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Here are some quotes from the CF Conventions document:<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">- The use of coordinate variables is required for all dimensions that correspond to one dimensional space or time coordinates.<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal">- Variables representing time must always explicitly include the units attribute; there is no default value<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">So, to fix one of your files, I did:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">1. ncdump ARAD4_test2.ncdf > ARAD4_test2_CF.txt<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">2. I added a new variable Time with the unit specified to the txt file.<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><br>
<br>
<span class="m_7649014869735212880m_-3292657193733308466gmail-apple-tab-span"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal">float Time(Time) ;<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Time:units = "minutes since 2016-07-29 12:00:00" ;<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">3.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">I generated a new nc file:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">ncgen ARAD4_test2_CF.txt -o ARAD4_test2_CF.nc<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">This new file is read correctly by ParaView.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Best,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Dan<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div><div><div class="m_7649014869735212880h5">
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Fri, Jun 16, 2017 at 6:31 PM, Aashish Chaudhary <<a href="mailto:aashish.chaudhary@kitware.com" target="_blank">aashish.chaudhary@kitware.com</a>> wrote:<u></u><u></u></p>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">Dear Cynthia, <u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Is there any chance that you can send me a sample dataset? <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">- Aashish<u></u><u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal">On Fri, Jun 16, 2017 at 6:15 PM Hart, Cynthia Joyce - (cyndyb) <<a href="mailto:cyndyb@email.arizona.edu" target="_blank">cyndyb@email.arizona.edu</a>> wrote:<u></u><u></u></p>
</div>
</div>
</div>
<blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">Folks –
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I am struggling with using the netCDF reader, with a bunch of timeseries of WRF data output (weather prediction modeling). I am a visualizer, NOT
a meteorologist, and so, my full understanding of the model’s output is shaky at best. The researcher I am trying to assist is currently on vacation, and I’m trying to animate the output of his WRF model.
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">So! I have been viewing some of the output using the built in reader in Paraview (the ‘generic and cf’ reader) and it works for some variables, not
for others. I am using NCO tools to extract some of the variables of interest into their own CDF file, to try and isolate the problem Paraview has, reading some variables.. but that doesn’t help. I am particularly trying to read a variable for radar data,
“REFL_10CM”. Here is the relevant bit of the file header from ncdump, below: </span>
<u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">netcdf ARAD4_test2 {</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">dimensions:</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> Time = UNLIMITED ; // (1 currently)</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> south_north = 393 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> west_east = 639 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> bottom_top = 37 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">variables:</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> float REFD_MAX(Time, south_north, west_east) ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:FieldType = 104 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:MemoryOrder = "XY " ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:description = "MAX DERIVED RADAR REFL" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:units = "dbZ" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:stagger = "" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFD_MAX:coordinates = "XLONG XLAT XTIME" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> float REFL_10CM(Time, bottom_top, south_north, west_east) ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:FieldType = 104 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:MemoryOrder = "XYZ" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:description = "Radar reflectivity (lamda = 10 cm)" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:units = "dBZ" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:stagger = "" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> REFL_10CM:coordinates = "XLONG XLAT XTIME" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> float XLAT(Time, south_north, west_east) ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:FieldType = 104 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:MemoryOrder = "XY " ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:units = "degree_north" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:stagger = "" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLAT:coordinates = "XLONG XLAT" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> float XLONG(Time, south_north, west_east) ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:FieldType = 104 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:MemoryOrder = "XY " ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:units = "degree_east" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:stagger = "" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XLONG:coordinates = "XLONG XLAT" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> float XTIME(Time) ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XTIME:FieldType = 104 ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XTIME:MemoryOrder = "0 " ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XTIME:description = "minutes since 2016-07-29 12:00:00" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XTIME:units = "minutes since 2016-07-29 12:00:00" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> XTIME:stagger = "" ;</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Paraview complains that it cannot read a variable with more than three dimensions without time. But there appears to ‘be’ time (I’m very unclear
what “UNLIMITED” time , as seen in the header, refers to in this context). Is there any way I can work around this? Most frustrating. I can actually read this variable into QGIS, but I’d rather use the strengths of Paraview to visualize this data.
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt">Thanks in advance.
</span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:11.0pt"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal">_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">
http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">
http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" target="_blank">
http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><u></u><u></u></p>
</blockquote>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">
http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">
http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" target="_blank">
http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
</blockquote></div><br></div>
</blockquote></div></div>