def METHOD_NAME(self):
    shape = 25, 180, 360
    dims = 'time', 'lat', 'lon'
    self.ts_a = np.linspace(1, 25, 25)
    self.ts_a_mean = np.linspace(1, 25, 25)
    self.ts_a_count = np.array(25 * [100])
    self.ts_a_std = np.array(25 * [0.0])
    self.ts_b = np.linspace(0, 1, 25)
    self.ts_b_mean = np.linspace(0, 1, 25)
    self.ts_b_count = np.array(25 * [100])
    self.ts_b_std = np.array(25 * [0.0])
    self.cube = new_cube(
        time_periods=25,
        variables=dict(
            A=xr.DataArray(
                np.broadcast_to(self.ts_a.reshape((25, 1, 1)), shape),
                dims=dims),
            B=xr.DataArray(
                np.broadcast_to(self.ts_b.reshape((25, 1, 1)), shape),
                dims=dims)
        )
    )
    self.cube = self.cube.chunk(chunks=dict(time=1, lat=180, lon=180))