#!/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(TRUE)
print(args)
date<-args[1]
ts_var<-args[2]
print("date:")
print(date)

print(format(as.Date(date), "%Y%m%d"))

FOI<-list.files(pattern="*.tif", full.names = TRUE)
print("files to be processed:")
print(FOI)

files_stack <- stack(FOI)

sum_files <- calc(files_stack, sum, na.rm=TRUE)

writeRaster(mean_files, paste0(ts_var,"_",gsub("-", "", date),'_monthly_sum.tif', sep = ''), format = 'GTiff', 
            options="COMPRESS=DEFLATE", overwrite = T)