[gradsusr] g2grb subregion

[Wesley Ebisuzaki]

Mike,

   The error message suggests that you haven't installed all the files
needed
for g2grb.gs to run. You need

 template.grb2 in the $GADDIR directory
 you need wgrib2 (v1.8.4+)   on your path


BTW The simpler way of making regional grib2 subset is to use wgrib2.

  For lat-lon -> subset lat-lon: use the -small_grib option

    wgrib2 IN.grb2 -set_grib_type j  -small_grib 100:170 -55:0 OUT.grb2

You can also interpolate to a lat-lon grid by the -new_grid option.

   Yours,
   Wesley



On Fri, May 18, 2012 at 3:19 AM, Mike Manning <michael at bsch.au.com> wrote:

> Hi all,
>
> I've got data downloaded from an archive now however I'd like to extract
> a subregion from the grib2 files using g2grb.  I've created a test.ctl
> files as follows:
>
> 'open gfs.00z.ctl'
> 'set lat -55 0'
> 'set lon 100 170'
> 'g2grb LFTXsfc extrap.grb UGRD:1000 mb:anl:packing=j:center=7'
> #'g2grb h(lev=925) +extrap.grb HGT:925 mb:anl:packing=j:center=7'
> #'g2grb h(lev=850) +extrap.grb HGT:850 mb:anl:packing=j:center=7'
> #'g2grb h(lev=700) +extrap.grb HGT:700 mb:anl:packing=j:center=7'
> #'g2grb h(lev=600) +extrap.grb HGT:600 mb:anl:packing=j:center=7'
> #'g2grb h(lev=500) +extrap.grb HGT:500 mb:anl:packing=j:center=7'
> quit
>
>
> I run it using grads -blc however it's throwing up some errors:
> Starting "/grads2/Linux/Versions/2.0.a9.oga.1/i686/grads   -blc test.ctl
> " ...
>
>
> Grid Analysis and Display System (GrADS) Version 2.0.a9.oga.1
> Copyright (c) 1988-2010 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> GrADS comes with ABSOLUTELY NO WARRANTY
> See file COPYRIGHT for more information
>
> Config: v2.0.a9.oga.1 little-endian readline printim grib2 netcdf
> hdf4-sds hdf5 opendap-grids,stn athena geotiff shapefile
> Issue 'q config' command for more detailed configuration information
> Loading User Defined Extensions table
> </grads2/Linux/Versions/2.0.a9.oga.1/i686/gex/udxt> ... ok.
> GX Package Initialization: Size = 11 8.5
> Running in Batch mode
> dlon=1 nx = 71 lon0 = 100
> Y is varying   Lat = -55 to 0   Y = 36 to 91
> dlat=1 ny = 56 lat0 = -55
> need to generate date
> new inv=d=2011121000:UGRD:1000 mb:anl:packing=j:center=7
> new grads found undef=-999000000.000000
> making template file
> ! /grads2/wgrib2 $GADDIR/template.grb2 -d 1 -lola 100:71:1 -55:56:1
> extrap.grb.template grib
> argument: -lola ????
> argument: 100:71:1 ????
> argument: -55:56:1 ????
> argument: extrap.grb.template ????
> argument: grib ????
> missing GRIB record(s)
> opening fwrite extrap.grb.tmp
> going to display LFTXsfc
> result is Wrote 3976 of 3976 elements to extrap.grb.tmp as Sequential
> Big_Endian
>
> metadata = 1:0:d=2011121000:UGRD:1000 mb:anl:packing=j:center=7:
> ! /grads2/wgrib2  extrap.grb.template -import_ieee extrap.grb.tmp
> -set_metadata extrap.grb.inv -undefine_val -999000000.000000 -grib_out
> extrap.grb
> argument: -import_ieee ????
> argument: extrap.grb.tmp ????
> argument: -set_grib_type ????
> argument: c3 ????
> argument: -set_metadata ????
> argument: extrap.grb.inv ????
> argument: -undefine_val ????
> argument: -999000000.000000 ????
> argument: -grib_out ????
> argument: extrap.grb ????
> could not open file: extrap.grb.template
>
> !rm extrap.grb.tmp extrap.grb.inv extrap.grb.template
> rm: cannot remove `extrap.grb.template': No such file or directory
>
>
> I'm not sure why it's not working :( All I'd like to do is extract the
> subregion from lat -55 to 0 and lon 100 to 170 for all variables in the
> grib data files I've downloaded.
>
> Cheers, Mike
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120518/233e23dc/attachment-0001.html