Skip to content
GitLab
Projects Groups Topics Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • s2dv s2dv
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributor statistics
    • Graph
    • Compare revisions
  • Issues 17
    • Issues 17
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 2
    • Merge requests 2
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Package Registry
    • Terraform modules
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • Earth SciencesEarth Sciences
  • s2dvs2dv
  • Issues
  • #64
Closed
Open
Issue created Jun 21, 2022 by aho@ahoMaintainer

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

Assignee
Assign to
Time tracking