diff --git a/R/RunModel_LRQ.R b/R/RunModel_LRQ.R
index 0d512d624a6de0127e5e637784ea9e88f7a6dc0e..46bcd6889487d8b50d8fb4e4da02271f2a69f19e 100644
--- a/R/RunModel_LRQ.R
+++ b/R/RunModel_LRQ.R
@@ -66,11 +66,10 @@ RunModel_LRQ <- function(InputsModel, RunOptions, Param, QcontribDown) {
 
   ## Parameters set up
   TParam <- Param[1L]
-  KParam <- Param[2L]
+  KParam <- Param[2L] # KParam en mm/pdt
 
   PT <-  InputsModel$LengthHydro * 1e3 / TParam / RunOptions$FeatFUN_MOD$TimeStep
-  PK <- KParam * PT
-
+  PK <- sqrt(KParam * InputsModel$BasinAreas[-length(InputsModel$BasinAreas)] * 1e-3) * PT #m3/2 pdt1/2
 
   ## Lag set up
 
@@ -120,8 +119,6 @@ RunModel_LRQ <- function(InputsModel, RunOptions, Param, QcontribDown) {
   QsimDown_input <- matrix(QsimDown *
                        InputsModel$BasinAreas[length(InputsModel$BasinAreas)] * 1e3, ncol = 1)
 
-
-
   for (upstream_basin in seq_len(NbUpBasins)) {
     Qupstream_m3 <- c(IniStates[[upstream_basin]],
                       InputsModel$Qupstream[IndPeriod1, upstream_basin])