as.data.frame.airGRt.R

as.data.frame.airGRt <- function(x, row.names = NULL, ...) {
  
  if (!(inherits(x, "PrepGR") | inherits(x, "CalGR") | inherits(x, "SimGR"))) {
    stop("'InputsCrit' must be of class 'PrepGR', 'CalGR', 'SimGR'")
  }
  
  TMGR <- .TypeModelGR(x)
  myGR <- list()
  
  myGR$FracSolid <- NA
  myGR$TempMean  <- NA
  
  if (inherits(x, "PrepGR")) {
    if (TMGR$CemaNeige) {
      PrecipSol <- rowMeans(as.data.frame(x$InputsModel$LayerPrecip) * as.data.frame(x$InputsModel$LayerFracSolidPrecip), na.rm = TRUE)
      PrecipSim <- rowMeans(as.data.frame(x$InputsModel$LayerPrecip), na.rm = TRUE)
      FracSolid <- PrecipSol / PrecipSim
      FracSolid <- ifelse(is.na(FracSolid) & PrecipSol == 0 & PrecipSim == 0, 0, FracSolid)
      myGR$FracSolid <- FracSolid
      TempMean  <- rowMeans(as.data.frame(x$InputsModel$LayerTempMean), na.rm = TRUE)
      myGR$TempMean  <- TempMean
    }
    myGR$DatesR  <- x$InputsModel$DatesR
    myGR$PotEvap <- x$InputsModel$PotEvap
    myGR$Precip  <- x$InputsModel$Precip
    myGR$Qobs    <- x$Qobs
    myGR$Qsim    <- NA
  }
  if (inherits(x, "CalGR") | inherits(x, "SimGR")) {
    if (TMGR$CemaNeige) {
      PrecipSol <- rowMeans(sapply(x$OutputsModel$CemaNeigeLayers, "[[", "Psol"), na.rm = TRUE)
      PrecipSim <- rowMeans(sapply(x$OutputsModel$CemaNeigeLayers, "[[", "Pliq"), na.rm = TRUE) + PrecipSol
      FracSolid <- PrecipSol / PrecipSim
      FracSolid <- ifelse(is.na(FracSolid) & PrecipSol == 0 & PrecipSim == 0, 0, FracSolid)
      myGR$FracSolid <- FracSolid
      TempMean  <- rowMeans(sapply(x$OutputsModel$CemaNeigeLayers, "[[", "Temp"), na.rm = TRUE)
      myGR$TempMean  <- TempMean
    }
    myGR$DatesR  <- x$OutputsModel$DatesR
    myGR$PotEvap <- x$OutputsModel$PotEvap
    myGR$Precip  <- x$OutputsModel$Precip
    myGR$Qobs    <- x$Qobs
    myGR$Qsim    <- x$OutputsModel$Qsim
  }    
TabSim <- data.frame(Dates                     = myGR$DatesR,
                     PotEvap                   = myGR$PotEvap,
                     PrecipObs                 = myGR$Precip,
                     PrecipFracSolid_CemaNeige = myGR$FracSolid,
                     TempMeanSim_CemaNeige     = myGR$TempMean,
                     Qobs                      = myGR$Qobs,
                     Qsim                      = myGR$Qsim)
as.data.frame(x = TabSim, row.names = row.names, ...)
}