diff --git a/tests/example_Calibration.R b/tests/example_Calibration.R deleted file mode 100644 index 694d336776d4443460e4e6164fe0b6f081d0cdff..0000000000000000000000000000000000000000 --- 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 efc267216aae50534cb72094a5926149edcb0494..0000000000000000000000000000000000000000 --- 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 c19c510563a54d66b956c418b592ea60c084ffa6..0000000000000000000000000000000000000000 --- 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 b240b8cf3c92053f08e856c1b317394dc3427f5a..0000000000000000000000000000000000000000 --- 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 134bd7d727d780f094194060e17ae74b8bcd0b82..0000000000000000000000000000000000000000 --- 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 6ed335ae3ace2b498373be27105adbc419d743a4..0000000000000000000000000000000000000000 --- 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 f543b5d560507e56f69cb769d57df76c18b07bba..0000000000000000000000000000000000000000 --- 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 a24e2098f3e90c84f23e883029f00c145c8efe11..0000000000000000000000000000000000000000 --- 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 35779d49c99d0a13f19ef448d100fa695440d24b..0000000000000000000000000000000000000000 --- 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 6ac7e2611e916a9968f4352030b2a0fb2141c61d..0000000000000000000000000000000000000000 --- 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 af99a43a54750385c18377c00737d388808407e8..0000000000000000000000000000000000000000 --- 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 192eee968468dadddc2dee7005b07bfae59d0105..0000000000000000000000000000000000000000 --- 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 4825f5ffc79190a77d25cf5092df5347fc50485f..0000000000000000000000000000000000000000 --- 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 387974853231877da469fbacacc3fb017ec5fd70..0000000000000000000000000000000000000000 --- 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 6018c5d6cb008f2ca81886f283ac4fbccd7e832f..0000000000000000000000000000000000000000 --- 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 4067affd733cba1932df1c91705cbc482575dc74..0000000000000000000000000000000000000000 --- 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 35bb6b1955942bc6de8c750f7ccdb21ff283109a..0000000000000000000000000000000000000000 --- 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 1727f6299fd8283e17171d9dbbf9e37fadb8e89f..0000000000000000000000000000000000000000 --- 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 3c9b3862efdbfe63cbd95b5c60d0575ed2e9f2c1..0000000000000000000000000000000000000000 --- 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 f711b8bf63b44bddc5e945d5f37f816c8501abd5..0000000000000000000000000000000000000000 --- 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 e3a6e12ff1f3964619c0784f34ea408ff3acc2cf..0000000000000000000000000000000000000000 --- 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 24ea140c8281229de97e2e8ea8e6ec65178335e6..0000000000000000000000000000000000000000 --- 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 5ba500625e140654e8d06a09efe3a4e34e989769..0000000000000000000000000000000000000000 --- 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 3cd8e8d2a5bf08990156919fbfb1f5a7239b3b9c..0000000000000000000000000000000000000000 --- 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]) - -