R/readwatbal.R

In RemkoDuursma/Maeswrap: Wrapper functions for MAESTRA/MAESPA

#' Reads the watbal.dat MAESPA output file
#'
#' @description Reads the hourly water balance output file ("watbal.dat").
#'
#'
#' @param filename Default name of the (half-)hourly water balance output file.
#' @return Returns a dataframe.
#' @author Remko Duursma
#' @keywords utilities
#' @examples
#'
#'
#' \dontrun{
#'
#' # Simple as this:
#' mywatbalresult <- readwatbal()
#'
#' # If you want to select the water balance file with a menu:
#' readwatbal(file.choose())
#'
#' }
#' @export
readwatbal <- function (filename = "watbal.dat"){
  watlines <- readLines(filename,100)
  colloc <- grep("Columns", watlines)
  namesline <- watlines[colloc]
  NAMES <- delempty(stringr::str_trim(strsplit(strsplit(namesline, ":")[[1]][2]," ")[[1]]))
  watbal <- data.table::fread(filename, header = FALSE, na.strings = "-999.0000", 
                              skip = colloc+1, data.table = FALSE)
  names(watbal) <- NAMES
  return(watbal)
}



#' Reads the watbalday.dat MAESPA output file
#'
#' @description Reads the daily water balance output file ("watbalday.dat").
#'
#'
#' @param filename Default name of the daily water balance output file.
#' @return Returns a dataframe.
#' @author Rémi Vezy
#' @keywords utilities
#' @examples
#'
#' \dontrun{
#'
#' # Simple as this:
#' mywatbalresult <- readwatbal()
#'
#' # If you want to select the water balance file with a menu:
#' readwatbalday(file.choose())
#'
#' }
#' @export
readwatbalday <- function (filename = "watbalday.dat") 
{
  watlines <- readLines(filename)
  colloc <- grep("Columns", watlines)
  namesline <- watlines[colloc]
  NAMES <- delempty(stringr::str_trim(strsplit(strsplit(namesline, ":")[[1]][2]," ")[[1]]))
  watbal <- data.table::fread(filename, header = FALSE, na.strings = "-999.0000",
                              skip = colloc, data.table = FALSE)
  names(watbal) <- NAMES
  watbal
}