[R-sig-Geo] Working with GRIB2 files

[Matthew Landis]

Ezra,

I've been messing around with GRIB2 myself, but as of July,  I hadn't 
had any luck finding ways to do it in R.  I ended up using NCL (NCAR 
command language) to convert the Grib2's to NetCDF, which various R 
packages can handle (notably raster).  I did figure out how to pass 
arguments from R to the NCL script (through Bash), which was a major 
achievement for me.   This only works on Linux though, as I haven't 
worked it out for the Cygwin version of NCL.

I would *love* to hear about ways of dealing with GRIB2 natively in R.  
Barring that, I'm happy to share NCL code.

Matt

On 9/30/2011 12:41 PM, Ezra Boyd wrote:
> Hello,
>
> I am trying to use the Real Time Ocean Forecasting System
> (http://polar.ncep.noaa.gov/ofs/download.shtml?) archival data which
> is in the GRIB2 format using a hybrid model-based coordinate system.
> The data can be downloaded from here:
> http://data.nodc.noaa.gov/opendap/ncep/rtofs/contents.html
>
> Has anyone had any experience using GRIB2 in R?  If so, please let me
> know how it worked out.  In particular, I need to get forecast data
> into a format where I can overlay a shapefile, but the wqgrib2 program
> won't work because hybrid coordinate system.
>
> Thanks,
> Ezra

-- 
~~~~~~~~~~~~~~~~~~~~~~~~~~
Matthew Landis, Ph.D.
Research Scientist
ISciences, LLC
61 Main St. Suite 200
Burlington VT 05401
802.864.2999
~~~~~~~~~~~~~~~~~~~~~~~~~~