diff --git a/DESCRIPTION b/DESCRIPTION
index 06032de184c1900ae4a1d4352c558600464be4f4..3ae179d55e76b3896ff0ea23958ac3af250cf372 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,8 +1,8 @@
 Package: airGR
 Type: Package
 Title: Suite of GR Hydrological Models for Precipitation-Runoff Modelling
-Version: 1.4.3.88
-Date: 2020-04-24
+Version: 1.6.0.0
+Date: 2020-04-05
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
   person("Olivier", "Delaigue", role = c("aut", "cre"), comment = c(ORCID = "0000-0002-7668-8468"), email = "airGR@inrae.fr"),
diff --git a/NEWS.md b/NEWS.md
index c6053c7d66ecd29c7733d0f247151e3041187bbb..776a9bfafe5e88c134e2d826c50e65ef3c2c33ef 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -2,6 +2,12 @@
 
 
 
+### 1.6.0.0 Release Notes (2020-04-05)
+
+#### New features
+- Semi-distributed GR models.
+
+____________________________________________________________________________________
 
 
 ### 1.4.3.88 Release Notes (2020-04-24)
diff --git a/R/CreateInputsModel.R b/R/CreateInputsModel.R
index 510ea517d40e075b206f1af096ec9ae24fc23829..0afb9f3ec2f9b76457248cf1e477062f062dc52a 100644
--- a/R/CreateInputsModel.R
+++ b/R/CreateInputsModel.R
@@ -4,6 +4,7 @@ CreateInputsModel <- function(FUN_MOD,
                               PotEvap = NULL,
                               TempMean = NULL, TempMin = NULL, TempMax = NULL,
                               ZInputs = NULL, HypsoData = NULL, NLayers = 5,
+                              Q_upstr = NULL, L_hyd = NULL, S_basin = NULL, 
                               verbose = TRUE) {
     
     
@@ -183,6 +184,35 @@ CreateInputsModel <- function(FUN_MOD,
       }
     }
     
+    if (!is.null(Q_upstr) & !is.null(L_hyd) & !is.null(S_basin)) {
+      ObjectClass <- c(ObjectClass, "SD")
+    }
+    # olivier: add a check to avoid any missing argument: 
+    # either we have the three of them (then it is SD) or none (then it is lumped)
+    
+    if("SD" %in% ObjectClass){
+      if (!("daily" %in% ObjectClass) & !("hourly" %in% ObjectClass)){
+        stop("Only daily and hourly time steps can be used in a semi-distributed mode")
+      }
+      if (!is.matrix(Q_upstr) | !is.matrix(L_hyd) | !is.matrix(S_basin)) {
+        stop("'Q_upstr', 'L_hyd' and 'S_basin' must be matrices of numeric values")
+      }
+      if (!is.numeric(Q_upstr) | !is.numeric(L_hyd) | !is.numeric(S_basin)) {
+        stop("'Q_upstr', 'L_hyd' and 'S_basin' must be matrices of numeric values")
+      }
+      if (ncol(Q_upstr) != ncol(L_hyd)){
+        stop("'Q_upstr' and 'L_hyd' must have the same number of columns")
+      }
+      if (ncol(Q_upstr)+1 != ncol(S_basin)){
+        stop("'S_basin' must have one column more than 'Q_upstr' and 'L_hyd'")
+      }
+      if (nrow(L_hyd) != 1 | nrow(S_basin) != 1){
+        stop("'L_hyd' and 'S_basin' must have only one row")
+      }
+      if(any(is.na(Q_upstr))){
+        stop("'Q_upstr' cannot contain any NA value")
+      }
+    }
     
     ##check_NA_values
     BOOL_NA <- rep(FALSE, length(DatesR))
@@ -296,6 +326,11 @@ CreateInputsModel <- function(FUN_MOD,
                                          LayerFracSolidPrecip = RESULT$LayerFracSolidPrecip,
                                          ZLayers              = RESULT$ZLayers))
     }
+    if("SD" %in% ObjectClass){
+      InputsModel <- c(InputsModel, list(Q_upstr = Q_upstr, 
+                                         L_hyd=L_hyd,
+                                         S_basin = S_basin)) 
+    }
     
     class(InputsModel) <- c("InputsModel", ObjectClass)