R/nc_date.R

In fengshy/sypaper: What the Package Does (One Line, Title Case)

#' nc_date
#'
#' @inheritParams ncread
#' @param fid An object of class ncdf4 retuned by [ncdf4::nc_open()] or nc file path
#' @param datastr Boolean. Whether convert date to string?
#'
#' @keywords internal
#'
#' @importFrom stringr str_extract
#' @import PCICt
#' @export
nc_date <- function(fid, unit = NULL, to_char = FALSE){
    if (class(fid)[1] != "ncdf4") {
        fid <- nc_open(fid)
        on.exit(nc_close(fid))
    }

    ctime    <- fid$dim$time # time class
    origin   <- ctime$units %>% str_extract("\\d{2,4}-\\d{1,2}-\\d{1,2}")
    calendar <- ctime$calendar

    if (is.null(calendar)) calendar <- "gregorian"


    if (is.null(unit)) unit = fid$dim$time$units %>% str_extract(".*(?= since)")
    secs <- switch(unit,
        "hours" = 3600,
        "days"  = 86400)
    date <- {as.PCICt(origin, cal=calendar) + timess*secs}
    if (to_char) date %<>% format(DATE_FORMAT)
    # date <- format(date, DATE_FORMAT)
    date
}