all three baby functions seem to work for single timeset datasets (bcc). now to check with multiple
issue with adding timesets to anomaly variable. will try with using an object that increases with each loop based on the previous one
all 3 baby functions work, put together in MEANS function with loop for all dx/dy data chunks. Now to be tested