Calibration_Michel.R

Calibration_Michel <- function(InputsModel, RunOptions, InputsCrit, CalibOptions,
                                    FUN_MOD, FUN_CRIT, FUN_TRANSFO = NULL, verbose = TRUE) {
  
  
  ##_____Arguments_check_____________________________________________________________________
  if (!inherits(InputsModel, "InputsModel")) {
    stop("InputsModel must be of class 'InputsModel' \n")
    return(NULL)
  }  
  if (!inherits(RunOptions, "RunOptions")) {
    stop("RunOptions must be of class 'RunOptions' \n")
    return(NULL)
  }  
  if (!inherits(InputsCrit, "InputsCrit")) {
    stop("InputsCrit must be of class 'InputsCrit' \n")
    return(NULL)
  }
  if (inherits(InputsCrit, "Multi")) {
    stop("InputsCrit must be of class 'Single' or 'Compo' \n")
    return(NULL)
  }
  if (!inherits(CalibOptions, "CalibOptions")) {
    stop("CalibOptions must be of class 'CalibOptions' \n")
    return(NULL)
  }  
  if (!inherits(CalibOptions, "HBAN")) {
    stop("CalibOptions must be of class 'HBAN' if Calibration_Michel is used \n")
    return(NULL)
  }
  if (!missing(FUN_CRIT)) {
    warning("argument 'FUN_CRIT' is deprecated. The error criterion function is now automatically get from the 'InputsCrit' object", call. = FALSE)
  }  

  
  
  ##_check_FUN_TRANSFO
  if (is.null(FUN_TRANSFO)) {
    if (identical(FUN_MOD, RunModel_GR4H         )) {
      FUN_TRANSFO <- TransfoParam_GR4H
    }
    if (identical(FUN_MOD, RunModel_GR4J         )) {
      FUN_TRANSFO <- TransfoParam_GR4J
    }
    if (identical(FUN_MOD, RunModel_GR5J         )) {
      FUN_TRANSFO <- TransfoParam_GR5J
    }
    if (identical(FUN_MOD, RunModel_GR6J         )) {
      FUN_TRANSFO <- TransfoParam_GR6J
    }
    if (identical(FUN_MOD, RunModel_GR2M         )) {
      FUN_TRANSFO <- TransfoParam_GR2M
    }
    if (identical(FUN_MOD, RunModel_GR1A         )) {
      FUN_TRANSFO <- TransfoParam_GR1A
    }
    if (identical(FUN_MOD, RunModel_CemaNeige    )) {
      FUN_TRANSFO <- TransfoParam_CemaNeige
    }
    if (identical(FUN_MOD, RunModel_CemaNeigeGR4J) | identical(FUN_MOD, RunModel_CemaNeigeGR5J) | identical(FUN_MOD, RunModel_CemaNeigeGR6J)) {
      if (identical(FUN_MOD, RunModel_CemaNeigeGR4J)) {
        FUN1 <- TransfoParam_GR4J
        FUN2 <- TransfoParam_CemaNeige
      }
      if (identical(FUN_MOD, RunModel_CemaNeigeGR5J)) {
        FUN1 <- TransfoParam_GR5J
        FUN2 <- TransfoParam_CemaNeige
      }
      if (identical(FUN_MOD,RunModel_CemaNeigeGR6J)) {
        FUN1 <- TransfoParam_GR6J
        FUN2 <- TransfoParam_CemaNeige