diff --git a/DESCRIPTION b/DESCRIPTION
index fdd7953ed95eaeca2b95f0c2cc6321e14c0741d9..73ffe97d2809398dedf2758cc77173dcfde5d435 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.4.8.0
+Version: 1.4.8.1
 Date: 2020-05-04
 Authors@R: c(
   person("Laurent", "Coron", role = c("aut", "trl"), comment = c(ORCID = "0000-0002-1503-6204")),
diff --git a/NEWS.md b/NEWS.md
index 4707b6d264e606fad316d75730047e944018a675..e2248d14e7ebf697bb13c9acb03ed29014ef654c 100644
--- a/NEWS.md
+++ b/NEWS.md
@@ -4,7 +4,7 @@
 
 
 
-### 1.4.8.0 Release Notes (2020-05-04)
+### 1.4.8.1 Release Notes (2020-05-04)
 
 #### New features
 
diff --git a/R/SeriesAggreg2.default.R b/R/SeriesAggreg2.default.R
index d50bc01bded95b16a16c1cb1c0f1bc21c88c14f1..fd4d44e3e15b78cb4891452ff4de8d29ec2473b4 100644
--- a/R/SeriesAggreg2.default.R
+++ b/R/SeriesAggreg2.default.R
@@ -28,6 +28,11 @@ SeriesAggreg2.default <- function(TabSeries,
   if (inherits(TabSeries, "CemaNeige")) {
     if (inherits(TabSeries, "InputsModel")) {
       CemaNeigeLayers <- TabSeries[grep("^Layer", names(TabSeries))]
+      CemaNeigeLayers <- lapply(seq_along(CemaNeigeLayers), function(iLayer) {
+        tmp <- CemaNeigeLayers[[iLayer]]
+        names(tmp) <- paste(names(CemaNeigeLayers)[iLayer], names(tmp), sep = ".")
+        tmp
+      })
     }
     if (inherits(TabSeries, "OutputsModel") ){
       CemaNeigeLayers <- TabSeries$CemaNeigeLayers
@@ -36,7 +41,7 @@ SeriesAggreg2.default <- function(TabSeries,
       tmp <- cbind(TabSeries$DatesR, as.data.frame(iLayer))
       res <- SeriesAggreg2(tmp, Format = Format,
                            TimeFormat = TimeFormat, NewTimeFormat = NewTimeFormat,
-                           ConvertFun = rep("sum", ncol(tmp)-1),
+                           ConvertFun = .AggregOutputs(gsub("[.].*", "", colnames(tmp)[-1L])),
                            YearFirstMonth = YearFirstMonth, TimeLag = TimeLag,
                            verbose = verbose)
       res <- res[, -1L]
@@ -49,7 +54,7 @@ SeriesAggreg2.default <- function(TabSeries,
   NewTabSeries <- SeriesAggreg2(TabSeries = TabSeries2,
                                 Format = Format,
                                 TimeFormat = TimeFormat, NewTimeFormat = NewTimeFormat,
-                                ConvertFun = rep("sum", ncol(TabSeries2)-1),
+                                ConvertFun = .AggregOutputs(colnames(TabSeries2)[-1L]),
                                 YearFirstMonth = YearFirstMonth, TimeLag = TimeLag,
                                 verbose = verbose)
   TabSeries$zzz <- NULL