[gradsusr] Problem with wgrib2 and RapidRefresh rotated grid

Ling,Alister [Edm] Alister.Ling at EC.gc.ca
Thu Oct 28 12:05:09 EDT 2010


Thank you for the clarification Wesley.

Outputting on the rotated grid is the decision of the RapidRefresh
group. They are continuing to export files for Alaska and CONUS on
Lambert Conformal. I will request a Lambert Conformal or Polar
Stereographic . The NCEP Grid 221 used by AWIP32 for NAM output would be
fine for us in Canada.

Regards,
Alister.

-----Original Message-----
From: gradsusr-bounces at gradsusr.org
[mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Wesley Ebisuzaki
Sent: Thursday, October 28, 2010 8:04 AM
To: GrADS Users Forum
Subject: Re: [gradsusr] Problem with wgrib2 and RapidRefresh rotated
grid

Alister,

Ling,Alister [Edm] wrote:
> Hello,
>
> I am having difficulty with displaying the new rotated grids for the 
> full RapidRefresh Domain and would appreciate any pointers to 
> files/tables that I need to update.
>   
Grid template  32769 is Rotated Latitude/Longitude (Arakawa Non-E
Staggered Grid) and is a new addition to the local NCEP grid templates.

Can't you find a nice Lambert Conformal grid to use?

So you really really want to see the data on the model grid.  Problem,
GrADS doesn't handle grid template 32769.

If you just want numbers, use the option "-g2clib 0" when running
wgrib2.  However, wgrib2 will not know the locations of the grid points
until someone writes some code.


      Wesley

> I have upgraded wgrib2 and compiled to  v0.1.8.6 7/2010.
>
> When I run wgrib2 wrfprs_rr_06.grib2 -domain I get:
> Getgridtemplate: GDT Template 3.32769 not defined.
> Fatal decode err=10 msg 1.1
> -s option spits out the fields.
>
> wgrib2 wrfprs_rr_06.grib2 -set_flag_table_3.3 48 -grib
> wrfprs_rr_06mod.grib2
>
> but it obviously does not understand the grid when I run g2ctl:
> *  produced by g2ctl v0.0.3h
> * griddef=1:0:(-1 x -1):grid_template=32769:no other grid info
>
> Similarly no luck with the GRIB1 files. But the verbose inventory
output
> (for 1 record) says
> rec 1:0:date 2010102015 HGT kpds5=7 kpds6=100 kpds7=50 levels=(0,50)
> grid=255 50 mb anl:
>   HGT=Geopotential height [gpm]
>   timerange 0 P1 0 P2 0 TimeU 1  nx 758 ny 567 GDS grid 205 num_in_ave
0
> missing 0
>   center 59 subcenter 0 process 125 Table 2
>   min/max data 20165.1 20739.5  num bits 15  BDS_Ref 20165.1  DecScale
0
> BinScale -5
>
> So there is grid information in the grib1 file, but then g2ctl  gives
me
> *** script needs to be modified ***
> unknown user-defined grid
>
> Thanks,
> Alister.
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>   

_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr




More information about the gradsusr mailing list