diff --git a/DESCRIPTION b/DESCRIPTION
index 9bcac62f718892d2de3c8dd9216214d5f5c90e71..6dcf9055bbeb1c7a6f0b2217429fd1205ab38b9b 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
Package: airGR
Type: Package
Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.0.10.0
+Version: 1.0.10.1
Date: 2018-05-22
Authors@R: c(
person("Laurent", "Coron", role = c("aut", "trl")),
diff --git a/NEWS.rmd b/NEWS.rmd
index de0315e23ee40357d8012e3b3ace6fe8be46a067..26be15f8176a8091977e978bc5d6a7722321df9d 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
-### 1.0.10.0 Release Notes (2018-05-22)
+### 1.0.10.1 Release Notes (2018-05-22)
____________________________________________________________________________________
diff --git a/R/CreateRunOptions.R b/R/CreateRunOptions.R
index 873921ba0ee0ee45dfd7837ab3e3f1c6c4264ad3..3ca652f54e785c3064029012f2de9b3aec450a46 100644
--- a/R/CreateRunOptions.R
+++ b/R/CreateRunOptions.R
@@ -206,7 +206,7 @@ CreateRunOptions <- function(FUN_MOD, InputsModel, IndPeriod_WarmUp = NULL, IndP
Outputs_all <- c(Outputs_all,"PotEvap", "Precip", "Prod", "Pn", "Ps", "AE", "Perc", "PR", "Q9", "Q1", "Rout", "Exch",
"AExch1", "AExch2", "AExch", "QR", "QRExp", "Exp", "QD", "Qsim"); }
if(identical(FUN_MOD,RunModel_GR2M)){
- Outputs_all <- c(Outputs_all,"PotEvap","Precip","AE","Perc","P3","Exch","Prod","Rout","Qsim"); }
+ Outputs_all <- c(Outputs_all,"PotEvap", "Precip", "AE", "Pn", "Perc", "PR", "Exch", "Prod", "Rout", "Qsim"); }
if(identical(FUN_MOD,RunModel_GR1A)){
Outputs_all <- c(Outputs_all,"PotEvap","Precip","Qsim"); }
if("CemaNeige" %in% ObjectClass){
diff --git a/R/RunModel_GR2M.R b/R/RunModel_GR2M.R
index f767bd211c42b6871784bc18a7e01e83b64b96f6..83fc84de2d346df457c0b10ba4e2628d78cb3416 100644
--- a/R/RunModel_GR2M.R
+++ b/R/RunModel_GR2M.R
@@ -1,7 +1,7 @@
RunModel_GR2M <- function(InputsModel,RunOptions,Param){
NParam <- 2;
- FortranOutputs <- c("PotEvap","Precip","AE","Perc","P3","Exch","Prod","Rout","Qsim");
+ FortranOutputs <- c("PotEvap", "Precip", "AE", "Pn", "Perc", "PR", "Exch", "Prod", "Rout", "Qsim")
##Arguments_check
if(inherits(InputsModel,"InputsModel")==FALSE){ stop("InputsModel must be of class 'InputsModel' \n"); return(NULL); }
diff --git a/man/RunModel_GR2M.Rd b/man/RunModel_GR2M.Rd
index bb32089517137bcbaf9ca48c94a98f96b272e4eb..c16a1b87995b4e1999ec8ffc464d01c5787ec393 100644
--- a/man/RunModel_GR2M.Rd
+++ b/man/RunModel_GR2M.Rd
@@ -20,8 +20,8 @@ RunModel_GR2M(InputsModel, RunOptions, Param)
\item{Param}{[numeric] vector of 2 parameters
\tabular{ll}{
-GR2M X1 \tab production store capacity [mm] \cr
-GR2M X2 \tab groundwater exchange coefficient [-] \cr
+GR2M X1 \tab production store capacity [mm] \cr
+GR2M X2 \tab groundwater exchange coefficient [-] \cr
}}
}
@@ -29,11 +29,31 @@ GR2M X2 \tab groundwater exchange coefficient [-] \cr
\value{
[list] list containing the function outputs organised as follows:
\tabular{ll}{
- \emph{$DatesR } \tab [POSIXlt] series of dates \cr
- \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/month] \cr
- \emph{$Precip } \tab [numeric] series of input total precipitation [mm/month] \cr
- \emph{$Qsim } \tab [numeric] series of simulated discharge [mm/month] \cr
- \emph{$StateEnd} \tab [numeric] states at the end of the run (production store level and routing store level) [mm], \cr\tab see \code{\link{CreateIniStates}} for more details \cr
+ \emph{$DatesR } \tab [POSIXlt] series of dates \cr
+ \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/month] \cr
+ \emph{$Precip } \tab [numeric] series of input total precipitation [mm/month] \cr
+ \emph{$Qsim } \tab [numeric] series of simulated discharge [mm/month] \cr
+ \emph{$StateEnd} \tab [numeric] states at the end of the run (production store level and routing store level) [mm], \cr\tab see \code{\link{CreateIniStates}} for more details \cr
+ }
+ (refer to the provided references or to the package source code for further details on these model outputs)
+}
+
+
+\value{
+[list] list containing the function outputs organised as follows:
+ \tabular{ll}{
+ \emph{$DatesR } \tab [POSIXlt] series of dates \cr
+ \emph{$PotEvap } \tab [numeric] series of input potential evapotranspiration [mm/month] \cr
+ \emph{$Precip } \tab [numeric] series of input total precipitation [mm/month] \cr
+ \emph{$AE } \tab [numeric] series of actual evapotranspiration [mm/month] \cr
+ \emph{$Pn } \tab [numeric] series of net rainfall (P1) [mm/month] \cr
+ \emph{$Perc } \tab [numeric] series of percolation (P2) [mm/month] \cr
+ \emph{$PR } \tab [numeric] series of PR=Pn+Perc (P3) [mm/month] \cr
+ \emph{$Exch } \tab [numeric] series of potential exchange between catchments [mm/month] \cr
+ \emph{$Prod } \tab [numeric] series of production store level [mm] \cr
+ \emph{$Rout } \tab [numeric] series of routing store level [mm] \cr
+ \emph{$Qsim } \tab [numeric] series of simulated discharge [mm/month] \cr
+ \emph{$StateEnd} \tab [numeric] states at the end of the run (res. levels, UH1 levels, UH2 levels) [mm], \cr\tab see \code{\link{CreateIniStates}} for more details \cr
}
(refer to the provided references or to the package source code for further details on these model outputs)
}
diff --git a/src/frun_GR2M.f b/src/frun_GR2M.f
index 99dd8920d4c64b5109d72d0488733f2b9315b24d..2cadd865f623dc8b023fd12dec411c0ebc8187e0 100644
--- a/src/frun_GR2M.f
+++ b/src/frun_GR2M.f
@@ -169,9 +169,9 @@ C Variables storage
MISC( 1)=E ! PE ! [numeric] observed potential evapotranspiration [mm/month]
MISC( 2)=P ! Precip ! [numeric] observed total precipitation [mm/month]
MISC( 3)=AE ! AE ! [numeric] actual evapotranspiration [mm/month]
- MISC( 4)=P1 ! P1 ! [numeric] net rainfall (P1) [mm/month]
- MISC( 5)=P2 ! P2 ! [numeric] percolation (P2) [mm/month]
- MISC( 6)=P3 ! P3 ! [numeric] P3=P1+P2 [mm/month]
+ MISC( 4)=P1 ! Pn ! [numeric] net rainfall (P1) [mm/month]
+ MISC( 5)=P2 ! Perc ! [numeric] percolation (P2) [mm/month]
+ MISC( 6)=P3 ! PR ! [numeric] P3=P1+P2 [mm/month]
MISC( 7)=EXCH ! EXCH ! [numeric] groundwater exchange (EXCH) [mm/month]
MISC( 8)=St(1) ! Prod ! [numeric] production store level (St(1)) [mm]
MISC( 9)=St(2) ! Rout ! [numeric] routing store level (St(2)) [mm]