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!