[ncl-talk] Issue reading grib in 6.3.0?
Wei Huang
huangwei at ucar.edu
Thu Apr 9 11:36:52 MDT 2015
We have fixed this problem in our source code.
The fix will be in the next release.
Thanks,
Wei
================================================
1850 Table Mesa Dr.
Boulder, CO 80307
Phone: 303-497-8924
On Thu, Apr 9, 2015 at 11:30 AM, David Brown <dbrown at ucar.edu> wrote:
> 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
>>
>>
>
> _______________________________________________
> 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/bdb70c3c/attachment.html
More information about the ncl-talk
mailing list