<div dir="ltr"><div><div>Initially, I did try the CDO on my local linux box but I got<br><br>%> cdo -f nc copy gonzalo08l.2014101412.hwrfsat.global.0p25.f024.grb2 gon_CDO.nc<br><br>Error (gribapiScanTimestep1) : GRIB_API support not compiled in!<br><br>---<br></div>I guess my local version needs an upgrade.<br><br>---<br></div><div>I did try the above on another machine and got<br><br>/glade/p/ncldev/data/grib2>ncdump -h gon_CDO.nc<br><br>netcdf gon_CDO {<br>dimensions:<br> lon = 1440 ;<br> lat = 721 ;<br> time = UNLIMITED ; // (18 currently)<br>variables:<br> double lon(lon) ;<br> lon:standard_name = "longitude" ;<br> lon:long_name = "longitude" ;<br> lon:units = "degrees_east" ;<br> lon:axis = "X" ;<br> double lat(lat) ;<br> lat:standard_name = "latitude" ;<br> lat:long_name = "latitude" ;<br> lat:units = "degrees_north" ;<br> lat:axis = "Y" ;<br> double time(time) ;<br> time:standard_name = "time" ;<br> time:units = "hours since 2014-10-14 12:00:00" ;<br> time:calendar = "proleptic_gregorian" ;<br> float param7.5.0(time, lat, lon) ;<br> param7.5.0:_FillValue = -9.e+33f ;<br> param7.5.0:missing_value = -9.e+33f ;<br><br>// global attributes:<br> :CDI = "Climate Data Interface version 1.6.3 (<a href="http://code.zmaw.de/projects/cdi">http://code.zmaw.de/projects/cdi</a>)" ;<br> :Conventions = "CF-1.4" ;<br> :history = "Wed Dec 28 12:04:14 2016: cdo -f nc copy gonzalo08l.2014101412.hwrfsat.global.0p25.f024.grb2 gon_CDO.nc" ;<br> :institution = "National Centers for Environmental Prediction" ;<br> :CDO = "Climate Data Operators version 1.6.3 (<a href="http://code.zmaw.de/projects/cdo">http://code.zmaw.de/projects/cdo</a>)" ;<br>}<br><br><br></div><div>The CDO returns the variable with a nonsensical name: <br><br> float param7.5.0(time, lat, lon) ;<br><br></div><div>but no units or long_name .... How do you know what it is?<br></div><div>I am sure the numbers are unpacked correctly ... but ... without using some external information, how do you know to what 'param7.5.0' refers?<br><br>===<br></div><div>Combining the CDO file with NCL naming (via an appropriate internal lookup table) you would get:<br><br> float BRTEMP_P32_GLL0 ( lat_0, lon_0 )<br> long_name : Brightness temperature<br> units : K<br><br></div><div>param7.5.0 ===> BRTEMP_P32_GLL0<br><br></div><div>and the attributes should be:<br><br> long_name : Brightness temperature<br> units : K<br></div><div><br></div>D<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 28, 2016 at 11:04 AM, Ligia Bernardet - NOAA Affiliate <span dir="ltr"><<a href="mailto:ligia.bernardet@noaa.gov" target="_blank">ligia.bernardet@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dennis,<div><br></div><div>Thanks for your willingness to help.</div><div><br></div><div>I was able to use CDO to convert the GRIB2 file to NetCDF. All 18 fields are in the output file (interestingly, it appears as one variable at 18 time levels, which is definitely not the case in the GRIB file ). I was able to open the resulting NetCDF file in NCL and manipulate it. So, I am on my way to get the results I need.</div><div><br></div><div>But if you guys get this working in NCL directly, let me know.</div><div><br></div><div>Thanks again,</div><div>Ligia</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 28, 2016 at 8:07 AM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Hi Ligea<br><br></div><div>FYI ... a tool unrelated to NCL also fails.<br></div><div><br></div>------<br></div>There is another tool available from NCEP: wgrib2<br><br><a href="http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/" target="_blank">http://www.cpc.ncep.noaa.gov/p<wbr>roducts/wesley/wgrib2/</a><br><br></div>There is an option to create netCDF<br><br><a href="http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/netcdf.html" target="_blank">http://www.cpc.ncep.noaa.gov/p<wbr>roducts/wesley/wgrib2/netcdf.h<wbr>tml</a><br><br>===<br></div>When I use this tool on the grib2 file you sent, it also fails:<br><br>%> wgrib2 gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024.grb2 -netcdf <a href="http://gon.nc" target="_blank">gon.nc</a><br><br>---<br>The wgrib2 tool prints many warning lines like:<br><br>code_table_4.5a: product definition template #32 not supported<br>code_table_4.5b: product definition template #32 not supported<br>....<br>code_table_4.5b: product definition template #32 not supported<br>1:0:d=<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>:var discipline=0 master_table=8 parmcat=5 parm=7:no_level:24 hour fcst:<br>code_table_4.5a: product definition template #32 not supported<br>...<br>code_table_4.5b: product definition template #32 not supported<br>2:279767:d=<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>:var discipline=0 master_table=8 parmcat=5 parm=7:no_level:24 hour fcst:<br>code_table_4.5a: product definition template #32 not supported<br>...<br><br></div>etc, etc<br><br></div>Like 'ncl_convert2nc', wgrib2 produces just one variable within the output netCDF. Specifically,<br><br><br>%> ncddump -h <a href="http://gon.nc" target="_blank">gon.nc</a> <br><br>netcdf gon {<br>dimensions:<br> latitude = 721 ;<br> longitude = 1440 ;<br> time = UNLIMITED ; // (1 currently)<br>variables:<br> double latitude(latitude) ;<br> latitude:units = "degrees_north" ;<br> latitude:long_name = "latitude" ;<br> double longitude(longitude) ;<br> longitude:units = "degrees_east" ;<br> longitude:long_name = "longitude" ;<br> double time(time) ;<br> time:units = "seconds since 1970-01-01 00:00:00.0 0:00" ;<br> time:long_name = "verification time generated by wgrib2 function verftime()" ;<br> time:reference_time = 1413288000. ;<br> time:reference_time_type = 3 ;<br> time:reference_date = "2014.10.14 12:00:00 UTC" ;<br> time:reference_time_descriptio<wbr>n = "forecast or accumulated, reference date is fixed" ;<br> time:time_step_setting = "auto" ;<br> time:time_step = 0. ;<br> float var0_5_7_no_level(time, latitude, longitude) ;<br> var0_5_7_no_level:_FillValue = 9.999e+20f ;<br> var0_5_7_no_level:short_name = "var0_5_7_no_level" ;<br> var0_5_7_no_level:long_name = "desc" ;<br> var0_5_7_no_level:level = "no_level" ;<br> var0_5_7_no_level:units = "unit" ;<br><br>// global attributes:<br> :Conventions = "COARDS" ;<br> :History = "created by wgrib2" ;<br> :GRIB2_grid_template = 0 ;<br>==============================<wbr>============<br><br></div><div><br><div><div><br><div><div><br></div></div></div></div></div></div><div class="m_3750456871960329632HOEnZb"><div class="m_3750456871960329632h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 27, 2016 at 5:01 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>6.3.0 will unpack one variable also. The name of the variable will be different because (! wild guess !) the needed table was not complete or available.<br><br></div>GRIB is a bit of a 'moving target' :-(<br><br></div>I have opened a JIRA ticket: NCL-2533<br><br></div>Unfortunately, this is not a great week to get things done.<br><br></div>Maybe, an IO developer will look?<span class="m_3750456871960329632m_6689671540169569900HOEnZb"><font color="#888888"><br><br></font></span></div><span class="m_3750456871960329632m_6689671540169569900HOEnZb"><font color="#888888">D<br><div><div><div><div><br><br><div><br><br></div></div></div></div></div></font></span></div><div class="m_3750456871960329632m_6689671540169569900HOEnZb"><div class="m_3750456871960329632m_6689671540169569900h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 27, 2016 at 4:48 PM, Ligia Bernardet - NOAA Affiliate <span dir="ltr"><<a href="mailto:ligia.bernardet@noaa.gov" target="_blank">ligia.bernardet@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Dennis,<div><br></div><div>There are 18 fields in the GRIB2 file. In the NetCDF file you sent back to me there is 1.</div><div>But more importantly, I need to find a way to do this myself, so I can process the multiple files I have.</div><div><br></div><div>The machine I am using, a NOAA's R&D platform, I was able to find NCL v6.3. Would this version be usable to read this GRIB2 file, until a newer version is available?</div><div><br></div><div>Thanks,</div><div>Ligia</div><div><br></div><div><br></div></div><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174HOEnZb"><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 27, 2016 at 4:11 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Attached is a gzip'd version of the netCDF file.<br><br>The 'soon-to-be-released' beta 6.4.0 returns a warning message. I will open a JIRA ticket.<br></div><div>Note: I believe the numbers are correct.<br></div><br><div>=========<br><br> NCAR Command Language Version 6.4.0 -27Dec2016_0148<br><br>warning:NclGRIB2: BRTEMP_P32 contains records that NCL cannot currently differentiate. One or more records will be ignored.<br><br>Variable: f<br>Type: file<br>filename: gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024<br>path: /Users/shea/Data/GRIB/gonzalo0<wbr>8l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.global.0<wbr>p25.f024.grb2<br> file global attributes:<br> dimensions:<br> lat_0 = 721<br> lon_0 = 1440<br> variables:<br> float BRTEMP_P32_GLL0 ( lat_0, lon_0 )<br> center : US National Weather Service - NCEP (WMC)<br> production_status : Operational products<br> long_name : Brightness temperature<br> units : K<br> _FillValue : 1e+20<br> grid_type : Latitude/longitude<br> parameter_discipline_and_categ<wbr>ory : Meteorological products, Long wave radiation<br> parameter_template_discipline_<wbr>category_number : ( 32, 0, 5, 7 )<br> forecast_time : 0<br> forecast_time_units : hours<br> initial_time : 10/14/2014 (12:00)<br><br> float lat_0 ( lat_0 )<br> long_name : latitude<br> grid_type : Latitude/Longitude<br> units : degrees_north<br> Dj : 0.25<br> Di : 0.25<br> Lo2 : 359.75<br> La2 : -90<br> Lo1 : 0<br> La1 : 90<br><br> float lon_0 ( lon_0 )<br> long_name : longitude<br> grid_type : Latitude/Longitude<br> units : degrees_east<br> Dj : 0.25<br> Di : 0.25<br> Lo2 : 359.75<br> La2 : -90<br> Lo1 : 0<br> La1 : 90<br><br>=========<br></div><div>The min and max are:<br></div><div><br> (0) Brightness temperature (K) : min=184.67 max=319.49<br><br></div><div><br></div></div><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777HOEnZb"><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 27, 2016 at 1:13 PM, Ligia Bernardet - NOAA Affiliate <span dir="ltr"><<a href="mailto:ligia.bernardet@noaa.gov" target="_blank">ligia.bernardet@noaa.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<div><br>I have transferred file gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024.grb2 to /ftp/incoming.<div><br></div><div>Thanks,</div><div>Ligia</div></div></div><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777m_7880590755414206150HOEnZb"><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777m_7880590755414206150h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Dec 27, 2016 at 1:06 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>[1]<br></div>6.1.2 (Feb. 2013) is almost 4 years old. Much has happened since.<br><div><br>[2]<br>To be able to determine the problem or if the current version of NCL has addressed the issue, we would need the file.<br><br></div><div>ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a><br></div><div>anonymous<br></div><div>your_email<br></div><div>cd incoming<br></div><div>prompt<br></div><div>put gonzalo08l.<a value="+12014101412" href="tel:%28201%29%20410-1412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024.grb2<br></div><div>quit<br><br></div><div>We can not look at /ftp/incoming, so ....<br><br>Let ncl-talk know when the file has been transferred.<br></div><div><br></div><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777m_7880590755414206150m_-3211259656709017677h5">On Tue, Dec 27, 2016 at 12:56 PM, Ligia Bernardet - NOAA Affiliate <span dir="ltr"><<a href="mailto:ligia.bernardet@noaa.gov" target="_blank">ligia.bernardet@noaa.gov</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_3750456871960329632m_6689671540169569900m_-3633709906522825174m_7277089200559379777m_7880590755414206150m_-3211259656709017677h5"><div dir="ltr">Hello,<div><br></div><div>I would like to use convert a GRIB2 file generated by NCEP's Hurricane WRF (HWRF) model to a NetCDF file. However, I am encountering the error below. Any suggestions?</div><div><br></div><div>Thanks,</div><div>Ligia<br><br><b>ncl_convert2nc gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024.grb2<br></b><br></div><div>Processing file: gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfsat.<wbr>global.0p25.f024.grb2...<br>Copyright (C) 1995-2013 - All Rights Reserved<br>University Corporation for Atmospheric Research<br>NCAR Command Language Version 6.1.2<br>The use of this software is governed by a License Agreement.<br>See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.</div><div>getpdstemplate: PDS Template 4.32 not defined.<br>warning:Error unpacking GRIB record section 4.<br>fatal:Could not open (./gonzalo08l.<a href="tel:(201)%20410-1412" value="+12014101412" target="_blank">2014101412</a>.hwrfs<wbr>at.global.0p25.f024.grb2)<br>Segmentation fault</div></div>
<br></div></div>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>