Skip to content
GitLab
  • Menu
Projects Groups Snippets
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • startR startR
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 24
    • Issues 24
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 7
    • Merge requests 7
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • 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 Sciences
  • startRstartR
  • Issues
  • #156
Closed
Open
Created Jun 16, 2022 by acarreri@acarreriDeveloper

Load reconstruction not working with new expid

Hi @aho,

@asanto1 uses one of my script which loads spatial mean timeseries computed with Earthdiagnostic.

It works well with a3t4 for instance but not anymore with new reconstructions, such as a403. I don't know why.

Could you have a look at it please?

The code is this one:

library(startR)

name_recon <- 'a3t4'  # 'a403'
if (name_recon == 'a3t4'){
    model <- 'EC-Earth3-LR'
} else if (name_recon == 'a403'){
    model <- 'EC-Earth3'
}
reg <- 'Nino3.4'
vari <- 'tosmean'

path_recon <- paste0('/esarchive/exp/ecearth/',name_recon,'/diags/CMIP/EC-Earth-Consortium/',
                     model, '/piControl/$memb$/Omon/$var$/gn/v*/$var$_Omon_',model,
                     '_piControl_$memb$_gn_$chunk$.nc')
data_recon_call <- startR::Start(dat = path_recon,
                                      var = vari,
                                      memb = 'r1i1p1f1',
                                      region = reg,
                                      region_var = 'region',
                                      time = 'all', 
                                      chunk = 'all',
                                      chunk_depends = 'sdate',
                                      time_across = 'chunk', 
                                      merge_across_dims = TRUE,
                                      return_vars = list(time=NULL, region=NULL),
                                      retrieve = T)

Thanks!

Assignee
Assign to
Time tracking