Season NA handling
The function Season
has currently no possibility to ignore NA values. In R terms that is na.omit, I think. For example it might be worth considering to calculate the annual mean of monthly anomalies if only 11 instead 12 values are present.
I think it would require some effort as currently the values for all months are added and the final array is simply divided by the number of months. This division would have to be specific for each element.
I will create a branch and give it a go.
Any comments or suggestions?