calc_climatology_stack.r

calc_climatology_stack=function(temp) {
                            require(raster)
                            res=list()
                            res$mean_WT=calc(temp,mean,na.rm=T);
                            res$sd_WT=calc(temp,sd,na.rm=T);
                            res$q33=calc(temp,fun=function(x) quantile(x,probs=c(0.33),na.rm=T));
                            res$median=calc(temp,fun=function(x) quantile(x,probs=c(0.5),na.rm=T));
                            res$q66=calc(temp,fun=function(x) quantile(x,probs=c(0.66),na.rm=T));
                            return(res)
                            }