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
  • #99
Closed
Open
Issue created May 26, 2023 by aclerici@aclerici

Problem when loading ECMWF seasonal forecast data

Hello, I am trying to load data from ECMWF (seasonal forecast daily and subdaily data on single levels). The data is from system 5. I am using the CSTools Load function, but I also tried to use the s2dv Load function directly, with the same issue. The data is not loaded correctly and data is missing when different number and time step is selected from the data cube.

The file I'm using: https://drive.google.com/file/d/17qiVUbH3TbK3_bfoMyuSTQQmI6aJr-XT/view?usp=sharing

The minimal loading code is the following:

exp <- list(list(
    name = "experiment_ecmwf",
    path = file.path("/experiment/$YEAR$/$VAR_NAME$_$START_DATE$.nc")
  ))

domain <- list(lonmin=-15, lonmax=60, latmin=30, latmax=70)

exp_dat <- CST_Load(
    var = "mx2t24",
    exp = exp,
    obs = NULL,
    sdates = c('19930201'),
    output = "lonlat",
    leadtimemin = 1,
    leadtimemax = 214,
    latmin = domain$latmin,
    latmax = domain$latmax,
    lonmin = domain$lonmin,
    lonmax = domain$lonmax,
    dimnames = list(lon='longitude', lat='latitude', member='number'))

test_expdat_equimap

Assignee
Assign to
Time tracking