diff --git a/DESCRIPTION b/DESCRIPTION
index 482c670545fd5de727843460876ce5836028eb59..02e1ad18f07a0263d73121b3f1de70abc49141d2 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.3.0.16
+Version: 1.3.0.17
 Date: 2019-05-21
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
diff --git a/NEWS.rmd b/NEWS.rmd
index 7f7b926dae33bbcf8d9236dc4390c5eac188927f..39370af9e5f091d015837de913f13897875158da 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
 
 
 
-### 1.3.0.16 Release Notes (2019-05-21)
+### 1.3.0.17 Release Notes (2019-05-21)
 
 
 #### New features
diff --git a/R/RunModel_CemaNeigeGR4H.R b/R/RunModel_CemaNeigeGR4H.R
index 8bba1294523c5a5690a63a57fd2f8f7650dec2ba..2247f5731251d0c7995b7f61655299eb62308310 100644
--- a/R/RunModel_CemaNeigeGR4H.R
+++ b/R/RunModel_CemaNeigeGR4H.R
@@ -59,9 +59,9 @@ RunModel_CemaNeigeGR4H <- function(InputsModel,RunOptions,Param){
     ##Call_DLL_CemaNeige_________________________
       for(iLayer in 1:NLayers){
         if (!IsHyst) {
-          StateStartCemaNeige <- RunOptions$IniStates[(7 + 20 + 40) + c(iLayer, iLayer+NLayers)]
+          StateStartCemaNeige <- RunOptions$IniStates[(7 + 20*24 + 40*24) + c(iLayer, iLayer+NLayers)]
         } else {
-          StateStartCemaNeige <- RunOptions$IniStates[(7 + 20 + 40) + c(iLayer, iLayer+NLayers, iLayer+2*NLayers, iLayer+3*NLayers)]
+          StateStartCemaNeige <- RunOptions$IniStates[(7 + 20*24 + 40*24) + c(iLayer, iLayer+NLayers, iLayer+2*NLayers, iLayer+3*NLayers)]
         }
         RESULTS <- .Fortran("frun_CemaNeige",PACKAGE="airGR",
                         ##inputs