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 3
    • Merge requests 3
  • 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
  • #46
Closed
Open
Issue created Oct 04, 2019 by aho@ahoMaintainer

Load(): problem when loading data from /esarchive/exp/ecmwf/system5_m1/monthly_mean/

The plot shows that when Load() retrieves the data under /esarchive/exp/ecmwf/system5_m1/monthly_mean/, the returned data is incorrect. No problem with other dataset (e.g., system4_m1).

library(startR)

  repos <- "/esarchive/exp/ecmwf/system5_m1/monthly_mean/$var$_f6h/$var$_$sdate$.nc"
  var <- 'tas'
  sdate <- c('20170101')#, '20170201')

  data_Start <- Start(dat = repos,
                var = var,
                sdate = sdate,
                ensemble = indices(1),
                time = 'all',
                latitude = 'all', #indices(1:100),
                longitude = 'all', #indices(1:100),
                return_vars = list(time = c('sdate'),
                                   longitude = NULL, latitude = NULL),
                retrieve = TRUE)

library(s2dverification)

#s2dv
   ecearth <- list(name = 'system5_m1',
                   path = file.path('/esarchive/exp/ecmwf',
                                    '$EXP_NAME$/$STORE_FREQ$_mean',  #month_mean',
                                    '$VAR_NAME$_f6h/$VAR_NAME$_$START_DATE$.nc'))
   startDates <- sdate
   time_length <- 7
   lon <- c(0, 359)
   lat <- c(-90, 90)

   data_Load <- Load(var = 'tas',
                     exp = list(ecearth),
                     sdates = startDates,
                     storefreq = "monthly",
                     leadtimemax = time_length,
                     output = 'lonlat',
                     nmember = 1)
                     #lonmin = lon[1], lonmax = lon[2],
                     #latmin = lat[1], latmax = lat[2])


#Plot
  lons <- data_Load$lon
  lats <- data_Load$lat

#plot Load
  PlotEquiMap(data_Load$mod[1, 1, 1, 1, , ], lons, lats)
#plot Start
  PlotEquiMap(data_Start[1, 1, 1, 1, 1, , ], lons, lats)

image image

Edited May 27, 2021 by aho
Assignee
Assign to
Time tracking