Changing level type

Hello,

i want to change the leveltype of a grib-file (DWD-COSMO-EU, 40 levels) from 109 to 110. This works with setltype, but the description of the level in the new Grib-file is like

hybrid 1-0 , 2-0, ... 41-0. It has to be 2-1, 3-2, ... 41-40

I then tried the setzaxis-command (setzaxis,vct1.txt infile outfile) but with this i get the error

Warning (varDefZaxis) : VCT missing
cdo setzaxis (Warning): No zaxis with 40 levels found!
Warning (cgribexDefLevel) : VCT missing. ( code = 40, zaxisID = 1 )

The file vct1.txt is attached.

Can someone give me a tipp how to solve that problem?

Thank you,

Ralf

vct1.txt (889 Bytes) vct1.txt

Replies (2)

RE: Changing level type - Added by Uwe Schulzweida over 13 years ago

Hi Ralf,

changing the GRIB leveltype from 109 to 110 does not work with setltype.
So you have to use setzaxis with a description of the lower and upper bounds of the layers.
The parameter name are lbounds and ubound. Here is an example:

zaxistype = hybrid
size      = 40
levels    = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
lbounds   = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
ubounds   = 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41
vctsize   = 82
vct       = 2000.00 4000.00 ... 99763.00 100000.00
            0.02000 0.04000 ... 0.99763 1.00000

But this is independent from the error message. You tried to set a zaxis with 40 layers but there is no 3D variable with
40 layers found in the datafile.

Regards,
Uwe

RE: Changing level type - Added by Ralf Thehos over 13 years ago

Hi Uwe,

that works with the right field.

Thank you!!!

(1-2/2)

Project

General

Profile

CDO