#!/usr/bin/env Rscript
################################################################
#script to aggregate time series to temporal composites
#1. create search intervals 
#2. search files to aggregate
#3. calculate mean for temporal steps
################################################################
library(raster)
library(lubridate)
rasterOptions(chunksize=1e+06, progress="text")

args<-commandArgs(trailingOnly=TRUE)
yr<-args[1]
ts_var<-args[2]

FOI<-list.files(path="./", pattern=paste0("swat_",".*.tif$"),full.names = TRUE,recursive = FALSE)
print(FOI)

files<-unique(grep(paste(paste0(ts_var,"_",yr),collapse="|"), FOI, value=TRUE))
files_stack<-stack(files)
annual<-calc(files_stack, sum, na.rm=TRUE)

writeRaster(annual, filename=paste0("output/",ts_var,"_",yr,".tif"), format = "GTiff", options="COMPRESS=DEFLATE")