diff --git a/DESCRIPTION b/DESCRIPTION
index bd1b6c00e838c1ffe5520613a1c3f9b6e545330a..482c670545fd5de727843460876ce5836028eb59 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.3.0.15
-Date: 2019-05-20
+Version: 1.3.0.16
+Date: 2019-05-21
 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@irstea.fr"),
diff --git a/NEWS.rmd b/NEWS.rmd
index 7ce587182b048847da38852944fef69a8bbba4e5..7f7b926dae33bbcf8d9236dc4390c5eac188927f 100644
--- a/NEWS.rmd
+++ b/NEWS.rmd
@@ -14,7 +14,7 @@ output:
 
 
 
-### 1.3.0.15 Release Notes (2019-05-20)
+### 1.3.0.16 Release Notes (2019-05-21)
 
 
 #### New features
diff --git a/src/frun_CEMANEIGE.f b/src/frun_CEMANEIGE.f
index ef866eb13014143476456acc5b5a4ea8d8787145..f3389245801fb3e86a16efdad7d0fe7727510c6e 100644
--- a/src/frun_CEMANEIGE.f
+++ b/src/frun_CEMANEIGE.f
@@ -11,7 +11,7 @@
      &                             MeanAnSolidPrecip    , ! [double]  value of annual mean solid precip [mm/y]
      &                             NParam               , ! [integer] number of model parameter
      &                             Param                , ! [double]  parameter set
-     &                             NStates              , ! [integer] number of state variables used for model initialising = 2
+     &                             NStates              , ! [integer] number of state variables used for model initialisation = 4
      &                             StateStart           , ! [double]  state variables used when the model run starts
      &                             IsHyst               , ! [logical] whether we should use the linear hysteresis or not
      &                             NOutputs             , ! [integer] number of output series