change in IDF.agg add year to output

R/IDF.R

@@ -82,7 +82,8 @@ IDF.agg <- function(data,ds,na.accept = 0,
     data.agg <- pbsapply(ds,agg.ts,simplify=TRUE,cl=cl)  
     
     df <- data.frame(xdat = as.vector(data.agg), ds = rep(ds,each=length(data.agg[,1])))
-    df$station <- rep(station,length(df$xdat))
+    df$station <- station
+    df$year <- rep(unique(as.POSIXlt(data.s[,names[1]])$year+1900),length(ds))
     return(df) # maxima for all durations at one station
   }
   # which stations should be used?