Project

General

Profile

remapnn destroys file structure?

Added by Matthias Demuzere almost 8 years ago

Hi all,

I am currently trying to perform a nearest neigbour remapping using the following command:
cdo remapnn,mod06_grid_1km cf_1step.rm.nc cf_1step_1km.nc

This works, but the values of the variable Band1 in the output file are gone.
Simultaneously, the variable crs has been removed (and this one seems linked as an attribute in Band1:grid_mapping = "crs"?).

SO it seems that the structure of my input file is not correct in order to properly work with the remapnn tool?

Any ideas on how to deal with this?

Thanks.
Matthias


Replies (3)

RE: remapnn destroys file structure? - Added by Ralf Mueller almost 8 years ago

Hi!

I could not reproduce the cleaning of Band1. Instead I got this (plotted with ncview)

In your input, the crs variable is empty, so I guess its just a placeholder to save the attributes for grid_mapping. CDO does not support the grid_mapping attribute, only for special cases like Lambert_Conformal, rotated_pole, laea and sinusoidal. In order to keep it in the file, you might add it later with NCO order use a global attribute instead. If the attributes of crs should be used by the remapping of CDO, I see no possibility to do that.

hth
ralf

Band1.png (41.4 KB) Band1.png

RE: remapnn destroys file structure? - Added by Ralf Mueller almost 8 years ago

BTW: I use cdo-1.7.2

cdo -V
Climate Data Operators version 1.7.2 (http://mpimet.mpg.de/cdo)
Compiled: by ram on luthien (x86_64-unknown-linux-gnu) Jun 28 2016 15:14:16
Compiler: gcc -g -O3 -std=gnu99 -Wall -fopenmp -march=native   
 version: gcc (GCC) 6.1.1 20160501
Features: DATA PTHREADS OpenMP45 HDF5 NC4/HDF5 OPeNDAP SZ UDUNITS2 PROJ.4 CURL FFTW3 AVX2
Libraries: HDF5/1.10.0 proj/4.92 curl/7.49.1
Filetypes: srv ext ieg grb grb2 nc nc2 nc4 nc4c 
     CDI library version : 1.7.2 of Jun 28 2016 15:13:56
 CGRIBEX library version : 1.7.5 of Jun  3 2016 14:44:00
GRIB_API library version : 1.15.0
  NetCDF library version : 4.4.1-rc1 of Apr 20 2016 12:06:30 $
    HDF5 library version : library undefined
 SERVICE library version : 1.4.0 of Jun 28 2016 15:13:53
   EXTRA library version : 1.4.0 of Jun 28 2016 15:13:52
     IEG library version : 1.4.0 of Jun 28 2016 15:13:52
    FILE library version : 1.8.2 of Jun 28 2016 15:13:52

RE: remapnn destroys file structure? - Added by Matthias Demuzere almost 8 years ago

Hi Ralf,

Thanks for your quick response.
I just upgraded my cdo version 1.7.0 to 1.7.2 and I also got rid of the error.

Perhaps I should do the same on our other system as well, since we had some issues there with the -t table statement (as in this post: https://code.zmaw.de/boards/1/topics/3834).
So perhaps a version upgrade can fix this as well?

Thanks for this!
Cheers,
Matthias

    (1-3/3)