[gradsusr] Problem matching GRIB data with descriptor file
Chidzambwa Simba
Simba.Chidzambwa at climatechange.qld.gov.au
Wed Sep 22 01:08:53 EDT 2010
I am using Unix. I want to display ACRE GRIB2 data in GRADS Version
2.0.a8 Copyright (c) 1988-2010. The data is for 3 hour
(sflxgrbfgensmean_1931??????_fhr03) and 6hour forecasts
(sflxgrbfgensmean_1931??????_fhr06) based on 6hour time-steps for each
day (00, 06, 12 and 18hrs). I used a template to create a control file
based on "grib2ctl.pl". My control file for 3 hour forecast based on 6
hour time-steps for 1-day looks like this;
"dset
//mnt/acre3/data/ensda_1890/18900101%h2/sflxgrbfgensmean_18900101%h2_fhr
03
index sflxgrbfgensmean_18900101_fhr03.idx
undef 9.999E+20
title sflxgrbfgensmean_18900101_fhr03
* produced by grib2ctl v0.9.12.5p45
dtype grib 98
options template
options yrev
xdef 192 linear 0.000000 1.875
ydef 94 levels
-88.542 -86.653 -84.753 -82.851 -80.947 -79.043 -77.139 -75.235 -73.331
-71.426...
tdef 4 linear 00Z1jan1890 6hr
zdef 1 linear 1 1
vars 88
ACONDsfc 0 179,1,0 ** surface Aerodynamic conductance [m/s]
ALBDOsfc 0 84,1,0 ** surface Albedo [%]
ENDVARS"
The output error message I get when I run gribmap is "Time out of
bounds" as below;
[chidzas at pollux simba]$ gribmap -v -i sflxgrbfgensmean_189001_fhr03.ctl
grib1map: opening GRIB file:
//mnt/acre3/data/ensda_1890/1890010100/sflxgrbfgensmean_1890010100_fhr03
----- Time out of bounds: 1 27173 0 98 1 0 124 1
0 96 0 btim: 1889123118:00 tau: 3 dtim: 1889123121:00
----- Time out of bounds: 2 54346 0 98 1 0 125 1
0 27269 0 btim: 1889123118:00 tau: 3 dtim: 1889123121:00
----- Time out of bounds: 3 77007 0 98 1 0 122 1
0 54442 0 btim: 1889123118:00 tau: 3 dtim: 1889123121:00
----- Time out of bounds: 4 99668 0 98 1 0 121 1
0 77103 0 btim: 1889123118:00 tau: 3 dtim: 1889123121:00
And "Time non-integral" as below;
grib1map: opening GRIB file:
//mnt/acre3/data/ensda_1890/1890010106/sflxgrbfgensmean_1890010106_fhr03
----- Time non-integral. 1 1.5: 1 27173 0 98 1 0
124 1 0 96 1724354 btim: 1890010100:00 tau:
3 dtim: 1890010103:00
----- Time non-integral. 1 1.5: 2 54346 0 98 1 0
125 1 0 27269 1724354 btim: 1890010100:00 tau:
3 dtim: 1890010103:00
----- Time non-integral. 1 1.5: 3 77007 0 98 1 0
122 1 0 54442 1724354 btim: 1890010100:00 tau:
3 dtim: 1890010103:00
----- Time non-integral. 1 1.5: 4 97412 0 98 1 0
121 1 0 77103 1724354 btim: 1890010100:00 tau:
3 dtim: 1890010103:00
----- Time non-integral. 1 1.5: 5 133609 0 98 1 0
11 1 0 97508 1724354 btim: 1890010100:00 tau: 3
dtim: 1890010103:00
I suspect it has something to do with definition of the 'TDEF" or the
options I am using or not using with "gribmap". I have no problem with
the 6 hour forecast control files which are matching well. When I
concatenate the files into a single file, it works ok with "gribmap -0"
option for the 3hr forecast, but no -0 is required for the 6hr forecast
file. Can anyone help?
Regards
Simba Chidzambwa
Coastal Unit, Dept of Environment, Queensland Australia
Simbarashe Chidzambwa
Senior Scientist, Coastal Unit, Queensland Climate Change Centre of
Excellence
Telephone: 07 3869 9507 Mob: 0420596025 Facsimile: 07 3869 9501
www.derm.qld.gov.au <http://www.derm.qld.gov.au/>
Department of Environment and Resource Management
27 Quinlan Street, Deagon QLD 4017
+----------------------------------------------------------------+
Think B4U Print
1 ream of paper = 6% of a tree and 5.4kg CO2 in the atmosphere
3 sheets of A4 paper = 1 litre of water
+----------------------------------------------------------------+
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100922/7859b08b/attachment-0003.html
More information about the gradsusr
mailing list