Interpolation of Decadal original files into a regular grids avoiding vertical empty lines and artifacts
Finally, we are able to load and regrid IPSL files avoiding artifacts with the code:
library(startR) library(s2dverification) # Path to original cmorfiles: path <- '/esarchive/exp/ipsl-cm6a-lr/cmip6-dcppA-hindcast_i1p1/original_files/cmorfiles/DCPP/IPSL/IPSL-CM6A-LR/dcppA-hindcast/r1i1p1f1/SImon/siconc/gn/v20200101/$var$_SImon_IPSL-CM6A-LR_dcppA-hindcast_s2016-r1i1p1f1_gn_201701-202612.nc' # Regrid: data <- Start(dat = path, var = 'siconc', x = indices(2:361), y = indices(2:331), time = 'all', return_vars = list(nav_lat = NULL, nav_lon = NULL), retrieve = TRUE) nav_lon <- attributes(data)$Variables$common$nav_lon[2:361, 2:331] nav_lat <- attributes(data)$Variables$common$nav_lat[2:361, 2:331] res <- CDORemap(data[1,1, , ,5], nav_lon, nav_lat, grid = 'r360x180', method = 'bicubic', crop = FALSE) dev.new() image(res$lons, res$lats, res$data_array, main = "Regrided without the 1st and last indices")
The figure shows the original data (top-left), the artifacts (top-right) and the results excluding some indices before the interpolation with 'bicubic' or 'bilinear' methods.
Please, Roberto, let me know about the case of CMCC data, so, we can report it here.
With @aho I would like to discuss in a chat the possibility of creating a FAQ and the possibility of using 'transform' parameters in Start, instead of using CDORemap.
Please, don't hesitate to add here any information you may found important.
Thanks, especially to @rfernand for his patient.