StartR and s2dv return different values
I tried two different ways of loading data, one with StartR, one with s2dv. And I think I ask for a similar interpolation but it doesn't return the same values. I use R/3.6.1-foss-2015a-bare on WS.
My codes:
library(startR)
rm(list=ls())
gc()
vari <- 'tos'
path_obs <- paste0('/esarchive/recon/ecmwf/era5/monthly_mean/', vari, '_f1h-r1440x721cds/$var$_$year$$month$.nc')
startdates_obs <- paste0(c(1980:1982))
lat_min <- -90
lat_max <- 90
lon_min <- 0
lon_max <- 360
data_obs <- startR::Start(dat = path_obs,
var = vari,
year = startdates_obs,
time = 'all',
lat = values(list(lat_min, lat_max)),
lon = values(list(lon_min, lon_max)),
month = 'all',
month_depends = 'year', #ca fait quoi si je mets ici 'time' ?
time_across = 'month',
merge_across_dims = TRUE,
merge_across_dims_narm = TRUE,
transform = startR::CDORemapper,
transform_extra_cells = 2,
transform_params = list(grid = 'r360x180',
method = 'conservative',
crop = c(lon_min, lon_max,
lat_min, lat_max)),
transform_vars = c('lat', 'lon'),
synonims = list(lon = c('lon', 'longitude'),
lat = c('lat', 'latitude')),
return_vars = list(lat='dat',
lon='dat',
time=c('year','month')),
#lon_reorder = CircularSort(-180,180),
lat_reorder = Sort(decreasing = T),
retrieve = TRUE)
and
library(maps)
library(s2dv)
# Clean memory
rm(list=ls())
gc()
vari <- 'tos'
startdates_obs <- paste0(c(1980:1982), '0101')
obs_era <- list(name = 'era5',
path = file.path('/esarchive/recon/ecmwf/era5/monthly_mean',
'$VAR_NAME$_f1h-r1440x721cds/$VAR_NAME$_$YEAR$$MONTH$.nc'))
data_obs <- s2dv::Load(var = vari,
exp = NULL,
obs = list(obs_era),
sdates= startdates_obs,
nmember=1,
leadtimemin = 1,
leadtimemax = 12,
output = 'lonlat',
grid = 'r360x180',
method = 'conservative',
storefreq = 'monthly',
nprocs= 1)
I checked that the retrieved lon and lat of both data are the same.
But I got different results:
- with StartR:
- with s2dv:
I'm doing maybe something wrong, but I can't see where.
Many thanks!
Aude