Templates
Douglas Clark
dbcl at CEH.AC.UK
Wed Feb 9 11:46:34 EST 2005
I should also have said that in fact I doubt a "simple time series" plot
will work.
eg Assuming you are looking at January data, as in your example, then
if you set t to cover 10 years, when GrADS looks for data in months
Feb-Dec it's not going to find them, since they're not in any file that
is described by the template. I haven't tried this, but I doubt it'll
work. Instead you'd have to explictly set the times to cover each
January in succession, eg. set t 1 248 for year 1, set t 365*8
(365+31)*8 for year 2 (noting that the arithmetic can't be used via the
command line, but you get the idea).
Doug
>>> dbcl at CEH.AC.UK 09/02/2005 16:32:48 >>>
Jean Pierre,
I don't think you have specified enough times in your tdef line.
tdef 1245 linear 00Z30jan1998 3hr means that the entire dataset
(across
all files) contains 1245 times. 1245*3hrs is much less than one year.
For N years of data, you will need approx N*365*8 in your tdef line.
Of
course most of these will be missing, since you only have one month of
data per year, so a simple time series plot will be rather patchy.
(I'm assuming that the fact the data are GRIB has no effect on
template, but I don't know, having rarely used GRIB.)
Doug
>>> arabonis at EGS.UCT.AC.ZA 09/02/2005 16:19:28 >>>
Hi all
I have several years worth of wavewatch 3 data in monthly files, I
would like to create a template ctl file to deal with specific months
over the years, ie nww3.hs.199801.grb nww3.hs.199901.grb
nww3.hs.200001.grb nww3.hs.200101.grb nww3.hs.200201.grb
so I need a ctl file to describe them all, note the sequencing is
unimportant as I am looking at frequencies of occurrences.
my monthly ctl file looks like:
dset /data/Wamhist/working/nww3.hs.199801.grb
index /data/Wamhist/working/nww3.hs.199801.grb.idx
undef 9.999E+20
title /data/Wamhist/working/nww3.hs.199801.grb
* produced by grib2ctl v0.9.12.5p32r
dtype grib 233
options yrev
ydef 157 linear -78.000000 1
xdef 288 linear 0.000000 1.250000
tdef 249 linear 00Z01jan1998 3hr
zdef 1 linear 1 1
vars 1
HTSGWsfc 0 100,1, 1 ** surface Sig height of wind waves and swell
[m]
ENDVARS
I created a new ctl file called nww3.hs.all01.grb.ctl, and made the
following changes. ie
altered the dset line to read
dset /data/Wamhist/working/nww3.hs.%y401.grb
the options line to
options template yrev
and the tdef line to
tdef 1245 linear 00Z30jan1998 3hr
Now I can see the first file but all following ones display "entire
grid
undefined"
Anyone see where I am going wrong.
Many Thanks
Jean Pierre
--
Jean Pierre Arabonis
arabonis at egs.uct.ac.za
Tel 021 780 1021 cell 084 401 1365
This email is 100% Microsoft Free
More information about the gradsusr
mailing list