From b98b572ac589264c0d1de4e740cbfd5d9c5d7389 Mon Sep 17 00:00:00 2001
From: Olivier Delaigue <olivier.delaigue@irstea.fr>
Date: Thu, 21 Apr 2016 10:13:59 +0200
Subject: [PATCH] =?UTF-8?q?Retrait=20du=20r=C3=A9pertoire=20des=20codes=20?=
 =?UTF-8?q?d'exemple=20(les=20codes=20sont=20dans=20les=20fichiers=20.Rd)?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 tests/example_Calibration.R            | 48 --------------------
 tests/example_Calibration_Michel.R     | 46 -------------------
 tests/example_ErrorCrit.R              | 61 --------------------------
 tests/example_RunModel.R               | 29 ------------
 tests/example_RunModel_CemaNeige.R     | 25 -----------
 tests/example_RunModel_CemaNeigeGR4J.R | 31 -------------
 tests/example_RunModel_CemaNeigeGR5J.R | 31 -------------
 tests/example_RunModel_CemaNeigeGR6J.R | 31 -------------
 tests/example_RunModel_GR1A.R          | 40 -----------------
 tests/example_RunModel_GR2M.R          | 38 ----------------
 tests/example_RunModel_GR4H.R          | 28 ------------
 tests/example_RunModel_GR4J.R          | 28 ------------
 tests/example_RunModel_GR5J.R          | 28 ------------
 tests/example_RunModel_GR6J.R          | 28 ------------
 tests/example_SeriesAggreg.R           | 18 --------
 tests/example_TransfoParam.R           | 15 -------
 tests/example_TransfoParam_CemaNeige.R | 15 -------
 tests/example_TransfoParam_GR1A.R      | 15 -------
 tests/example_TransfoParam_GR2M.R      | 15 -------
 tests/example_TransfoParam_GR4H.R      | 15 -------
 tests/example_TransfoParam_GR4J.R      | 15 -------
 tests/example_TransfoParam_GR5J.R      | 15 -------
 tests/example_TransfoParam_GR6J.R      | 15 -------
 tests/example_plot_OutputsModel.R      | 54 -----------------------
 24 files changed, 684 deletions(-)
 delete mode 100644 tests/example_Calibration.R
 delete mode 100644 tests/example_Calibration_Michel.R
 delete mode 100644 tests/example_ErrorCrit.R
 delete mode 100644 tests/example_RunModel.R
 delete mode 100644 tests/example_RunModel_CemaNeige.R
 delete mode 100644 tests/example_RunModel_CemaNeigeGR4J.R
 delete mode 100644 tests/example_RunModel_CemaNeigeGR5J.R
 delete mode 100644 tests/example_RunModel_CemaNeigeGR6J.R
 delete mode 100644 tests/example_RunModel_GR1A.R
 delete mode 100644 tests/example_RunModel_GR2M.R
 delete mode 100644 tests/example_RunModel_GR4H.R
 delete mode 100644 tests/example_RunModel_GR4J.R
 delete mode 100644 tests/example_RunModel_GR5J.R
 delete mode 100644 tests/example_RunModel_GR6J.R
 delete mode 100644 tests/example_SeriesAggreg.R
 delete mode 100644 tests/example_TransfoParam.R
 delete mode 100644 tests/example_TransfoParam_CemaNeige.R
 delete mode 100644 tests/example_TransfoParam_GR1A.R
 delete mode 100644 tests/example_TransfoParam_GR2M.R
 delete mode 100644 tests/example_TransfoParam_GR4H.R
 delete mode 100644 tests/example_TransfoParam_GR4J.R
 delete mode 100644 tests/example_TransfoParam_GR5J.R
 delete mode 100644 tests/example_TransfoParam_GR6J.R
 delete mode 100644 tests/example_plot_OutputsModel.R

