Problem in writing .ctl file for global GRIB file
Wesley Ebisuzaki
Wesley.Ebisuzaki at NOAA.GOV
Tue Feb 5 10:29:55 EST 2008
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
>>
>>> 69 69 70 70 71
>>> 71 71 72 72 72 73 73 73 73
>>>
>> 73
>>
>>> 73 73 73
>>> min/max data 1038 1571 num bits 10 BDS_Ref
>>>
>> 1038
>>
>>> DecScale 0 BinScale 0
>>>
>>> rec 3:9207:date 2008011500 HGT kpds5=7 kpds6=100
>>> kpds7=850 levels=(3,82) grid=43 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 150.000000 to -120.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 929.963 1632.96 num bits 10
>>>
>> BDS_Ref
>>
>>> 929.963 DecScale 0 BinScale 0
>>>
>>> rec 4:13790:date 2008011500 HGT kpds5=7 kpds6=100
>>> kpds7=850 levels=(3,82) grid=44 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 -120.000000 to -30.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 1023.33 1632.33 num bits 10
>>>
>> BDS_Ref
>>
>>> 1023.33 DecScale 0 BinScale 0
>>>
>>>
>>> Above there are 4 first records : lat changes from
>>>
>> -90
>>
>>> to 0 by 1.25, lon changes in each record (for
>>>
>> example,
>>
>>> in rec 1 : lon -30 60, in rec 2 : lon 60 150,
>>>
>> etc),
>>
>>> and nx -1, ny 73, nxny 3447 grid pts.
>>>
>>> Then from rec 5 to rec 8 the lat changes from 0 to
>>>
>> 90
>>
>>> by 1.25, and the lon changes as in each rec above.
>>>
>>> So how to write the xdef and ydef in this case ???
>>> Please help me.
>>> Thank you very much.
>>> Phuong
>>>
>>>
>>> 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
>>>
> === 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
>
More information about the gradsusr
mailing list