[gradsusr] Problem matching GRIB data with descriptor file
Chidzambwa Simba
Simba.Chidzambwa at climatechange.qld.gov.au
Wed Sep 8 19:22:14 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
+----------------------------------------------------------------+
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/20100909/3d4c81ae/attachment-0003.html
More information about the gradsusr
mailing list