Resolve "Adjust TransfoParam_GR5H (production store)"

Closes #50 (closed)

R/TransfoParam_GR5H.R

@@ -2,13 +2,13 @@ TransfoParam_GR5H <- function(ParamIn, Direction) {
 
   ## number of model parameters
   NParam <- 5L
-  
-  
+
+
   ## check arguments
   isVecParamIn <- is.vector(ParamIn)
   if (isVecParamIn) {
     ParamIn <- matrix(ParamIn, nrow = 1)
-  }  
+  }
   if (!inherits(ParamIn, "matrix")) {
     stop("'ParamIn' must be of class 'matrix'")
   }
@@ -18,32 +18,32 @@ TransfoParam_GR5H <- function(ParamIn, Direction) {
   if (ncol(ParamIn) != NParam) {
     stop(sprintf("the GR5H model requires %i parameters", NParam))
   }
-  
-  
+
+
   ## transformation
   if (Direction == "TR") {
     ParamOut <- ParamIn
-    ParamOut[, 1] <- exp(1.5 * ParamIn[, 1])                       ### GR5H X1 (production store capacity)    
-    ParamOut[, 2] <- sinh(ParamIn[, 2])                            ### GR5H X2 (groundwater exchange coefficient)      
-    ParamOut[, 3] <- exp(ParamIn[, 3])                             ### GR5H X3 (routing store capacity)       
+    ParamOut[, 1] <- exp(ParamIn[, 1])                             ### GR5H X1 (production store capacity)
+    ParamOut[, 2] <- sinh(ParamIn[, 2])                            ### GR5H X2 (groundwater exchange coefficient)
+    ParamOut[, 3] <- exp(ParamIn[, 3])                             ### GR5H X3 (routing store capacity)
     ParamOut[, 4] <- 480 + (480 - 0.01) * (ParamIn[, 4] - 10) / 20 ### GR5H X4 (unit hydrograph time constant)
     ParamOut[, 5] <- (ParamIn[, 5] + 10) / 20                      ### GR5H X5 (groundwater exchange coefficient 2)
   }
   if (Direction == "RT") {
     ParamOut <- ParamIn
-    ParamOut[, 1] <- log(ParamIn[, 1]) / 1.5                       ### GR5H X1 (production store capacity)    
-    ParamOut[, 2] <- asinh(ParamIn[, 2])                           ### GR5H X2 (groundwater exchange coefficient)      
-    ParamOut[, 3] <- log(ParamIn[, 3])                             ### GR5H X3 (routing store capacity)       
+    ParamOut[, 1] <- log(ParamIn[, 1])                             ### GR5H X1 (production store capacity)
+    ParamOut[, 2] <- asinh(ParamIn[, 2])                           ### GR5H X2 (groundwater exchange coefficient)
+    ParamOut[, 3] <- log(ParamIn[, 3])                             ### GR5H X3 (routing store capacity)
     ParamOut[, 4] <- (ParamIn[, 4] - 480) * 20 / (480 - 0.01) + 10 ### GR5H X4 (unit hydrograph time constant)
     ParamOut[, 5] <- ParamIn[, 5] * 20 - 10                        ### GR5H X5 (groundwater exchange coefficient 2)
   }
 
-  
+
   ## check output
   if (isVecParamIn) {
     ParamOut <- as.vector(ParamOut)
   }
-  
+
   return(ParamOut)
 
 }