diff --git a/tests/example_Calibration.R b/tests/example_Calibration.R
deleted file mode 100644
index 694d3367..00000000
--- a/tests/example_Calibration.R
+++ /dev/null
@@ -1,48 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## calibration period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## calibration criterion: preparation of the InputsCrit object
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-
-## preparation of CalibOptions object
-CalibOptions <- CreateCalibOptions(FUN_MOD=RunModel_GR4J,FUN_CALIB=Calibration_Michel)
-
-## calibration
-OutputsCalib <- Calibration(InputsModel=InputsModel,RunOptions=RunOptions,InputsCrit=InputsCrit,
-                            CalibOptions=CalibOptions,FUN_MOD=RunModel_GR4J,FUN_CRIT=ErrorCrit_NSE,
-                            FUN_CALIB=Calibration_Michel)
-
-## simulation
-Param <- OutputsCalib$ParamFinalR
-OutputsModel <- RunModel(InputsModel=InputsModel,RunOptions=RunOptions,
-                         Param=Param,FUN=RunModel_GR4J)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-## efficiency criterion: Kling-Gupta Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_KGE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_KGE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-
diff --git a/tests/example_Calibration_Michel.R b/tests/example_Calibration_Michel.R
deleted file mode 100644
index efc26721..00000000
--- a/tests/example_Calibration_Michel.R
+++ /dev/null
@@ -1,46 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## calibration period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## calibration criterion: preparation of the InputsCrit object
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-
-## preparation of CalibOptions object
-CalibOptions <- CreateCalibOptions(FUN_MOD=RunModel_GR4J,FUN_CALIB=Calibration_Michel)
-
-## calibration
-OutputsCalib <- Calibration_Michel(InputsModel=InputsModel,RunOptions=RunOptions,
-                                 InputsCrit=InputsCrit,CalibOptions=CalibOptions,
-                                 FUN_MOD=RunModel_GR4J,FUN_CRIT=ErrorCrit_NSE)
-
-## simulation
-Param <- OutputsCalib$ParamFinalR
-OutputsModel <- RunModel_GR4J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-## efficiency criterion: Kling-Gupta Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_KGE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_KGE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_ErrorCrit.R b/tests/example_ErrorCrit.R
deleted file mode 100644
index c19c5105..00000000
--- a/tests/example_ErrorCrit.R
+++ /dev/null
@@ -1,61 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(734.568,-0.840,109.809,1.971)
-OutputsModel <- RunModel(InputsModel=InputsModel,RunOptions=RunOptions,
-                         Param=Param,FUN=RunModel_GR4J)
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-## efficiency criterion: Nash-Sutcliffe Efficiency on log-transformed flows
-transfo <- "log"
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run],transfo=transfo)
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-## efficiency criterion: Nash-Sutcliffe Efficiency above a threshold (q75%)
-BoolCrit <- rep(TRUE,length(BasinObs$Qmm[Ind_Run])); 
-BoolCrit[BasinObs$Qmm[Ind_Run]<quantile(BasinObs$Qmm[Ind_Run],0.75,na.rm=TRUE)] <- FALSE;
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run],BoolCrit=BoolCrit)
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
-## efficiency criterion: Kling-Gupta Efficiency
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_KGE,InputsModel=InputsModel,
-                               RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_KGE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-cat(paste("SubCrit  ",OutputsCrit$SubCritNames,"  ",round(OutputsCrit$SubCritValues,4),"\n",sep=""))
-
-## efficiency criterion: Kling-Gupta Efficiency below a threshold (q10%) on log-trqansformed flows
-transfo  <- "log"
-BoolCrit <- rep(TRUE,length(BasinObs$Qmm[Ind_Run])); 
-BoolCrit[BasinObs$Qmm[Ind_Run]>quantile(BasinObs$Qmm[Ind_Run],0.10,na.rm=TRUE)] <- FALSE;
-InputsCrit <- CreateInputsCrit(FUN_CRIT=ErrorCrit_KGE,InputsModel=InputsModel,RunOptions=RunOptions,
-                               Qobs=BasinObs$Qmm[Ind_Run],BoolCrit=BoolCrit,transfo=transfo)
-OutputsCrit <- ErrorCrit_KGE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-cat(paste("SubCrit  ",OutputsCrit$SubCritNames,"  ",round(OutputsCrit$SubCritValues,4),"\n",sep=""))
-
-
-
-
diff --git a/tests/example_RunModel.R b/tests/example_RunModel.R
deleted file mode 100644
index b240b8cf..00000000
--- a/tests/example_RunModel.R
+++ /dev/null
@@ -1,29 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(734.568,-0.840,109.809,1.971)
-OutputsModel <- RunModel(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param,
-                         FUN_MOD=RunModel_GR4J)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_CemaNeige.R b/tests/example_RunModel_CemaNeige.R
deleted file mode 100644
index 134bd7d7..00000000
--- a/tests/example_RunModel_CemaNeige.R
+++ /dev/null
@@ -1,25 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123002)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_CemaNeige,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,TempMean=BasinObs$T,
-                                 ZInputs=BasinInfo$HypsoData[51],HypsoData=BasinInfo$HypsoData,
-                                 NLayers=5)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_CemaNeige,InputsModel=InputsModel,
-                               IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(0.962,2.249)
-OutputsModel <- RunModel_CemaNeige(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel)
-
diff --git a/tests/example_RunModel_CemaNeigeGR4J.R b/tests/example_RunModel_CemaNeigeGR4J.R
deleted file mode 100644
index 6ed335ae..00000000
--- a/tests/example_RunModel_CemaNeigeGR4J.R
+++ /dev/null
@@ -1,31 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123002)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_CemaNeigeGR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E,TempMean=BasinObs$T,
-                                 ZInputs=BasinInfo$HypsoData[51],HypsoData=BasinInfo$HypsoData,
-                                 NLayers=5)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_CemaNeigeGR4J,InputsModel=InputsModel,
-                               IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(408.774,2.646,131.264,1.174,0.962,2.249)
-OutputsModel <- RunModel_CemaNeigeGR4J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_CemaNeigeGR5J.R b/tests/example_RunModel_CemaNeigeGR5J.R
deleted file mode 100644
index f543b5d5..00000000
--- a/tests/example_RunModel_CemaNeigeGR5J.R
+++ /dev/null
@@ -1,31 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123002)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_CemaNeigeGR5J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E,TempMean=BasinObs$T,
-                                 ZInputs=BasinInfo$HypsoData[51],HypsoData=BasinInfo$HypsoData,
-                                 NLayers=5)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_CemaNeigeGR5J,InputsModel=InputsModel,
-                               IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(179.139,-0.100,203.815,1.174,2.478,0.977,2.774)
-OutputsModel <- RunModel_CemaNeigeGR5J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_CemaNeigeGR6J.R b/tests/example_RunModel_CemaNeigeGR6J.R
deleted file mode 100644
index a24e2098..00000000
--- a/tests/example_RunModel_CemaNeigeGR6J.R
+++ /dev/null
@@ -1,31 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123002)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_CemaNeigeGR6J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E,TempMean=BasinObs$T,
-                                 ZInputs=BasinInfo$HypsoData[51],HypsoData=BasinInfo$HypsoData,
-                                 NLayers=5)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_CemaNeigeGR6J,InputsModel=InputsModel,
-                               IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(116.482,0.500,72.733,1.224,0.278,30.333,0.977,2.776)
-OutputsModel <- RunModel_CemaNeigeGR6J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR1A.R b/tests/example_RunModel_GR1A.R
deleted file mode 100644
index 35779d49..00000000
--- a/tests/example_RunModel_GR1A.R
+++ /dev/null
@@ -1,40 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## conversion of example data from daily to yearly time-step
-TabSeries       <- data.frame(BasinObs$DatesR,BasinObs$P,BasinObs$E,BasinObs$T,BasinObs$Qmm)
-TimeFormat      <- "daily"
-NewTimeFormat   <- "yearly"
-ConvertFun      <- c("sum","sum","mean","sum")
-YearFirstMonth  <- 09;
-NewTabSeries    <- SeriesAggreg(TabSeries=TabSeries,TimeFormat=TimeFormat,
-                                NewTimeFormat=NewTimeFormat,ConvertFun=ConvertFun,
-                                YearFirstMonth=YearFirstMonth)
-BasinObs        <- NewTabSeries
-names(BasinObs) <- c("DatesR","P","E","T","Qmm")
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR1A,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%Y")=="1990"),
-               which(format(BasinObs$DatesR,format="%Y")=="1999"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR1A,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(0.840)
-OutputsModel <- RunModel_GR1A(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR2M.R b/tests/example_RunModel_GR2M.R
deleted file mode 100644
index 6ac7e261..00000000
--- a/tests/example_RunModel_GR2M.R
+++ /dev/null
@@ -1,38 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## conversion of example data from daily to monthly time-step
-TabSeries       <- data.frame(BasinObs$DatesR,BasinObs$P,BasinObs$E,BasinObs$T,BasinObs$Qmm)
-TimeFormat      <- "daily"
-NewTimeFormat   <- "monthly"
-ConvertFun      <- c("sum","sum","mean","sum")
-NewTabSeries    <- SeriesAggreg(TabSeries=TabSeries,TimeFormat=TimeFormat,
-                                NewTimeFormat=NewTimeFormat,ConvertFun=ConvertFun)
-BasinObs        <- NewTabSeries
-names(BasinObs) <- c("DatesR","P","E","T","Qmm")
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR2M,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%m/%Y")=="01/1990"),
-               which(format(BasinObs$DatesR,format="%m/%Y")=="12/1999"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR2M,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(265.072,1.040)
-OutputsModel <- RunModel_GR2M(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR4H.R b/tests/example_RunModel_GR4H.R
deleted file mode 100644
index af99a43a..00000000
--- a/tests/example_RunModel_GR4H.R
+++ /dev/null
@@ -1,28 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123003)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4H,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/03/2004 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/2008 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4H,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(521.113,-2.918,218.009,4.124)
-OutputsModel <- RunModel_GR4H(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR4J.R b/tests/example_RunModel_GR4J.R
deleted file mode 100644
index 192eee96..00000000
--- a/tests/example_RunModel_GR4J.R
+++ /dev/null
@@ -1,28 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(257.238,1.012,88.235,2.208)
-OutputsModel <- RunModel_GR4J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR5J.R b/tests/example_RunModel_GR5J.R
deleted file mode 100644
index 4825f5ff..00000000
--- a/tests/example_RunModel_GR5J.R
+++ /dev/null
@@ -1,28 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR5J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR5J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(245.918,1.027,90.017,2.198,0.434)
-OutputsModel <- RunModel_GR5J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_RunModel_GR6J.R b/tests/example_RunModel_GR6J.R
deleted file mode 100644
index 38797485..00000000
--- a/tests/example_RunModel_GR6J.R
+++ /dev/null
@@ -1,28 +0,0 @@
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR6J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR6J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(242.257,0.637,53.517,2.218,0.424,4.759)
-OutputsModel <- RunModel_GR6J(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-## efficiency criterion: Nash-Sutcliffe Efficiency
-InputsCrit  <- CreateInputsCrit(FUN_CRIT=ErrorCrit_NSE,InputsModel=InputsModel,
-                                RunOptions=RunOptions,Qobs=BasinObs$Qmm[Ind_Run])
-OutputsCrit <- ErrorCrit_NSE(InputsCrit=InputsCrit,OutputsModel=OutputsModel)
-cat(paste("  Crit  ",OutputsCrit$CritName,"  ",round(OutputsCrit$CritValue,4),"\n",sep=""))
-
diff --git a/tests/example_SeriesAggreg.R b/tests/example_SeriesAggreg.R
deleted file mode 100644
index 6018c5d6..00000000
--- a/tests/example_SeriesAggreg.R
+++ /dev/null
@@ -1,18 +0,0 @@
-## load of catchment data
-require(airGR)
-
-## preparation of the initial time series data frame at the daily time-step
-data(L0123002)
-TabSeries  <- data.frame(BasinObs$DatesR,BasinObs$P,BasinObs$E,BasinObs$T,BasinObs$Qmm)
-TimeFormat <- "daily"
-
-## conversion at the monthly time-step
-NewTimeFormat <- "monthly"
-ConvertFun    <- c("sum","sum","mean","sum")
-NewTabSeries  <- SeriesAggreg(TabSeries=TabSeries,TimeFormat,NewTimeFormat,ConvertFun)
-
-## conversion at the yearly time-step
-NewTimeFormat <- "yearly"
-ConvertFun <- c("sum","sum","mean","sum")
-NewTabSeries <- SeriesAggreg(TabSeries=TabSeries,TimeFormat,NewTimeFormat,ConvertFun)
-
diff --git a/tests/example_TransfoParam.R b/tests/example_TransfoParam.R
deleted file mode 100644
index 4067affd..00000000
--- a/tests/example_TransfoParam.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR4J model
-    Xraw  <- matrix( c( 221.41, -3.63,  30.00, 1.37,
-                        347.23, -1.03,  60.34, 1.76,
-                        854.06, -0.10, 148.41, 2.34),
-                        ncol=4,byrow=TRUE)
-    Xtran <- TransfoParam(ParamIn=Xraw ,Direction="RT",FUN_TRANSFO=TransfoParam_GR4J)
-
-## transformation Transformed->Raw for the GR4J model
-    Xtran <- matrix( c( +3.60, -2.00, +3.40, -9.10,
-                        +3.90, -0.90, +4.10, -8.70,
-                        +4.50, -0.10, +5.00, -8.10),
-                        ncol=4,byrow=TRUE)
-    Xraw  <- TransfoParam(ParamIn=Xtran,Direction="TR",FUN_TRANSFO=TransfoParam_GR4J)
diff --git a/tests/example_TransfoParam_CemaNeige.R b/tests/example_TransfoParam_CemaNeige.R
deleted file mode 100644
index 35bb6b19..00000000
--- a/tests/example_TransfoParam_CemaNeige.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the CemaNeige module
-    Xraw  <- matrix( c( 0.19, 1.73,
-                        0.39, 2.51,
-                        0.74, 4.06),
-                        ncol=2,byrow=TRUE)
-    Xtran <- TransfoParam_CemaNeige(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the CemaNeige module
-    Xtran <- matrix( c( -6.26, +0.55,
-                        -2.13, +0.92,
-                        +4.86, +1.40)
-                        ,ncol=2,byrow=TRUE)
-    Xraw  <- TransfoParam_CemaNeige(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR1A.R b/tests/example_TransfoParam_GR1A.R
deleted file mode 100644
index 1727f629..00000000
--- a/tests/example_TransfoParam_GR1A.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR1A model
-    Xraw  <- matrix( c(  0.625,
-                         1.250,
-                         1.875),
-                        ncol=1,byrow=TRUE)
-    Xtran <- TransfoParam_GR1A(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR1A model
-    Xtran <- matrix( c( -5.00,
-                        +0.00,
-                        +5.00),
-                        ncol=1,byrow=TRUE)
-    Xraw  <- TransfoParam_GR1A(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR2M.R b/tests/example_TransfoParam_GR2M.R
deleted file mode 100644
index 3c9b3862..00000000
--- a/tests/example_TransfoParam_GR2M.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR2M model
-    Xraw  <- matrix( c( 221.41, -1.00,
-                        347.23, +0.00,
-                        854.01, +1.00),
-                        ncol=2,byrow=TRUE)
-    Xtran <- TransfoParam_GR2M(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR2M model
-    Xtran <- matrix( c( +3.60, -5.00,
-                        +3.90, +0.00,
-                        +4.50, +5.00),
-                        ncol=2,byrow=TRUE)
-    Xraw  <- TransfoParam_GR2M(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR4H.R b/tests/example_TransfoParam_GR4H.R
deleted file mode 100644
index f711b8bf..00000000
--- a/tests/example_TransfoParam_GR4H.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR4H model
-    Xraw  <- matrix( c( 221.41, -3.63,  30.00, 1.37,
-                        347.23, -1.03,  60.34, 1.76,
-                        854.06, -0.10, 148.41, 2.34),
-                        ncol=4,byrow=TRUE)
-    Xtran <- TransfoParam_GR4H(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR4H model
-    Xtran <- matrix( c( +3.60, -2.00, +3.40, -9.10,
-                        +3.90, -0.90, +4.10, -8.70,
-                        +4.50, -0.10, +5.00, -8.10),
-                        ncol=4,byrow=TRUE)
-    Xraw  <- TransfoParam_GR4H(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR4J.R b/tests/example_TransfoParam_GR4J.R
deleted file mode 100644
index e3a6e12f..00000000
--- a/tests/example_TransfoParam_GR4J.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR4J model
-    Xraw  <- matrix( c( 221.41, -3.63,  30.00, 1.37,
-                        347.23, -1.03,  60.34, 1.76,
-                        854.06, -0.10, 148.41, 2.34),
-                        ncol=4,byrow=TRUE)
-    Xtran <- TransfoParam_GR4J(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR4J model
-    Xtran <- matrix( c( +3.60, -2.00, +3.40, -9.10,
-                        +3.90, -0.90, +4.10, -8.70,
-                        +4.50, -0.10, +5.00, -8.10),
-                        ncol=4,byrow=TRUE)
-    Xraw  <- TransfoParam_GR4J(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR5J.R b/tests/example_TransfoParam_GR5J.R
deleted file mode 100644
index 24ea140c..00000000
--- a/tests/example_TransfoParam_GR5J.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR5J model
-    Xraw  <- matrix( c( 221.41, -2.65,  27.11, 1.37, -0.76,
-                        347.23, -0.64,  60.34, 1.76,  0.30,
-                        854.01, -0.10, 148.41, 2.34,  0.52),
-                        ncol=5,byrow=TRUE)
-    Xtran <- TransfoParam_GR5J(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR5J model
-    Xtran <- matrix( c( +3.60, -1.70, +3.30, -9.10, -0.70,
-                        +3.90, -0.60, +4.10, -8.70, +0.30,
-                        +4.50, -0.10, +5.00, -8.10, +0.50),
-                        ncol=5,byrow=TRUE)
-    Xraw  <- TransfoParam_GR5J(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_TransfoParam_GR6J.R b/tests/example_TransfoParam_GR6J.R
deleted file mode 100644
index 5ba50062..00000000
--- a/tests/example_TransfoParam_GR6J.R
+++ /dev/null
@@ -1,15 +0,0 @@
-require(airGR)
-
-## transformation Raw->Transformed for the GR6J model
-    Xraw  <- matrix( c( 221.41, -1.18,  27.11, 1.37, -0.18,  20.09,
-                        347.23, -0.52,  60.34, 1.76,  0.02,  54.60,
-                        854.06,  0.52, 148.41, 2.34,  0.22, 148.41),
-                        ncol=6,byrow=TRUE)
-    Xtran <- TransfoParam_GR6J(ParamIn=Xraw ,Direction="RT")
-
-## transformation Transformed->Raw for the GR6J model
-    Xtran <- matrix( c( +3.60, -1.00, +3.30, -9.10, -0.90, +3.00,
-                        +3.90, -0.50, +4.10, -8.70, +0.10, +4.00,
-                        +4.50, +0.50, +5.00, -8.10, +1.10, +5.00),
-                        ncol=6,byrow=TRUE)
-    Xraw  <- TransfoParam_GR6J(ParamIn=Xtran,Direction="TR")
diff --git a/tests/example_plot_OutputsModel.R b/tests/example_plot_OutputsModel.R
deleted file mode 100644
index 3cd8e8d2..00000000
--- a/tests/example_plot_OutputsModel.R
+++ /dev/null
@@ -1,54 +0,0 @@
-#### example 1 without snow module
-
-## load of catchment data
-require(airGR)
-data(L0123001)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_GR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_GR4J,InputsModel=InputsModel,IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(734.568,-0.840,109.809,1.971)
-OutputsModel <- RunModel(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param,
-                         FUN_MOD=RunModel_GR4J)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-
-#### example 2 with snow module
-
-## load of catchment data
-require(airGR)
-data(L0123002)
-
-## preparation of the InputsModel object
-InputsModel <- CreateInputsModel(FUN_MOD=RunModel_CemaNeigeGR4J,DatesR=BasinObs$DatesR,
-                                 Precip=BasinObs$P,PotEvap=BasinObs$E,TempMean=BasinObs$T,
-                                 HypsoData=BasinInfo$HypsoData,NLayers=5)
-
-## run period selection
-Ind_Run <- seq(which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="01/01/1990 00:00"),
-               which(format(BasinObs$DatesR,format="%d/%m/%Y %H:%M")=="31/12/1999 00:00"))
-
-## preparation of the RunOptions object
-RunOptions <- CreateRunOptions(FUN_MOD=RunModel_CemaNeigeGR4J,InputsModel=InputsModel,
-                               IndPeriod_Run=Ind_Run)
-
-## simulation
-Param <- c(408.774,2.646,131.264,1.174,0.962,2.249)
-OutputsModel <- RunModel(InputsModel=InputsModel,RunOptions=RunOptions,Param=Param,
-                         FUN_MOD=RunModel_CemaNeigeGR4J)
-
-## results preview
-plot_OutputsModel(OutputsModel=OutputsModel,Qobs=BasinObs$Qmm[Ind_Run])
-
-
-- 
GitLab