R/retrieve.dsensemble.R

In metno/esd: Climate analysis and empirical-statistical downscaling (ESD) package for monthly and daily data

#' @exportS3Method
#' @export 
retrieve.dsensemble <- function(ncfile,param="auto",
                             path=NULL,verbose=FALSE,...) {
  if(verbose) print("retrieve.dsensemble")
  nc <- nc_open(file.path(path,ncfile))
  #browser()
  dimnames <- names(nc$dim)
  lon <- ncvar_get(nc,dimnames[grep("lon|x",tolower(dimnames))])
  lat <- ncvar_get(nc,dimnames[grep("lat|y",tolower(dimnames))])
  nc_close(nc)
  if ( (length(dim(lon))==1) & (length(dim(lat))==1) )  {
    if (verbose) print('Regular grid field found')
    X <- retrieve.ncdf4(ncfile,path=path,param=param,verbose=verbose,...)
  } else {
    if (verbose) print('Irregular grid field found')
    X <- retrieve.rcm(ncfile,path=path,param=param,verbose=verbose,...) 
  }
}