Skip to content
GitLab
Projects Groups Topics 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
    • Contributor statistics
    • Graph
    • Compare revisions
  • Issues 29
    • Issues 29
    • 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 SciencesEarth Sciences
  • startRstartR
  • Issues
  • #115
Closed
Open
Issue created Sep 22, 2021 by Nuria Pérez-Zanón@nperezMaintainer

Interpolation of Decadal original files into a regular grids avoiding vertical empty lines and artifacts

Hi @rfernand

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")

Captura_de_pantalla_2021-09-22_a_las_12.25.56

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.

This issue is related to s2dverification#259 (closed) but at that time the files were processed by @pabretonniere and Marga. It may also be interesting for @cdelgado.

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.

Núria

Assignee
Assign to
Time tracking