intlevel3d and intlevelx3d: getting slightly wrong results
Added by Andrea Kaiser-Weiss over 4 years ago
Hi,
I'd like to interpolate methane from ECMWF IFS model levels (source) onto ICON model levels (target), using geopotential as vertical coordinate. I have not used intlevel3d before, and as results looked funny, I tested with synthetic fields (attached) with unexpected results:
When I set for test reasons the target geopotential equivalent to source geopotential, I would expect to get as resulting output the same ch4 as I fed in as input:
cdo -v intlevel3d,source_geopot_cut.nc CH4_float.nc source_geopot_cut.nc out_soll_gleich.nc
However, the resulting output file is shifted vertically for one level, and first level gets missing values - also when I try intlevelx3d.
Is something wrong with my input files? Or is there some issue/bug with intlevel3d and intlevelx3d? Or is my usage incorrect?
Many thanks for help,
Andrea Kaiser-Weiss
| CH4_float.nc.gz (35.3 KB) CH4_float.nc.gz | data file to be interpolated | ||
| source_geopot_cut.nc.gz (34.7 KB) source_geopot_cut.nc.gz | source 3d vertical coordinate, also test target | ||
| out_soll_gleich.nc.gz (65.7 KB) out_soll_gleich.nc.gz | resulting output |
Replies (1)
RE: intlevel3d and intlevelx3d: getting slightly wrong results - Added by Uwe Schulzweida over 4 years ago
The operator intlevel3d obviously does not work correctly since version 1.9.4.
In version 1.9.4 and 1.9.5 there is a memory access error. Unfortunately, this was fixed incorrectly in 1.9.6.
We will have solved the problem in the next version 1.9.9.
Thanks for this example!
Cheers,
Uwe