[ncl-talk] Issue reading grib in 6.3.0?
David Brown
dbrown at ucar.edu
Thu Apr 9 11:30:35 MDT 2015
I will confirm that this change was not intended and is not part of any
change to the actual GRIB handling code. The code that determines whether a
file is GRIB1 or GRIB2 is unchanged. The issue that needs fixing is simply
the recognition of the alternate suffixes.
-dave
On Thu, Apr 9, 2015 at 11:03 AM, Dennis Shea <shea at ucar.edu> wrote:
> One of the NCL core developers will likely respond.
>
> To my knowledge, all pre-6.3.0 GRIB extensions *should still be
> recognized*. Backward compatibility is important.
>
> --
> That said .... .grib, gr1, gr2 do produces a failure in 6.3.0. I will open
> a JIRA ticket.
>
> NCL-2197 - GRIB extensions not recognized in 6.3.0
> <https://vets.development.ucar.edu/jira/browse/NCL-2197>
>
> --
> Also, I see in my original reply to CS that I should have included the
> following,
>
> The following illustrate how to force a 'time dimension' when importing
> variables from a file with just on time step. By default, NCL does not
> import singleton (aka, degenerate) dimensions. Hence, no 'time' dimension.
> Users can force one or more 'time' dimensions via:
>
> %> ncl_filedump -itime nam.t06z.awip1203.grb
>
> [snip]
> path: nam.t06z.awip1203
> file global attributes:
> dimensions:
> initial_time0_hours = 1
> [snip]
> variables:
> float TMP_P0_L1_GLC0 ( initial_time0_hours, ygrid_0, xgrid_0 )
> [SNIP]
>
> ===
> %> ncl_filedump -ftime nam.t06z.awip1203.grb
>
> [snip]
> path: nam.t06z.awip1203
> file global attributes:
> dimensions:
> forecast_time0 = 1
> [snip]
> variables:
> float TMP_P0_L1_GLC0 ( forecast_time0, ygrid_0, xgrid_0 )
> [SNIP]
>
> ===
> %> ncl_filedump -ftime -itime nam.t06z.awip1203.grb
>
> [snip]
> dimensions:
> initial_time0_hours = 1
> forecast_time0 = 1
> [snip]
>
> float TMP_P0_L1_GLC0 ( initial_time0_hours, forecast_time0, ygrid_0,
> xgrid_0 )
> [SNIP]
>
> *************
> From a script
>
> setfileoption("grb","SingleElementDimensions","Initial_time") ;
> initial_time0_hours
> and/or
> setfileoption("grb","SingleElementDimensions","Forecast_time")
>
>
>
> On Thu, Apr 9, 2015 at 8:40 AM, Alan Brammer <abrammer at albany.edu> wrote:
>
>> As an ignorant user of various types of grib without worrying about which
>> type it is. Does this sentence no longer apply then, from the addfile doc
>> page.
>> " *Also note, in the case of GRIB, that you do not need to indicate
>> whether the file is GRIB 1 or GRIB 2. NCL treats all the acceptable GRIB
>> extensions as meaning the same thing. The GRIB version is determined
>> internally.*"
>>
>> Is there a way to determine which extension is needed then? If this isn't
>> a bug, this seems like a step backwards in my uneducated opinion.
>>
>> On 9 April 2015 at 10:27, Carl Schreck <cjschrec at ncsu.edu> wrote:
>>
>>> Aha! I was adding ".grib", which worked in the previous version, but not
>>> now. The addfile documentation
>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>lists
>>> these as possible grib extensions:
>>> ".gr", ".gr1", ".grb", ".grib", ".grb1", ".grib1", ".gr2", "grb2",
>>> ".grib2"
>>>
>>> I get errors in 6.3.0 if I use ".gr1:, ".gr2", or ".grib" for this
>>> particular file. I see there were significant changes to the handling of
>>> grib files, but I can't tell whether this change was either intentional or
>>> documented.
>>>
>>> Carl
>>>
>>> On Thu, Apr 9, 2015 at 10:11 AM, Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> Hi Carl,
>>>>
>>>> In the following, the file name is: nam.t06z.awip1203
>>>>
>>>> You must add the '.grb' extension
>>>>
>>>> %> ncl_filedump nam.t06z.awip1203.grb
>>>>
>>>> NCL will look for a file named "nam.t06z.awip1203.grb" on the disk.
>>>> If it does not find it, NCL will look for "nam.t06z.awip1203" and
>>>> 'remember' to treat it as a grib file.
>>>>
>>>> One file
>>>>
>>>> f = addfile("nam.t06z.awip1203.grb","r")
>>>>
>>>>
>>>> For multiple files
>>>>
>>>>
>>>> diri = "./"
>>>> fili = systemfunc("cd "+diri+" ; ls nam.t06z.awip*")
>>>> f = addfiles(diri+fili+".grb")
>>>> ListSetType (f, "join")
>>>>
>>>> HTH
>>>> D
>>>>
>>>> Copyright (C) 1995-2015 - All Rights Reserved
>>>> University Corporation for Atmospheric Research
>>>> NCAR Command Language Version 6.3.0
>>>> The use of this software is governed by a License Agreement.
>>>> See http://www.ncl.ucar.edu/ for more details.
>>>>
>>>> Variable: f
>>>> Type: file
>>>> filename: nam.t06z.awip1203
>>>> path: nam.t06z.awip1203
>>>> file global attributes:
>>>> dimensions:
>>>> ygrid_0 = 428
>>>> xgrid_0 = 614
>>>> lv_ISBL0 = 4
>>>> lv_SIGL1 = 5
>>>> lv_SPDL2 = 6
>>>> lv_HTGL3 = 2
>>>> lv_SPDL4 = 5
>>>> lv_ISBL5 = 2
>>>> lv_ISBL6 = 5
>>>> lv_SPDL7 = 3
>>>> lv_HTGL8 = 2
>>>> lv_HTGL9 = 2
>>>> lv_DBLL10 = 4
>>>> variables:
>>>> float TMP_P0_L1_GLC0 ( ygrid_0, xgrid_0 )
>>>> center : US National Weather Service - NCEP (WMC)
>>>> production_status : Operational products
>>>> long_name : Temperature
>>>> units : K
>>>> _FillValue : 1e+20
>>>> coordinates : gridlat_0 gridlon_0
>>>> grid_type : Lambert Conformal can be secant or tangent,
>>>> conical or bipolar
>>>> parameter_discipline_and_category : Meteorological
>>>> products, Temperature
>>>> parameter_template_discipline_category_number : ( 0, 0,
>>>> 0, 0 )
>>>> level_type : Ground or water surface
>>>> level : 0
>>>> forecast_time : 3
>>>> forecast_time_units : hours
>>>> initial_time : 04/09/2015 (06:00)
>>>>
>>>> [SNIPu
>>>>
>>>>
>>>>
>>>>
>>>> On Thu, Apr 9, 2015 at 7:54 AM, Carl Schreck <cjschrec at ncsu.edu> wrote:
>>>>
>>>>> I'm having a strange issue where NAM forecast grib files from ldm are
>>>>> not able to be read in 6.3.0. A sample file is:
>>>>> ftp://filsrv.cicsnc.org/carl/temp/nam.t06z.awip1203
>>>>>
>>>>> The error I get is:
>>>>> --
>>>>> fatal:["FileSupport.c":4159]:(/home/carl/ftp/temp/nam.t06z.awip1203.grib)
>>>>> has no file extension, can't determine type of file to open
>>>>>
>>>>> Variable: f
>>>>> Type: file
>>>>> (0) File Missing Value : -1
>>>>> --
>>>>> Renaming the file to have a .grib extension makes no difference. It
>>>>> does work if I roll back to 6.2.1, however.
>>>>>
>>>>> The binaries I downloaded are:
>>>>> ncl_ncarg-6.2.1.Linux_SL6.5_x86_64_gcc447.tar.gz
>>>>> ncl_ncarg-6.3.0.Linux_RHEL6.4_x86_64_gcc447.tar.gz
>>>>>
>>>>> Thanks!
>>>>>
>>>>> Carl
>>>>> --
>>>>>
>>>>> <http://www.cicsnc.org/>
>>>>> <http://www.researcherid.com/rid/B-8711-2011>
>>>>> * Cyclone <http://www.cyclonecenter.org/>*Center.org
>>>>> <http://www.cyclonecenter.org/> *Carl J. Schreck III, PhD*
>>>>> *Research Associate*
>>>>> North Carolina Institute for Climate Studies (NCICS) /
>>>>> Cooperative Institute for Climate and Satellites NC (CICS-NC)
>>>>> <http://www.cicsnc.org/>
>>>>> North Carolina State University <http://ncsu.edu/>
>>>>> NOAA National Centers for Environmental Information (NCEI)
>>>>> <http://ncdc.noaa.gov/>
>>>>> 151 Patton Ave, Asheville, NC 28801
>>>>> e: cjschrec at ncsu.edu
>>>>> o: +1 828 257 3140
>>>>> Publications
>>>>> <http://scholar.google.com/citations?hl=en&user=th8ONEcAAAAJ&view_op=list_works&sortby=pubdate>
>>>>> monitor.cicsnc.org/mjo
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>
>>>
>>>
>>> --
>>>
>>> <http://www.cicsnc.org/>
>>> <http://www.researcherid.com/rid/B-8711-2011>
>>> * Cyclone <http://www.cyclonecenter.org/>*Center.org
>>> <http://www.cyclonecenter.org/> *Carl J. Schreck III, PhD*
>>> *Research Associate*
>>> North Carolina Institute for Climate Studies (NCICS) /
>>> Cooperative Institute for Climate and Satellites NC (CICS-NC)
>>> <http://www.cicsnc.org/>
>>> North Carolina State University <http://ncsu.edu/>
>>> NOAA National Centers for Environmental Information (NCEI)
>>> <http://ncdc.noaa.gov/>
>>> 151 Patton Ave, Asheville, NC 28801
>>> e: cjschrec at ncsu.edu
>>> o: +1 828 257 3140
>>> Publications
>>> <http://scholar.google.com/citations?hl=en&user=th8ONEcAAAAJ&view_op=list_works&sortby=pubdate>
>>> monitor.cicsnc.org/mjo
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150409/f376e2cc/attachment.html
More information about the ncl-talk
mailing list