Estou trabalhando com a observação diária de dados climáticos organizados em arquivos .nc. Eu os li usando o comando stack do pacote raster. Cada arquivo (correspondente a um ano) é um elemento RasterStack com as seguintes características:

class       : RasterStack 
dimensions  : 360, 720, 259200, 365  (nrow, ncol, ncell, nlayers)
resolution  : 0.5, 0.5  (x, y)
extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax) 

cada camada é a varredura de valores de um dia.
Gostaria de somar as camadas para calcular os valores mensais. Acredito que a solução deve estar usando calc ou stackApply {raster}, mas não consegui encontrar uma maneira de somar da camada x para a camada y ou uma maneira de subconjunto do RasterStack antes da soma.

Eu preparei umarquivo de exemplo com apenas 12 camadas (para reduzir o tamanho).

Não sei exatamente como propor um código, desculpe, mas deve ser algo como:

library(raster)
setwd("myfolder")
data<-stack(mydata.nc)

datasum<- stackApply(data, ??? ,fun=sum)

Obrigado