Skip to content

CDORemap: CDO warning "Coordinates variable xxx can't be assigned!"

CDORemap() returns the warnings like Warning (cdfInqContents): Coordinates variable time can't be assigned! sometimes. The following example can produce this problem. If grid = 'r100x50, the warnings disappear.

# Load a data array

  hcst <- Start(dat ="/esarchive/exp/ecmwf/system5c3s/daily_mean/$var$_f6h/$var$_$file_date$.nc",
                var = "tas",
                file_date = "19931101",
                time = 1,
                latitude = values(list(10, 20)),
                latitude_reorder = Sort(),
                longitude = values(list(0, 20)),
                longitude_reorder = CircularSort(0, 360),
                synonims = list(latitude = c('lat', 'latitude'),
                                longitude = c('lon', 'longitude'),
                                ensemble = c('member', 'ensemble')),
                ensemble = indices(1),
                return_vars = list(latitude = 'dat', longitude = 'dat',
                                   time = 'file_date'),
                split_multiselected_dims = F,
                retrieve = TRUE)
data <- drop(hcst)
lon <- attr(data, 'Variables')$dat1$longitude
lat <- attr(data, 'Variables')$dat1$latitude

# Interpolation
res <- CDORemap(data, lons = lon, lats = lat, grid = "/esarchive/recon/ecmwf/era5/monthly_mean/tas_f1h-r1440x721cds/tas_201805.nc", method = 'bil')

Warning (cdfInqContents): Coordinates variable time can't be assigned!
Warning (cdfInqContents): Coordinates variable height can't be assigned!

A similar issue (has been solved) is https://earth.bsc.es/gitlab/es/s2dv/-/merge_requests/92, but the warning there is Warning (find_time_vars): Time variable >time< not found!

The result is correct, just the warnings are annoying, especially when startR::CDORemapper is used, CDORemap() is called several times and creates a sea of warnings.

FYI @vagudets