[gradsusr] g2grb subregion
Mike Manning
michael at bsch.au.com
Sun May 20 19:03:52 EDT 2012
Thanks Wesley,
I'll give that a try this week :)
Cheers, Mike
On 19/05/2012 1:07 AM, Wesley Ebisuzaki wrote:
> Mike,
>
> The error message suggests that you haven't installed all the files
> needed
> for g2grb.gs <http://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
> <mailto: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 <mailto:gradsusr at gradsusr.org>
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
>
> _______________________________________________
> 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/20120521/21f56adc/attachment-0003.html
More information about the gradsusr
mailing list