Minor changes to allow simple broadcasting (add and sub)

Acceptance criteria: Given a cube of phenomenon 'ph' extensive in x,y,t collapse over t and calculate mean to give cube 'ph_av' calculate anomaly cube 'ph_anomaly' as ph_anomaly = ph - ph_av (analogous to array broadcasting in numpy)