Problems interpolating datasets in startR
Hi @vagudets ,
I have a problem when interpolating the observed tos data from esarchive/obs/ukmo/hadisst_v1.1/monthly_mean/tos/ to the same grid in /esarchive/exp/CMIP6/dcppA-hindcast/CanESM5/DCPP/CCCma/CanESM5/dcppA-hindcast/r1i1p2f1/Omon/tos/gr/v20190429/tos_Omon_CanESM5_dcppA-hindcast_s2019-r1i1p2f1_gr_202001-202912.nc
The dimensions of the observed tos data after interpolating should be 64 (in latitude) and 128 (in longitude), however, after interpolating the dimension in longitude is 127 instead of 128.
#### Script example:
library(maps)
library(RColorBrewer)
library(s2dv)
library(abind)
library(multiApply)
library(ncdf4)
library(startR)
library(plyr)
library(psych)
library(boot)
library(ClimProjDiags)
syear=paste(seq(1971,1990,by=1))
smonth=c('01','02','03','04','05','06','07','08','09','10','11','12')
ss <- apply(expand.grid(syear, smonth), 1, paste0, collapse = "")
ss
path2= paste0('/esarchive/obs/ukmo/hadisst_v1.1/monthly_mean/tos/$var$_$sdates$.nc')
# Reading tos data from hadisstv1.1
data_obs <- Start(dat = path2,
var = 'tos',
sdates = ss,
synonims = list(longitude = c('lon', 'longitude'),latitude = c('lat', 'latitude')),
time = 'all',
latitude = 'all',
longitude = 'all',
transform = CDORemapper,
transform_params = list(grid = '/esarchive/exp/CMIP6/dcppA-hindcast/CanESM5/DCPP/CCCma/CanESM5/dcppA-hindcast/r1i1p2f1/Omon/tos/gr/v20190429/tos_Omon_CanESM5_dcppA-hindcast_s2019-r1i1p2f1_gr_202001-202912.nc',
method = 'bil', force_remap = TRUE),
transform_vars = c('latitude', 'longitude'),
latitude_reorder = Sort(),
longitude_reorder = CircularSort(0, 360),
# Para que nos devuelva los valores numericos de lon, lat y time
return_vars = list(latitude = NULL,
longitude = NULL,
time = 'sdates'),
retrieve = TRUE)
# Para tener los valores de lon y lat de las observaciones
lat_obs <- as.vector(attr(data_obs, 'Variables')$common$latitude)
lon_obs <- as.vector(attr(data_obs, 'Variables')$common$longitude)
Module and Package Version
Which R version are you using? R version 3.6.1 (2019-07-05) Which R package versions are you using? Check with sessionInfo(). startR_2.2.0-2 Which machine are you using? WS
Many thanks in advance, Best, Veronica