displaying problem

[Stegert]

Dear Leva,

sorry for answering late. now returning to GrADS things again...
So, I tried out GrADS' behaviour on different ctls.
Using two grids:
1. cn3.ctl
XDEF 88 LINEAR  -15.0833 0.33333333
YDEF 82 LINEAR  47.6833  0.2
ZDEF  5 LEVELS 10 15 20 25 30
TDEF 13 LINEAR 01JAN2001 1DY

2. cn4.dtl
XDEF 58 LINEAR  -5.08333333 0.33333333
YDEF 65 LINEAR   48.8833 0.2
ZDEF  5 LEVELS 5 10 15 20
TDEF 13 LINEAR 01JAN1986 1DY

I opened cn3 first, afterwards the cn4.
Using q dims GrADS returned:
ga-> q dims
Default file number is: 1
X is varying   Lon = -5.08333 to 13.9167   X = 0.999999 to 58
Y is varying   Lat = 48.8833 to 61.6833   Y = 1 to 65
Z is fixed     Lev = 5  Z = 1
T is fixed     Time = 00Z01JAN1986  T = 1

then:
ga-> set dfile 2
Default file set to: e:\ecoham3\res\t05\hcn4.grads
ga-> q dims
Default file number is: 2
X is varying   Lon = -5.08333 to 13.9167   X = 30.9999 to 87.9999
Y is varying   Lat = 48.8833 to 61.6833   Y = 7 to 71
Z is fixed     Lev = 5  Z = 0
T is fixed     Time = 00Z01JAN1986  T = -5478

He memorized the LON, LAT, LEV and T data and transponed them
to the X,Y,Z values depending on the second grid.

X,Y are easy as the dx and dy are identical.
T set to -5478 is of course out of bounds, there's nothing to
plot for the cn4.

Additionally: As there is no 5m level in the second ctl, z is set
to zero and GrADS told me:
ga-> set t 1
Time values set: 2001:1:1:0 2001:1:1:0
ga-> d zoc
Data Request Warning:  Request beyond file limits
   Entire grid contents are set to missing data
   Dimension ranges are:  X = 31 88  Y = 7 82  Z = 0 24  T = 1 13
   Warning issued for variable = zoc
Cannot contour grid - all undefined values

Maybe this is where your problem lies!?
If you firstly open the mean maybe with Z=1/LEV=10m
and switch to file1 with LEV(1)=5m,LEV(2)=10m it works
(z will be 2 then!)
But vice versa you change from z=1/LEV=5m to where there is no
5m level, GrADS can't diplay.

Hope to help with this...

Greetings,
Christoph


Leva Djoneidi schrieb:
> Dear Christoph,
> I have one file (file1) that I use first on its own and display
> something with it (display1).
> Then I have 3 files (meanfiles) that have the same size and I use to
> display a mean value with them (displaymean).
> (file1) and (meanfiles) do not have the same size.
> Each display works on its own.
> When I first display the mean (displaymean), then do a set dfile
> #(file1) and display (display1), it also works.
> But if I first open file1 and display (display1), and then want to do a
> mean with the meanfiles (which are now number 2,3,4), it doesnt
> understand the mean.
> I tried to open file 1 then close it and then open the 3 meanfiles and
> do a mean (var.2+var.3+var.4)/but it says 'request beyond file limits'.
> If I didnt need to do a mean, I could set dfile #file display, then do
> another set dfile #file and display, but when I get to the meanfiles, I
> cannot do set #file because there are 3 files, and it doesnt understand
> the var.#meanfile1, var.#meanfile2...
>
> Do even my meanfiles and file1 have to have the same size, do you think
> its another problem...?
>
> Thank you for your help
>
> Leva

--
Christoph Stegert
Institute for Oceanography - Ecological Modelling group
University of Hamburg, GER - (ZMK/ZMAW alliance member)
Bundesstr.53 20146 Hamburg - +49-40/42838-7486 Room 348