R/forecast-ensemble.R

Defines functions print.Forecast.ensemble ensemble.forecast

Documented in ensemble.forecast print.Forecast.ensemble

##' @rdname forecast
##' @description 
##' \command{ensemble.forecast} crée un objet de classe 
##' \sQuote{\var{Forecast.ens}}, à partir d'un jeu de prévisions
##' d'ensemble (GEFS ou CEP).
##'  
##' @param data le jeu de données sous la forme d'une liste de
##' \sQuote{\var{rasterBricks}}.
##' @param model le mot-clef désignant le modèle météorologique
##' à l'origine de la simulation.  
##' @param run la date de la simulation au format 
##' \sQuote{\var{POSIXc}}.
##' @param variable le mot-clef désigant la variable contenue 
##' dans la simulation.
##' @param echeances le vecteur des échéances de prévisions en 
##' heures. 
##' 
##' @export
ensemble.forecast <- 
  function(data, model, run, variable, echeances, units)
  {
    if (missing(variable)) variable="unknown"
    if (missing(units)) units="unknown"
<<<<<<< HEAD

    variable <- fcst.var(variable)    
    units <- switch(tolower(substr(variable,1,2)), 
                    pr="mm", te="°C", sn="mm", 
                    gu="m/s", ff="m/s", units)

    names(data) <- sprintf("H%d", echeances)
=======
    
    units <- switch(tolower(substr(variable,1,2)), 
                    pr="mm", te="°C", sn="mm", 
                    gu="m/s", ff="m/s", units)
    variable <- switch(tolower(variable), 
                       precip="Precip", temp="Tempe", snow="Snow", 
                       ff250="FF250m", ff500="FF500m",
                       gust="FFRaf", variable)
    names(data) <- sprintf("H%d",echeances)
>>>>>>> 53ae716ce8775c9fd36538db2ce00c6b918e4778
    structure(
      list( model=model, run=run, variable=variable, units=units, 
            echeances=echeances, dates=run + 3600 * echeances, 
            data=data), 
      class="Forecast.ensemble")
  }

##' @rdname forecast
##' @method print Forecast
##' @export 
print.Forecast.ensemble <- 
  function(x,...) 
  {
    # cat("Jeu de prévisions météorologique :\n")
    cat("\nModèle      :", x$model)
    cat("\nRéseau      :", format(x$run,"%A %d %B %Y %Hh %Z"))
    cat("\nParamètre   :", x$variable)
    cat("\nUnités      :", x$units)
    dech <- sort(unique(diff(x$echeances)))
    cat(sprintf("\nEcheances   : %d (%d -> %d [%d])", 
                length(x$echeances), min(x$echeances), max(x$echeances),
                dech))
    cat(sprintf("\nScenarios   : %d", length(x$data)))
    cat("\n-- Raster [1]--\n")
    print(x$data[[1]])
  }

# ##' @rdname forecast
# ##' @method plot Forecast
# ##' @export 
# plot.Forecast.ensemble <- 
#   function(x, layer=1, ...) 
#   {
#     
#   }
coolTot/DtgRecup documentation built on May 12, 2017, 9:45 a.m.