Problem in writing .ctl file for global GRIB file

Wed Feb 27 05:24:39 EST 2008

Dear Wesley, Elif and GRADS users

Following your guidance I tried " grib2ctl.pl
v0.9.12.5p39c " and got the ctl file and pdef file and
could go on to plot the pictures.
Thank you very much.
Best regards,
Phuong
--- Wesley Ebisuzaki <Wesley.Ebisuzaki at NOAA.GOV>
wrote:

> Hi,
> > *  produced by grib2ctl v0.9.12.5p33k
> There have been 21 upgrades since 1/2005.  Try a
> newer version.
>
>          Wesley
>
>
>
>
> PHUONG NGUYEN wrote:
> > Dear Elif and GRADS users,
> > Following your guidance I tried
> >
> > grib2ctl -verf a01150000 > a01150000.ctl
> >
> > and received a .ctl file without xdef and ydef as
> > following :
> >
> >
> > dset ^a01150000
> > index ^a01150000.idx
> > undef 9.999E+20
> > title a01150000
> > *  produced by grib2ctl v0.9.12.5p33k
> > dtype grib 41
> > tdef 6 linear 06Z15jan2008 6hr
> > zdef 10 levels
> > 850 700 600 500 400 300 250 200 150 100
> > vars 10
> > HGTprs 10 7,100,0 ** Geopotential height [gpm]
> > HGTmwl  0 7,6,0 ** max wind level Geopotential
> height
> > [gpm]
> > HGTtrp  0 7,7,0 ** tropopause Geopotential height
> > [gpm]
> > RHprs 4 52,100,0 ** Relative humidity [%]
> > TMPprs 10 11,100,0 ** Temp. [K]
> > TMPtrp  0 11,7,0 ** tropopause Temp. [K]
> > UGRDprs 10 33,100,0 ** u wind [m/s]
> > UGRDmwl  0 33,6,0 ** max wind level u wind [m/s]
> > VGRDprs 10 34,100,0 ** v wind [m/s]
> > VGRDmwl  0 34,6,0 ** max wind level v wind [m/s]
> > ENDVARS
> >
> > This .ctl file still doesn't contain xdef and
> ydef.
> >
> > Please advise what could I try more ???
> >
> > Thank you very much.
> > Best regards
> >
> >
> > --- Elif SERTEL <esertel at ENVSCI.RUTGERS.EDU>
> wrote:
> >
> >
> >> Did you try grib2ctl?
> >>
> >>
> >>
> >
>
http://www.cpc.ncep.noaa.gov/products/wesley/grib2ctl.html
> >
> >> Elif.
> >>
> >>
> >> On Mon, February 4, 2008 3:54 am, PHUONG NGUYEN
> >> wrote:
> >>
> >>> Dear GRADS users,
> >>> I have a global GRIB file. I would like to write
> >>>
> >> the
> >>
> >>> .ctl file for plotting it.
> >>>
> >>> I use wgrib to see it structure. Based on the
> >>>
> >> results
> >>
> >>> I  can write a . ctl file as the following :
> >>>
> >>> dset ^a01150000
> >>> index ^a01150000.idx
> >>> undef 9.999E+20
> >>> title a01150000
> >>> dtype grib 41
> >>> xdef
> >>> ydef
> >>> tdef 1 linear 00Z15jan2008 1mo
> >>> zdef 10 levels
> >>> 850 700 600 500 400 300 250 200 150 100
> >>> vars 10
> >>> HGTprs 10 7,100,0 ** Geopotential height [gpm]
> >>> HGTmwl  0 7,6,0 ** max wind level Geopotential
> >>>
> >> height
> >>
> >>> [gpm]
> >>> HGTtrp  0 7,7,0 ** tropopause Geopotential
> height
> >>> [gpm]
> >>> RHprs 4 52,100,0 ** Relative humidity [%]
> >>> TMPprs 10 11,100,0 ** Temp. [K]
> >>> TMPtrp  0 11,7,0 ** tropopause Temp. [K]
> >>> UGRDprs 10 33,100,0 ** u wind [m/s]
> >>> UGRDmwl  0 33,6,0 ** max wind level u wind [m/s]
> >>> VGRDprs 10 34,100,0 ** v wind [m/s]
> >>> VGRDmwl  0 34,6,0 ** max wind level v wind [m/s]
> >>> ENDVARS
> >>>
> >>> The problem is in writing xdef and ydef : it is
> >>>
> >> global
> >>
> >>> file :
> >>>
> >>> rec 1:41:date 2008011500 HGT kpds5=7 kpds6=100
> >>> kpds7=850 levels=(3,82) grid=41 850 mb 6hr fcst:
> >>>   HGT=Geopotential height [gpm]
> >>>   timerange 0 P1 6 P2 0 TimeU 1  nx -1 ny 73 GDS
> >>>
> >> grid
> >>
> >>> 0 num_in_ave 0 missing 0
> >>>   center 74 subcenter 0 process 45 Table 2
> >>>   thinned latlon: lat  -90.000000 to 0.000000 by
> >>> 1.250000  nxny 3447
> >>>           long -30.000000 to 60.000000, 3447
> grid
> >>>
> >> pts
> >>
> >>>  (-1 x 73) scan 64 mode 128 bdsgrid 1
> >>>        2    3    5    6    8    9   11   12   14
> >>>
> >> 16
> >>
> >>>  17   19   20   22   23
> >>>       25   26   28   29   30   32   33   35   36
> >>>
> >> 38
> >>
> >>>  39   40   42   43   44
> >>>       45   47   48   49   50   51   52   54   55
> >>>
> >> 56
> >>
> >>>  57   58   59   60   60
> >>>       61   62   63   64   65   65   66   67   67
> >>>
> >> 68
> >>
> >>>  69   69   70   70   71
> >>>       71   71   72   72   72   73   73   73   73
> >>>
> >> 73
> >>
> >>>  73   73   73
> >>>   min/max data 959.999 1580  num bits 10
> BDS_Ref
> >>> 959.999  DecScale 0 BinScale 0
> >>>
> >>> rec 2:4624:date 2008011500 HGT kpds5=7 kpds6=100
> >>> kpds7=850 levels=(3,82) grid=42 850 mb 6hr fcst:
> >>>   HGT=Geopotential height [gpm]
> >>>   timerange 0 P1 6 P2 0 TimeU 1  nx -1 ny 73 GDS
> >>>
> >> grid
> >>
> >>> 0 num_in_ave 0 missing 0
> >>>   center 74 subcenter 0 process 45 Table 2
> >>>   thinned latlon: lat  -90.000000 to 0.000000 by
> >>> 1.250000  nxny 3447
> >>>           long 60.000000 to 150.000000, 3447
> grid
> >>>
> >> pts
> >>
> >>>  (-1 x 73) scan 64 mode 128 bdsgrid 1
> >>>        2    3    5    6    8    9   11   12   14
> >>>
> >> 16
> >>
> >>>  17   19   20   22   23
> >>>       25   26   28   29   30   32   33   35   36
> >>>
> >> 38
> >>
> >>>  39   40   42   43   44
> >>>       45   47   48   49   50   51   52   54   55
> >>>
> >> 56
> >>
> >>>  57   58   59   60   60
> >>>       61   62   63   64   65   65   66   67   67
> >>>
> >> 68
> >>
>
=== message truncated ===


MS NGUYEN THI MINH PHUONG
National Center for Hydrometeorological Forecasts
Address:     4 DANG THAI THAN, HANOI, VIETNAM,
Email :        ntmphuong at hotmail.com, ntmphuong at yahoo.com.au
======


      Get the name you always wanted with the new y7mail email address.
www.yahoo7.com.au/y7mail