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 !92 (merged), 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