UPDATE: adjusted what is being returned by EotCycle so that we also get pred variables. added anomalise.R to create anomaly stacks bsed on either a reference raster or the mean of the stack