Changing level type
Added by Ralf Thehos over 13 years ago
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
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!!!